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.

Magmatic–hydrothermal molybdenum isotope fractionation and its relevance to the igneous crustal signature

We analysed the Mo isotope composition of a comprehensive series of molybdenite samples from the porphyry- type Questa deposit (NM, USA), as well as one rhyolite and one granite sample, directly associated with the Mo mineralization. The δ98Mo of the molybdenites ranges between −0.48‰ and +0.40‰, with a median at −0.05‰. The median Mo isotope composition increases from early magmatic (−0.29‰) to hydrothermal (−0.05‰) breccia mineralization (median bulk breccia = −0.17‰) to late stockwork veining (+0.22‰). Moreover, variations of up to 0.34‰ are found between different molybdenite crystals within an individual hand specimen. The rhyolite sample with 0.12 μg g−1 Mo has δ98Mo = −0.57‰ and is lighter than all molybde- nites from the Questa deposit, interpreted to represent the igneous leftover after aqueous ore fluid exsolution. We recognize three Mo isotope fractionation processes that occur between about 700 and 350 °C, affecting the Mo iso- tope composition of magmatic–hydrothermal molybdenites. Δ1Mo: Minerals preferentially incorporate light Mo isotopes during progressive fractional crystallization in subvolcanic magma reservoirs, leaving behind a melt enriched in heavy Mo isotopes. Δ2Mo: Magmatic–hydrothermal fluids preferentially incorporate heavy Mo iso- topes upon fluid exsolution. Δ3Mo: Light Mo isotopes get preferentially incorporated in molybdenite during crys- tallization from an aqueous fluid, leaving behind a hydrothermal fluid that gets heavier with progressive molybdenite crystallization. The sum of all three fractionation processes produces molybdenites that record heavier δ98Mo compositions than their source magmas. This implies that the mean δ98Mo of molybdenites published so far (~0.4‰) likely represents a maximum value for the Mo isotope composition of Phanerozoic igneous upper crust.

Weathering geochemistry and Sr-Nd fingerprints of equatorial upper Nile and Congo muds

This study investigates processes of sediment generation in equatorial central Africa. An original, complete and integrated mineralogical-geochemical database on silt-sized sediments derived from different parent rocks (basalt, granite, gneiss, metapsammite, sandstone) along the East African Rift from 5°S in Tanzania to 5°N in Sudan is presented and used to assess the incidence of diverse factors controlling sediment composition (source-rock lithology, geomorphology, hydraulic sorting, grain size, recycling), with particular emphasis on chemical weathering.

Solute transport in crystalline rocks at Äspö — I: Geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted A¨ spo¨ Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175]. The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments. While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours–days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption Kds are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.

Solute transport in crystalline rocks at Äspö — II: Blind predictions, inverse modelling and lessons learnt from test STT1

Based on the results from detailed structural and petrological characterisation and on up-scaled laboratory values for sorption and diffusion, blind predictions were made for the STT1 dipole tracer test performed in the Swedish A¨ spo¨ Hard Rock Laboratory. The tracers used were nonsorbing, such as uranine and tritiated water, weakly sorbing 22Na+, 85Sr2 +, 47Ca2 +and more strongly sorbing 86Rb+, 133Ba2 +, 137Cs+. Our model consists of two parts: (1) a flow part based on a 2D-streamtube formalism accounting for the natural background flow field and with an underlying homogeneous and isotropic transmissivity field and (2) a transport part in terms of the dual porosity medium approach which is linked to the flow part by the flow porosity. The calibration of the model was done using the data from one single uranine breakthrough (PDT3). The study clearly showed that matrix diffusion into a highly porous material, fault gouge, had to be included in our model evidenced by the characteristic shape of the breakthrough curve and in line with geological observations. After the disclosure of the measurements, it turned out that, in spite of the simplicity of our model, the prediction for the nonsorbing and weakly sorbing tracers was fairly good. The blind prediction for the more strongly sorbing tracers was in general less accurate. The reason for the good predictions is deemed to be the result of the choice of a model structure strongly based on geological observation. The breakthrough curves were inversely modelled to determine in situ values for the transport parameters and to draw consequences on the model structure applied. For good fits, only one additional fracture family in contact with cataclasite had to be taken into account, but no new transport mechanisms had to be invoked. The in situ values for the effective diffusion coefficient for fault gouge are a factor of 2–15 larger than the laboratory data. For cataclasite, both data sets have values comparable to laboratory data. The extracted Kd values for the weakly sorbing tracers are larger than Swedish laboratory data by a factor of 25–60, but agree within a factor of 3–5 for the more strongly sorbing nuclides. The reason for the inconsistency concerning Kds is the use of fresh granite in the laboratory studies, whereas tracers in the field experiments interact only with fracture fault gouge and to a lesser extent with cataclasite both being mineralogically very different (e.g. clay-bearing) from the intact wall rock.

Field-scale water transport in unsaturated crystalline rock

Safe disposal of toxic wastes in geologic formations requires minimal water and gas movement in the vicinity of storage areas, Ventilation of repository tunnels or caverns built in solid rock can desaturate the near field up to a distance of meters from the rock surface, even when the surrounding geological formation is saturated and under hydrostatic pressures. A tunnel segment at the Grimsel test site located in the Aare granite of the Bernese Alps (central Switzerland) has been subjected to a resaturation and, subsequently, to a controlled desaturation, Using thermocouple psychrometers (TP) and time domain reflectometry (TDR), the water potentials psi and water contents theta were measured within the unsaturated granodiorite matrix near the tunnel wall at depths between 0 and 160 cm. During the resaturation the water potentials in the first 30 cm from the rock surface changed within weeks from values of less than -1.5 MPa to near saturation. They returned to the negative initial values during desaturation, The dynamics of this saturation-desaturation regime could be monitored very sensitively using the thermocouple psychrometers, The TDR measurements indicated that water contents changed dose to the surface, but at deeper installation depths the observed changes were within the experimental noise. The field-measured data of the desaturation cycle were used to test the predictive capabilities of the hydraulic parameter functions that were derived from the water retention characteristics psi(theta) determined in the laboratory. A depth-invariant saturated hydraulic conductivity k(s) = 3.0 x 10(-11) m s(-1) was estimated from the psi(t) data at all measurement depths, using the one-dimensional, unsaturated water flow and transport model HYDRUS Vogel er al., 1996, For individual measurement depths, the estimated k(s) varied between 9.8 x 10(-12) and 6.1 x 10(-11) The fitted k(s) values fell within the range of previously estimated k(s) for this location and led to a satisfactory description of the data, even though the model did not include transport of water vapor.

Exhumation rates in the Archean from pressure-ime paths: Examplefrom the Skjoldungen Orogen (SE Greenland)

The discussions on the orogenic evolution during Earth's history converge to the question of a different thermal structure in the Archean compared to the Phanerozoic and the applicability of the plate tectonic paradigm. However, geothermal structures are transient in orogens and are difficult to translate into large-scale tectonics and exhumation rates. Therefore, we propose depth–time data in the Archean Skjoldungen Orogen (SE Greenland, North Atlantic Craton) that allow for reconstruction of an exhumation rate independent of geothermal gradients. The resulting exhumation rate of ca. 0.4 km/Ma is similar to exhumation rates during erosion-controlled processes in modern orogens. These exhumation rates can only be established by erosion time constants similar to modern orogens. The occurrence of erosion-controlled exhumation is best explained by a stiff foreland promoting localized deformation in the orogen. Therefore, a switch from magmatic-dominated processes to localized deformation is proposed in the Skjoldungen Orogen area. This is supported by a change in magma composition and volume, from widespread granodiorite to localized alkaline intrusions. In addition, the involved metasedimentary rocks include detrital zircons of the only 50 Ma older foreland, which also correspond to erosion and tectonics as in modern orogens, i.e. flysh-type sediments. Relatively fast exhumation rates and the structural-magmatic evolution of the Neoarchean Skjoldungen Orogen thus indicate modern-style tectonic processes where stiff Mesoarchean continental crust forms a foreland to a collisional orogen instead of typical accretionary tectonics of weak island arc-like terranes in granite-greenstone terranes.

Bedload composition, transport and modification in rivers of Westland, New Zealand, with implications for the distribution of alluvial pounamu (jade)

River bedload surveyed at 50 sites in Westland is dominated by Alpine Schist or Torlesse Greywacke from the Alpine Fault hanging wall, with subordinate Pounamu Ultramafics or footwall-derived Western Province rocks. Tumbling experiments found ultramafics to have the lowest attrition rates, compared with greywacke sandstone and granite (which abrade to produce silt to medium-sand), or incompetent schist (which fragments). Arahura has greater total concentrations (103–105 t/km2) and proportions (5–40%) of ultramafic bedload compared with Hokitika and Taramakau catchments (101–104 t/km2, mostly <10%), matching relative areas of mapped Pounamu Ultramafic bedrock, but enriched relative to absolute areal proportions. Western Province rocks downthrown by the Alpine Fault are under-represented in the bedload. Enriched concentrations of ultramafic bedload decrease rapidly with distance downstream from source rock outcrops, changing near prominent ice-limit moraines. Bedload evolution with transport involves both downstream fining and dilution from tributaries, in a sediment supply regime more strongly influenced by tectonics and the imprint of past glaciation. Treasured New Zealand pounamu (jade) is associated with ultramafic rocks. Chances of discovery vary between catchments, are increased near glacial moraines, and are highest near source-rock outcrops in remote mountain headwaters.

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.

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.