The tectonically controlled emplacement of a vertically sheeted gabbro-pyroxenite intrusion: Feeder-zone of an ocean-island volcano (Fuerteventura, Canary Islands)

The Miocene PX1 gabbro-pyroxenite pluton, Fuerteventura, Canary Islands, is a 3.5 x 5.5 km shallow-level intrusion (0.15-0.2 GPa and 1100-1120 degrees C), interpreted as the feeder-zone to an ocean-island volcano. It displays a vertical magmatic banding expressed in five 50 to 100 metre-wide NNE-SSW trending alkaline gabbro sequences alternating with pyroxenites. This emplacement geometry was controlled by brittle to ductile shear zones, generated by a regional E-W extensional tectonic setting that affected Fuerteventura during the Miocene. At a smaller scale, the PX1 gabbro and pyroxenite bands consist of metre-thick differentiation units, which suggest emplacement by periodic injection of magma pulses as vertical dykes that amalgamated, similarly to a sub-volcanic sheeted dyke complex. Individual dykes underwent internal differentiation following a solidification front parallel to the dyke edges. This solidification front may have been favoured by a significant lateral/horizontal thermal gradient, expressed by the vertical banding in the gabbros, the fractionation asymmetry within individual dykes and the migmatisation of the wall rocks. Pyroxenitic layers result from the fractionation and accumulation of clinopyroxene +/- olivine +/- plagioclase crystals from a mildly alkaline basaltic liquid. They are interpreted as truncated differentiation sequences, from which residual melts were extracted at various stages of their chemical evolution by subsequent dyke intrusions, either next to or within the crystallising unit. Compaction and squeezing of the crystal mush is ascribed to the incoming and inflating magma pulses. The expelled interstitial liquid was likely collected and erupted along with the magma flowing through the newly injected dykes. Clinopyroxene mineral orientation - as evidenced by EBSD and micro X-ray tomography investigations - displays a marked pure-shear component, supporting the interpretation of the role of compaction in the generation of the pyroxenites. Conversely, gabbro sequences underwent minor melt extraction and are believed to represent crystallised coalesced magma batches emplaced at lower rates at the end of eruptive cycles. Clinopyroxene orientations in gabbros record a simple shear component suggesting syn-magmatic deformation parallel to observed NNE-SSW trending shear zones induced by the regional tensional stress field. This emplacement model implies a crystallisation time of 1 to 5 years for individual dykes, consistent with PX1 emplacement over less than 0.5 My. A minimum amount of approximately 150 km(3) of magma is needed to generate the pluton, part of it having been erupted through the Central Volcanic Centre of Fuerteventura. If the regional extensional tectonic regime controls the PX1 feeder-zone initiation and overall geometry, rates and volumes of magma depend on other, source-related factors. High injection rates are likely to induce intrusion growth rates larger than could be accommodated by the regional extension. In this case, dyke intrusion by propagation of a weak tip, combined with the inability of magma to circulate through previously emplaced and crystallised dykes could result in an increase of non-lithostatic pressure on previously emplaced mushy dyke walls; thus generating strong pure-shear compaction within the pluton feeder-zone and interstitial melt expulsion. These compaction-dominated processes are recorded by the cumulitic pyroxenite bands. (C) 2010 Elsevier B.V. All rights reserved.

Growth mechanisms and timing of emplacement of a vertically layered mafic intrusion in the feeder zone of an ocean island volcano (Fuerteventura, Canary islands)

Résumé pour le grand public L'île de Fuerteventura (Canaries) offre l'occasion rare d'observer les racines d'un volcan océanique édifié il y a 25 à 30 millions d'années et complètement érodé. On y voit de nombreux petits plutons de forme et composition variées, témoignant d'autant d'épisodes de l'activité magmatique. L'un de ces plutons, appelé PX1, présente une structure inhabituelle formée d'une alternance de bandes verticales d'épaisseur métrique à hectométrique de roches sombres de composition pyroxénilique ou gabbroïque. Les pyroxénites résultent clairement de l'accumulation de cristaux de pyroxènes et non de la simple solidification d'un magma? Se pose dès lors la question de la nature du processus qui a conduit à l'accumulation verticale de niveaux concentrés en pyroxènes. En effet, les litages pyroxénitiques classiques sont subhorizontaux, car ils résultent de l'accumulation gravitaire des cristaux séparés du magma dont ils cristalli¬sent par sédimentation. Cette étude vise à identifier et comprendre les mécanismes qui ont engendré ce Iitage minéralogique vertical et l'im¬portant volume de ces faciès cumulatifs. Nous nous sommes également intéressés aux conditions de pression et de température régnant au moment de la mise en place du pluton, ainsi qu'à sa durée de vie et à sa vitesse de refroidis¬sement. Enfin une approche géochimique nous a permis de préciser la nature de la source mantellique des magmas liés à cette activité magmatique. PX1 est en réalité un complexe filonien formé à des conditions de pression et de température de 1-2 kbar et 1050- 1100°C; sa construction a nécessité au moins 150 km3 de magma. L'alternance d'horizons gabbroïques et pyroxéniti¬ques représente des injections successives de magma sous la forme de filons verticaux, mis en place dans un contexte régional en extension. L'étude des orientations des minéraux dans ces faciès révèle que les horizons gabbroïques enregistrent l'extension régionale, alors que les pyroxénites sont générées par une compaction au sein du pluton. Ceci suggère que le régime des contraintes, qui était extensif lors de l'initiation de la mise en place de PX1, est pério¬diquement devenu compressif au sein même du pluton. Cette compression serait liée à des cycles de mise en place où la vitesse de croissance du pluton dépassait celle de l'extension régionale. La différenciation observée au sein de chaque horizon, depuis des pyroxénites riches en olivine jusqu'à des pyroxé¬nites à plagioclase interstitiel et des gabbros, ainsi que la composition géochimique des minéraux qui les constituent suggèrent que chaque filon vertical s'est mis en place à partir d'un magma de composition identique, puis a évolué indépendamment des autres en fonction du régime thermique et du régime des contraintes local. Lorsque le magma en train de cristalliser s'est trouvé en compression, le liquide résiduel a été séparé des cristaux déjà formés et extrait du système, laissant derrière lui une accumulation de cristaux dont la nature et les proportions dépendaient du stade de cristallisation atteint par le magma au moment de l'extraction. Ainsi, les niveaux de pyroxénites à olivine (premier minéral à cristalliser) ont été formés lorsque le magma correspondant était encore peu cristallisé; à l'inverse, les py¬roxénites riches en plagioclase (minéral plus tardif dans la séquence de cristallisation) et certains gabbros à caractère cumulatif résultent d'une compression tardive dans le processus de cristallisation du filon concerné. Les liquides résiduels extraits des niveaux pyroxénitiques sont rarement observés dans PX1, certaines poches et filonets de com¬position anorthositique pourraient en être les témoins. L'essentiel de ces liquides a probablement gagné des niveaux supérieurs du pluton, voire la surface du volcan. L'origine du régime compressif périodique affectant les filons en voie de cristallisation est attribuée aux injections suivantes de magma au sein du pluton, qui se sont succédées à un rythme plus rapide que la vitesse de consolidation des filons. Des datations U/Pb de haute précision sur des cristaux de zircon et de baddeleyite ainsi que40Ar/39Ar sur des cris¬taux d'amphibole révèlent une initiation de la mise en place de PX1 il y a 22.1 ± 0,7 Ma; celle-ci a duré quelque 0,48 ± 0,22 à 0,52 ± 0,29 Ma. Ce laps de temps est compatible avec celui nécessaire à la cristallisation des filons individuels, qui va de moins d'une année lors de l'initiation du magmatisme à 5 ans lors du maximum d'activité de PX1. La présence de cristaux résorbés enregistrant une cristallisation complexe suggère l'existence d'une chambre mag¬matique convective sous-jacente à PX1 et périodiquement rechargée. Les compositions isotopiques des roches étu¬diées révèlent une source mantellique profonde de type point chaud avec une contribution du manteau lithosphéri- que métasomatisé présent sous les îles Canaries. Résumé L'intrusion mafique Miocène PX1 fait partie du soubassement superficiel (0.15-0.2 GPa, 1100 °Q d'un volcan d'île océanique. La particularité de ce pluton est l'existence d'alternances d'unités de gabbros et de pyroxénites qui met¬tent en évidence un litage magmatique vertical (NNE-SSW). Les horizons gabbroiques et pyroxénitiques sont constitués d'unités de différenciation métriques qui suggèrent tine mise en place par injections périodiques de filons verticaux de magma formant un complexe filonien. Chaque filon vertical a subi une différenciation parallèle à un front de solidification sub-vertical parallèle aux bords du filon. Les pyroxénites résultent du fractionnement et de l'accumulation d'olivine ± clinopyroxene ± plagioclase à partir d'un magma basaltique faiblement alcalin et sont interprétées comme étant des imités de différenciation tronquées dont le liquide interstitiel a été extrait par compaction. L'orientation préférentielle des clinopyroxènes dans ces pyroxe- nites (obtenues par analyse EBSD et micro-tomographique) révèle une composante de cisaillement simple dans la genèse de ces roches, ce qui confirme cette interprétation. La compaction des pyroxénites est probablement causée par a mise en place de filons de magma suivants. Le liquide interstitiel expulsé est probablement par ces derniers. Les clinopyroxènes des gabbros, montrent une composante de cisaillement pure suggérant qu'ils sont affectés par une déformation syn-magmatique parallèle aux zones de cisaillement NNE-SSW observées autour de PX1 et liées au contexte tectonique Miocène d'extension régionale. Ceci suggère que les gabbros sont liés à des taux de mise en place faibles à la fin de cycles d'activité magmatique et sont peu ou pas affectés par la compaction. L'initiation et la géométrie de PX1 sont donc contrôlées par le contexte tectonique régional d'extension alors que les taux et les volumes de magma dépendent de facteurs liés à la source. Des taux d'injection élevés résultent probable¬ment en une croissance du pluton supérieure à la place crée par cette extension. Dans ce cas de figure, la propagation des nouveaux dykes et l'inaptitude du magma à circuler à travers les anciens dykes cristallisés pourrait causer une augmentation de la pression non-lithostatique sur ces derniers, exprimée par un cisaillement simple et l'expulsion du liquide interstitiel qu'ils contiennent (documenté par les zones de collecte anorthositiques). Les compositions en éléments majeurs et traces des gabbros et pyroxenites de PX1 sont globalement homogènes et dépendent de la nature cumulative des échantillons. Cependant, de petites variations des concentrations en éléments traces ainsi que les teneurs en éléments traces des bordures de clinopyroxenes suggèrent que ces derniers ont subi un processus de rééquilibrage et de cristallisation in situ. L'homogénéité des compositions chimiques des échantillons, ainsi que la présence de grains de clinopyroxene résorbés suggère que le complexe filonien PX1 s'est mis en place au dessus d'une chambre magmatique périodiquement rechargée dans laquelle la convection est efficace. Chaque filon est donc issu d'un même magma, mais a subi une différenciation par cristallisation in situ (jusqu'à 70% de fraction¬nement) indépendamment des autres. Dans ces filons cristallisés, les minéraux cumulatifs subissent un rééquilibrage partiel avec les liquide interstitiel avant que ce dernier ne soit expulsé lors de la compaction (mettant ainsi un terme à la différenciation). Ce modèle de mise en place signifie qu'un minimum de 150Km3 de magma est nécessaire à la genèse de PX1, une partie de ce volume ayant été émis par le 'Central Volcanic Complex' de Fuerteventura. Les rapports isotopiques radiogéniques mesurés révèlent la contribution de trois pôles mantelliques dans la genèse du magma formant PX1. Le mélange de ces pôles HIMU, DMM et EM1 refléterai l'interaction du point chaud Cana¬rien avec un manteau lithosphérique hétérogène métasomatisé. Les petites variations de ces rapports et des teneurs en éléments traces au sein des faciès pourrait refléter des taux de fusion partielle variable de la source, résultant en un échantillonnage variable du manteau lithosphérique métasomatisé lors de son interaction avec le point chaud. Des datations U/Pb de haute précision (TIMS) sur des cristaux de zircon et de baddeleyite extraits de gabbros de PX1 révèlent que l'initiation de la cristallisation du magma a eu lieu il y a 22.10±0.07 Ma et que l'activité magmatique a duré un minimum de 0.48 à 0.52 Ma. Des âges 40Ar/39Ar obtenus sur amphibole sont de 21.9 ± 0.6 à 21.8 ± 0.3 Ma, identiques aux âges U/Pb. La combinaison de ces méthodes de datations, suggère que le temps maximum nécessaire à PX1 pour se refroidir en dessous de la température de fermeture de l'amphibole est de 0.8Ma. Ceci signifie que la durée de vie de PX1 est de 520 000 à 800 000 ans. La coexistence de cristaux de baddeleyite et de zircon dans un gabbro est attribuée à son interaction avec un fluide riche en C02 relâché par les carbonatites encaissantes lors du métamorphisme de contact généré par la mise en place de PX1 environ 160 000 ans après le début de sa mise en place. Les durées de vie obtenue sont en accord avec le modèle de mise en place suggérant une durée de cristallisation poux chaque filon allant de 1 an à 5 ans. Abstract The Miocene PX1 gabbro-pyroxenite intrusion (Fuerteventura, Canary Islands), is interpreted as the shallow-level feeder-zone (0.15-0.2 GPa and 1100-1120°C), to an ocean island volcano. The particularity of PX1 is that it displays a NNE-SSW trending vertical magmatic banding expressed by alternating gabbro and pyroxeriite sequences. The gabbro and pyroxenite sequences consist of metre-thick differentiation units, which suggest emplacement by pe¬riodic injection of magma pulses as vertical dykes that amalgamated, similarly to a sub-volcanic sheeted dyke com¬plex. Individual dykes underwent internal differentiation following a solidification front (favoured by a significant lateral/horizontal thermal gradient) parallel to the dyke edges. Pyroxenitic layers result from the fractionation and accumulation of clinopyroxene ± olivine ± plagioclase crystals from a mildly alkaline basaltic liquid and are interpre¬ted as truncated differentiation sequences, from which residual melts were extracted by compaction. Clinopyroxene mineral orientation in pyroxenites (evidenced by EBSD and micro X-ray tomography analysis) display a marked pure shear component, supporting this interpretation. Compaction and squeezing of the crystal mush is ascribed to the incoming and inflating magma pulses. The resulting expelled interstitial liquid was likely collected and erupted along with the magma flowing through the newly injected dykes. Gabbro sequences represent crystallised coalesced magma batches, emplaced at lower rates at the end of eruptive cycles, and underwent minor melt extraction as evi¬denced by clinopyroxene orientations that record a simple shear component suggesting syn-magmatic deformation parallel to observed NNF.-SSW trending shear-zones induced by the regional tensional Miocene stress-field. The initiation and geometry of PX1 is controlled by the regional extensional tectonic regime whereas rates and vo¬lumes of magma depend on source-related factors. High injection rates are likely to induce intrusion growth rates larger than could be accommodated by the regional extension. In this case, dyke tip geometry and the inability of magma to circulate through previously emplaced and crystallised dykes could result in an increase of non-lithostatic pressure on previously emplaced mushy dyke walls; generating strong pure-shear compaction and interstitial melt expulsion within the feeder-zone as recorded by the cumulitic pyroxenite bands and anorthositic collection zones. The whole-rock major and trace-element chemistry of PX1 gabbros and pyroxenites is globally homogeneous and controlled by the cumulate nature of the samples (i.e. on the modal proportions of olivine, pyroxene, plagioclase and oxides). However, small variations of whole-rock trace-element contents as well as trace-element contents of clinopyroxene rims suggest that in-situ re-equilibration and crystallisation has occurred. Additionally, the global homogeneity and presence of complex zoning of rare resorbed clinopyroxene crystals suggest that the PX1 feeder- zone overlies a periodically replenished and efficiently mixed magma chamber. Each individual dyke of magma thus originated from a compositionally constant mildly alkaline magma and differentiated independently from the others reaching up to 70% fractionation. Following dyke arrest these are affected by interaction with the trapped interstitial liquid prior to its compaction-linked expulsion (thus stopping the differentiation process). This emplacement model implies that minimum amount of approximately 150 km3 of magma is needed to generate PX1, part of it having been erupted through the overlying Central Volcanic Complex of Fuerteventura. The radiogenic isotope ratios of PX1 samples reveal the contribution on three end-members during magma genesis. This mixing of the H1MU, EMI and DMM end-members could reflect the interaction of the deep-seated Canarian mantle plume with a heterogeneous metasomatic and sepentininsed lithospheric mantle. Additionally, the observed trace-element and isotopic variations within the same fades groups could reflect varying degrees of partial melting of the source region, thus tapping more or less large areas of the metasomatised lithospheric mantle during interac¬tion with the plume. High precision ID-TIMS U/Pb zircon and baddeleyite ages from the PX1 gabbro samples, indicate initiation of magma crystallisation at 22.10 ± 0.07 Ma. The magmatic activity lasted a minimum of 0.48 to 0.52 Ma. 40Ar/39Ar amphibole ages are of 21.9 ± 0.6 to 21.8 ± 0.3, identical within errors to the U/Pb ages. The combination of the 40Ar/39Ar and U/Pb datasets imply that the maximum amount of time PX1 took to cool below amphibole Tc is 0.8 Ma, suggesting PX1 lifetime of 520 000 to 800 000 years. On top of this, the coexistence of baddeleyite and zircon in a single sample is ascribed to the interaction of PX1 with C02-rich carbonatite-derived fluids released from the host-rock carbonatites during contact metamorphism 160 000 years after PX1 initiation. These ages are in agreement with the emplacement model, implying a crystallisation time of less than 1 to 5 years for individual dykes.

Complex, 3D strain patterns in a synkinematic tonalite batholith from the Aracuai Neoproterozoic orogen (Eastern Brazil): Evidence from combined magnetic and isotopic chronology studies

This work combines structural and geochronological data to improve our understanding of the mechanical behaviour of continental crust involving large amount of magma or partially melted material in an abnormally hot collisional belt. We performed a magnetic and geochronological (U/Pb) study on a huge tonalitic batholith from the Neoproterozoic Aracual belt of East Brazil to determine the strain distribution through space and time. Anisotropy of magnetic susceptibility, combined with rock magnetism investigations, supports that the magnetic fabric is a good proxy of the structural fabric. Field measurements together with the magnetic fabrics highlight the presence in the batholith of four domains characterized by contrasted magmatic flow patterns. The western part is characterized by a gently dipping, orogen-parallel (similar to NS) magmatic foliation that bears down-dip lineations, in agreement with westward thrusting onto the Sao Francisco craton. Eastward, the magmatic foliation progressively turns sub-vertical with a lineation that flips from sub-horizontal to sub-vertical over short distances. This latter domain involves an elongated corridor in which the magmatic foliation is sub-horizontal and bears an orogen-parallel lineation. Finally the fourth, narrow domain displays sub-horizontal lineations on a sub-vertical magmatic foliation oblique (similar to N150 degrees E) to the trend of the belt. U/Pb dating of zircons from the various domains revealed homogeneity in age for all samples. This, together with the lack of solid-state deformation suggests that: 1) the whole batholith emplaced during a magmatic event at similar to 580 Ma, 2) the deformation occurred before complete solidification. and 3) the various fabrics are roughly contemporaneous. The complex structural pattern mapped in the studied tonalitic batholith suggests a 3D deformation of a slowly cooling, large magmatic body and its country rock. We suggest that the development of the observed 3D flow field was promoted by the low viscosity of the middle crust that turned gravitational force as an active tectonic force combining with the East-West convergence between the Sao Francisco and Congo cratons. (C) 2012 Elsevier Ltd. All rights reserved.

Intrusion of Lamprophyre dyke and related deformation effects in the host rock salt: A case study from the Loulé Diapir, Portugal

A rock salt-lamprophyre dyke contact zone (sub-vertical, NE-SW strike) was investigated for its petrographic, mechanic and physical properties by means of anisotropy of magnetic susceptibility CAMS) and rock magnetic properties, coupled with quantitative microstructural analysis and thermal mathematical modelling. The quantitative microstructural analysis of halite texture and solid inclusions revealed good spatial correlation with AMS and halite fabrics. The fabrics of both lamprophyre and rock salt record the magmatic intrusion, "plastic" flow and regional deformation (characterized by a NW-SE trending steep foliation). AMS and microstructural analysis revealed two deformation fabrics in the rock salt: (1) the deformation fabrics in rock salt on the NW side of the dyke are associated with high temperature and high fluid activity attributed to the dyke emplacement; (2) On the opposite side of the dyke, the emplacement-related fabric is reworked by localized tectonic deformation. The paleomagnetic results suggest significant rotation of the whole dyke, probably during the diapir ascent and/or the regional Tertiary to Quaternary deformation. (C) 2014 Elsevier B.V. All rights reserved.

Petrology and chemistry of the Chalten Plutonic Complex and implications on the magmatic and tectonic evolution of the Southernmost Andes (Patagonia) during the Miocene

Résumé Scientific:Pétrologie et Géochimie du Complexe Plutonique de Chaltén et les conséquences pour l'évolution magmatique et tectonique du Andes du Sud (Patagonia) pendant le MiocèneLe sujet de cette thèse est le Complexe Plutonique de Chaltén (CHPC), situé à la frontière entre le Chili et l'Argentine, en Patagonie (49°15'S). Ce complexe s'est mis en place au début du Miocène, dans un contexte de changements tectoniques importants. La géométrie et la vitesse de migration des plaques en Patagonie a été modifiée suite l'ouverture de la plaque Farallon il y a 25Ma (Pardo-Casas and Molnar 1987) et la subduction de la ride active du Chili sous la plaque sud-américaine il y a 14Ma (Cande and Leslie 1986). Les effets de cette reconfiguration tectonique sur la morphologie et le magmatisme de la plaque supérieure sont encore sujets à discussion. Dans ce contexte, un groupe d'intrusions miocènes - telle que le CHPC - est particulièrement intriguant, car en position transitionnelle entre le batholithe patagonien et l'arc volcanique cénozoïque et récent à l'ouest, et les laves de plateau de Patagonie à l'est (Fig. 1). A cause de leur position tectonique transitoire, ces plutons isolés hors du batholithe représentent un endroit clé pour comprendre les interactions entre la tectonique à large échelle et le magmatisme en Patagonie. Ici, je présente de nouvelles données de terrain, petrologiques, géochimiques et géochronologiques dans le but de caractériser la nature du CHPC, qui était largement inconnu avant cette étude, dans le but de tester l'hypothèse de migration de l'arc et erosion par subduction.Les résultats de l'investigation géochimique (chapitre 2) montrent que le CHPC n'est qu'un exemple parmi les plutons isolés d'arrière arc ave une composition calco-alcaline caractéristique, c-à-d une signature d'arc. La plupart de ces plutons isolés ont une composition alcaline. Le CHPC, contrairement, a une signature calco-alcaline avec Κ intermédiaire, tel que le batholithe patagonien et la plupart des roches volcaniques quaternaires liées à l'arc le long des Andes.De nouvelles données géochronologiques U-Pb de haute précision sur des zircons, acquis par TIMS, sur le CHPC donnent des âges entre 17.0 et 16.4Ma. Les âges absolus sont en accord avec la séquence intrusive déduite des relations de terrain (chapitre 1). Ces données sont les premières contraintes d'âge U-Pb sur le CHPC. Elles montrent clairement que l'histoire magmatique du CHPC n'a pas de lien direct avec la subduction de la ride à cette latitude (Cande and Leslie 1986), car le complexe est au moins 6Ma plus ancien.Une comparaison en profondeur avec les autres intrusions d'âge Miocène en Patagonie révèlent - pour la première fois - une évolution temporelle intéressante. Il y a une tendance E-W distincte au magmatisme calco-alcalin entre 20-16Ma avec une diminution de l'âge vers l'est - le CHPC est l'expression la plus orientale de cette tendance. Je suggère que la relation espace-temps reflète une migration vers l'est (vers le continent) de l'arc magmatique. Je propose que le facteur principal contrôlant cette migration est la subduction rapide suite à la reconfiguration de la vitesse des plaques tectoniques après l'ouverture la plaque Farallon (à ~26Ma) qui résulterait en une déformation importante ainsi qu'à des taux élevés d'érosion dans la fosse de subduction.Les rapports d'isotopes radiogéniques (Pb, Sr, Nd) élevés, une signature 6018 basse et un rapport Th/La élevé sont des paramètres distinctifs pour les roches mafiques du CHPC. Le modèle isotopique présenté (chapitre 2) suggère que cette signature reflète une contamination de la source, dans le coin de manteau, plutôt qu'une contamination crustale. La signature des éléments en trace du CHPC indiquent que le coin de manteau a été contaminé par des composés terrigènes, le plus vraisemblablement par des sédiments paléozoïques.Les travaux de terrain, la pétrographie et la géothermobarométrie ont été utilisés dans le but de comprendre l'histoire interne du CHPC (chapitre 3). Ces données suggèrent deux niveaux distincts de cristallisation : l'un dans la croûte moyenne (6 à 4.5kbar) et l'autre à un niveau peu profond (3.5 à 2kbar). La modélisation isotopique AFC de la contamination crustale indique des taux variables d'assimilation, qui ne sont pas corrélés avec le degré de différenciation. Cela suggère que différents volumes de magma se sont différenciés en profondeur, de façon indépendante. Cela implique que le CHPC se serait formés en plusieurs puises de magmas provenant d'au moins trois sources différentes. Les textures des granodiorites et des granites indiquent des teneurs élevées en cristaux avant la mise en place et, par conséquent, des températures d'emplacement faibles. Les observations de terrain montrent que les roches mafiques sont déformées, alors que ce n'est pas le cas pour les granodiorites et granites (plus jeunes). La déformation des roches mafiques est encore sujet de recherche, afin de savoir si elle est liée à la déformation régionale en régime compressif ou à l'emplacement lui-même. Cependant, la mise en place de grand volume de magma felsique riche en cristaux suggère un régime d'extension.Scientific Abstract:Petrology and chemistry of the Chaltén Plutonic Complex and implications on the magmatic and tectonic evolution of the Southernmost Andes (Patagonia) during the MioceneThe subject of this thesis is the Chaltén Plutonic Complex (CHPC) located at the frontier between Chile and Argentina in Patagonia (at 49° 15 'Southern latitude). This complex intruded during early Miocene in a context of major tectonics changes. The plate geometry of Patagonia has been modified by changes in the plate motions after the break up of the Farallôn plate at 25Ma (Pardo-Casas and Molnar 1987) and by the subduction of the Chile spreading Ridge beneath South-America at 14 Ma (Cande and Leslie 1986). The effects of this tectonic setting on the morphology and the magmatism of the overriding plate are a matter of on-going discussion. Particularly intriguing in this context is a group of isolated Miocene intrusions - like the CHPC - which are located in a transitional position between the Patagonian Batholith and the Cenozoic and Recent volcanic arc in the West, and the Patagonian plateau lavas in the East (Fig. 1). Due to their transient tectonic position these isolated plutons outside the batholith represent a key to understanding the interaction between global-scale tectonics and magmatism in Patagonia. Here, I present new field, penological, geochemical and geochronological data to characterize the nature of the CHPC, which was largely unknown before this study, in order to test the hypothesis of time- transgressive magmatism.The results of the geochemical investigation (Chapter 2) show that the CHPC is only one among these isolated back-arc plutons with a characteristic calc-alkaline composition, i.e. arc signature. Most of these isolated intrusives have an alkaline character. The CHPC, in contrast, has a medium Κ calc-alkaline signature, like the Patagonian batholith and most of the Quaternary arc-related volcanic rocks along the Andes.New high precision TIMS U-Pb zircon dating of the CHPC yield ages between 17.0 to 16.4 Ma. The absolute ages support the sequence of intrusion relations established in the field (Chapter 1). These data are the first U-Pb age constraints on the CHPC, and clearly show that the magmatic history of CHPC has no direct link to the subduction of the ridge, since this complex is at least 6 Ma older than the time of collision of the Chile ridge at this latitude (Cande and Leslie 1986).An in-depth comparison with other intrusion of Miocene age in Patagonia reveals - for the first time - an interesting temporal pattern. There is a distinct E-W trend of calc-alkaline magmatism between 20-16 Ma with the younging of ages in the East - the CHPC is the easternmost expression of this trend. I suggest that this time-space relation reflects an eastward (landward) migration of the magmatic arc. I propose that main factor controlling this migration is the fast rates of subduction after the major reconfigurations of the plate tectonic motions after the break up of the Farallôn Plate (at -26 ) resulting in strong deformation and high rates of subduction erosion.High radiogenic isotope ratios (Pb, Sr, Nd) ratios, low 5018 signature and high Th/La ratios in mafic rocks are distinctive features of the CHPC. The presented isotopic models (Chapter 2) suggest that this signature reflects source contamination of the mantle wedge rather than crustal contamination. The trace element signature of the CHPC indicates that the mantle wedge was contaminated with a terrigenous component, most likely from Paleozoic sediments.Fieldwork, petrography and geothermobarometry were used to further unravel the internal history of the CHPC (Chapter 3). These data suggest two main levels of crystallization: one a mid crustal levels (6 to 4.5 kbar) and other a shallow level (3.5 to 2 kbar). Isotopic AFC modeling of crustal contamination indicate variable rates of assimilation, which are not correlated with the degree of differentiation. This suggests that different batches of magma differentiate independently at depths. This implies that the CHPC would have formed by several pulses of magmas from at least 3 different sources. Textures of granodiorites and granites indicate a high content of crystals previous to the emplacement and consequently low emplacement temperatures. Field observations show that the mafic rocks are deformed, whereas the (younger) granodiorites and granites are not. It is still subject of investigation whether the deformation of the mafic rocks is related to regional deformation during a compressional regime or to the emplacement it self. However, the emplacement of huge amount of crystal rich felsic magmas suggests an extensional regime.Résumé Grand PublicPétrologie et Géochimie du Complexe Plutonique de Chaltén et les conséquences pour l'évolution magmatique et tectonique du Andes du Sud (Patagonia) pendant le MiocèneLe Complexe Plutonique de Chaltén (CHPC) est un massif montagneux situé à 49°S à la frontière entre le Chili et l'Argentine, en Patagonie (région la plus au sud de l'Amérique du Sud). Il est composé de montagnes qui peuvent atteindre plus de 3000 mètres d'altitude, telles que le Cerro Fitz Roy (3400m) et le Cerro Torre (3100m). Ces montagnes sont composées de roches plutoniques, c.-à-d. des magmas qui se sont refroidis et ont cristallisés sous la surface terrestre.La composition chimique de ces roches montre que les magmas, qui ont formé ce complexe plutonique, font partie d'un volcanisme d'arc. Celui-ci se forme lorsqu'une plaque océanique plonge sous une plaque continentale. Les géologues appellent ce processus « subduction ». Dans un tel scénario, le manteau terrestre, qui se fait prendre entre ces deux plaques, fond pour former ainsi du magma. Ce magma remonte à travers la plaque continentale vers la surface. Si celui-ci atteint la surface, il forme les roches volcaniques, comme par exemple des laves. S'il n'atteint pas la surface, le magma se refroidit pour former finalement les roches plutoniques.Le long de la marge ouest d'Amérique du Sud, la plaque Nazca - qui se situe au sud-est de la plaque océanique pacifique - passe en dessous de la plaque d'Amérique du Sud. La bordure ouest du sud de la plaque sud-américaine a également été affectée par d'autres processus tectoniques, tels que des changements dramatiques dans les déplacements de plaques (il y a 25Ma) et la collision de la ride du Chili (depuis 15 Ma jusqu'à aujourd'hui). Ces caractéristiques tectoniques et magmatiques font de cette région un haut lieu pour les géologues. La plaque Nazca, s'est formée suite à l'ouverture d'une plaque océanique plus ancienne, il y a 25Ma. Cette ouverture est liée aux vitesses de subduction les plus rapides jamais connues. La ride du Chili est l'endroit où le sol de l'Océan Pacifique s'ouvre, formant deux plaques océaniques : les plaques Nazca et Antarctique. La ride du Chili subducte sous la plaque sud-américaine depuis 15Ma, en association avec la formation de grands volumes de magma ainsi que des changements morphologiques importants. La question de savoir lequel de ces changements tectoniques globaux affecte la géologie et la géographie de Patagonie a été, et est encore, discutée pendant de nombreuses années. De nombreux chercheurs suggèrent que la plupart des caractéristiques morphologiques et magmatiques en Patagonie sont liés à la subduction de la ride du Chili, mais cette suggestion est encore débattue comme le montre notre étude.Le batholithe de Patagonie du sud (SPB) est un énorme massif composé de roches plutoniques et il s'étend tout au long de la côte ouest de Patagonie (au sud de 47°S). Ces roches correspondent certainement aux racines d'un ancien arc volcanique, qui a été soulevé et érodé. Le CHPC, ainsi que d'autres petites intrusions dans la région, se situe dans une position exotique, à 100km à l'est du SPB. Certains chercheurs suggèrent que ces intrusions pourraient être liées à la subduction de la ride du Chili.Afin de débattre de cette problématique, nous avons utilisé différentes méthodes géochronologiques pour déterminer l'âge du CHPC et le comparer (a) à l'âge des roches intrusives similaires du SPB et (b) à l'âge de la collision de la ride du Chili. Dans ce travail, nous prouvons que le CHPC s'est formé au moins 7Ma avant la collision avec la ride du Chili. Sur la base des âges du CHPC et de la composition chimique de ses roches et minéraux, nous proposons que le CHPC fait partie d'un arc volcanique ancien. La migration de l'arc volcanique plus profondément dans le continent résulte de la grande vitesse de subduction entre 25 et lOMa. Des caractéristiques évidentes pour un tel processus - telles qu'une déformation importante et une vitesse d'érosion élevée - peuvent être rencontrées tout au long de la bordure ouest de l'Amérique du sud.

Alpine tectonic evolution of a Jurassic subduction-accretionary complex: Deformation, kinematics and 40Ar/39Ar age constraints on the Mesozoic low-grade schists of the Circum-Rhodope Belt in the eastern Rhodope-Thrace region, Bulgaria-Greece

Deformation of the Circum-Rhodope Belt Mesozoic (Middle Triassic to earliest Lower Cretaceous) low-grade schists underneath an arc-related ophiolitic magmatic suite and associated sedimentary successions in the eastern Rhodope-Thrace region occurred as a two-episode tectonic process: (i) Late Jurassic deformation of arc to margin units resulting from the eastern Rhodope-Evros arc-Rhodope terrane continental margin collision and accretion to that margin, and (ii) Middle Eocene deformation related to the Tertiary crustal extension and final collision resulting in the closure of the Vardar ocean south of the Rhodope terrane. The first deformational event D-1 is expressed by Late Jurassic NW-N vergent fold generations and the main and subsidiary planar-linear structures. Although overprinting, these structural elements depict uniform bulk north-directed thrust kinematics and are geometrically compatible with the increments of progressive deformation that develops in same greenschist-facies metamorphic grade. It followed the Early-Middle Jurassic magmatic evolution of the eastern Rhodope-Evros arc established on the upper plate of the southward subducting Maliac-Meliata oceanic lithosphere that established the Vardar Ocean in a supra-subduction back-arc setting. This first event resulted in the thrust-related tectonic emplacement of the Mesozoic schists in a supra-crustal level onto the Rhodope continental margin. This Late Jurassic-Early Cretaceous tectonic event related to N-vergent Balkan orogeny is well-constrained by geochronological data and traced at a regional-scale within distinct units of the Carpatho-Balkan Belt. Following subduction reversal towards the north whereby the Vardar Ocean was subducted beneath the Rhodope margin by latest Cretaceous times, the low-grade schists aquired a new position in the upper plate, and hence, the Mesozoic schists are lacking the Cretaceous S-directed tectono-metamorphic episode whose effects are widespread in the underlying high-grade basement. The subduction of the remnant Vardar Ocean located behind the colliding arc since the middle Cretaceous was responsible for its ultimate closure, Early Tertiary collision with the Pelagonian block and extension in the region caused the extensional collapse related to the second deformational event D-2. This extensional episode was experienced passively by the Mesozoic schists located in the hanging wall of the extensional detachments in Eocene times. It resulted in NE-SW oriented open folds representing corrugation antiforms of the extensional detachment surfaces, brittle faulting and burial history beneath thick Eocene sediments as indicated by 42.1-39.7 Ma Ar-40/Ar-39 mica plateau ages obtained in the study. The results provide structural constraints for the involvement components of Jurassic paleo-subduction zone in a Late Jurassic arc-continental margin collisional history that contributed to accretion-related crustal growth of the Rhodope terrane. (C) 2011 Elsevier Ltd. All rights reserved.

Alpine deformation and 40Ar/39Ar geochronology of synkinematic white mica in the Siviez-Mischabel Nappe, western Pennine Alps, Switzerland

We explore the timing of deformation and exhumation of the Siviez-Mischabel Nappe (western Swiss Alps), which has been considered a classic example of a midcrustal crystalline nappe since the studies of Argand [1916]. This study presents Ar-40/Ar-39 ages obtained on both synkinematic white mica from Permo-Triassic cover sediments and more complex white mica populations from basement gneisses of the Siviez-Mischabel and middle Pennine Nappes. Primary foliation developed in cover units by nucleation, growth, and rigid rotation of mica grains during noncoaxial Alpine deformation. Although some samples show a crenulation of this primary foliation, mica growth appears to have occurred only during the development of primary foliation, the main phase of greenschist facies deformation related to imbrication of the Siviez-Mischabel Nappe and other middle Pennine Nappes. Good agreement exists between independent estimates of the timing of deformation and reported Ar-40/Ar-39, white mica ages from cover units of the central and southern Siviez-Mischabel Nappe. In cover units from the central and southern Siviez-Mischabel regions of the study area, Ar-40/Ar-39 ages appear to date synkinematic white mica growth. Results suggest that the Siviez-Mischabel :Nappe was emplaced and developed foliation during a 5 m.y. period from 41 to 36 Ma. In cover units from the eastern Siviez-Mischabel, however, Ar-40/Ar-39 white mica ages appear to date postkinematic thermal events. These thermal events may be related to Oligocene magmatic activity in the lower Pennine Nappes or to Miocene development of the Simplon fault zone. Variations in the relation between Alpine age and grain size for cover samples from the central, eastern, and southern Siviez-Mischabel correlate well with the regional variations in temperature inferred from quartz microfabrics and the pattern of regional metamorphism. When considered in concert with other recent isotopic studies on the timing of major tectonic and thermal events in the western Swiss Alps, these data support arguments that the relative timing of events such as thrusting and back thrusting of crystalline nappes in hinterland units and exhumation of high-pressure units in the suture zone of the western Alps are intimately related and synchronous on the scale of a few million years. Copyright 1998 by the American Geophysical Union.

Palaeozoic deformation and magmatism in the northern area of the Anatolide block (Konya), witness of the Palaeotethys active margin

The study area. located north of Konva (Central Turkey), is composed of Silurian to Cretaceous metamorphosed rocks. The lower unit of the oldest formation (Silurian-Early Permian) is mostly made up of Silurian-Early Carboniferous metacarbonates. These rocks pass laterally and vertically to Devonian-Early Permian series having continental margin, shallow water and pelagic characteristics. They are intruded or juxtaposed to different kinds of metamagmatic rocks. which show MORB. continental arc and within plate characteristics. The Palaeozoic units are covered unconformably by Triassic-Cretaceous metasedimentary units. All these rocks are overthrusted by Mesozoic ophiolites. The Palaeozoic sequence can be seen as a northern Palaeotethys passive, then active margin. The northward subduction of the Palaeotethys ocean during the Carboniferous-Triassic times, induced the development of a magmatic arc and fore-arc sequence (Carboniferous-Permian). Before the Early Triassic (?Late Permian) time. the fore-arc sequence was uplifted above sea level and eroded. The Triassic sequences are regarded as marking the onset of back-arc opening and detachment of the Anatolian Konya block from the active Eurasian margin. Finally. a suture zone formed during the Carman between the Konya region and the Menderes-Tauride Cimmerian block due to the closing of Palaeotethvs. This geodynamic evolution can be correlated with the evolution of the Karaburun sequence in western Turkey.

Deformation and geochronology of syntectonic granitoids emplaced in the Major Gercino Shear Zone, southeastern South America

The Major Gercino Shear Zone is one of the NE-SW lineaments that separate the Neoproterozoic Dom Feliciano Belt, of Brazil and Uruguay, into two different domains: a northwestern supracrustal domain from a southeastern granitoid domain. The shear zone, striking NE, is composed of protomylonites to ultramylonites with mainly dextral kinematic indicators. In Santa Catarina State, southern Brazil, the shear zone is composed of two mylonite belts. The mylonites have mineral orientations produced under greenschist fades conditions at a high strain rate. Strong flattening and coaxial deformation indicate the transpressive character, while the role of pure shear is emphasized by the orientation of the mylonite belts in relation to the inferred stress field component. The quartz microstructures point out that different dynamic recrystallization regimes and crystal plasticity were the dominant mechanisms of deformation during the mylonitization process. Additionally, the fabrics suggest that the glide systems are activated for deformation conditions compatible with the metamorphism in the middle greenschist facies. Elongated granitoid intrusions belonging to two petrographically, geochemically and isotopically distinct rock associations occur between the two mylonite belts. The structures observed in the granites result from a deformation range from magmatic to solid-state conditions points to a continuum of magma straining during and just after its crystallization. Conventional U-Pb analysis of multi-crystal zircon fractions yielded essentially identical ages of 609 +/- 16 Ma and 614 +/- 2 Ma for the two granitic associations, and constrain the transpressive phase of the shear zone. K-Ar ages of biotites between 585 and 560 Ma record the slow cooling and uplift of the intrusions. Some K-Ar ages of micas in regional mylonites are similar, suggesting that thermo-tectonic activity was intense up to this time, probably related to the agglutination of the granite belt to the supracrustal belt NW of the MGSZ. (C) 2009 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

Contrasting magmatic signatures in the Rairakhol and Koraput alkaline complexes, Eastern Ghats belt, India

The relation between alkaline magmatism and tectonism has been a contentious issue, particularly for the Precambrian continental regions. Alkaline complexes at the southwestern margin of Eastern Ghats belt, India, have been interpreted as rift-valley magmatism. However, those complexes occurring in granulite ensemble in the interior segments of the Eastern Ghats belt could not possibly be related to the rift-system, assumed for the western margin of the Eastern Ghats belt. Koraput complex was emplaced in a pull-apart structure, dominated by magmatic fabrics and geochemically similar to a fractionated alkaline complex, compatible with an alkalibasalt series. Rairakhol complex, on the other hand, shows dominantly solid-state deformation fabrics and geochemically similar to a fractionated calc-alkaline suite. Isotopic data for the Koraput complex indicate ca. 917 Ma alkaline magmatism from a depleted mantle source and postcrystalline thermal overprint at ca. 745 Ma, also recorded from sheared metapelitic country rocks. The calc-alkaline magmatism of the Rairakhol complex occurred around 938 Ma, from an enriched mantle source, closely following Grenvillian granulite facies imprint in the charnockitic country rocks.

Magnetic fabric and zircon U-Pb geochronology of the Itaoca pluton: Implications for the Brasiliano deformation of the southern Ribeira Belt (SE Brazil)

The Itaoca pluton consists of porphyritic monzogranite that intruded the upper crust into low-grade metasedimentary rocks of the Apiai Dornain (Ribeira Belt). Anisotropy of magnetic susceptibility and zircon U-Pb (Shrimp) geochronology were combined to determine pluton emplacement mechanisms and its chronology relative to the collision structures of the Paranapiacaba (Brasiliano II) orogenic system. Magnetic susceptibility ranges between 4 and 38 x 10(-3) SI, and thermomagnetic measurements indicate multidomain magnetite is the main carrier of anisotropy. The pluton shows an ""onion-skin"" structure roughly elongated to the northeast with its hinge zone including kilometer-wide roof-pendants. Magnetic lineations are variable in orientation in consistency with the dominant oblate symmetry of the magnetic fabric. A distinct NE-trending point-maxima, however, indicates the mean lineation is parallel to the stretching direction of the transpressive deformation that affected the regional host rocks. Prismatic zircon from the monzogranite, both in the core and in the finely-zoned margins, yielded an age of 623 +/- 10 Ma. These results suggest the magmatic fabric recorded the earlier strain increments of the regional shear deformation. It may correspond to the transition from continental arc to collision tectonics of the southern Ribeira Belt. (C) 2008 Elsevier Ltd. All rights reserved.

Neoproterozoic tectonic and magmatic episodes in the NW sector of Borborema Province, NE Brazil, during assembly of Western Gondwana

The Medio Coreau domain of NE Brazil is located along the northwest margin of Borborema Province, the western branch of a Brasiliano/Pan-African collisional belt that formed during the assembly of Western Gondwana. The early Paleoproterozoic basement of the Medio Coreau domain is composed of migmatitic gneisses and juvenile granulites, overlain by late Paleoproterozoic and Neoproterozoic rocks intruded by syn- to post-tectonic Brasiliano granitoids. According to integrated structural and geochronological data (U-Pb zircon and monazite ages), the Neoproterozoic tectonic evolution of the Medio Coreau is characterized by low-angle thrusting and transcurrent deformation. U-Pb geochronological data from plutons intruded during this compressional regime indicate the collisional evolution began at approximately 622 Ma and continued until about 591 Ma. The continuation of convergence until approximately 560 Ma resulted in the formation of NE-SW and E-W shear zones within the Borborema Province and adjoining West African provinces. The final stage of the ductile tectonism was characterized by uplift and high-angle fault generation between approximately 560 and 545 Ma. The last tectonic event was an extensional phase, resulting in the formation of the Jaibaras graben and intrusion of post-orogenic granites at around 532 Ma. (c) 2007 Elsevier Ltd. All rights reserved.

Geology, petrology, U-Pb (shrimp) geochronology of the Morrinhos granite -Paraguá terrane, SW Amazonian craton: implications for the magmatic evolution of the San Ignácio orogeny

O Granito Morrinhos é um corpo batolítico levemente alongado segundo a direção NNW, de aproximadamente 1.140 km², localizado no município de Vila Bela da Santíssima Trindade, estado de Mato Grosso. Situa-se no Terreno Paraguá, Província Rondoniana-San Ignácio, na porção SW do Cráton Amazônico. Essa intrusão exibe uma variação composicional entre tonalito a monzogranito, textura inequigranular média a grossa, localmente, porfirítica, tendo biotita como máfico predominante em uma das fácies e hornblenda na outra, ambas metamorfizadas na fácies xisto verde. As rochas estudadas caracterizam uma sequência intermediária a ácida formada por um magmatismo subalcali no, do tipo álcali-cálcico, metaluminoso a levemente peraluminoso evoluído por meio de mecanismos de cristalização fracionada. Dados estruturais exibem registros de duas fases deformacionais, representadas pela foliação penetrativa (S₁) e dobras abertas (D₂) ambas, provavelmente, relacionadas à Orogenia San Ignácio. A investigação geocronológica (U-Pb SHRIMP) e geoquímica isotópica (Sm-Nd) dessas rochas indicaram, respectivamente, idade de cristalização 1.350 ± 12 Ma, T_DM em torno de 1,77 Ga e valor negativo para εNd(_1,35) de -2,57, sugerindo uma geração relacionada com processo de fusão parcial de uma crosta continental paleoproterozoica (estateriana). Os resultados aqui obtidos indicam que o Granito Morrinhos foi gerado em arco magmático continental, em estágio tardi a pós-orogênico, da Orogenia San Ignácio e permite reconhecê-lo como pertencente à Suíte Intrusiva Pensamiento.

Cerro Porã Batholith: post-orogenic A-type granite from the Amoguijá Magmatic Arc - Rio ApaTerrane - South of the Amazonian Craton

O Batólito Cerro Porã é um corpo de aproximadamente 30 por 4 km de extensão, localizado na região de Porto Murtinho, Mato Grosso do Sul. Situa-se nos domínios do Terreno Rio Apa, porção sul do Cráton Amazônico. Constitui-se pela Fácies sienogranítica rosa e Fácies monzogranítica cinza. A primeira é caracterizada por textura equi a, essencialmente, inequigranular xenomórfica e pela presença constante de intercrescimentos gráfico e granofíric; constitui-se por feldspatos alcalinos, quartzo e plagioclásio, tendo biotita como único máfico primário. A Fácies monzogranítica cinza apresenta textura porfirítica, com uma matriz de granulação fina gráfica a granofírica e consiste de quartzo, plagioclásio, feldspatos alcalinos e agregados máficos (biotita e anfibólio). Ambas foram metamorfizadas na fácies xisto verde e a Fácies sienogranítica rosa mostra-se milonitizada quando em zonas de cisalhamento. Foi identificado um evento deformacional dúctil-rúptil originado em regime compressivo, responsável pela geração de xistosidade e lineação de estiramento mineral. A Zona de Cisalhamento Esperança relaciona-se a esta fase e reflete a história cinemática convergente, reversa a de cavalgamento, com transporte de topo para NWW. Quimicamente, esses litotipos classificam-se como granitoides do tipo A2 da série alcalina potássica saturada em sílica. Determinação geocronológica obtida pelo método U-Pb (SHRIMP) em zircão, forneceu idade de 1749 ±45 Ma para sua cristalização. Do ponto vista geotectônico, admite-se que o Granito Cerro Porã corresponda a um magmatismo associado a um arco vulcânico desenvolvido no Estateriano e que sua colocação se deu no estágio tardi a pós-orogênico.

Numerical simulations of magma chamber dynamics at Campi Flegrei, and associated seismicity, deformation and gravity changes

Understanding the complex relationships between quantities measured by volcanic monitoring network and shallow magma processes is a crucial headway for the comprehension of volcanic processes and a more realistic evaluation of the associated hazard. This question is very relevant at Campi Flegrei, a volcanic quiescent caldera immediately north-west of Napoli (Italy). The system activity shows a high fumarole release and periodic ground slow movement (bradyseism) with high seismicity. This activity, with the high people density and the presence of military and industrial buildings, makes Campi Flegrei one of the areas with higher volcanic hazard in the world. In such a context my thesis has been focused on magma dynamics due to the refilling of shallow magma chambers, and on the geophysical signals detectable by seismic, deformative and gravimetric monitoring networks that are associated with this phenomenologies. Indeed, the refilling of magma chambers is a process frequently occurring just before a volcanic eruption; therefore, the faculty of identifying this dynamics by means of recorded signal analysis is important to evaluate the short term volcanic hazard. The space-time evolution of dynamics due to injection of new magma in the magma chamber has been studied performing numerical simulations with, and implementing additional features in, the code GALES (Longo et al., 2006), recently developed and still on the upgrade at the Istituto Nazionale di Geofisica e Vulcanologia in Pisa (Italy). GALES is a finite element code based on a physico-mathematical two dimensional, transient model able to treat fluids as multiphase homogeneous mixtures, compressible to incompressible. The fundamental equations of mass, momentum and energy balance are discretised both in time and space using the Galerkin Least-Squares and discontinuity-capturing stabilisation technique. The physical properties of the mixture are computed as a function of local conditions of magma composition, pressure and temperature.The model features enable to study a broad range of phenomenologies characterizing pre and sin-eruptive magma dynamics in a wide domain from the volcanic crater to deep magma feeding zones. The study of displacement field associated with the simulated fluid dynamics has been carried out with a numerical code developed by the Geophysical group at the University College Dublin (O’Brien and Bean, 2004b), with whom we started a very profitable collaboration. In this code, the seismic wave propagation in heterogeneous media with free surface (e.g. the Earth’s surface) is simulated using a discrete elastic lattice where particle interactions are controlled by the Hooke’s law. This method allows to consider medium heterogeneities and complex topography. The initial and boundary conditions for the simulations have been defined within a coordinate project (INGV-DPC 2004-06 V3_2 “Research on active volcanoes, precursors, scenarios, hazard and risk - Campi Flegrei”), to which this thesis contributes, and many researchers experienced on Campi Flegrei in volcanological, seismic, petrological, geochemical fields, etc. collaborate. Numerical simulations of magma and rock dynamis have been coupled as described in the thesis. The first part of the thesis consists of a parametric study aimed at understanding the eect of the presence in magma of carbon dioxide in magma in the convection dynamics. Indeed, the presence of this volatile was relevant in many Campi Flegrei eruptions, including some eruptions commonly considered as reference for a future activity of this volcano. A set of simulations considering an elliptical magma chamber, compositionally uniform, refilled from below by a magma with volatile content equal or dierent from that of the resident magma has been performed. To do this, a multicomponent non-ideal magma saturation model (Papale et al., 2006) that considers the simultaneous presence of CO2 and H2O, has been implemented in GALES. Results show that the presence of CO2 in the incoming magma increases its buoyancy force promoting convection ad mixing. The simulated dynamics produce pressure transients with frequency and amplitude in the sensitivity range of modern geophysical monitoring networks such as the one installed at Campi Flegrei . In the second part, simulations more related with the Campi Flegrei volcanic system have been performed. The simulated system has been defined on the basis of conditions consistent with the bulk of knowledge of Campi Flegrei and in particular of the Agnano-Monte Spina eruption (4100 B.P.), commonly considered as reference for a future high intensity eruption in this area. The magmatic system has been modelled as a long dyke refilling a small shallow magma chamber; magmas with trachytic and phonolitic composition and variable volatile content of H2O and CO2 have been considered. The simulations have been carried out changing the condition of magma injection, the system configuration (magma chamber geometry, dyke size) and the resident and refilling magma composition and volatile content, in order to study the influence of these factors on the simulated dynamics. Simulation results allow to follow each step of the gas-rich magma ascent in the denser magma, highlighting the details of magma convection and mixing. In particular, the presence of more CO2 in the deep magma results in more ecient and faster dynamics. Through this simulations the variation of the gravimetric field has been determined. Afterward, the space-time distribution of stress resulting from numerical simulations have been used as boundary conditions for the simulations of the displacement field imposed by the magmatic dynamics on rocks. The properties of the simulated domain (rock density, P and S wave velocities) have been based on data from literature on active and passive tomographic experiments, obtained through a collaboration with A. Zollo at the Dept. of Physics of the Federici II Univeristy in Napoli. The elasto-dynamics simulations allow to determine the variations of the space-time distribution of deformation and the seismic signal associated with the studied magmatic dynamics. In particular, results show that these dynamics induce deformations similar to those measured at Campi Flegrei and seismic signals with energies concentrated on the typical frequency bands observed in volcanic areas. The present work shows that an approach based on the solution of equations describing the physics of processes within a magmatic fluid and the surrounding rock system is able to recognise and describe the relationships between geophysical signals detectable on the surface and deep magma dynamics. Therefore, the results suggest that the combined study of geophysical data and informations from numerical simulations can allow in a near future a more ecient evaluation of the short term volcanic hazard.