Characterization of granite matrix porosity and pore-space geometry by in situ and laboratory methods

Most available studies of interconnected matrix porosity of crystalline rocks are based on laboratory investigations; that is, work on samples that have undergone stress relaxation and were affected by drilling and sample preparation. The extrapolation of the results to in situ conditions is therefore associated with considerable uncertainty, and this was the motivation to conduct the ‘in situ Connected Porosity’ experiment at the Grimsel Test Site (Central Swiss Alps). An acrylic resin doped with fluorescent agents was used to impregnate the microporous granitic matrix in situ around an injection borehole, and samples were obtained by overcoring. The 3-D structure of the porespace, represented by microcracks, was studied by U-stage fluorescence microscopy. Petrophysical methods, including the determination of porosity, permeability and P -wave velocity, were also applied. Investigations were conducted both on samples that were impregnated in situ and on non-impregnated samples, so that natural features could be distinguished from artefacts. The investigated deformed granites display complex microcrack populations representing a polyphase deformation at varying conditions. The crack population is dominated by open cleavage cracks in mica and grain boundary cracks. The porosity of non-impregnated samples lies slightly above 1 per cent, which is 2–2.5 times higher than the in situ porosity obtained for impregnated samples. Measurements of seismic velocities (Vp ) on spherical rock samples as a function of confining pressure, spatial direction and water saturation for both non-impregnated and impregnated samples provide further constraints on the distinction between natural and induced crack types. The main conclusions are that (1) an interconnected network of microcracks exists in the whole granitic matrix, irrespective of the distance to ductile and brittle shear zones, and (2) conventional laboratory methods overestimate the matrix porosity. Calculations of contaminant transport through fractured media often rely on matrix diffusion as a retardation mechanism.

Dating prograde fluid pulses during subduction by in situ U–Pb and oxygen isotope analysis

Keywords High-pressure fluids · Whiteschists · U–Pb dating · Oxygen isotopes · Ion microprobe · Metasomatism Introduction The subduction of crustal material to mantle depths and its chemical modification during burial and exhumation contribute to element recycling in the mantle and the formation of new crust through arc magmatism. Crustal rocks that Abstract The Dora-Maira whiteschists derive from metasomatically altered granites that experienced ultrahighpressure metamorphism at ~750 °C and 40 kbar during the Alpine orogeny. In order to investigate the P–T–time– fluid evolution of the whiteschists, we obtained U–Pb ages from zircon and monazite and combined those with trace element composition and oxygen isotopes of the accessory minerals and coexisting garnet. Zircon cores are the only remnants of the granitic protolith and still preserve a Permian age, magmatic trace element compositions and δ18O of ~10 ‰. Thermodynamic modelling of Si-rich and Si-poor whiteschist compositions shows that there are two main fluid pulses during prograde subduction between 20 and 40 kbar. In Si-poor samples, the breakdown of chlorite to garnet + fluid occurs at ~22 kbar. A first zircon rim directly overgrowing the cores has inclusions of prograde phlogopite and HREE-enriched patterns indicating zircon growth at the onset of garnet formation. A second main fluid pulse is documented close to peak metamorphic conditions in both Si-rich and Si-poor whiteschist when talc + kyanite react to garnet + coesite + fluid. A second metamorphic overgrowth on zircon with HREE depletion was observed in the Si-poor whiteschists, whereas a single metamorphic overgrowth capturing phengite and talc inclusions was observed in the Si-rich whiteschists. Garnet rims, zircon rims and monazite are in chemical and isotopic equilibrium for oxygen, demonstrating that they all formed at peak metamorphism at 35 Ma as constrained by the age of monazite (34.7 ± 0.4 Ma) and zircon rims (35.1 ± 0.8 Ma). The prograde zircon rim in Si-poor whiteschists has an age that is within error indistinguishable from the age of peak metamorphic conditions, consistent with a minimum rate of subduction of 2 cm/year for the Dora-Maira unit. Oxygen isotope values for zircon rims, monazite and garnet are equal within error at 6.4 ± 0.4 ‰, which is in line with closed-system equilibrium fractionation during prograde to peak temperatures. The resulting equilibrium Δ18Ozircon-monazite at 700 ± 20 °C is 0.1 ± 0.7 ‰. The in situ oxygen isotope data argue against an externally derived input of fluids into the whiteschists. Instead, fluidassisted zircon and monazite recrystallisation can be linked to internal dehydration reactions during prograde subduction. We propose that the major metasomatic event affecting the granite protolith was related to hydrothermal seafloor alteration post-dating Jurassic rifting, well before the onset of Alpine subduction.

Tectonic evolution of the Proterozoic Itremo Group metasediments in central Madagascar

The major geologic units of the Itremo region in central Madagascar include: (1) upper amphibolite to granulite facies (higher grade) Precambrian rocks, mainly para- and orthogneisses, and migmatites; (2) the newly defined Itremo Nappes, a fold-and-thrust belt containing the Proterozoic Itremo Group sediments, metamorphosed at greenschist to lower amphibolite facies (lower grade) conditions: (3) Middle Neoproterozoic and Late Neoproterozoic-Cambrian intrusives. The stratigraphic succession of the Itremo Group in the eastern part of the Itremo region is, from bottom to top: quartzites, metapelites, metacarbonates and metapelites overlain by metacarbonates. During D1 the Itremo Group sediments were detached from their continental substratum, deformed into a fold-and-thrust nappe (Itremo Nappes), and transported on top of higher grade rocks that are intruded by Middle Neoproterozoic (c. 797–780 Ma) granites and gabbros. A second phase of deformation shortening (D2) affected both the Itremo Sedimentary Nappes and structurally underlying higher-grade rocksunits, and formed large-scale N-S-trending F2 folds. S1 axial plane foliations in Itremo Group sediments are truncated by Late Neoproterozoic-Cambrian granites (c. 570–540 Ma). The age of the formation of the Itremo Nappes is not well constrained: they formed in Neoproterozoic times between 780 and 570 Ma.

Age constraints on the tectonic evolution of the Itremo region in Central Madagascar

The Itremo region in Central Madagascar comprises a deformed metasedimentary sequence (Itremo Group) that has undergone greenschist to lower amphibolite facies metamorphism. During a first phase of deformation (D1) Itremo Group sediments were deformed into a fold-and-thrust belt and transported toward the E to NE on top of migmatitic gneisses rocks of Anatananarivo block. A second phase of deformation (D2) affected both the fold-and-thrust belt and structurally underlying units, and formed large-scale N-S trending folds with steeply dipping axial planes. A Late Neoproterozoic Th–U–Pb XRF monazite age (565±17 Ma) dates the emplacement of a granite that truncates first-phase structures in the Itremo Group, and indicates that the fold-and-thrust belt formed prior to ≈565 Ma. Th–U–Pb electron microprobe dating was applied to elongated monazites that lie within the first-phase foliation of Itremo Group metapelites. The detrital cores of zoned monazites reveal two distinct age populations at ∼2000 and 1700 Ma, the latter age giving a maximum depositional age for the Itremo Group. Statistical analysis of ages determined from the rims of zoned monazites and from unzoned monazites indicates three Late Proterozoic–Early Paleozoic monazite growth events at about 565–540, 500 and 430 Ma. The oldest age population is contemporaneous within error, with the intrusion of the dated granite. The two younger age populations are found both in the Th–U–Pb and Ar–Ar data; together with the perturbation of the Rb–Sr system we interpret both ages as due to alteration related to fluid circulation events, possibly connected to the emplacement of pegmatite fields in Central Madagascar. Syn-D1 tectonic growth of contact metamorphism minerals such as andalusite has been observed locally in metapelites along the margin of Middle Neoproterozoic (≈800 Ma) granites, suggesting that D1 in the Itremo Group is contemporaneous with the intrusion of granites at ≈800 Ma. The N-S trending D2 folds are associated with ≈E-W shortening during the final assembly of Gondwana in Late Neoproterozoic–Early Cambrian times.