The tectonometamorphic evolution of the Sesia–Dent Blanche nappes (internal Western Alps): review and synthesis

This study reviews and synthesizes the present knowledge on the Sesia–Dent Blanche nappes, the highest tectonic elements in the Western Alps (Switzerland and Italy), which comprise pieces of pre-Alpine basement and Mesozoic cover. All of the available data are integrated in a crustal-scale kinematic model with the aim to reconstruct the Alpine tectono-metamorphic evolution of the Sesia–Dent Blanche nappes. Although major uncertainties remain in the pre-Alpine geometry, the basement and cover sequences of the Sesia–Dent Blanche nappes are seen as part of a thinned continental crust derived from the Adriatic margin. The earliest stages of the Alpine evolution are interpreted as recording late Cretaceous subduction of the Adria-derived Sesia–Dent Blanche nappes below the South-Alpine domain. During this subduction, several sheets of crustal material were stacked and separated by shear zones that rework remnants of their Mesozoic cover. The recently described Roisan-Cignana Shear Zone of the Dent Blanche Tectonic System represents such a shear zone, indicating that the Sesia–Dent Blanche nappes represent a stack of several individual nappes. During the subsequent subduction of the Piemonte–Liguria Ocean large-scale folding of the nappe stack (including the Roisan-Cignana Shear Zone) took place under greenschist facies conditions, which indicates partial exhumation of the Dent Blanche Tectonic System. The entrance of the Briançonnais micro-continent within the subduction zone led to a drastic change in the deformation pattern of the Alpine belt, with rapid exhumation of the eclogite-facies ophiolite bearing units and thrust propagation towards the foreland. Slab breakoff probably was responsible for allowing partial melting in the mantle and Oligocene intrusions into the most internal parts of the Sesia–Dent Blanche nappes. Finally, indentation of the Adriatic plate into the orogenic wedge resulted in the formation of the Vanzone back-fold, which marks the end of the pervasive ductile deformation within the Sesia–Dent Blanche nappes during the earliest Miocene.

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.

Valle d’Aosta section of the Sesia Zone: multi-stage HP metamorphism and assembly of a rifted continental margin

The spectrum characteristic of the EMC ranges from eclogites (containing omphacite and/or jadeite, garnet, phengite, glaucophane, zoisite, chloritoid, rutile) to phengite schists, calcschists, and marbles, as well as a variety of orthogneisses. Despite the intense polyphase deformation and HP-metamorphic recrystallization, it is possible in some locations to recognize pre-Alpine characteristics in some of the protoliths. For instance, two types of felsic orthogneiss can be distinguished in the Aosta Valley, one derived from Permian granitoids (with local preservation of intrusive contacts, magmatic inclusions, leucocratic veins and other magmatic structures; Stop 3), the other derived from pre-Variscan leuco-monzogranite, such as the building stone mined at the “Argentera” quarry near Settimo Vittone / Montestrutto (Stop 2; so-called “Verde Argento” contains jadeite, phengite, K-feldspar, quartz). Polycyclic and more rarely monocyclic metasediments contain evidence of a complex Alpine PTDt-evolution, locally including relics of their prograde history from blueschist, one or more stages at eclogite facies. Recent petrochronological studies have dated this HP-evolution of the Sesia Zone in some detail. In the area visited, clear evidence of HP-cycling has been identified in one km-size tectonic slice (Stop 1), but not in adjacent parts of the EMC, indicating “yo-yo tectonics”. Partial retrogression and attendant ductile to brittle deformation of the HP-rocks is evident in one of the outcrops (Stop 4). Apart from the four localities in the Sesia Zone, a final outcrop introduces HP-rocks of the adjacent Piemonte oceanic unit, specifically calc-schists and ophiolite members of the “Zermatt-Saas” zone. The hilltop outcrop (Stop 5) displays foliated antigorite schist with peridotite relics (clinopyroxene, spinel) containing lenses derived from doleritic dykes. These fine-grained metarodingites and the folded veins containing Mg-chlorite and titanoclinohumite within serpentinite once again indicate equilibration under low-temperature eclogite facies conditions. However, these units reached that HP stage more than 20 Ma after the youngest eclogite facies imprint recognized in the Sesia Zone. Despite nearly half a century of intense study in the Sesia Zone, the complex assembly of its HP-terranes and their relation to more external parts of the Western Alps remains incompletely understood. This field guide merely introduces a few of the classic outcrops and discusses some of the critical evidence they contain, but it could not incorporate details on each stage of the evolution recognized so far.

Multiple metamorphic stages within an eclogite-facies terrane (Sesia Zone, Western Alps) revealed by Th-U-Pb petrochronology

Convergent plate margins typically experience a transition from subduction to collision dynamics as massive continental blocks enter the subduction channel. Studies of high-pressure rocks indicate that tectonic fragments are rapidly exhumed from eclogite facies to midcrustal levels, but the details of such dynamics are controversial.To understand the dynamics of a subduction channel we report the results of a petrochronological study from the central Sesia Zone, a key element of the internalWestern Alps.This comprises two polymetamorphic basement complexes (Eclogitic Micaschist Complex and Gneiss Minuti Complex) and a thin, dismembered cover sequence (Scalaro Unit) associated with pre-Alpine metagabbros and metasediments (Bonze Unit). Structurally controlled samples from three of these units (Eclogitic Micaschist Complex and Scalaro-Bonze Units) yield unequivocal petrological and geochronological evidence of two distinct high-pressure stages. Ages (U-Th-Pb) of growth zones in accessory allanite and zircon, combined with inclusion and textural relationships, can be tied to the multi-stage evolution of single samples.Two independent tectono-metamorphic ‘slices’ showing a coherent metamorphic evolution during a given time interval have been recognized: the Fondo slice (which includes Scalaro and Bonze rocks) and the Druer slice (belonging to the Eclogitic Micaschist Complex).The new data indicate separate stages of deformation at eclogite-facies conditions for each recognized independent kilometer-sized tectono-metamorphic slice, between ~85 and 60 Ma, with evidence of intermittent decompression (∆P~0.5 GPa) within only the Fondo slice. The evolution path of the Druer slice indicates a different P-T-time evolution with prolonged eclogite-facies metamorphism between ~85 and 75Ma. Our approach, combining structural, petrological and geochronological techniques, yields field-based constraints on the duration and rates of dynamics within a subduction channel.

Geological map of the Scalaro Valley – Sesia Zone (Italian Western Alps)

In the Sesia Zone (Italian Western Alps), slivers of continental crust characterised by an Alpine high-pressure imprint are intermingled with abundant mafic rocks and Mesozoic metasediments. An extensive study of the central Sesia Zone was undertaken to identify and reconstruct the lithological setting of the mono-cyclic sediments of the Scalaro Unit. A new geological map (1:5000) and schematic cross sections across the Scalaro Unit and the adjoining Eclogitic Micaschist Complex are presented here. In order to delimit the size and shape of the mono-metamorphic unit and understand its internal geometry with respect to the poly-metamorphic basement, an integrated approach was used. Linking observations and data across a range of scales, from kilometres in the field down to petrological and chronological data obtained at micrometre scale, we define for the first time the real size and internal geometry of the Scalaro Unit, as well as its large-scale structural context.

Quaternary uplift rates of the Central Anatolian Plateau, Turkey: insights from cosmogenic isochron-burial nuclide dating of the Kızılırmak River terraces

The Central Anatolian Plateau (CAP) in Turkey is a relatively small plateau (300 × 400 km) with moderate average elevations of ∼1 km situated between the Pontide and Tauride orogenic mountain belts. Kızılırmak, which is the longest river (1355 km) within the borders of Turkey, flows within the CAP and slowly incises into lacustrine and volcaniclastic units before finally reaching the Black Sea. We dated the Cappadocia section of the Kızılırmak terraces in the CAP by using cosmogenic burial and isochron-burial dating methods with 10Be and 26Al as their absolute dating can provide insight into long-term incision rates, uplift and climatic changes. Terraces at 13, 20, 75 and 100 m above the current river indicate an average incision rate of 0.051 ± 0.01 mm/yr (51 ± 1 m/Ma) since ∼1.9 Ma. Using the base of a basalt fill above the modern course of the Kızılırmak, we also calculated 0.05–0.06 mm/yr mean incision and hence rock uplift rate for the last 2 Ma. Although this rate might be underestimated due to normal faulting along the valley sides, it perfectly matches our results obtained from the Kızılırmak terraces. Although up to 5–10 times slower, the Quaternary uplift of the CAP is closely related to the uplift of the northern and southern plateau margins respectively.

Benchmark problems for reactive transport modeling of the generation and attenuation of acid rock drainage

Acid rock drainage (ARD) is a problem of international relevance with substantial environmental and economic implications. Reactive transport modeling has proven a powerful tool for the process-based assessment of metal release and attenuation at ARD sites. Although a variety of models has been used to investigate ARD, a systematic model intercomparison has not been conducted to date. This contribution presents such a model intercomparison involving three synthetic benchmark problems designed to evaluate model results for the most relevant processes at ARD sites. The first benchmark (ARD-B1) focuses on the oxidation of sulfide minerals in an unsaturated tailing impoundment, affected by the ingress of atmospheric oxygen. ARD-B2 extends the first problem to include pH buffering by primary mineral dissolution and secondary mineral precipitation. The third problem (ARD-B3) in addition considers the kinetic and pH-dependent dissolution of silicate minerals under low pH conditions. The set of benchmarks was solved by four reactive transport codes, namely CrunchFlow, Flotran, HP1, and MIN3P. The results comparison focused on spatial profiles of dissolved concentrations, pH and pE, pore gas composition, and mineral assemblages. In addition, results of transient profiles for selected elements and cumulative mass loadings were considered in the intercomparison. Despite substantial differences in model formulations, very good agreement was obtained between the various codes. Residual deviations between the results are analyzed and discussed in terms of their implications for capturing system evolution and long-term mass loading predictions.

Implementation and evaluation of permeability-porosity and tortuosity-porosity relationships linked to mineral dissolution-precipitation

Changes of porosity, permeability, and tortuosity due to physical and geochemical processes are of vital importance for a variety of hydrogeological systems, including passive treatment facilities for contaminated groundwater, engineered barrier systems (EBS), and host rocks for high-level nuclear waste (HLW) repositories. Due to the nonlinear nature and chemical complexity of the problem, in most cases, it is impossible to verify reactive transport codes analytically, and code intercomparisons are the most suitable method to assess code capabilities and model performance. This paper summarizes model intercomparisons for six hypothetical scenarios with generally increasing geochemical or physical complexity using the reactive transport codes CrunchFlow, HP1, MIN3P, PFlotran, and TOUGHREACT. Benchmark problems include the enhancement of porosity and permeability through mineral dissolution, as well as near complete clogging due to localized mineral precipitation, leading to reduction of permeability and tortuosity. Processes considered in the benchmark simulations are advective-dispersive transport in saturated media, kinetically controlled mineral dissolution-precipitation, and aqueous complexation. Porosity changes are induced by mineral dissolution-precipitation reactions, and the Carman-Kozeny relationship is used to describe changes in permeability as a function of porosity. Archie’s law is used to update the tortuosity and the pore diffusion coefficient as a function of porosity. Results demonstrate that, generally, good agreement is reached amongst the computer models despite significant differences in model formulations. Some differences are observed, in particular for the more complex scenarios involving clogging; however, these differences do not affect the interpretation of system behavior and evolution.

Time-lapse pressure tomography for characterizing CO2 plume evolution in a deep saline aquifer

A time-lapse pressure tomography inversion approach is applied to characterize the CO2 plume development in a virtual deep saline aquifer. Deep CO2 injection leads to flow properties of the mixed-phase, which vary depending on the CO2 saturation. Analogous to the crossed ray paths of a seismic tomographic experiment, pressure tomography creates streamline patterns by injecting brine prior to CO2 injection or by injecting small amounts of CO2 into the two-phase (brine and CO2) system at different depths. In a first step, the introduced pressure responses at observation locations are utilized for a computationally rapid and efficient eikonal equation based inversion to reconstruct the heterogeneity of the subsurface with diffusivity (D) tomograms. Information about the plume shape can be derived by comparing D-tomograms of the aquifer at different times. In a second step, the aquifer is subdivided into two zones of constant values of hydraulic conductivity (K) and specific storage (Ss) through a clustering approach. For the CO2 plume, mixed-phase K and Ss values are estimated by minimizing the difference between calculated and “true” pressure responses using a single-phase flow simulator to reduce the computing complexity. Finally, the estimated flow property is converted to gas saturation by a single-phase proxy, which represents an integrated value of the plume. This novel approach is tested first with a doublet well configuration, and it reveals a great potential of pressure tomography based concepts for characterizing and monitoring deep aquifers, as well as the evolution of a CO2 plume. Still, field-testing will be required for better assessing the applicability of this approach.

Experimental Characterization and Quantification of Cement-Bentonite Interaction using Core Infiltration Techniques Coupled with X-ray Tomography

Deep geological storage of radioactive waste foresees cementitious materials as reinforcement of tunnels and as backfill. Bentonite is proposed to enclose spent fuel canisters and as drift seals. Sand/bentonite (s/b) is foreseen as backfill material of access galleries or as drift seals. The emplacement of cementitious material next to clay material generates an enormous chemical gradient in pore-water composition that drives diffusive solute transport. Laboratory studies and reactive transport modeling predicted significant mineral alteration at and near interfaces, mainly resulting in a decrease of porosity in bentonite. The goal of this thesis was to characterize and quantify the cement/bentonite interactions both spatially and temporally in laboratory experiments. A newly developed mobile X-ray transparent core infiltration device was used to perform X-ray computed tomography (CT) scans without interruption of running experiments. CT scans allowed tracking the evolution of the reaction plume and changes in core volume/diameter/density during the experiments. In total 4 core infiltration experiments were carried out for this study with the compacted and saturated cores consisting of MX-80 bentonite and sand/MX-80 bentonite mixture (s/b; 65/35%). Two different high-pH cementitious pore-fluids were infiltrated: a young (early) ordinary Portland cement pore-fluid (APWOPC; K+–Na+–OH-; pH 13.4; ionic strength 0.28 mol/kg) and a young ‘low-pH’ ESDRED shotcrete pore-fluid (APWESDRED; Ca2+–Na+–K+–formate; pH 11.4; ionic strength 0.11 mol/kg). The experiments lasted between 1 and 2 years. In both bentonite experiments, the hydraulic conductivity was strongly reduced after switching to high-pH fluids, changing eventually from an advective to a diffusion-dominated transport regime. The reduction was mainly induced by mineral precipitation and possibly partly also by high ionic strength pore-fluids. Both bentonite cores showed a volume reduction and a resulting transient flow in which pore-water was squeezed out during high-pH infiltration. The outflow chemistry was characterized by a high ionic strength, while chloride in the initial pore water got replaced as main anionic charge carrier by sulfate, originating from gypsum dissolution. The chemistry of the high-pH fluids got strongly buffered by the bentonite, consuming hydroxide and in case of APWESDRED also formate. Hydroxide got consumed by mineral reactions (saponite and possibly talc and brucite precipitation), while formate being affected by bacterial degradation. Post-mortem analysis showed reaction zones near the inlet of the bentonite core, characterized by calcium and magnesium enrichment, consisting predominately of calcite and saponite, respectively. Silica got enriched in the outflow, indicating dissolution of silicate-minerals, identified as preferentially cristobalite. In s/b, infiltration of APWOPC reduced the hydraulic conductivity strongly, while APWESDRED infiltration had no effect. The reduction was mainly induced by mineral precipitation and probably partly also by high ionic strength pore-fluids. Not clear is why the observed mineral precipitates in the APWESDRED experiment had no effect on the fluid flow. Both s/b cores showed a volume expansion along with decreasing ionic strengths of the outflow, due to mineral reactions or in case of APWESDRED infiltration also mediated by microbiological activity, consuming hydroxide and formate, respectively. The chemistry of the high-pH fluids got strongly buffered by the s/b. In the case of APWESDRED infiltration, formate reached the outflow only for a short time, followed by enrichment in acetate, indicating most likely biological activity. This was in agreement to post-mortem analysis of the core, observing black spots on the inflow surface, while the sample had a rotten-egg smell indicative of some sulfate reduction. Post-mortem analysis showed further in both cores a Ca-enrichment in the first 10 mm of the core due to calcite precipitation. Mg-enrichment was only observed in the APWOPC experiment, originating from newly formed saponite. Silica got enriched in the outflow of both experiments, indicating dissolution of silicate-minerals, identified in the OPC experiment as cristobalite. The experiments attested an effective buffering capacity for bentonite and s/b, a progressing coupled hydraulic-chemical sealing process and also the preservation of the physical integrity of the interface region in this setup with a total pressure boundary condition on the core sample. No complete pore-clogging was observed but the hydraulic conductivity got rather strongly reduced in 3 experiments, explained by clogging of the intergranular porosity (macroporosity). Such a drop in hydraulic conductivity may impact the saturation time of the buffer in a nuclear waste repository, although the processes and geometry will be more complex in repository situation.