Topoclimatological case-study of Alpine pastures near the Albula Pass in the eastern Swiss Alps

Alpine grasslands are an important source of fodder for the cattle of Alpine farmers. Only during the short summer season can these pastures be used for grazing. With the anticipated climate change, it is likely that plant production – and thus the fodder basis for the cattle – will be influenced. Investigating the dependence of biomass production on topoclimatic factors will allow us to better understand how anticipated climate change may influence this traditional Alpine farming system. Because small-scale topoclimatological variations of the main meteorological variables: temperature, humidity, precipitation, shortwave incoming radiation and wind speed are not easily derived from available long-term climate stations in mountainous terrain, it was our goal to investigate the topoclimatic variations over the pastures belonging to the Alp Weissenstein research station north of the Albula Pass in the eastern Swiss Alps. We present a basic assessment of current topoclimatic conditions as a site characterization for ongoing ecological climate change studies. To be able to link short-term studies with long-term climate records, we related agrometeorological measurements with those of surrounding long-term sites run by MeteoSwiss, both on valley bottoms (Davos, Samedan), and on mountain tops (Weissfluhjoch, Piz Corvatsch). We found that the Davos climate station north of the study area is most closely correlated with the local climate of Alp Weissenstein, although a much closer site (Samedan) exists on the other side of the Albula Pass. Mountain top stations, however, did not provide a convincing approximation for the climate at Alp Weissenstein. Direct comparisons of near-surface measurements from a set of 11 small weather stations distributed over the domain where cattle and sheep are grazed indicate that nocturnal minimum air temperature and minimum vapor pressure deficit are mostly governed by the altitudinal gradient, whereas daily maxima – including also wind speed – are more strongly depending on vegetation cover and less on the altitude.

Observation errors in early historical upper-air observations

Upper-air observations are a fundamental data source for global atmospheric data products, but uncertainties, particularly in the early years, are not well known. Most of the early observations, which have now been digitized, are prone to a large variety of undocumented uncertainties (errors) that need to be quantified, e.g., for their assimilation in reanalysis projects. We apply a novel approach to estimate errors in upper-air temperature, geopotential height, and wind observations from the Comprehensive Historical Upper-Air Network for the time period from 1923 to 1966. We distinguish between random errors, biases, and a term that quantifies the representativity of the observations. The method is based on a comparison of neighboring observations and is hence independent of metadata, making it applicable to a wide scope of observational data sets. The estimated mean random errors for all observations within the study period are 1.5 K for air temperature, 1.3 hPa for pressure, 3.0 ms−1for wind speed, and 21.4° for wind direction. The estimates are compared to results of previous studies and analyzed with respect to their spatial and temporal variability.

Signatures of slope failures and river-delta collapses in a perialpine lake (Lake Lucerne, Switzerland)

Historical reports from the 17th Century document two destructive tsunamis with runups exceeding 5 m, affecting proximal basins of Lake Lucerne (Switzerland). One event in AD 1601 is coeval with a strong nearby earthquake (MW ca 5.9), which caused extensive slope failures in many parts of the lake. The second event in AD 1687 is associated with an apparently spontaneous partial collapse of the Muota river delta. This study combines high-resolution bathymetry, reflection seismic and lithological data in order to document the sedimentary and morphological signatures of the two subaqueous mass movements that probably generated the observed tsunamis. Such mass movements are significant as a common sedimentation process and as natural hazard in fjord-type lakes and similar environments. The deposits, covering large parts of the basins with thicknesses reaching >10 m, consist of two subunits: A lower ‘massflow deposit’ contains variably deformed sediments from the source areas. Its emplacement affected pre-existing sediments, incorporating thin sediment slices into the deposit and increasing its volume. Deep-reaching deformation near This is an Accepted Article that has been peer-reviewed and approved for publication in the Sedimentology, but has yet to undergo copy-editing and proof correction.

The history of subaquatic volcanism recorded in the sediments of Lake Kivu

Subaquatic volcanic activity has been ongoing in Lake Kivu since the early Holocene and has a dynamic effect on the biological productivity in the surface water, and the preservation of carbonate in the deep anoxic water. Groundwater discharge into the lake’s deepwater propels the upward advection of the water column that ultimately supplies nutrients to the surface water for biological production. The amount of nutrients supplied from the deepwater can be increased suddenly by (1) a cold meteorological event that drives deep seasonal mixing resulting in increased nutrients from below and oxygen from above, or (2) subaquatic volcanic activity that induces a buoyant hydrothermal plume, which entrains nutrients from the deepwater and results in anoxia or suboxic conditions in the surface water. Previous sedimentological studies in Lake Kivu have hypothesized that regional climatic changes are responsible for sudden changes in the preservation of carbonates in the Main Basin. Here we reveal that sublacustrine volcanic events most likely induce the abrupt changes to the geochemistry in the sediment in Lake Kivu. An unprecedented look into the sediment stratigraphy and geochemistry from high-resolution seismic-reflection, and 15N-isotope analyses was conducted in the Main Basin. The results reveal that buoyant hydrothermal plumes caused by subaquatic volcanic activity are a possible trigger for increased biological productivity and organic matter preservation, and that ongoing hydrothermal activity increases the alkalinity in the deepwater, leading to carbonate preservation. The onset of carbonate preservation since the 1970s that is currently observed in the sediment could indicate that hydrothermal discharge has recently increased in the lake.

Lake-level rise in the late Pleistocene and active subaquatic volcanism since the Holocene in Lake Kivu; East African Rift

The history of Lake Kivu is strongly linked to the activity of the Virunga volcanoes. Subaerial and subaquatic volcanoes, in addition to lake-level changes, shape the subaquatic morphologic and structural features in Lake Kivu's Main Basin. Previous studies revealed that volcanic eruptions blocked the former outlet of the lake to the north in the late Pleistocene, leading to a substantial rise in the lake level and subsequently the present- day thermohaline stratification. Additional studies have speculated that volcanic and seismic activities threaten to trigger a catastrophic release of the large amount of gases dissolved in the lake. The current study presents a bathymetric mapping and seismic profiling survey that covers the volcanically active area of the Main Basin at a resolution that is unprecedented for Lake Kivu. New geomorphologic features identified on the lake floor can accurately describe related lake-floor processes for the first time. The late Pleistocene lowstand is observed at 425 m depth, and volcanic cones, tuff rings, and lava flows observed above this level indicate both subaerial and subaquatic volcanic activities during the Holocene. The geomorphologic analysis yields new implications on the geologic processes that have shaped Lake Kivu's basin, and the presence of young volcanic features can be linked to the possibility of a lake overturn.

Integrated reconstruction of Holocene millennial-scale environmental changes in Tierra del Fuego, southernmost South America

This study presents new paleoenvironmental data obtained from sedimentary cores from Lago Fagnano, an elon- gated lake located at 54°S in southernmost South America. Data from palynomorphs (pollen, spores and algae) and associated palynofacies as well as from diatom taxa retrieved from these cores compared with other regional proxies contribute to evaluate the similarities and differences in the climate patterns based on different proxies from southernmost Patagonia. The pollen analysis reveals that a grass steppe environment existed during the early Holocene (11,300–~8000 cal a BP) followed by a major vegetation change characterized by development of forest-steppe ecotone communities between ~8000 and ~6500 cal a BP, under more humid conditions. Between ~ 6500 and ~ 4000 cal a BP, expansion and colonization by Nothofagus forests reflect an increase in effec- tive moisture levels, while openness in the forest communities characterizes the region after ~ 1100 cal a BP. The palynological organic matter combined with the algal content reflects hydrological changes occurring in the lake and its nutrient status, probably in close relation with past climate oscillations. All these past ecological changes are closely related to oscillations in precipitation and temperature as a response to the variations in the latitudinal position and/or strength of the Southern Westerlies wind belt during the Holocene.

Sediment dynamics in the subaquatic channel of the Rhone delta (Lake Geneva, France/Switzerland)

bstract With its smaller size, well-known boundary conditions, and the availability of detailed bathymetric data, Lake Geneva’s subaquatic canyon in the Rhone Delta is an excellent analogue to understand sedimentary pro- cesses in deep-water submarine channels. A multidisciplinary research effort was undertaken to unravel the sediment dynamics in the active canyon. This approach included innovative coring using the Russian MIR sub- mersibles, in situ geotechnical tests, and geophysical, sedimentological, geochemical and radiometric analysis techniques. The canyon floor/levee complex is character- ized by a classic turbiditic system with frequent spillover events. Sedimentary evolution in the active canyon is controlled by a complex interplay between erosion and sedimentation processes. In situ profiling of sediment strength in the upper layer was tested using a dynamic penetrometer and suggests that erosion is the governing mechanism in the proximal canyon floor while sedimen- tation dominates in the levee structure. Sedimentation rates progressively decrease down-channel along the levee structure, with accumulation exceeding 2.6 cm/year in the proximal levee. A decrease in the frequency of turbidites upwards along the canyon wall suggests a progressive confinement of the flow through time. The multi-proxy methodology has also enabled a qualitative slope-stability assessment in the levee structure. The rapid sediment loading, slope undercutting and over-steepening, and increased pore pressure due to high methane concentrations hint at a potential instability of the proximal levees. Fur- thermore, discrete sandy intervals show very high methane concentrations and low shear strength and thus could cor- respond to potentially weak layers prone to scarp failures.

Paleozoic-aged brine and authigenic helium preserved in an Ordovician shale aquiclude

Consideration of the geosphere for isolation of nuclear waste has generated substantial interest in the origin, age, and movement of fl uids and gases in low-permeability rock formations. Here, we present profi les of isotopes, solutes, and helium in porewaters recovered from 860 m of Cambrian to Devonian strata on the eastern fl ank of the Michigan Basin. Of particular interest is a 240-m-thick, halite-mineralized, Ordovician shale and carbonate aquiclude, which hosts Br–-enriched, post-dolomitic brine (5.8 molal Cl) originating as evaporated Silurian seawater. Authigenic helium that has been accumulating in the aquiclude for more than 260 m.y. is found to be isolated from underlying allochthonous, 3He-enriched helium that originated from the rifted base of the Michigan Basin and the Canadian Shield. The Paleozoic age and immobility of the pore fl uids in this Ordovician aquiclude considerably strengthen the safety case for deep geological repositories, but also provide new insights into the origin of deep crustal brines and opportunities for research on other components of a preserved Paleozoic porewater system.

Geochemical synthesis for the Effingen Member in boreholes at Oftringen, Gösgen and Küttigen

The Effingen Member is a low-permeability rock unit of Oxfordian age (ca. 160 Ma) that occurs across northern Switzerland. It comprises sandy calcareous marls and (argillaceous) limestones. This report describes the hydrogeochemistry, mineralogy and supporting physical properties of the Effingen Member in three boreholes in the Jura-Südfuss area: Oftringen, Gösgen and Küttigen, where it is 220–240 m thick. The top of the Effingen Member is at 420, 66 and 32 m depths at the three sites. Core materials are available from Oftringen and Gösgen, whereas information from Küttigen is limited to cuttings, in-situ hydrogeological testing and geophysical logging. Hydrogeological boundaries of the Effingen Member vary between locations. Ground-water flows were identified during drilling at the top (Geissberg Member), but not at the base, of the Effingen Member at Oftringen, at the base (Hauptrogenstein Formation) of the Effingen Member at Gösgen, and in a limestone layer (Gerstenhübel unit) within the Effingen Member at Küttigen. The marls and limestones of the Effingen Member have carbonate contents of 46–91 wt.-% and clay-mineral contents of 5–37 wt.-%. Pyrite contents are up to 1.6 wt.-%, but no sulphate minerals were detected by routine analyses. Clay minerals are predominantly mixed-layer illite-smectite, illite and kaolinite, with sporadic traces of chlorite and smectite. Veins filled with calcite ± celestite occur through the Effingen Member at Oftringen but not at Gösgen or Küttigen. They formed at 50–70 ºC from externally derived fluids, probably of Miocene age. Water contents are 0.7–4.2 wt.-%, corresponding to a water-loss porosity range of 1.9–10.8 vol.-%. Specific surface areas, measured by the BET method, are 2–30 m2/g, correlating with clay-mineral contents. Water activity has been measured and yielded surprisingly low values down to 0.8. These cannot be explained by pore-water salinity alone and include other effects, such as changes in the fabric due to stress release or partial saturation. Observed variations in measurements are not fully understood. Cation exchange capacity (CEC) and exchangeable cation populations have been studied by the Ni-en method. CEC, derived from the consumption of the index cation Ni, is 9–99 meq/kgrock at a solid:liquid ratio of 1, correlating with the clay-mineral content. Cation concentrations in Ni-en extract solutions are in the order Na+≥Ca2+>Mg2+>K+>Sr2+. However, the analytical results from the Ni-en extractions have additional contributions from cations originating from pore water and from mineral dissolution reactions that occurred during extraction, and it was not possible to reliably quantify these contributions. Therefore, in-situ cation populations and selectivity coefficients could not be derived. A suite of methods have been used for characterising the chemical compositions of pore waters in the Effingen Member. Advective displacement was used on one sample from each Oftringen and Gösgen and is the only method that produces results that approach complete hydrochemical compositions. Aqueous extraction was used on core samples from these two boreholes and gives data only for Cl- and, in some cases, Br-. Out-diffusion was used on core samples from Oftringen and similarly gives data for Cl- and Br- only. For both aqueous extraction and out-diffusion, reaction of the experimental water with rock affected concentrations of cations, SO42 and alkalinity in experimental solutions. Another method, centrifugation, failed to extract pore water. Stable isotope ratios (δ18O and δ2H) of pore waters in core samples from Oftringen were analysed by the diffusive exchange method and helium contents of pore water in Oftringen samples were extracted for mass spectrometric analysis by quantitative outgassing of preserved core samples. Several lines of evidence indicate that drillcore samples might not have been fully saturated when opened and subsampled in the laboratory. These include comparisons of water-loss porosities with physical porosities, water-activity measurements, and high contents of dissolved gas as inferred from ground-water samples. There is no clear proof of partial saturation and it is unclear whether this might represent in-situ conditions or is due to exsolution of gas due to the pressure release since drilling. Partial saturation would have no impact on the recalculation of pore-water compositions from aqueous extraction experiments using water-loss porosity data. The largest uncertainty in the pore-water Cl- concentrations recalculated from aqueous extraction and out-diffusion experiments is the magnitude of the anion-accessible fraction of water-loss porosity. General experience of clay-mineral rich formations suggests that the anion-accessible porosity fraction is very often about 0.5 and generally in a range of 0.3 to 0.6 and tends to be inversely correlated with clay-mineral contents. Comparisons of the Cl- concentration in pore water obtained by advective displacement with that recalculated from aqueous extraction of an adjacent core sample suggests a fraction of 0.27 for an Oftringen sample, whereas the same procedure for a Gösgen sample suggests a value of 0.64. The former value for anion-accessible porosity fraction is presumed to be unrepresentative given the local mineralogical heterogeneity at that depth. Through-diffusion experiments with HTO and 36Cl- suggest that the anion-accessible porosity fraction in the Effingen Member at Oftringen and Gösgen is around 0.5. This value is proposed as a typical average for rocks of the Effingen Member, bearing in mind that it varies on a local scale in response to the heterogeneity of lithology and pore-space architecture. The substantial uncertainties associated with the approaches to estimating anion-accessible porosity propagate into the calculated values of in-situ pore-water Cl- concentrations. On the basis of aqueous extraction experiments, and using an anion-accessible porosity fraction of 0.5, Cl- concentrations in the Effingen Member at Oftringen reach a maximum of about 14 g/L in the centre. Cl- decreases upwards and downwards from that, forming a curved depth profile. Cl- contents in the Effingen Member at Gösgen increase with depth from about 3.5 g/L to about 14 g/L at the base of the cored profile (which corresponds to the centre of the formation). Out-diffusion experiments were carried out on four samples from Oftringen, distributed through the Effingen Member. Recalculated Cl- concentrations are similar to those from aqueous extraction for 3 out of the 4 samples, and somewhat lower for one sample. Concentrations of other components, i.e. Na+, K+, Ca2+, Mg2+, Sr2+, SO42- and HCO3- cannot be obtained from the aqueous extraction and out-diffusion experimental data because of mineral dissolution and cation exchange reactions during the experiments. Pore-water pH also is not constrained by those extraction experiments. The only experimental approach to obtain complete pore-water compositions for samples from Oftringen and Gösgen is advective displacement of pore water. The sample from Oftringen used for this experiment is from 445 m depth in the upper part of the Effingen Member and gave eluate with 16.5 g/L Cl- whereas aqueous extraction from a nearby sample indicated about 9 g/L Cl-. The sample from Gösgen used for advective displacement is from 123 m depth in the centre of the Effingen Member sequence and gave eluate with about 9 g/L Cl- whereas aqueous extraction gave 11.5 g/L Cl-. In both cases the pore waters have Na-(Ca)-Cl compositions and SO42- concentrations of about 1.1 g/L. The Gösgen sample has a Br/Cl ratio similar to that of sea water, whereas this ratio is lower for the Oftringen sample. Taking account of uncertainties in the applied experimental approaches, it is reasonable to place an upper limit of ca. 20 g/L on Cl- concentration for pore water in the Effingen Member in this area. There are major discrepancies between pore-water SO42- concentrations inferred from aqueous extraction or out-diffusion experiments and those obtained from advective displacement in both the Oftringen and Gösgen cases. A general conclusion is that all or at least part of the discrepancies are attributable to perturbation of the sulphur system and enhancement of SO42- by sulphate mineral dissolution and possibly minor pyrite oxidation during aqueous extraction and out-diffusion. Therefore, data for SO42- calculated from those pore-water sampling methods are considered not to be representative of in-situ conditions. A reference pore-water composition was defined for the Effingen Member in the Jura Südfuss area. It represents the probable upper limits of Cl- contents and corresponding anion and cation concentrations that are reasonably constrained by experimental data. Except for Cl- and possibly Na+ concentrations, this composition is poorly constrained especially with respect to SO42- and Ca2+ concentrations, and pH and alkalinity. Stable isotope compositions, δ18O and δ2H, of pore waters in the Effingen Member at Oftringen plot to the right of the meteoric water line, suggesting that 18O has been enriched by water-rock exchange, which indicates that the pore waters have a long residence time. A long residence time of pore water is supported by the level of dissolved 4He that has accumulated in pore water of the Effingen Member at Oftringen. This is comparable with, or slightly higher than, the amounts of 4He in the Opalinus Clay at Benken. Ground waters were sampled from flowing zones intersected by boreholes at the three locations. The general interpretation is that pore waters and ground-water solutes may have similar origins in Mesozoic and Cenozoic brackish-marine formations waters, but ground-water solutes have been diluted rather more than pore waters by ingress of Tertiary and Quaternary meteoric waters. The available hydrochemical data for pore waters from the Effingen Member at these three locations in the Jura-Südfuss area suggest that the geochemical system evolved slowly over geological periods of time, in which diffusion was an important mechanism of solute transport. The irregularity of Cl- and δ18O profiles and spatial variability of advective ground-water flows in the Malm-Dogger system suggests that palaeohydrogeological and hydrochemical responses to changing tectonic and surface environmental conditions were complex.

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.

Characterization of granite matrix porosity and pore-space geometry by in situ and laboratory methods

Most available studies of interconnected matrix porosity of crystalline rocks are based on laboratory investigations; that is, work on samples that have undergone stress relaxation and were affected by drilling and sample preparation. The extrapolation of the results to in situ conditions is therefore associated with considerable uncertainty, and this was the motivation to conduct the ‘in situ Connected Porosity’ experiment at the Grimsel Test Site (Central Swiss Alps). An acrylic resin doped with fluorescent agents was used to impregnate the microporous granitic matrix in situ around an injection borehole, and samples were obtained by overcoring. The 3-D structure of the porespace, represented by microcracks, was studied by U-stage fluorescence microscopy. Petrophysical methods, including the determination of porosity, permeability and P -wave velocity, were also applied. Investigations were conducted both on samples that were impregnated in situ and on non-impregnated samples, so that natural features could be distinguished from artefacts. The investigated deformed granites display complex microcrack populations representing a polyphase deformation at varying conditions. The crack population is dominated by open cleavage cracks in mica and grain boundary cracks. The porosity of non-impregnated samples lies slightly above 1 per cent, which is 2–2.5 times higher than the in situ porosity obtained for impregnated samples. Measurements of seismic velocities (Vp ) on spherical rock samples as a function of confining pressure, spatial direction and water saturation for both non-impregnated and impregnated samples provide further constraints on the distinction between natural and induced crack types. The main conclusions are that (1) an interconnected network of microcracks exists in the whole granitic matrix, irrespective of the distance to ductile and brittle shear zones, and (2) conventional laboratory methods overestimate the matrix porosity. Calculations of contaminant transport through fractured media often rely on matrix diffusion as a retardation mechanism.