Assessing Initial Conditions for Chloride Transport Across Low- permeability Argillaceous Rocks, Wellenberg, Switzerland

Information about fluid evolution and solute transport in a low-permeability metamorphic rock sequence has been obtained by comparing chloride concentrations and chlorine isotope ratios of pore water, groundwater, and fluid inclusions. The similarity of d37Cl values in fluid inclusions and groundwater suggests a closed-system evolution during the metamorphic overprint, and signatures established at this time appear to form the initial conditions for chloride transport after exhumation of the rock sequence.

Reconstruction of in-situ porosity and porewater compositions of low-permeability crystalline rocks: Magnitude of artefacts induced by drilling and sample recovery

Geological site characterisation programmes typically rely on drill cores for direct information on subsurface rocks. However, porosity, transport properties and porewater composition measured on drill cores can deviate from in-situ values due to two main artefacts caused by drilling and sample recovery: (1) mechanical disruption that increases porosity and (2) contamination of the porewater by drilling fluid. We investigated the effect and magnitude of these perturbations on large drill core samples (12–20 cm long, 5 cmdiameter) of high-grade, granitic gneisses obtained from 350 to 600 m depth in a borehole on Olkiluoto Island (SW Finland). The drilling fluid was traced with sodium–iodide. By combining out-diffusion experiments, gravimetry, UV-microscopy and iodide mass balance calculations, we successfully quantified the magnitudes of the artefacts: 2–6% increase in porosity relative to the bulk connected porosity and 0.9 to 8.9 vol.% contamination by drilling fluid. The spatial distribution of the drilling-induced perturbations was revealed by numerical simulations of 2D diffusion matched to the experimental data. This showed that the rims of the samples have a mechanically disrupted zone 0.04 to 0.22 cm wide, characterised by faster transport properties compared to the undisturbed centre (1.8 to 7.7 times higher pore diffusion coefficient). Chemical contamination was shown to affect an even wider zone in all samples, ranging from 0.15 to 0.60 cm, inwhich iodide enrichmentwas up to 180 mg/kgwater, compared to 0.5 mg/kgwater in the uncontaminated centre. For all samples in the present case study, it turned out that the magnitude of the artefacts caused by drilling and sample recovery is so small that no correction is required for their effects. Therefore, the standard laboratory measurements of porosity, transport properties and porewater composition can be taken as valid in-situ estimates. However, it is clear that the magnitudes strongly depend on site- and drilling-specific factors and therefore our results cannot be transferred simply to other locations. We recommend the approach presented in this study as a route to obtain reliable values in future drilling campaigns aimed at characterising in-situ bedrock properties.

Characterisation of Gas Transport Properties of the Opalinus Clay, a Potential Host Rock Formation for Radioactive Waste Disposal

The Opalinus Clay in Northern Switzerland has been identified as a potential host rock formation for the disposal of radioactive waste. Comprehensive understanding of gas transport processes through this low-permeability formation forms a key issue in the assessment of repository performance. Field investigations and laboratory experiments suggest an intrinsic permeability of the Opalinus Clay in the order of 10(-20) to 10(-21) m(2) and a moderate anisotropy ratio < 10. Porosity depends on clay content and burial depth; values of similar to 0.12 are reported for the region of interest. Porosimetry indicates that about 10-30 of voids can be classed as macropores, corresponding to an equivalent pore radius > 25 nm. The determined entry pressures are in the range of 0.4-10 MPa and exhibit a marked dependence on intrinsic permeability. Both in situ gas tests and gas permeameter tests on drillcores demonstrate that gas transport through the rock is accompanied by porewater displacement, suggesting that classical flow concepts of immiscible displacement in porous media can be applied when the gas entry pressure (i.e. capillary threshold pressure) is less than the minimum principal stress acting within the rock. Essentially, the pore space accessible to gas flow is restricted to the network of connected macropores, which implies a very low degree of desaturation of the rock during the gas imbibition process. At elevated gas pressures (i.e. when gas pressure approaches the level of total stress that acts on the rock body), evidence was seen for dilatancy controlled gas transport mechanisms. Further field experiments were aimed at creating extended tensile fractures with high fracture transmissivity (hydro- or gasfracs). The test results lead to the conclusion that gas fracturing can be largely ruled out as a risk for post-closure repository performance.

Modeling immune functions of the mouse blood-cerebrospinal fluid barrier in vitro: primary rather than immortalized mouse choroid plexus epithelial cells are suited to study immune cell migration across this brain barrier.

BACKGROUND The blood-cerebrospinal fluid barrier (BCSFB) established by the choroid plexus (CP) epithelium has been recognized as a potential entry site of immune cells into the central nervous system during immunosurveillance and neuroinflammation. The location of the choroid plexus impedes in vivo analysis of immune cell trafficking across the BCSFB. Thus, research on cellular and molecular mechanisms of immune cell migration across the BCSFB is largely limited to in vitro models. In addition to forming contact-inhibited epithelial monolayers that express adhesion molecules, the optimal in vitro model must establish a tight permeability barrier as this influences immune cell diapedesis. METHODS We compared cell line models of the mouse BCSFB derived from the Immortomouse(®) and the ECPC4 line to primary mouse choroid plexus epithelial cell (pmCPEC) cultures for their ability to establish differentiated and tight in vitro models of the BCSFB. RESULTS We found that inducible cell line models established from the Immortomouse(®) or the ECPC4 tumor cell line did not express characteristic epithelial proteins such as cytokeratin and E-cadherin and failed to reproducibly establish contact-inhibited epithelial monolayers that formed a tight permeability barrier. In contrast, cultures of highly-purified pmCPECs expressed cytokeratin and displayed mature BCSFB characteristic junctional complexes as visualized by the junctional localization of E-cadherin, β-catenin and claudins-1, -2, -3 and -11. pmCPECs formed a tight barrier with low permeability and high electrical resistance. When grown in inverted filter cultures, pmCPECs were suitable to study T cell migration from the basolateral to the apical side of the BCSFB, thus correctly modelling in vivo migration of immune cells from the blood to the CSF. CONCLUSIONS Our study excludes inducible and tumor cell line mouse models as suitable to study immune functions of the BCSFB in vitro. Rather, we introduce here an in vitro inverted filter model of the primary mouse BCSFB suited to study the cellular and molecular mechanisms mediating immune cell migration across the BCSFB during immunosurveillance and neuroinflammation.

Comparison of immortalized bEnd5 and primary mouse brain microvascular endothelial cells as in vitro blood-brain barrier models for the study of T cell extravasation

Important insights into the molecular mechanism of T cell extravasation across the blood-brain barrier (BBB) have already been obtained using immortalized mouse brain endothelioma cell lines (bEnd). However, compared with bEnd, primary brain endothelial cells have been shown to establish better barrier characteristics, including complex tight junctions and low permeability. In this study, we asked whether bEnd5 and primary mouse brain microvascular endothelial cells (pMBMECs) were equally suited as in vitro models with which to study the cellular and molecular mechanisms of T cell extravasation across the BBB. We found that both in vitro BBB models equally supported both T cell adhesion under static and physiologic flow conditions, and T cell crawling on the endothelial surface against the direction of flow. In contrast, distances of T cell crawling on pMBMECs were strikingly longer than on bEnd5, whereas diapedesis of T cells across pMBMECs was dramatically reduced compared with bEnd5. Thus, both in vitro BBB models are suited to study T cell adhesion. However, because pMBMECs better reflect endothelial BBB specialization in vivo, we propose that more reliable information about the cellular and molecular mechanisms of T cell diapedesis across the BBB can be attained using pMBMECs.

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.

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.

(234)U/(238)U activity ratio disequilibrium technique for studying uranium mobility in the Opalinus Clay at Mont Terri, Switzerland

Mobility of naturally occurring U-238 and U-234 radionuclides was studied in a low permeability, reducing claystone formation (Opalinus Clay) near its contact with an overlying oxidising aquifer (Dogger Limestones) at Mont Terri, Switzerland. Our data point to a limited redistribution of U in some of the studied samples. Observed centimetre-scale U mobility is explained by slow diffusive transport of U-234 in the pore waters of the Opalinus Clay driven by spatially variable in situ supply (by alpha-recoil) of U-234 from the rock matrix. Metre-scale mobility is interpreted as a result of infiltration of meteoric water into the overlying aquifer which developed gradients of U concentration across the two rock formations. This triggered a slow in-diffusion of U with (U-234/U-238) > 1 into the Opalinus Clay as attested by a clear-cut pattern of decreasing bulk rock (U-234/U-238) inwards the Opalinus Clay, away from the Dogger Limestones.

From transient to steady state deformation and grain size: A thermodynamic approach using elasto-visco-plastic numerical modeling

Numerical simulation experiments give insight into the evolving energy partitioning during high-strain torsion experiments of calcite. Our numerical experiments are designed to derive a generic macroscopic grain size sensitive flow law capable of describing the full evolution from the transient regime to steady state. The transient regime is crucial for understanding the importance of micro structural processes that may lead to strain localization phenomena in deforming materials. This is particularly important in geological and geodynamic applications where the phenomenon of strain localization happens outside the time frame that can be observed under controlled laboratory conditions. Ourmethod is based on an extension of the paleowattmeter approach to the transient regime. We add an empirical hardening law using the Ramberg-Osgood approximation and assess the experiments by an evolution test function of stored over dissipated energy (lambda factor). Parameter studies of, strain hardening, dislocation creep parameter, strain rates, temperature, and lambda factor as well asmesh sensitivity are presented to explore the sensitivity of the newly derived transient/steady state flow law. Our analysis can be seen as one of the first steps in a hybrid computational-laboratory-field modeling workflow. The analysis could be improved through independent verifications by thermographic analysis in physical laboratory experiments to independently assess lambda factor evolution under laboratory conditions.

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.

Matrix Porewater In Crystalline Rocks: Extraction and Analysis

The chemical and isotopic characterization of porewater residing in the inter- and intragranular pore space of the low-permeability rock matrix is an important component with respect to the site characterization and safety assessment of potential host rocks for a radioactive waste disposal. The chemical and isotopic composition of porewater in such low permeability rocks has to be derived by indirect extraction techniques applied to naturally saturated rock material. In most of such indirect extraction techniques – especially in case of rocks of a porosity below about 2 vol.% – the original porewater concentrations are diluted and need to be back-calculated to in-situ concentrations. This requires a well-defined value for the connected porosity – accessible to different solutes under in-situ conditions. The derivation of such porosity values, as well as solute concentrations, is subject to various perturbations during drilling, core sampling, storage and experiments in the laboratory. The present study aims to demonstrate the feasibility of a variety of these techniques to charac-terize porewater and solute transport in crystalline rocks. The methods, which have been de-veloped during multiple porewater studies in crystalline environments, were applied on four core samples from the deep borehole DH-GAP04, drilled in the Kangerlussuaq area, Southwest Greenland, as part of the joint NWMO–Posiva–SKB Greenland Analogue Project (GAP). Potential artefacts that can influence the estimation of in situ porewater chemistry and isotopes, as well as their controls, are described in detail in this report, using specific examples from borehole DH-GAP04

Time variable Earth’s gravity field from SLR satellites

The time variable Earth’s gravity field contains information about the mass transport within the system Earth, i.e., the relationship between mass variations in the atmosphere, oceans, land hydrology, and ice sheets. For many years, satellite laser ranging (SLR) observations to geodetic satellites have provided valuable information of the low-degree coefficients of the Earth’s gravity field. Today, the Gravity Recovery and Climate Experiment (GRACE) mission is the major source of information for the time variable field of a high spatial resolution. We recover the low-degree coefficients of the time variable Earth’s gravity field using SLR observations up to nine geodetic satellites: LAGEOS-1, LAGEOS-2, Starlette, Stella, AJISAI, LARES, Larets, BLITS, and Beacon-C. We estimate monthly gravity field coefficients up to degree and order 10/10 for the time span 2003–2013 and we compare the results with the GRACE-derived gravity field coefficients. We show that not only degree-2 gravity field coefficients can be well determined from SLR, but also other coefficients up to degree 10 using the combination of short 1-day arcs for low orbiting satellites and 10-day arcs for LAGEOS-1/2. In this way, LAGEOS-1/2 allow recovering zonal terms, which are associated with long-term satellite orbit perturbations, whereas the tesseral and sectorial terms benefit most from low orbiting satellites, whose orbit modeling deficiencies are minimized due to short 1-day arcs. The amplitudes of the annual signal in the low-degree gravity field coefficients derived from SLR agree with GRACE K-band results at a level of 77 %. This implies that SLR has a great potential to fill the gap between the current GRACE and the future GRACE Follow-On mission for recovering of the seasonal variations and secular trends of the longest wavelengths in gravity field, which are associated with the large-scale mass transport in the system Earth.

Constraining porewater chemistry in a 250m thick argillaceous rock sequence

The geochemistry of an argillaceous rock sequence from a deep borehole in NE-Switzerland was investigated. The focus was to constrain the porewater chemistry in low permeability Jurassic rocks comprising the Liassic, the Opalinus Clay formation, the 'Brown Dogger' unit and the Effingen Member (Malm). A multi-method approach including mineralogical analysis, aqueous and Ni-ethylenediamine extraction, squeezing tests and pCO(2) measurements as well as geochemical modelling was applied for this purpose. A consistent dataset was obtained with regard to the main solutes in the porewaters. A fairly constant anion-accessible porosity of similar to 50% of the total porosity was deduced for all analysed samples which displayed variable clay-mineral contents. Sulphate concentrations were shown to be constrained by a sulphate-bearing phase, presumably by celestite or a Sr-Ba sulphate. Application of a simple equilibrium model, including cation exchange reactions, calcite and celestite equilibrium showed good agreement with squeezing data, indicating the suitability of the modelling approach to simulate porewater chemistry in the studied argillaceous rocks. The modelling highlighted the importance of correct determination of the exchangeable cation population. The analysis corroborates that squeezing of the studied rocks is a viable and efficient way to sample porewater.