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.

Geometry and kinematics of the Roisan-Cignana Shear Zone, and the orogenic evolution of the Dent Blanche Tectonic System (Western Alps)

The Dent Blanche Tectonic System (DBTS) is a composite thrust sheet derived from the previously thinned passive Adriatic continental margin. A kilometric high-strain zone, the Roisan-Cignana Shear Zone (RCSZ) defines the major tectonic boundary within the DBTS and separates it into two subunits, the Dent Blanche s.s. nappe to the northwest and the Mont Mary nappe to the southeast. Within this shear zone, tectonic slices of Mesozoic and pre-Alpine meta-sediments became amalgamated with continental basement rocks of the Adriatic margin. The occurrence of high pressure assemblages along the contact between these tectonic slices indicates that the amalgamation occurred prior to or during the subduction process, at an early stage of the Alpine orogenic cycle. Detailed mapping, petrographic and structural analysis show that the Roisan-Cignana Shear Zone results from several superimposed Alpine structural and metamorphic stages. Subduction of the continental fragments is recorded by blueschist-facies deformation, whereas the Alpine collision is reflected by a greenschist facies overprint associated with the development of large-scale open folds. The postnappe evolution comprises the development of low-angle brittle faults, followed by large-scale folding (Vanzone phase) and finally brittle extensional faults. The RCSZ shows that fragments of continental crust had been torn off the passive continental margin prior to continental collision, thus recording the entire history of the orogenic cycle. The role of preceding Permo-Triassic lithospheric thinning, Jurassic rifting, and ablative subduction processes in controlling the removal of crustal fragments from the reactivated passive continental margin is discussed. Results of this study constrain the temporal sequence of the tectono-metamorphic processes involved in the assembly of the DBTS, but they also show limits on the interpretation. In particular it remains difficult to judge to what extent precollisional rifting at the Adriatic continental margin preconditioned the efficiency of convergent processes, i.e. accretion, subduction, and orogenic exhumation.

Distributed transpressive continental deformation: The Varto Fault Zone, eastern Turkey

The convergence between the Eurasian and Arabian plates has created a complicated structural setting in the Eastern Turkish high plateau (ETHP), particularly around the Karlıova Triple Junction (KTJ) where the Eurasian, Arabian, and Anatolian plates intersect. This region of interest includes the junction of the North Anatolian Shear Zone (NASZ) and the East Anatolian Shear Zone (EASZ), which forms the northern border of the westwardly extruding Anatolian Scholle and the western boundary of the ETHP, respectively. In this study, we focused on a poorly studied component of the KTJ, the Varto Fault Zone (VFZ), and the adjacent secondary structures, which have complex structural settings. Through integrated analyses of remote sensing and field observations, we identified a widely distributed transpressional zone where the Varto segment of the VFZ forms the most northern boundary. The other segments, namely, the Leylekdağ and Çayçatı segments, are oblique-reverse faults that are significantly defined by uplifted topography along their strikes. The measured 515 and 265 m of cumulative uplifts for Mt. Leylek and Mt. Dodan, respectively, yield a minimum uplift rate of 0.35 mm/a for the last 2.2 Ma. The multi-oriented secondary structures were mostly correlated with “the distributed strike-slip” and “the distributed transpressional” in analogue experiments. The misfits in strike of some of secondary faults between our observations and the experimental results were justified by about 20° to 25° clockwise restoration of all relevant structures that were palaeomagnetically measured to have happened since ~ 2.8 Ma ago. Our detected fault patterns and their true nature are well aligned as being part of a transpressional tectonic setting that supports previously suggested stationary triple junction models.

Constraints on the thermal evolution of the Adriatic margin during Jurassic continental break-up: U–Pb dating of rutile from the Ivrea–Verbano Zone, Italy

The Ivrea–Verbano Zone (IVZ), northern Italy, exposes an attenuated section through the Permian lower crust that records high-temperature metamorphism under lower crustal conditions and a protracted history of extension and exhumation associated partly with the Jurassic opening of the Alpine Tethys ocean. This study presents SHRIMP U–Pb geochronology of rutile from seven granulite facies metapelites from the base of the IVZ, collected from locations spanning ~35 km along the strike of Paleozoic fabrics. Rutile crystallised during Permian high-temperature metamorphism and anatexis, yet all samples give Jurassic rutile U–Pb ages that record cooling through 650–550 °C. Rutile age distributions are dominated by a peak at ~160 Ma, with a subordinate peak at ~175 Ma. Both ~160 and ~175 Ma age populations show excellent agreement between samples, indicating that the two distinctive cooling stages they record were synchronous on a regional scale. The ~175 Ma population is interpreted to record cooling in the footwall of rift-related faults and shear zones, for which widespread activity in the Lower Jurassic has been documented along the western margin of the Adriatic plate. The ~160 Ma age population postdates the activity of all known rift-related structures within the Adriatic margin, but coincides with extensive gabbroic magmatism and exhumation of sub-continental mantle to the floor of the Alpine Tethys, west of the Ivrea Zone. We propose that this ~160 Ma early post-rift age population records regional cooling following episodic heating of the distal Adriatic margin, likely related to extreme lithospheric thinning and associated advection of the asthenosphere to shallow levels. The partial preservation of the ~175 Ma age cluster suggests that the post-rift (~160 Ma) heating pulse was of short duration. The regional consistency of the data presented here, which is in contrast to many other thermochronometers in the IVZ, demonstrates the value of the rutile U–Pb technique for probing the thermal evolution of high-grade metamorphic terrains. In the IVZ, a significant decoupling between Zr-in-rutile temperatures and U–Pb ages of rutile is observed, with the two systems recording events ~120 Ma apart.

Inherited structural controls on fault geometry, architecture and hydrothermal activity: an example from Grimsel Pass, Switzerland

Exhumed faults hosting hydrothermal systems provide direct insight into relationships between faulting and fluid flow, which in turn are valuable for making hydrogeological predictions in blind settings. The Grimsel Breccia Fault (Aar massif, Central Swiss Alps) is a late Neogene, exhumed dextral strike-slip fault with a maximum displacement of 25–45 m, and is associated with both fossil and active hydrothermal circulation. We mapped the fault system and modelled it in three dimensions, using the distinctive hydrothermal mineralisation as well as active thermal fluid discharge (the highest elevation documented in the Alps) to reveal the structural controls on fluid pathway extent and morphology. With progressive uplift and cooling, brittle deformation inherited the mylonitic shear zone network at Grimsel Pass; preconditioning fault geometry into segmented brittle reactivations of ductile shear zones and brittle inter-shear zone linkages. We describe ‘pipe’-like, vertically oriented fluid pathways: (1) within brittle fault linkage zones and (2) through alongstrike- restricted segments of formerly ductile shear zones reactivated by brittle deformation. In both cases, low-permeability mylonitic shear zones that escaped brittle reactivation provide important hydraulic seals. These observations show that fluid flow along brittle fault planes is not planar, but rather highly channelised into sub-vertical flow domains, with important implications for the exploration and exploitation of geothermal energy.

Sandbox experiments on basin inversion: testing the influence of basin orientation and basin fill

Analogue modelling experiments using brittle materials are performed to study the inversion of extensional structures. Asymmetric grabens of two different orientations are first created during a phase of extension and progressively filled. They are subsequently shortened in the same direction. The aim of our experiments is to determine factors affecting the style of deformation during inversion. We specifically investigate variations in thickness and distribution of strong and weak layers constituting the graben fill and in initial basin orientation. The main advantage of our experimental set-up is that we have a complete control on graben location, width, infill and orientation before inversion. The experiments show that shortening results only in limited reactivation of pre-existing normal faults. In general, forward thrusts and backthrusts cut across normal faults into the footwall of the graben. The forward thrusts either propagate parallel to the enveloping surface of faulted blocks or they cut across basin-limiting normal faults at various angles. The graben fill is mechanically extruded by displacement along forward thrusts that accommodate most of the shortening. Both pre-existing faults and weak graben fill act as zones of weakness during inversion and determine the orientation and location of both backthrusts and forward thrusts. The results of our experiments conform well to natural examples of inverted graben structures.

Comparisons between analogue and numerical models of thrust wedge development

Analogue and finite element numerical models with frictional and viscous properties are used to model thrust wedge development. Comparison between model types yields valuable information about analogue model evolution, scaling laws and the relative strengths and limitations of the techniques. Both model types show a marked contrast in structural style between ‘frictional-viscous domains’ underlain by a thin viscous layer and purely ‘frictional domains’. Closely spaced thrusts form a narrow and highly asymmetric fold-and-thrust belt in the frictional domain, characterized by in-sequence propagation of forward thrusts. In contrast, the frictional-viscous domain shows a wide and low taper wedge and a thrust belt with a more symmetrical vergence, with both forward and back thrusts. The frictional-viscous domain numerical models show that the viscous layer initially simple shears as deformation propagates along it, while localized deformation resulting in the formation of a pop-up structure occurs in the overlying frictional layers. In both domains, thrust shear zones in the numerical model are generally steeper than the equivalent faults in the analogue model, because the finite element code uses a non-associated plasticity flow law. Nevertheless, the qualitative agreement between analogue and numerical models is encouraging. It shows that the continuum approximation used in numerical models can be used to model frictional materials, such as sand, provided caution is taken to properly scale the experiments, and some of the limitations are taken into account.

An insight into the breakup of Gondwana: Identifying events through low-temperature thermochronology from the basement rocks of Madagascar

Fission track analysis was applied to the Precambrian suites of Madagascar in order to identify the lower-temperature cooling histories and their relationships to the Phanerozoic events that affected the island. Apatite ages range from 431 to 68 Ma, and zircon ages range from 452 to 238 Ma. Thermochronologically, the island can be divided into a southern, central, and northern region each with a subdivision on an east-west basis. The southern region is sharply separated from the central region by strongly contrasting apparent apatite ages over the northwest-southeast striking Ranotsara Shear Zone (RSZ). The change in apparent ages over the RSZ is indicative of later reactivation along younger brittle faults. The central region has the oldest ages of the island and has a diffuse contact to the third region northward. Along the entire western margin of the Precambrian basement initial Paleozoic exhumation was followed by heating (burial by sediments) during Jurassic and Cretaceous times. A decrease in ages along the eastern margin from 119 to 68 Ma coincides with the predicted positions of the Marion hot spot after effects of erosion are considered. On the other hand, these ages may represent progressive opening of the margin in a southward direction together with associated denudation of the rift shoulder. The eastern part of the central region has remained very stable since at least Devonian times, undergoing only long-term very slow exhumation at rates of 1–5 m/Myr.

Fault development and interaction in distributed strike-slip shear zones: an experimental approach

Analogue model experiments using both brittle and viscous materials were performed to investigate the development and interaction of strike-slip faults in zones of distributed shear deformation. At low strain, bulk dextral shear deformation of an initial rectangular model is dominantly accommodated by left-stepping, en echelon strike-slip faults (Riedel shears, R) that form in response to the regional (bulk) stress field. Push-up zones form in the area of interaction between adjacent left-stepping Riedel shears. In cross sections, faults bounding push-up zones have an arcuate shape or merge at depth. Adjacent left-stepping R shears merge by sideways propagation or link by short synthetic shears that strike subparallel to the bulk shear direction. Coalescence of en echelon R shears results in major, through-going faults zones (master faults). Several parallel master faults develop due to the distributed nature of deformation. Spacing between master faults is related to the thickness of the brittle layers overlying the basal viscous layer. Master faults control to a large extent the subsequent fault pattern. With increasing strain, relatively short antithetic and synthetic faults develop mostly between old, but still active master faults. The orientation and evolution of the new faults indicate local modifications of the stress field. In experiments lacking lateral borders, closely spaced parallel antithetic faults (cross faults) define blocks that undergo clockwise rotation about a vertical axis with continuing deformation. Fault development and fault interaction at different stages of shear strain in our models show similarities with natural examples that have undergone distributed shear.

Tectono-metamorphic evolution of a nappe stack: A case study of the Swiss Alps

Fold-and-thrust belts are prominent structures that occur at the front of compressional orogens. To unravel the tectonic and metamorphic evolution of such complexes, kinematic investigations, quantitative microstructural analysis and geothermometry (calcite–graphite, calcite–dolomite) were performed on carbonate mylonites from thrust faults of the Helvetic nappe stack in Central Switzerland. Paleo-isotherms of peak temperature conditions and cooling stages (fission track) of the nappe pile were reconstructed in a vertical section and linked with the microstructural and kinematic evolution. Mylonitic microstructures suggest that under metamorphic conditions close to peak temperature, strain was highly localized within thrust faults where deformation temperatures spatially continuously increased in both directions, from N to S within each nappe and from top–down in the nappe stack, covering a temperature range of 180–380 °C. Due to the higher metamorphic conditions, thrusting of the lowermost nappe, the Doldenhorn nappe, was accompanied by a much more pronounced nappe internal ductile deformation of carbonaceous rock types than was the case for the overlying Wildhorn- and Gellihorn nappes. Ongoing thrusting brought the Doldenhorn nappe closer to the surface. The associated cooling resulted in a freezing in of the paleo-isotherms of peak metamorphic conditions. Contemporaneous shearing localized in the basal thrust, initially still in the ductile deformation regime and finally as brittle faulting and cataclasis inducing ultimately an inverse metamorphic zonation. With ongoing exhumation and the formation of the Helvetic antiformal nappe stack, a bending of large-scale tectonic structures (thrusts, folds), peak temperature isotherms and cooling isotherms occurred. While this local bending can directly be attributed to active deformation underneath the section investigated up to times of 2–3 ma, a more homogeneous uplift of the entire region is suggested for the very late and still active exhumation stage.