Effects of drilling and stress release on hydraulic properties and porewater chemistry of crystalline rocks

The experimental verification of matrix diffusion in crystalline rocks largely relies on indirect methods performed in the laboratory. Such methods are prone to perturbations of the rock samples by collection and preparation and therefore the laboratory-derived transport properties and fluid composition might not represent in situ conditions. We investigated the effects induced by the drilling process and natural rock stress release by mass balance considerations and sensitivity analysis of analytical out-diffusion data obtained from originally saturated, large-sized drillcore material from two locations drilled using traced drilling fluid. For in situ stress-released drillcores of quartz-monzodiorite composition from the Aspo HRL, Sweden, tracer mass balance considerations and 1D and 2D diffusion modelling consistently indicated a contamination of <1% of the original pore water. This chemically disturbed zone extends to a maximum of 0.1 mm into the drillcore (61.8 mm x 180.1 mm) corresponding to about 0.66% of the total pore volume (0.77 vol.%). In contrast, the combined effects of stress release and the drilling process, which have influenced granodioritic drillcore material from 560 m below surface at Forsmark. Sweden, resulted in a maximum contamination of the derived porewater Cl(-) concentration of about 8%. The mechanically disturbed zone with modified diffusion properties covers the outermost similar to 6 mm of the drillcore (50 mm x 189 mm), whereas the chemically disturbed zone extends to a maximum of 0.3 mm based on mass balance considerations, and to 0.15 mm to 0.2 mm into the drillcore based on fitting the observed tracer data. This corresponds to a maximum of 2.4% of the total pore volume (0.62 vol.%) being affected by the drilling-fluid contamination. The proportion of rock volume affected initially by drilling fluid or subsequently with experiment water during the laboratory diffusion and re-saturation experiments depends on the size of the drillcore material and will become larger the smaller the sample used for the experiment. The results are further in support of matrix diffusion taking place in the undisturbed matrix of crystalline rocks at least in the cm range.

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.

Matrix Porewater In Crystalline Rocks: Extraction and Analysis

The chemical and isotopic characterization of porewater residing in the inter- and intragranular pore space of the low-permeability rock matrix is an important component with respect to the site characterization and safety assessment of potential host rocks for a radioactive waste disposal. The chemical and isotopic composition of porewater in such low permeability rocks has to be derived by indirect extraction techniques applied to naturally saturated rock material. In most of such indirect extraction techniques – especially in case of rocks of a porosity below about 2 vol.% – the original porewater concentrations are diluted and need to be back-calculated to in-situ concentrations. This requires a well-defined value for the connected porosity – accessible to different solutes under in-situ conditions. The derivation of such porosity values, as well as solute concentrations, is subject to various perturbations during drilling, core sampling, storage and experiments in the laboratory. The present study aims to demonstrate the feasibility of a variety of these techniques to charac-terize porewater and solute transport in crystalline rocks. The methods, which have been de-veloped during multiple porewater studies in crystalline environments, were applied on four core samples from the deep borehole DH-GAP04, drilled in the Kangerlussuaq area, Southwest Greenland, as part of the joint NWMO–Posiva–SKB Greenland Analogue Project (GAP). Potential artefacts that can influence the estimation of in situ porewater chemistry and isotopes, as well as their controls, are described in detail in this report, using specific examples from borehole DH-GAP04

Reconstruction of in-situ porosity and porewater compositions of low-permeability crystalline rocks: Magnitude of artefacts induced by drilling and sample recovery

Geological site characterisation programmes typically rely on drill cores for direct information on subsurface rocks. However, porosity, transport properties and porewater composition measured on drill cores can deviate from in-situ values due to two main artefacts caused by drilling and sample recovery: (1) mechanical disruption that increases porosity and (2) contamination of the porewater by drilling fluid. We investigated the effect and magnitude of these perturbations on large drill core samples (12–20 cm long, 5 cmdiameter) of high-grade, granitic gneisses obtained from 350 to 600 m depth in a borehole on Olkiluoto Island (SW Finland). The drilling fluid was traced with sodium–iodide. By combining out-diffusion experiments, gravimetry, UV-microscopy and iodide mass balance calculations, we successfully quantified the magnitudes of the artefacts: 2–6% increase in porosity relative to the bulk connected porosity and 0.9 to 8.9 vol.% contamination by drilling fluid. The spatial distribution of the drilling-induced perturbations was revealed by numerical simulations of 2D diffusion matched to the experimental data. This showed that the rims of the samples have a mechanically disrupted zone 0.04 to 0.22 cm wide, characterised by faster transport properties compared to the undisturbed centre (1.8 to 7.7 times higher pore diffusion coefficient). Chemical contamination was shown to affect an even wider zone in all samples, ranging from 0.15 to 0.60 cm, inwhich iodide enrichmentwas up to 180 mg/kgwater, compared to 0.5 mg/kgwater in the uncontaminated centre. For all samples in the present case study, it turned out that the magnitude of the artefacts caused by drilling and sample recovery is so small that no correction is required for their effects. Therefore, the standard laboratory measurements of porosity, transport properties and porewater composition can be taken as valid in-situ estimates. However, it is clear that the magnitudes strongly depend on site- and drilling-specific factors and therefore our results cannot be transferred simply to other locations. We recommend the approach presented in this study as a route to obtain reliable values in future drilling campaigns aimed at characterising in-situ bedrock properties.

Structural, metamorphic and geochronological relations between the Zanskar Shear Zone and the Miyar Shear Zone (NW Indian Himalaya): Evidence for two distinct tectonic structures and implications for the evolution of the High Himalayan Crystalline of Zanskar

The geologic structures and metamorphic zonation of the northwestern Indian Himalaya contrast significantly with those in the central and eastern parts of the range, where the high-grade metamorphic rocks of the High Himalayan Crystalline (HHC) thrust southward over the weakly metamorphosed sediments of the Lesser Himalaya along the Main Central Thrust (MCT). Indeed, the hanging wall of the MCT in the NW Himalaya mainly consists of the greenschist facies metasediments of the Chamba zone, whereas HHC high-grade rocks are exposed more internally in the range as a large-scale dome called the Gianbul dome. This Gianbul dome is bounded by two oppositely directed shear zones, the NE-dipping Zanskar Shear Zone (ZSZ) on the northern flank and the SW-dipping Miyar Shear Zone (MSZ) on the southern limb. Current models for the emplacement of the HHC in NW India as a dome structure differ mainly in terms of the roles played by both the ZSZ and the MSZ during the tectonothermal evolution of the HHC. In both the channel flow model and wedge extrusion model, the ZSZ acts as a backstop normal fault along which the high-grade metamorphic rocks of the HHC of Zanskar are exhumed. In contrast, the recently proposed tectonic wedging model argues that the ZSZ and the MSZ correspond to one single detachment system that operates as a subhorizontal backthrust off of the MCT. Thus, the kinematic evolution of the two shear zones, the ZSZ and the MSZ, and their structural, metamorphic and chronological relations appear to be diagnostic features for discriminating the different models. In this paper, structural, metamorphic and geochronological data demonstrate that the MSZ and the ZSZ experienced two distinct kinematic evolutions. As such, the data presented in this paper rule out the hypothesis that the MSZ and the ZSZ constitute one single detachment system, as postulated by the tectonic wedging model. Structural, metamorphic and geochronological data are used to present an alternative tectonic model for the large-scale doming in the NW Indian Himalaya involving early NE-directed tectonics, weakness in the upper crust, reduced erosion at the orogenic front and rapid exhumation along both the ZSZ and the MSZ.

39Ar-40Ar geochronology of mono- and polymetamorphic basement rocks

The 39Ar-40Ar technique is often used to date the metamorphic evolution of basement rocks. The present review article examines systematic aspects of the K-Ar decay system in different mineral chronometers frequently found in mono- and polymetamorphic basements (amphibole, biotite, muscovite/phengite, K-feldspar). A key observation is that the measured dissolution rate of silicates in aqueous fluids is many orders of magnitude faster, and has a much lower activation energy, than the rate of Fickian diffusion of Ar. The effects of this inequality are patchy age zonations, very much like those observed in many U-Pb chronometers, unaccompanied by intra-crystalline bell¬shaped Ar loss profiles. Recognizing the importance of the respective rate constants in field situations leads to re-evaluating the ages and the interpretive paradigms in classic examples such as the Central Alpine "Lepontine" amphibolite event and the Western Alpine eclogitic event.