U-Pb SHRIMP and 40Ar/39Ar constraints on the timing of mineralization in the Paleoproterozoic Caxias orogenic gold deposit, São Luís cratonic fragment, Brazil

Caxias é um depósito de ouro orogênico do fragmento cratônico São Luís, que é correlacionável aos terrenos Riacianos do Cráton Oeste-Africano. O depósito se formou após o metamorfismo regional (estimado em 2100 ± 15 Ma) e está hospedado em zona de cisalhamento que secciona xistos do Grupo Aurizona (2240 ± 5 Ma) e o Microtonalito Caxias. O microtonalito foi aqui datado em 2009 ± 11 Ma, e representa um estágio magmático tardio na evolução do fragmento cratônico São Luís. Cristais de zircão com idades de 2139 ± 10 Ma foram herdados da fonte magmática ou são produto de contaminação durante a intrusão. A composição dos isótopos de chumbo sugere que granitoides de arco de ilhas de ca. 2160 Ma são a fonte provável para o Pb incorporado na pirita relacionada com o minério. Sericita hidrotermal mostra idade ⁴⁰Ar/³⁹Ar de 1990 ± 30 Ma, que, combinada com a idade de posicionamento do microtonalito hospedeiro, limita o evento mineralizador ao intervalo 2020-1960 Ma.

The mineralizing fluid in the Piaba gold deposit, São Luís cratonic fragment (NW-Maranhão, Brazil) based on fluid inclusion studies in quartz veins

O depósito aurífero de Piaba tornou-se a primeira mina em operação no fragmento cratônico São Luís, noroeste do Maranhão. Seu ambiente geológico compreende rochas metavulcanossedimentares do Grupo Aurizona e granitoides da Suíte Tromaí, entre outras unidades menores, formadas em ambiente de arcos de ilhas entre 2240 e 2150 Ma, juntamente com outras unidades menores. A mineralização é caracterizada por uma trama stockwork de veios e vênulas de quartzo com seus halos de alteração (clorita + muscovita + carbonato + pirita + calcopirita e ouro) hospedada em um granodiorito granofírico fino (Granófiro Piaba) e em rocha subvulcânica andesítica do Grupo Aurizona. O corpo mineralizado é espacialmente limitado à zona de cisalhamento ENE-WSW rúptil-dúctil (Falha Piaba). Estudos petrográficos, microtermométricos e por espectroscopia microRaman no quartzo definiram inclusões aquo-carbônicas bifásicas e trifásicas, produzidas por aprisionamento heterogêneo durante separação de fases, e fluidos aquosos tardios. A solução mineralizadora corresponde a um fluido aquo-carbônico composto por CO₂ (5 - 24 mol%, densidade de 0,96 - 0,99 g/cm³), H₂O (74 - 93 mol%), N₂ (< 1 mol%), CH₄ (<1mol%) e 5,5 % em peso NaCl equivalente. O minério depositou a 267 - 302ºC e 1,25 - 2,08 kbar, correspondendo a profundidades de 4 a 7 km, em consonância com o regime estrutural. A composição e o intervalo de P-T do fluido mineralizador, combinados com o caráter redutor (log ƒ_O2 -31,3 a -34,3) e a sulfetação da rocha hospedeira, sugerem que o ouro foi transportado como um complexo sulfetado. O minério foi depositado em consequência da separação de fase, redução da atividade de enxofre e da ƒO₂ pela interação fluido-rocha.

Spectral-Element and Adjoint 3D Full-Wave Tomography for the Lithosphere of Central Italy

The primary objective of this thesis is to obtain a better understanding of the 3D velocity structure of the lithosphere in central Italy. To this end, I adopted the Spectral-Element Method to perform accurate numerical simulations of the complex wavefields generated by the 2009 Mw 6.3 L’Aquila event and by its foreshocks and aftershocks together with some additional events within our target region. For the mainshock, the source was represented by a finite fault and different models for central Italy, both 1D and 3D, were tested. Surface topography, attenuation and Moho discontinuity were also accounted for. Three-component synthetic waveforms were compared to the corresponding recorded data. The results of these analyses show that 3D models, including all the known structural heterogeneities in the region, are essential to accurately reproduce waveform propagation. They allow to capture features of the seismograms, mainly related to topography or to low wavespeed areas, and, combined with a finite fault model, result into a favorable match between data and synthetics for frequencies up to ~0.5 Hz. We also obtained peak ground velocity maps, that provide valuable information for seismic hazard assessment. The remaining differences between data and synthetics led us to take advantage of SEM combined with an adjoint method to iteratively improve the available 3D structure model for central Italy. A total of 63 events and 52 stations in the region were considered. We performed five iterations of the tomographic inversion, by calculating the misfit function gradient - necessary for the model update - from adjoint sensitivity kernels, constructed using only two simulations for each event. Our last updated model features a reduced traveltime misfit function and improved agreement between data and synthetics, although further iterations, as well as refined source solutions, are necessary to obtain a new reference 3D model for central Italy tomography.

Seismic tomographic full-waveform inversion for the Vrancea sinking lithosphere structure using the adjoint method.

The Vrancea region, at the south-eastern bend of the Carpathian Mountains in Romania, represents one of the most puzzling seismically active zones of Europe. Beside some shallow seismicity spread across the whole Romanian territory, Vrancea is the place of an intense seismicity with the presence of a cluster of intermediate-depth foci placed in a narrow nearly vertical volume. Although large-scale mantle seismic tomographic studies have revealed the presence of a narrow, almost vertical, high-velocity body in the upper mantle, the nature and the geodynamic of this deep intra-continental seismicity is still questioned. High-resolution seismic tomography could help to reveal more details in the subcrustal structure of Vrancea. Recent developments in computational seismology as well as the availability of parallel computing now allow to potentially retrieve more information out of seismic waveforms and to reach such high-resolution models. This study was aimed to evaluate the application of a full waveform inversion tomography at regional scale for the Vrancea lithosphere using data from the 1999 six months temporary local network CALIXTO. Starting from a detailed 3D Vp, Vs and density model, built on classical travel-time tomography together with gravity data, I evaluated the improvements obtained with the full waveform inversion approach. The latter proved to be highly problem dependent and highly computational expensive. The model retrieved after the first two iterations does not show large variations with respect to the initial model but remains in agreement with previous tomographic models. It presents a well-defined downgoing slab shape high velocity anomaly, composed of a N-S horizontal anomaly in the depths between 40 and 70km linked to a nearly vertical NE-SW anomaly from 70 to 180km.

Geochemistry of the gabbroic lithosphere in IODP Holes 304-U1309D and 305-U1309D

IODP Hole U1309D (Atlantis Massif, Mid-Atlantic Ridge 30°N) is the second deepest hole drilled into slow spread gabbroic lithosphere. It comprises 5.4% of olivine-rich troctolites (~ > 70% olivine), possibly the most primitive gabbroic rocks ever drilled at mid-ocean ridges. We present the result of an in situ trace element study carried out on a series of olivine-rich troctolites, and neighbouring troctolites and gabbros, from olivine-rich intervals in Hole U1309D. Olivine-rich troctolites display poikilitic textures; coarse-grained subhedral to medium-grained rounded olivine crystals are included into large undeformed clinopyroxene and plagioclase poikiloblasts. In contrast, gabbros and troctolites have irregularly seriate textures, with highly variable grain sizes, and locally poikilitic clinopyroxene oikocrysts in troctolites. Clinopyroxene is high Mg# augite (Mg# 87 in olivine-rich troctolites to 82 in gabbros), and plagioclase has anorthite contents ranging from 77 in olivine-rich troctolites to 68 in gabbros. Olivine has high forsterite contents (82-88 in olivine-rich troctolites, to 78-83 in gabbros) and is in Mg-Fe equilibrium with clinopyroxene. Clinopyroxene cores and plagioclase are depleted in trace elements (e.g., Ybcpx ~ 5-11 * Chondrite), they are in equilibrium with the same MORB-type melt in all studied rock-types. These compositions are not consistent with the progressively more trace element enriched (evolved) compositions expected from olivine rich primitive products to gabbros in a MORB cumulate sequence. They indicate that clinopyroxene and plagioclase crystallized concurrently, after melts having the same trace element composition, consistent with crystallization in an open system with a buffered magma composition. The slight trace element enrichments and lower Cr contents observed in clinopyroxene rims and interstitial grains results from crystallization of late-stage differentiated melts, probably indicating the closure of the magmatic system. In contrast to clinopyroxene and plagioclase, olivine is not in equilibrium with MORB, but with a highly fractionated depleted melt, similar to that in equilibrium with refractory oceanic peridotites, thus possibly indicating a mantle origin. In addition, textural relationships suggest that olivine was in part assimilated by the basaltic melts after which clinopyroxene and plagioclase crystallized (impregnation). These observations suggest a complex crystallization history in an open system involving impregnation by MORB-type melt(s) of an olivine-rich rock or mush. The documented magmatic processes suggest that olivine-rich troctolites were formed in a zone with large magmatic transfer and accumulation, similar to the mantle-crust transition zone documented in ophiolites and at fast spreading ridges.

Strength of the Iberian intraplate lithosphere: Cenozoic deformations and seismicity

Previous studies about the strength of the lithosphere in the center of Iberia fail to resolve the depth of earthquakes because of the rheological uncertainties. Therefore, new contributions are considered (the crustal structure from a density model) and several parameters (tectonic regime, mantle rheology, strain rate) are checked in this paper to properly examine the role of lithospheric strength in the intraplate seismicity and the Cenozoic evolution. The strength distribution with depth, the integrated strength, the effective elastic thickness and the seismogenic thickness have been calculated by a finite element modelling of the lithosphere across the Central System mountain range and the bordering Duero and Madrid sedimentary basins. Only a dry mantle under strike-slip/extension and a strain rate of 10-15 s-1, or under extension and 10-16 s-1, causes a strong lithosphere. The integrated strength and the elastic thickness are lower in the mountain chain than in the basins. This heterogeneity has been maintained since the Cenozoic and determine the mountain uplift and the biharmonic folding of the Iberian lithosphere during the Alpine deformations. The seismogenic thickness bounds the seismic activity in the upper–middle crust, and the decreasing crustal strength from the Duero Basin towards the Madrid Basin is related to a parallel increase in Plio–Quaternary deformations and seismicity. However, elasto–plastic modelling shows that current African–Eurasian convergence is resolved elastically or ductilely, which accounts for the low seismicity recorded in this region.