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.

Ar-40/Ar-39 Evidence for Middle Proterozoic (1300-1500 Ma) Slow Cooling of the Southern Black Hills, South Dakota, Midcontinent, North America: Implications for Early Proterozoic P-T Evolution and Posttectonic Magmatism

Ar-40/Ar-39 total gas and plateau dates from muscovite and biotite in the southern Black Hills, South Dakota, provide evidence for a period of Middle Proterozoic slow cooling. Early Proterozoic (1600-1650 Ma) mica dates were obtained from metasedimentary rocks located in a synformal structure between the Harney Peak and Bear Mountain domes and also south of Bear Mountain. Metamorphic rocks from the dome areas and undeformed samples of the similar to 1710 Ma Harney Peak Granite (HPG) yield Middle Proterozoic mica dates (similar to 1270-1500 Ma). Two samples collected between the synform and Bear Mountain dome yield intermediate total gas mica dates of similar to 1550 Ma. We suggest two end-member interpretations to explain the map pattern of cooling ages: (1) subhorizontal slow cooling of an area which exhibits variation in mica Ar retention intervals or (2) mild folding of a Middle Proterozoic (similar to 1500 Ma) similar to 300 degrees C isotherm. According to the second interpretation, the preservation of older dates between the domes may reflect reactivation of a preexisting synformal structure (and downwarping of relatively cold rocks) during a period of approximately east-west contraction and slow uplift during the Middle Proterozoic. The mica data, together with hornblende data from the Black Hills published elsewhere, indicate that the ambient country-rock temperature at the 3-4 kbar depth of emplacement of the HPG was between 350 degrees C and 500 degrees C, suggesting that the average upper crustal geothermal gradient was 25 degrees-40 degrees C/km prior to intrusion. The thermochronologic data suggest HPG emplacement was followed by a similar to 200 m.y. period of stability and tectonic quiescence with little uplift. We propose that crust thickened during the Early Proterozoic was uplifted and erosionally(?) thinned prior to similar to 1710 Ma and that the HPG magma was emplaced into isostatically stable crust of relatively normal thickness. We speculate that uplift and crustal thinning prior to HPG intrusion was the result of differential thinning of the subcrustal lithosphere beneath the Black Hills. If so, this process would have also caused an increase in mantle heat flux across the Moho and triggered vapor-absent melting of biotite to produce the HPG magma. This scenario for posttectonic granite generation is supported, in part, by the fact that in the whole of the Black Hills, the HPG is spatially associated with the deepest exposed Early Proterozoic country rock.

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.

Geochemistry of bulk samples and various minerals of sediment core CRP-1 from the Ross Sea, Antarctica

Thirty-nine medium and fine grained sandstones from between 19,26 and 147,23 mbsf in the Cape Roberts-l core (CRP-1) were analysed for 10 major and 16 trace elements. Using whole-lock compositions, 9 samples were selected for analyses of mineral and glass grains by energy dispersive electron microscope. Laser-Ablation Mass-Spectrometry was used to determine rare earth elements and 14 additional trace elements in glass shards, pyroxenes and feldspars in order to examine their contribution to the bulk rock chemistry. Geochemical data reveal the major contribution played by the Granite Harbour Intrusives to the whole rock composition, even if a significant input is supplied by McMurdo volcanics and Ferrar dolerite pyroxenes McMurdo volcanics were studied in detail; they appeal to derive from a variety of litologies, and a dominant role of wind transpoitation from exposures of volcanic rocks may be inferred from the contemporary occurrence of different compositions at all depths. Only at 116.55 mbsf was a thin layer of tephra found, linked to an explosive eruption McMurdo volcanic rocks exhibit larger abundances at depths above 62 mbsf, in correspondence with the onset of volcanic activity in the McMurdo Sound area. From 62 mbsf to the bottom of the core, McMurdo volcanics are less abundant and probably issued from some centres in the McMurdo Sound region. However, available data do not allow the exclusion of wind transport from some eruptive centres active in north Victoria Land at the beginning of the Miocene Epoch.

Geochemistry of basement clasts from sediment core CRP-3 in the ross Sea, Antarctica

Distribution patterns and petrographical and mineral chemistry data are described for the most representative basement lithologies occuring as clast in the c. 824 m thick Tertiary sedimentary sequence at the CRP-3 drillsite. These are granule to bolder grain size clasts of igneous and metamorphic rocks. Within the basement clast assemblage, granitoid pebbles are the predominant lithology. They consist of dominant grey biotic-bearing monzogranite, pink biotite-hornblende monzogranite, and biotite-bearing leucomomonzgranite. Minor lithologies include: actinolite-bearing leucotonalite, microgranite, biotite-hornblende quartz-monzonitic porphyr, and foliated biotic leucomonzogranite. Metamorphic clasts include rocks of both granitic and sedimentary derivation. They include mylonitic biotic orthogneiss, with or without garnet, muscovite-bearing quartzite, sillimanite-biotite paragneiss, biotite meta-sandstone, biotite-spotted schist, biotite-clacite-clinoamphibole meta-feldspathic arenite, biotite-calcite-clinozoisite meta-siltstone, biotite±clinoamphibole meta-marl, and graphite-bearing marble. As in previous CRP drillcores, the ubiquitous occurence of biotite±hornblende monzogranite pebbles is indicative of a local provenance, closely mirroring the dominance of these lithologies in the on-shore basement, where the Cambro-Ordovician Granite Harbour Intrusive Complex forms the most extensively exposed rock unit.

Sm-Nd isotopic investigation for rocks of Prince Charles Mountains, East Antarctica

The Prince Charles Mountains have been subject to extensive geological and geophysical investigations by former Soviet, Russian and Australian scientists from the early 1970s. In this paper we summarise, and review available geological and isotopic data, and report results of new isotopic studies (Sm-Nd, Pb-Pb, and U-Pb SHRIMP analyses); field geological data obtained during the PCMEGA 2002/2003 are utilised. The structure of the region is described in terms of four tectonic terranes. Those include Archaean Ruker, Palaeoproterozoic Lambert, Mesoproterozoic Fisher, and Meso- to Neoproterozoic Beaver Terranes. Pan-African activities (granite emplacement and probably tectonics) in the Lambert Terrane are reported. We present a summary of the composition of these terranes, discuss their origin and relationships. We also outline the most striking geological features, and problems, and try to draw attention to those rocks and regional geological features which are important in understanding the composition and evolution of the PCM and might suggest targets for further investigations.

Tab. 2: U-Pb ion-microprobe ages from Taymyr peninsula

We present the initial results of a U-Th-Pb zircon ion-microprobe investigation on samples from the Central Belt of Taimyr, in order to constrain its tectono-magmatic evolution. The zircon samples are from a deformed twomica granite (Faddey Massif), deformed metamorphosed gabbroic dike entrained as pods and lenses within metamorphosed tholeiitic basalts of the Kunar-Mod volcanic suite (Klyaz'ma River region), a metamorphosed rhyolite of the same volcanic suite overlying the basic metavolcanic rocks, as well as an undeformed dolerite dike which intrudes the metamorphosed Kunar-Mod basic volcanic rocks. Preliminary results on zircons from the two-mica granite suggest a crystallization age of ~630 Ma for this rock, with inheritance from assimilated crust 840 Ma to 1.1 Ga in age. In the Klyaz'ma River region, zircons from the meta-rhyolite yield a concordant age of -630 Ma. Zircons from the entrained metagabbroic dikes have so far yielded an age of -615 Ma (1 grain), as well as Archean ages (5 grains, concordant at 2.6-2.8 Ga). It seems likely that the Archean grains represent assimilation of older crustal material. Zircons from the post-tectonic dolerite dike have a bimodal age distribution. A well-defined younger age of 281 ±9 Ma is interpreted to represent the crystallization age of the dike, while older, concordant ages of 2.6-2.9 Ga likely represent assimilation of Archean crust (Siberian craton at depth). Several important conclusions can be drawn from the data. (1) The mafic and felsic lithologies of the Kunar-Mod volcanic suite are genetically related and should be the same age. Ages of-630 Ma (meta-rhyolite) and -615 Ma (metagabbroic dikes representing the latest stage of mafic magmatism associated the Kunar-Mod suite) suggest that these lithologies may be the same age, but more data are required to confirm this hypothesis. (2) The 630 Ma two-mica granite is similar in age to the time of high-grade metamorphism, suggesting that syntectonic granite emplacement accompanied obduction of the accretionary Central Belt to the Siberian craton. (3) An Early Permian age is well defined for the undeformed dolerite dike. Dolerite dikes occur across the whole of Taimyr, but are deformed to the south. If related, this single magmatic event pre-dates Permo-Triassic Siberian trap magmatism. Furthermore, it suggests that deformation was localized to southeastern Taimyr.