Trace element partitioning in HP-LT metamorphic assemblages during subduction-related metamorphism, Ile de Groix, France: a detailed LA-ICPMS study

Devolatilization reactions and subsequent transfer of fluid from subducted oceanic crust into the overlying mantle wedge are important processes, which are responsible for the specific geochemical characteristics of subduction-related metamorphic rocks, as well as those of arc magmatism. To better understand the geochemical fingerprint induced by fluid mobilization during dehydration and rehydration processes related to subduction zone metamorphism, the trace element and rare earth element (REE) distribution patterns in HP-LT metamorphic assemblages in eclogite-, blueschist- and greenschist-facies rocks of the Ile de Groix were obtained by laser ablation inductively coupled plasma mass spectrometry (LA-ICPMS) analysis. This study focuses on 10 massive basic rocks representing former hydrothermally altered mid-ocean ridge basalts (MORB), four banded basic rocks of volcano-sedimentary origin and one micaschist. The main hosts for incompatible trace elements are epidote (REE, Th, U, Pb, Sr), garnet [Y, heavy REE (HREE)], phengite (Cs, Rb, Ba, B), titanite [Ti, Nb, Ta, REE; HREE > LREE (light REE)], rutile (Ti, Nb, Ta) and apatite (REE, Sr). The trace element contents of omphacite, amphibole, albite and chlorite are low. The incompatible trace element contents of minerals are controlled by the stable metamorphic mineral assemblage and directly related to the appearance, disappearance and reappearance of minerals, especially epidote, garnet, titanite, rutile and phengite, during subduction zone metamorphism. Epidote is a key mineral in the trace element exchange process because of its large stability field, ranging from lower greenschist- to blueschist- and eclogite-facies conditions. Different generations of epidote are generally observed and related to the coexisting phases at different stages of the metamorphic cycle (e.g. lawsonite, garnet, titanite). Epidote thus controls most of the REE budget during the changing P-T conditions along the prograde and retrograde path. Phengite also plays an important role in determining the large ion lithophile element (LILE) budget, as it is stable to high P-T conditions. The breakdown of phengite causes the release of LILE during retrogression. A comparison of trace element abundances in whole-rocks and minerals shows that the HP-LT metamorphic rocks largely retain the geochemical characteristics of their basic, volcano-sedimentary and pelitic protoliths, including a hydrothermal alteration overprint before the subduction process. A large part of the incompatible trace elements remained trapped in the rocks and was recycled within the various metamorphic assemblages stable under changing metamorphic conditions during the subduction process, indicating that devolatilization reactions in massive basic rocks do not necessarily imply significant simultaneous trace element and REE release.

A global marine sedimentary response to the end-Permian mass extinction: Examples from southern Turkey and the western United States

The end-Permian mass extinction greatly diminished marine diversity and brought about a whole-scale restructuring of marine ecosystems; these ecosystem changes also profoundly affected the sedimentary record. Data presented here, attained through facies analyses of strata deposited during the immediate aftermath of the end-Permian mass extinction (southern Turkey) and at the close of the Early Triassic (southwestern United States), in combination with a literature review, show that sedimentary systems were profoundly affected by: (1) a reduction in biotic diversity and abundance and (2) long-term environmental fluctuations that resulted from the end-Permian crisis. Lower Triassic strata display widespread microbialite and carbonate seafloor fan development and contain indicators of suppressed infaunal bioturbation such as flat-pebble conglomerates and wrinkle structures (facies considered unusual in post-Cambrian subtidal deposits). Our observations suggest that depositional systems, too, respond to biotic crises, and that certain facies may act as barometers of ecologic and environmental change independent of fossil assemblage analyses. Close investigation of facies changes during other critical times in Earth history may serve as an important tool in interpreting the ecology of metazoans and their environment.

Petrology, geochemistry and geodynamic implications of Jurassic island arc magmatism as revealed by mafic volcanic rocks in the Mesozoic low-grade sequence, eastern Rhodope, Bulgaria

In the eastern Bulgarian Rhodope, mafic extrusive rocks and underlying greenschists are found in the Mesozoic low-grade unit, which represents the northern extension of similar sequences including the Evros ophiolites in Thrace (Greece). Both rock types define a suite of low-Ti tholeiitic basalts to transitional boninitic basaltic andesites and andesites and associated metapyroclastites (greenschists), intruded at its base by diorite dikes of a boninitic affinity. Mafic lavas and greenschists display large ion lithophile element (LILE) enrichment relative to high-field strength elements (HFSE), flat REE patterns of a slight light REE depletion, a strong island arc tholeiite (IAT) and weak MORB-like signature. All these rocks are characterized by negative Nb anomalies ascribed to arc lavas. They have positive epsilon Nd(i) values in the range of +4.87 to +6.09, approaching the lower limit of MORB-like source, and relatively high ((207)Pb/(204)Pb)(i) (15.57-15.663) at low ((206)Pb/(204)Pb)(i) (18.13-18.54) ratios. The Nd isotopic compositions coupled with trace element data imply a dominantly depleted MORB-like mantle source and a contribution of subduction modified LILE-enriched component derived from the mantle wedge. The diorite dike has a low eNdi value of -2.61 and is slightly more Pb radiogenic ((207)Pb/(204)Pb)(i) (15.64) and ((206)Pb/(204)Pb)(i) (18.56), respectively, reflecting crustal contamination. Petrologic and geochemical data indicate that the greenschists and mafic extrusive rocks represent a magmatic assemblage formed in an island arc setting. The magmatic suite is interpreted as representing an island arc-accretionary complex related to the southward subduction of the Meliata-Maliac ocean under the supra-subduction back-arc Vardar ocean/island arc system. Magmatic activity appears to have initiated in the north during the inception of the island arc system by the Early-Middle Jurassic time in the eastern Rhodope that most likely graded to back-arc spreading southwards as represented by the Late Jurassic MORB-type Samothraki Island ophiolites. This tectonic scenario is further constrained by paleotectonic reconstructions. The arc-trench system collided with the Rhodope in the Late Jurassic times. (c) 2007 Elsevier B.V. All rights reserved.

Fluid migration through thrust faults in the Helvetic nappes (Western Swiss Alps)

The delta(18)O, delta(13)C and Sr-87/Sr-86 values of calcite and organic matter were measured for carbonate mylonites from numerous thrusts in the Helvetic Alps. Carbonate mylonites in most of the thrusts retain essentially unaltered protolith delta(18)O and delta(13)C values, consistent with there having been little to no advection of isotopically distinct fluid through these faults. Only carbonate mylonites from the basal thrusts of the largest nappes have delta(18)O and/or delta(13)C values that differ from those of their protoliths. The zone of isotopic alteration/exchange is confined to c. 10 to 20 meters of these fault contacts. We propose the fluids that migrated through these faults contained variable amounts of organically derived carbon and radiogenic strontium, and were probably derived from dewatering of the sedimentary rocks and prograde metamorphic reactions in the nappes' root zones. Apart from the basal thrusts of the largest nappes that behaved as narrow, laterally extensive conduits for fluids, there is little isotopic evidence that large quantities of fluids passed through most of the carbonate-hosted thrusts in the Helvetic Alps.

Elemental (C/N ratios) and isotopic (δ15N(org), δ13C(org)) compositions of sedimentary organic matter from a high-altitude mountain lake (Meidsee, 2661 m a.s.l., Switzerland): Implications for Lateglacial and Holocene Alpine landscape evolution

The deposition of Late Pleistocene and Holocene sediments in the high-altitude lake Meidsee (located at an altitude of 2661 m a.s.l. in the Southwestern Alps) strikingly coincided with global ice-sheet and mountain-glacier decay in the Alpine forelands and the formation of perialpine lakes. Radiocarbon ages of bottom-core sediments point out (pre-) Holocene ice retreat below 2700 m a.s.l., at about 16, 13, 10, and 9 cal. kyr BP. The Meidsee sedimentary record therefore provides information about the high-altitude Alpine landscape evolution since the Late Pleistocene/Holocene deglaciation in the Swiss Southwestern Alps. Prior to 5 cal. kyr BP, the C/N ratio and the isotopic composition of sedimentary organic matter (delta N-15(org), delta C-13(org)) indicate the deposition of algal-derived organic matter with limited input of terrestrial organic matter. The early Holocene and the Holocene climatic optimum (between 7.0 and 5.5 cal. kyr BP) were characterized by low erosion (decreasing magnetic susceptibility, chi) and high content of organic matter (C-org > 13 wt.%), enriched in C-13(org) (>-18 parts per thousand) with a low C/N (similar to 10) ratio, typical of modern algal matter derived from in situ production. During the late Holocene, there was a long-term increasing contribution of terrestrial organic matter into the lake (C/N > 11), with maxima between 2.4 and 0.9 cal. kyr BP. A major environmental change took place 800 years ago, with an abrupt decrease in the relative contribution of terrestrial organic material into the lake compared with aquatic organic material which subsequently largely dominated (C/N drop from 16 to 10). Nonetheless, this event was marked by a rise in soil erosion (chi), in nutrients input (N and P contents) and in anthropogenic lead deposition, suggesting a human disturbance of Alpine ecosystems 800 years ago. Indeed, this time period coincided with the migration of the Walser Alemannic people in the region, who settled at relatively high altitude in the Southwestern Alps for farming and maintaining Alpine passes.

Settlement evolution, field systems and sedimentary processes in a Mediterranean mid-mountain (Montseny Massif, Catalunya)

This paper presents the first results of a current research project about human – environmental interactions in the Montseny Massif. Our work sets out to integrate two research lines in the studied area: - Archaeological and archaeo-morphological surveys in a lower part of the mountains in order to characterize the evolution of the settlements and field systems. - The geological and geomorphological characterization of the slope and terrace deposits in relation with field systems and archaeological data. First results point out the intensive occupation of these inland areas during the Iberian and the Roman periods. Post-Roman sediments show different processes of erosion.

High-pressure metamorphism in the Adula nappe - Guide to the excursion of the Swiss Society of Mineralogy and Petrology (September 29 - October 5, 1991)

1st day: Lithology and structure of the northern Adula nappe around Zervreila 2nd day: High-pressure rocks of the Suretta nappe, the middle Adula nappe and its Mesozoic cover - Eclogite near Innerferrara, Suretta nappe - Crossite-bearing prasinite from schistes lustres near Nufenen - Eclogites south of Hinterrhein, Adula nappe - Blueschists and eclogites from Neu-Wahli, Misox zone 3 days: Eclogite boudin and associated whiteschists in the uppermost Calanca valley, middle Adula nappe 4th day: Eclogites, associated metapelites and granitoid gneisses of Trescolmen, middle Adula nappe 5th day: The ultramafic-mafic suite of the Cima Lunga nappe around Cima di Gagnone 6th day: Garnet peridotites and eclogites from Alpe Arami, Cima Lunga nappe

Stratigraphy and petrology of the Miocene Montjuic delta: (Barcelona, Spain).

The Neogene rift in the Catalan Coastal Ranges, which is located in the NE part of the Eastern Iberian Margin, corresponds to a system of grabens formed at the nort h - we s t e rn edge of the Valencia Trough. In the central part of the Catalan Coastal Ranges are the Valls - Peneds half-graben in the onshore and the Barcelona half-graben in the offshore, which are separated by the Garraf and the Collserola-Montnegre horsts. Montjuc hill is a tilted block, which is located to the S of the Barcelona city, between the Collserola-Montn egre horst and the Barcelona half-graben . The Middle Miocene section of Montjuc is constituted by an alternation of conglomerate, sandstone, mudstone and marlstone beds. The Montjuc section was divided into four lithostratigraphic units from base to top: (1) The Morrot conglomerate and sandstone Unit, interpreted as delta plain deposits; (2) the Castell conglomerate, sandstone and mudstone Unit considered as proximal delta front deposits; (3) the Miramar marlstone Unit attributed

Macrofossils in Raraku Lake (Easter Island) integrated with sedimentary and geochemical records towards a palaeoecological synthesis for the last 34,000 years

Macrofossil analysis of a composite 19 m long sediment core from Rano Raraku Lake (Easter Island)was related to litho-sedimentary and geochemical features of the sediment. Strong stratigraphical patterns are shown by indirect gradient analyses of the data. The good correspondence between the stratigraphical patterns derived from macrofossil (Correspondence Analysis) and sedimentary and geochemical data (Principal Component Analysis) shows that macrofossil associations provide sound palaeolimnological information in conjunction with sedimentary data. The main taphonomic factors in fluencing the macrofossil assemblages are run-off from the catchment, the littoral plant belt, and the depositional environment within the basin. Five main stages during the last 34,000 calibrated years BP (cal yr BP) are characterised from the lithological, geochemical, and macrofossil data. From 34 to 14.6 cal kyr BP (last glacial period) the sediments were largely derived from the catchment, indicating a high energy lake environment with much erosion and run-off bringing abundant plant trichomes, lichens, and mosses into the centre of Raraku Lake.

The El Masnou Infralittoral sedimentary environment (Barcelona province, NW Mediterranean Sea): morphology and Holocene seismic stratigraphy

A detailed analysis of the morphology and the Holocene seismic and sequence stratigraphy and architecture of the infralittoral sedimentary environment of the El Masnou coast (Catalonia, NW Mediterranean Sea) was carried out using multibeam bathymetry and GeoPulse seismic data. This environment extends down to 26-30 m water depth, and is defined morphologically by two depositional wedges whose seafloor is affected by erosive furrows, slides, fields of large- and small-scale wavy bedforms, and dredging trenches and pits. Erosive terraces are also identified in the transition domain toward the inner continental shelf. The Holocene stratigraphy of the infralittoral environment is defined by two major seismic sequences (lower and upper), each one formed by internal seismic units. The sequences and units are characterised by downlapping surfaces made up of deposits formed by progradation of coastal lithosomes. The stratigraphy and stratal architecture, displaying a retrogradational arrangement with progradational patterns of minor order, were controlled by different sea-level positions. The stratigraphic division represents the coastal response to the last fourth-order transgressive and highstand conditions, modulated by small-scale sea-level oscillations (≈1-2 m) of fith to sixth order. This study also highlights the advantage of an integrated analysis using acoustic/seismic methods for practical assessment of the anthropogenic effects on infralittoral domains based on the association of marine geological observations.

Seismic characterization of heterogeneous sedimentary and fractured rocks based on Biot's theory of poroelasticity

Understanding and quantifying seismic energy dissipation, which manifests itself in terms of velocity dispersion and attenuation, in fluid-saturated porous rocks is of considerable interest, since it offers the perspective of extracting information with regard to the elastic and hydraulic rock properties. There is increasing evidence to suggest that wave-induced fluid flow, or simply WIFF, is the dominant underlying physical mechanism governing these phenomena throughout the seismic, sonic, and ultrasonic frequency ranges. This mechanism, which can prevail at the microscopic, mesoscopic, and macroscopic scale ranges, operates through viscous energy dissipation in response to fluid pressure gradients and inertial effects induced by the passing wavefield. In the first part of this thesis, we present an analysis of broad-band multi-frequency sonic log data from a borehole penetrating water-saturated unconsolidated glacio-fluvial sediments. An inherent complication arising in the interpretation of the observed P-wave attenuation and velocity dispersion is, however, that the relative importance of WIFF at the various scales is unknown and difficult to unravel. An important generic result of our work is that the levels of attenuation and velocity dispersion due to the presence of mesoscopic heterogeneities in water-saturated unconsolidated clastic sediments are expected to be largely negligible. Conversely, WIFF at the macroscopic scale allows for explaining most of the considered data while refinements provided by including WIFF at the microscopic scale in the analysis are locally meaningful. Using a Monte-Carlo-type inversion approach, we compare the capability of the different models describing WIFF at the macroscopic and microscopic scales with regard to their ability to constrain the dry frame elastic moduli and the permeability as well as their local probability distribution. In the second part of this thesis, we explore the issue of determining the size of a representative elementary volume (REV) arising in the numerical upscaling procedures of effective seismic velocity dispersion and attenuation of heterogeneous media. To this end, we focus on a set of idealized synthetic rock samples characterized by the presence of layers, fractures or patchy saturation in the mesocopic scale range. These scenarios are highly pertinent because they tend to be associated with very high levels of velocity dispersion and attenuation caused by WIFF in the mesoscopic scale range. The problem of determining the REV size for generic heterogeneous rocks is extremely complex and entirely unexplored in the given context. In this pilot study, we have therefore focused on periodic media, which assures the inherent self- similarity of the considered samples regardless of their size and thus simplifies the problem to a systematic analysis of the dependence of the REV size on the applied boundary conditions in the numerical simulations. Our results demonstrate that boundary condition effects are absent for layered media and negligible in the presence of patchy saturation, thus resulting in minimum REV sizes. Conversely, strong boundary condition effects arise in the presence of a periodic distribution of finite-length fractures, thus leading to large REV sizes. In the third part of the thesis, we propose a novel effective poroelastic model for periodic media characterized by mesoscopic layering, which accounts for WIFF at both the macroscopic and mesoscopic scales as well as for the anisotropy associated with the layering. Correspondingly, this model correctly predicts the existence of the fast and slow P-waves as well as quasi and pure S-waves for any direction of wave propagation as long as the corresponding wavelengths are much larger than the layer thicknesses. The primary motivation for this work is that, for formations of intermediate to high permeability, such as, for example, unconsolidated sediments, clean sandstones, or fractured rocks, these two WIFF mechanisms may prevail at similar frequencies. This scenario, which can be expected rather common, cannot be accounted for by existing models for layered porous media. Comparisons of analytical solutions of the P- and S-wave phase velocities and inverse quality factors for wave propagation perpendicular to the layering with those obtained from numerical simulations based on a ID finite-element solution of the poroelastic equations of motion show very good agreement as long as the assumption of long wavelengths remains valid. A limitation of the proposed model is its inability to account for inertial effects in mesoscopic WIFF when both WIFF mechanisms prevail at similar frequencies. Our results do, however, also indicate that the associated error is likely to be relatively small, as, even at frequencies at which both inertial and scattering effects are expected to be at play, the proposed model provides a solution that is remarkably close to its numerical benchmark. -- Comprendre et pouvoir quantifier la dissipation d'énergie sismique qui se traduit par la dispersion et l'atténuation des vitesses dans les roches poreuses et saturées en fluide est un intérêt primordial pour obtenir des informations à propos des propriétés élastique et hydraulique des roches en question. De plus en plus d'études montrent que le déplacement relatif du fluide par rapport au solide induit par le passage de l'onde (wave induced fluid flow en anglais, dont on gardera ici l'abréviation largement utilisée, WIFF), représente le principal mécanisme physique qui régit ces phénomènes, pour la gamme des fréquences sismiques, sonique et jusqu'à l'ultrasonique. Ce mécanisme, qui prédomine aux échelles microscopique, mésoscopique et macroscopique, est lié à la dissipation d'énergie visqueuse résultant des gradients de pression de fluide et des effets inertiels induits par le passage du champ d'onde. Dans la première partie de cette thèse, nous présentons une analyse de données de diagraphie acoustique à large bande et multifréquences, issues d'un forage réalisé dans des sédiments glaciaux-fluviaux, non-consolidés et saturés en eau. La difficulté inhérente à l'interprétation de l'atténuation et de la dispersion des vitesses des ondes P observées, est que l'importance des WIFF aux différentes échelles est inconnue et difficile à quantifier. Notre étude montre que l'on peut négliger le taux d'atténuation et de dispersion des vitesses dû à la présence d'hétérogénéités à l'échelle mésoscopique dans des sédiments clastiques, non- consolidés et saturés en eau. A l'inverse, les WIFF à l'échelle macroscopique expliquent la plupart des données, tandis que les précisions apportées par les WIFF à l'échelle microscopique sont localement significatives. En utilisant une méthode d'inversion du type Monte-Carlo, nous avons comparé, pour les deux modèles WIFF aux échelles macroscopique et microscopique, leur capacité à contraindre les modules élastiques de la matrice sèche et la perméabilité ainsi que leur distribution de probabilité locale. Dans une seconde partie de cette thèse, nous cherchons une solution pour déterminer la dimension d'un volume élémentaire représentatif (noté VER). Cette problématique se pose dans les procédures numériques de changement d'échelle pour déterminer l'atténuation effective et la dispersion effective de la vitesse sismique dans un milieu hétérogène. Pour ce faire, nous nous concentrons sur un ensemble d'échantillons de roches synthétiques idéalisés incluant des strates, des fissures, ou une saturation partielle à l'échelle mésoscopique. Ces scénarios sont hautement pertinents, car ils sont associés à un taux très élevé d'atténuation et de dispersion des vitesses causé par les WIFF à l'échelle mésoscopique. L'enjeu de déterminer la dimension d'un VER pour une roche hétérogène est très complexe et encore inexploré dans le contexte actuel. Dans cette étude-pilote, nous nous focalisons sur des milieux périodiques, qui assurent l'autosimilarité des échantillons considérés indépendamment de leur taille. Ainsi, nous simplifions le problème à une analyse systématique de la dépendance de la dimension des VER aux conditions aux limites appliquées. Nos résultats indiquent que les effets des conditions aux limites sont absents pour un milieu stratifié, et négligeables pour un milieu à saturation partielle : cela résultant à des dimensions petites des VER. Au contraire, de forts effets des conditions aux limites apparaissent dans les milieux présentant une distribution périodique de fissures de taille finie : cela conduisant à de grandes dimensions des VER. Dans la troisième partie de cette thèse, nous proposons un nouveau modèle poro- élastique effectif, pour les milieux périodiques caractérisés par une stratification mésoscopique, qui prendra en compte les WIFF à la fois aux échelles mésoscopique et macroscopique, ainsi que l'anisotropie associée à ces strates. Ce modèle prédit alors avec exactitude l'existence des ondes P rapides et lentes ainsi que les quasis et pures ondes S, pour toutes les directions de propagation de l'onde, tant que la longueur d'onde correspondante est bien plus grande que l'épaisseur de la strate. L'intérêt principal de ce travail est que, pour les formations à perméabilité moyenne à élevée, comme, par exemple, les sédiments non- consolidés, les grès ou encore les roches fissurées, ces deux mécanismes d'WIFF peuvent avoir lieu à des fréquences similaires. Or, ce scénario, qui est assez commun, n'est pas décrit par les modèles existants pour les milieux poreux stratifiés. Les comparaisons des solutions analytiques des vitesses des ondes P et S et de l'atténuation de la propagation des ondes perpendiculaires à la stratification, avec les solutions obtenues à partir de simulations numériques en éléments finis, fondées sur une solution obtenue en 1D des équations poro- élastiques, montrent un très bon accord, tant que l'hypothèse des grandes longueurs d'onde reste valable. Il y a cependant une limitation de ce modèle qui est liée à son incapacité à prendre en compte les effets inertiels dans les WIFF mésoscopiques quand les deux mécanismes d'WIFF prédominent à des fréquences similaires. Néanmoins, nos résultats montrent aussi que l'erreur associée est relativement faible, même à des fréquences à laquelle sont attendus les deux effets d'inertie et de diffusion, indiquant que le modèle proposé fournit une solution qui est remarquablement proche de sa référence numérique.