Lake sediments as natural seismographs: A compiled record of Late Quaternary earthquakes in Central Switzerland and its implication for Alpine deformation

Central Switzerland lies tectonically in an intraplate area and recurrence rates of strong earthquakes exceed the time span covered by historic chronicles. However, many lakes are present in the area that act as natural seismographs: their continuous, datable and high-resolution sediment succession allows extension of the earthquake catalogue to pre-historic times. This study reviews and compiles available data sets and results from more than 10 years of lacustrine palaeoseismological research in lakes of northern and Central Switzerland. The concept of using lacustrine mass-movement event stratigraphy to identify palaeo-earthquakes is showcased by presenting new data and results from Lake Zurich. The Late Glacial to Holocene mass-movement units in this lake document a complex history of varying tectonic and environmental impacts. Results include sedimentary evidence of three major and three minor, simultaneously triggered basin-wide lateral slope failure events interpreted as the fingerprints of palaeoseismic activity. A refined earthquake catalogue, which includes results from previous lake studies, reveals a non-uniform temporal distribution of earthquakes in northern and Central Switzerland. A higher frequency of earthquakes in the Late Glacial and Late Holocene period documents two different phases of neotectonic activity; they are interpreted to be related to isostatic post-glacial rebound and relatively recent (re-)activation of seismogenic zones, respectively. Magnitudes and epicentre reconstructions for the largest identified earthquakes provide evidence for two possible earthquake sources: (i) a source area in the region of the Alpine or Sub-Alpine Front due to release of accumulated north-west/south-east compressional stress related to an active basal thrust beneath the Aar massif; and (ii) a source area beneath the Alpine foreland due to reactivation of deep-seated strike-slip faults. Such activity has been repeatedly observed instrumentally, for example, during the most recent magnitude 4.2 and 3.5 earthquakes of February 2012, near Zug. The combined lacustrine record from northern and Central Switzerland indicates that at least one of these potential sources has been capable of producing magnitude 6.2 to 6.7 events in the past.

Temporal variations of perfluoroalkyl substances and polybrominated diphenyl ethers in alpine snow

The occurrence and temporal variation of 18 perfluoroalkyl substances (PFASs) and 8 polybrominated diphenyl ethers (PBDEs) in the European Alps was investigated in a 10 m shallow firn core from Colle Gnifetti in the Monte Rosa Massif (4455 m above sea level). The firn core encompasses the years 1997-2007. Firn core sections were analyzed by liquid chromatography-tandem mass spectrometry (PFASs) and gas chromatography-mass spectrometry (PBDEs). We detected 12 PFASs and 8 PBDEs in the firn samples. Perfluorobutanoic acid (PFBA; 0.3-1.8 ng L(-1)) and perfluorooctanoic acid (PFOA; 0.2-0.6 ng L(-1)) were the major PFASs while BDE 99 (

Neotethys closure history of Anatolia: insights from 40Ar–39Ar geochronology and P–T estimation in high pressure metasedimentary rocks

The multiple high-pressure (HP), low-temperature (LT) metamorphic units of Western and Central Anatolia offer a great opportunity to investigate the subduction- and continental accretion-related evolution of the eastern limb of the long-lived Aegean subduction system. Recent reports of the HP–LT index mineral Fe-Mg-carpholite in three metasedimentary units of the Gondwana-derived Anatolide–Tauride continental block (namely the Afyon Zone, the Ören Unit and the southern Menderes Massif) suggest a more complicated scenario than the single-continental accretion model generally put forward in previous studies. This study presents the first isotopic dates (white mica 40Ar–39Ar geochronology), and where possible are combined with P–T estimates (chlorite thermometry, phengite barometry, multi-equilibrium thermobarometry), on carpholite-bearing rocks from these three HP–LT metasedimentary units. It is shown that, in the Afyon Zone, carpholite-bearing assemblages were retrogressed through greenschist-facies conditions at c. 67–62 Ma. Early retrograde stages in the Ören Unit are dated to 63–59 Ma. In the Kurudere–Nebiler Unit (HP Mesozoic cover of the southern Menderes Massif), HP retrograde stages are dated to c. 45 Ma, and post-collisional cooling to c. 26 Ma. These new results support that the Ören Unit represents the westernmost continuation of the Afyon Zone, whereas the Kurudere–Nebiler Unit correlates with the Cycladic Blueschist Unit of the Aegean Domain. In Western Anatolia, three successive HP–LT metamorphic belts thus formed: the northernmost Tavşanlı Zone (c. 88–82 Ma), the Ören–Afyon Zone (between 70 and 65 Ma), and the Kurudere–Nebiler Unit (c. 52–45 Ma). The southward younging trend of the HP–LT metamorphism from the upper and internal to the deeper and more external structural units, as in the Aegean Domain, points to the persistence of subduction in Western Anatolia between 93–90 and c. 35 Ma. After the accretion of the Menderes–Tauride terrane, in Eocene times, subduction stopped, leading to continental collision and associated Barrovian-type metamorphism. Because, by contrast, the Aegean subduction did remain active due to slab roll-back and trench migration, the eastern limb (below Southwestern Anatolia) of the Hellenic slab was dramatically curved and consequently teared. It therefore is suggested that the possibility for subduction to continue after the accretion of buoyant (e.g. continental) terranes probably depends much on palaeogeography.

Mid- to late Holocene Indian Ocean Monsoon variability recorded in four speleothems from Socotra Island, Yemen

Four stalagmites covering the last 7.0 ka were sampled on Socotra, an island in the northern Indian Ocean to investigate the evolution of the northeast Indian Ocean Monsoon (IOM) since the mid Holocene. On Socotra, rain is delivered at the start of the southwest IOM in May–June and at the start of the northeast IOM from September to December. The Haggeher Mountains act as a barrier forcing precipitation brought by the northeast winds to fall preferentially on the eastern side of the island, where the studied caves are located. δ18O and δ13C and Mg/Ca and Sr/Ca signals in the stalagmites reflect precipitation amounts brought by the northeast winds. For stalagmite STM6, this amount effect is amplified by kinetic effects during calcite deposition. Combined interpretation of the stalagmites' signals suggest a weakening of the northeast precipitation between 6.0 and 3.8 ka. After 3.8 ka precipitation intensities remain constant with two superimposed drier periods, between 0 and 0.6 ka and from 2.2 to 3.8 ka. No link can be established with Greenland ice cores and with the summer IOM variability. In contrast to the stable northeast rainy season suggested by the records in this study, speleothem records from western Socotra indicate a wettening of the southwest rainy season on Socotra after 4.4 ka. The local wettening of western Socotra could relate to a more southerly path (more over the Indian Ocean) taken by the southwest winds. Stalagmite STM5, sampled at the fringe between both rain areas displays intermediate δ18O values. After 6.2 ka, similar precipitation changes are seen between eastern Socotra and northern Oman indicating that both regions are affected similarly by the monsoon. Different palaeoclimatologic records from the Arabian Peninsula currently located outside the ITCZ migration pathway display an abrupt drying around 6 ka due to their disconnection from the southwest rain influence. Records that are nowadays still receiving rain by the southwest winds, suggest a more gradual drying reflecting the weakening of the southwest monsoon.

A new sensitive method for the quantification of glyoxal and methylglyoxal in snow and ice by stir bar sorptive extraction and liquid desorption-HPLC-ESI-MS

In this study, the development of a new sensitive method for the analysis of alpha-dicarbonyls glyoxal (G) and methylglyoxal (MG) in environmental ice and snow is presented. Stir bar sorptive extraction with in situ derivatization and liquid desorption (SBSE-LD) was used for sample extraction, enrichment, and derivatization. Measurements were carried out using high-performance liquid chromatography coupled to electrospray ionization tandem mass spectrometry (HPLC-ESI-MS/MS). As part of the method development, SBSE-LD parameters such as extraction time, derivatization reagent, desorption time and solvent, and the effect of NaCl addition on the SBSE efficiency as well as measurement parameters of HPLC-ESI-MS/MS were evaluated. Calibration was performed in the range of 1–60 ng/mL using spiked ultrapure water samples, thus incorporating the complete SBSE and derivatization process. 4-Fluorobenzaldehyde was applied as internal standard. Inter-batch precision was <12 % RSD. Recoveries were determined by means of spiked snow samples and were 78.9 ± 5.6 % for G and 82.7 ± 7.5 % for MG, respectively. Instrumental detection limits of 0.242 and 0.213 ng/mL for G and MG were achieved using the multiple reaction monitoring mode. Relative detection limits referred to a sample volume of 15 mL were 0.016 ng/mL for G and 0.014 ng/mL for MG. The optimized method was applied for the analysis of snow samples from Mount Hohenpeissenberg (close to the Meteorological Observatory Hohenpeissenberg, Germany) and samples from an ice core from Upper Grenzgletscher (Monte Rosa massif, Switzerland). Resulting concentrations were 0.085–16.3 ng/mL for G and 0.126–3.6 ng/mL for MG. Concentrations of G and MG in snow were 1–2 orders of magnitude higher than in ice core samples. The described method represents a simple, green, and sensitive analytical approach to measure G and MG in aqueous environmental samples.

Natural Tracers in recent groundwaters from different Alpine aquifers

Groundwater with underground residence times between days and a few years have been investigated over more than 20 years from 487 remote sites located in different aquifer types in the Alpine belt. Analysis of the data reveals that groundwaters evolved in crystalline, evaporite, carbonate, molasse, and flysch aquifers can be clearly distinguished based on their major and trace element composition and degree of mineralisation. A further subdivision can be made even within one aquifer type based on the trace element compositions, which are characteristic for the lithologic environment. Major and trace element concentrations can be quantitatively described by interaction of the groundwater with the aquifer- specific mineralogy along the flow path. Because all investigated sites show minimal anthropogenic influences, the observed concentration ranges represent the natural background concentrations and can thus serve as a “geo-reference” for recent groundwaters from these five aquifer types. This “geo-reference” is particularly useful for the identification of groundwater contamination. It further shows that drinking water standards can be grossly exceeded for critical elements by purely natural processes

Recharge Areas and Geochemical Evolution of Groundwater in an Alluvial Aquifer System in the Sultanate of Oman

A regional hydrogeochemical model was developed to evaluate the geochemical evolution of different groundwaters in an alluvial aquifer system in the Interior of Oman. In combination with environmental isotopes the model is able to extract qualitative and quantitative information about recharge, groundwater flow paths and hydraulic connections between different aquifers. The main source of water to the alluvial aquifer along the flow paths ofWadi Abyadh andWadi M’uaydin in the piedmont is groundwater from the high-altitude areas of the Jabal Akhdar and local infiltration along the wadi channels. In contrast, the piedmont alluvial aquifer alongWadi Halfayn is primarily replenished by lateral recharge from the ophiolite foothills to the east besides smaller contributions from the Jabal Akhdar and local infiltration. Further down gradient in the Southern Alluvial Plain aquifer a significant source of recharge is direct infiltration of rain and surface runoff, originating from a moisture source that approaches Oman from the south. The model shows that the main geochemical evolution of the alluvial groundwaters occurs along the flow path from the piedmont to the Southern Alluvial Plain, where dedolomitization is responsible for the observed changes in the chemical and carbon isotope composition in these waters.

Transient sediment supply in a high-altitude Alpine environment evidenced through a 10Be budget of the Etages catchment (French Western Alps)

Although beryllium-10 (10Be) concentrations in stream sediments provide useful synoptic views of catchment-wide erosion rates, little is known on the relative contributions of different sediment supply mechanisms to the acquisition of their initial signature in the headwaters. Here we address this issue by conducting a 10Be-budget of detrital materials that characterize the morphogenetic domains representative of high-altitude environments of the European Alps. We focus on the Etages catchment, located in the Ecrins-Pelvoux massif (southeast France), and illustrate how in situ 10Be concentrations can be used for tracing the origin of the sand fraction from the bedload in the trunk stream. The landscape of the Etages catchment is characterized by a geomorphic transient state, high topographic gradients, and a large variety of modern geomorphic domains ranging from glacial environments to scarcely vegetated alluvial plains. Beryllium-10 concentrations measured in the Etages catchment vary from similar to 1 x 104 to 4.5 x 105 atoms per gram quartz, while displaying consistent 10Be signatures within each representative morphogenetic unit. We show that the basic requirements for inferring catchment-wide denudation from 10Be concentration measurements are not satisfied in this small, dynamic catchment. However, the distinct 10Be signature observed for the geomorphic domains can be used as a tracer. We suggest that a terrestrial cosmogenic nuclide (TCN) budget approach provides a valuable tool for the tracing of material origin in basins where the let nature do the averaging' principles may be violated.

Snow shielding factors for cosmogenic nuclide dating inferred from long-term neutron detector monitoring

The depth-dependent attenuation of the secondary cosmic-ray particle flux due to snow cover and its effects on production rates of cosmogenic nuclides constitutes a potential source of uncertainty for studies conducted in regions characterized by frequent seasonal snow burial. Recent experimental and numerical modelling studies have yielded new constraints on the effect of hydrogen-rich media on the production rates of cosmogenic nuclides by low- and high-energy neutrons (<10(-3) MeV and >10(2) MeV, respectively). Here we present long-term neutron-detector monitoring data from a natural setting that we use to quantify the effect of snow cover on the attenuation of fast neutrons (0.1-10 MeV), which are responsible for the production of Ne-21 from Mg and Cl-36 from K. We use data measured between July 2001 and May 2008 at seven stations located throughout the Ecrins-Pelvoux massif (French Western Alps) and its surroundings, at elevations ranging from 200 to 2500 m a.s.l. From the cosmic-ray fluxes recorded during summer, when snow is absent, we infer an apparent attenuation length of 148 g cm(-2) in the atmosphere at a latitude of similar to 45 degrees N and for altitudes ranging from similar to 200 to 2500 m a.s.l. Using snow water-equivalent (SWE) values obtained through snow-coring campaigns that overlap in time the neutron monitoring for five stations, we show that fast neutrons are much more strongly attenuated in snow than predicted by a conventional mass-shielding formulation and the attenuation length estimated in the atmosphere. We suggest that such strong attenuation results from boundary effects at the atmosphere/snow interface induced by the high efficiency of water as a neutron moderator. Finally, we propose an empirical model that allows calculating snow-shielding correction factors as a function of SWE for studies using Ne-21 and Cl-36 analyses in Mg- and K-rich minerals, respectively. This empirical model is of interest for studies with a focus on cosmic-ray exposure dating, particularly if the target rocks are made up of mafic to ultramafic units where seasonal snow-cover is a common phenomenon.

XMapTools: a MATLAB©-based program for electron microprobe X-ray image processing and geothermobarometry

XMapTools is a MATLAB©-based graphical user interface program for electron microprobe X-ray image processing, which can be used to estimate the pressure–temperature conditions of crystallization of minerals in metamorphic rocks. This program (available online at http://www.xmaptools.com) provides a method to standardize raw electron microprobe data and includes functions to calculate the oxide weight percent compositions for various minerals. A set of external functions is provided to calculate structural formulae from the standardized analyses as well as to estimate pressure–temperature conditions of crystallization, using empirical and semi-empirical thermobarometers from the literature. Two graphical user interface modules, Chem2D and Triplot3D, are used to plot mineral compositions into binary and ternary diagrams. As an example, the software is used to study a high-pressure Himalayan eclogite sample from the Stak massif in Pakistan. The high-pressure paragenesis consisting of omphacite and garnet has been retrogressed to a symplectitic assemblage of amphibole, plagioclase and clinopyroxene. Mineral compositions corresponding to ~165,000 analyses yield estimates for the eclogitic pressure–temperature retrograde path from 25 kbar to 9 kbar. Corresponding pressure–temperature maps were plotted and used to interpret the link between the equilibrium conditions of crystallization and the symplectitic microstructures. This example illustrates the usefulness of XMapTools for studying variations of the chemical composition of minerals and for retrieving information on metamorphic conditions on a microscale, towards computation of continuous pressure–temperature-and relative time path in zoned metamorphic minerals not affected by post-crystallization diffusion.

Interpretation of fluid inclusions in quartz deformed by weak ductile shearing: reconstruction of differential stress magnitudes and pre-deformation fluid properties

A well developed theoretical framework is available in which paleofluid properties, such as chemical composition and density, can be reconstructed from fluid inclusions in minerals that have undergone no ductile deformation. The present study extends this framework to encompass fluid inclusions hosted by quartz that has undergone weak ductile deformation following fluid entrapment. Recent experiments have shown that such deformation causes inclusions to become dismembered into clusters of irregularly shaped relict inclusions surrounded by planar arrays of tiny, new-formed (neonate) inclusions. Comparison of the experimental samples with a naturally sheared quartz vein from Grimsel Pass, Aar Massif, Central Alps, Switzerland, reveals striking similarities. This strong concordance justifies applying the experimentally derived rules of fluid inclusion behaviour to nature. Thus, planar arrays of dismembered inclusions defining cleavage planes in quartz may be taken as diagnostic of small amounts of intracrystalline strain. Deformed inclusions preserve their pre-deformation concentration ratios of gases to electrolytes, but their H2O contents typically have changed. Morphologically intact inclusions, in contrast, preserve the pre-deformation composition and density of their originally trapped fluid. The orientation of the maximum principal compressive stress (σ1σ1) at the time of shear deformation can be derived from the pole to the cleavage plane within which the dismembered inclusions are aligned. Finally, the density of neonate inclusions is commensurate with the pressure value of σ1σ1 at the temperature and time of deformation. This last rule offers a means to estimate magnitudes of shear stresses from fluid inclusion studies. Application of this new paleopiezometer approach to the Grimsel vein yields a differential stress (σ1–σ3σ1–σ3) of ∼300 MPa∼300 MPa at View the MathML source390±30°C during late Miocene NNW–SSE orogenic shortening and regional uplift of the Aar Massif. This differential stress resulted in strain-hardening of the quartz at very low total strain (<5%<5%) while nearby shear zones were accommodating significant displacements. Further implementation of these experimentally derived rules should provide new insight into processes of fluid–rock interaction in the ductile regime within the Earth's crust.

P-T-t estimation of syn-kinematic strain in low-grade quartz-feldspar bearing rocks using thermodynamic modeling and 40Ar/39Ar dating techniques: example of the Plan-de-Phasy shear zone unit (Briançonnais Zone, Western Alps)

Pressure–Temperature–time (P–T–t) estimates of the syn-kinematic strain at the peak-pressure conditions reached during shallow underthrusting of the Briançonnais Zone in the Alpine subduction zone was made by thermodynamic modelling and 40Ar/39Ar dating in the Plan-de-Phasy unit (SE of the Pelvoux Massif, Western Alps). The dated phengite minerals crystallized syn-kinematically in a shear zone indicating top-to-the-N motion. By combining X-ray mapping with multi-equilibrium calculations, we estimate the phengite crystallization conditions at 270 ± 50 °C and 8.1 ± 2 kbar at an age of 45.9 ± 1.1 Ma. Combining this P–T–t estimate with data from the literature allows us to constrain the timing and geometry of Alpine continental subduction. We propose that the Briançonnais units were scalped on top of the slab during ongoing continental subduction and exhumed continuously until collision.

Minor inheritance inhibits the calibration of the 10Be production rate from the AD 1717 Val Ferret rock avalanche, European Alps

To calibrate the in situ 10Be production rate, we collected surface samples from nine large granitic boulders within the deposits of a rock avalanche that occurred in AD 1717 in the upper Ferret Valley, Mont Blanc Massif, Italy. The 10Be concentrations were extremely low and successfully measured within 10% analytical uncertainty or less. The concentrations vary from 4829 ± 448 to 5917 ± 476 at g−1. Using the historical age exposure time, we calculated the local and sea level-high latitude (i.e. ≥60°) cosmogenic 10Be spallogenic production rates. Depending on the scaling schemes, these vary between 4.60 ± 0.38 and 5.26 ± 0.43 at g−1 a−1. Although they correlate well with global values, our production rates are clearly higher than those from more recent calibration sites. We conclude that our 10Be production rate is a mean and an upper bound for production rates in the Massif region over the past 300 years. This rate is probably influenced by inheritance and will yield inaccurate (e.g. too young) exposure ages when applied to surface-exposure studies in the area. Other independently dated rock-avalanche deposits in the region that are approximately 103 years old could be considered as possible calibration sites.

The triplite–triploidite supergroup: structural modulation in wagnerite, discreditation of magniotriplite, and the new mineral hydroxylwagnerite

Electron-microprobe analysis, single-crystal X-ray diffraction with an area detector, and high-resolution transmission electron microscopy show that minerals related to wagnerite, triplite and triploidite, which are monoclinic Mg, Fe and Mn phosphates with the formula Me2+ 2PO4(F,OH), constitute a modulated series based on the average triplite structure. Modulation occurs along b and may be commensurate with (2b periodicity) or incommensurate but generally close to integer values (∼3b, ∼5b, ∼7b, ∼9b), i.e. close to polytypic behaviour. As a result, the Mg- and F-dominant minerals magniotriplite and wagnerite can no longer be considered polymorphs of Mg2PO4F, i.e., there is no basis for recognizing them as distinct species. Given that wagnerite has priority (1821 vs. 1951), the name magniotriplite should be discarded in favour of wagnerite. Hydroxylwagnerite, end-member Mg2PO4OH, occurs in pyrope megablasts along with talc, clinochlore, kyanite, rutile and secondary apatite in two samples from lenses of pyrope–kyanite–phengite–quartz-schist within metagranite in the coesite-bearing ultrahigh-pressure metamorphic unit of the Dora-Maira Massif, western Alps, Vallone di Gilba, Val Varaita, Piemonte, Italy. Electron microprobe analyses of holotype hydroxylwagnerite and of the crystal with the lowest F content gave in wt%: P2O5 44.14, 43.99; SiO2 0.28, 0.02; SO3 –, 0.01; TiO2 0.20, 0.16; Al2O3 0.06, 0.03; MgO 48.82, 49.12; FeO 0.33, 0.48; MnO 0.01, 0.02; CaO 0.12, 0.10; Na2O 0.01, –; F 5.58, 4.67; H2O (calc) 2.94, 3.36; –O = F 2.35, 1.97; Sum 100.14, 99.98, corresponding to (Mg1.954Fe0.007Ca0.003Ti0.004Al0.002Na0.001)Σ=1.971(P1.003Si0.008)Σ=1.011O4(OH0.526F0.474)Σ=1 and (Mg1.971Fe0.011Ca0.003Ti0.003Al0.001)Σ=1.989(P1.002Si0.001)Σ=1.003O4(OH0.603F0.397)Σ=1, respectively. Due to the paucity of material, H2O could not be measured, so OH was calculated from the deficit in F assuming stoichiometry, i.e., by assuming F + OH = 1 per formula unit. Holotype hydroxylwagnerite is optically biaxial (+), α 1.584(1), β 1.586(1), γ 1.587(1) (589 nm); 2V Z(meas.) = 43(2)°; orientation Y = b. Single-crystal X-ray diffraction gives monoclinic symmetry, space group P21/c, a = 9.646(3) Å, b = 12.7314(16) Å, c = 11.980(4) Å, β = 108.38(4) , V = 1396.2(8) Å3, Z = 16, i.e., hydroxylwagnerite is the OH-dominant analogue of wagnerite [β-Mg2PO4(OH)] and a high-pressure polymorph of althausite, holtedahlite, and α- and ε-Mg2PO4(OH). We suggest that the group of minerals related to wagnerite, triplite and triploidite constitutes a triplite–triploidite super-group that can be divided into F-dominant phosphates (triplite group), OH-dominant phosphates (triploidite group), O-dominant phosphates (staněkite group) and an OH-dominant arsenate (sarkinite). The distinction among the three groups and a potential fourth group is based only on chemical features, i.e., occupancy of anion or cation sites. The structures of these minerals are all based on the average triplite structure, with a modulation controlled by the ratio of Mg, Fe2+, Fe3+ and Mn2+ ionic radii to (O,OH,F) ionic radii.

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.