Abstract. Estimating the porosity of slates is of great interest for the industries dealing with sub-surface areas such as CO2 sequestration, nuclear waste disposal and shale gas but also for engineering purposes in terms of mechanical stability for underground or surface constructions. In this study, we aim at understanding estimates of the porosity of slates from the Infrahelvetic flysch units (IFUs) in the Glarus Alps (eastern Switzerland). Surface and sub-surface samples were collected along a temperature gradient from 200 to 320 °C and therefore give the opportunity to link pore types along this temperature and deformation path. In addition, we indicate which porosity is the effect of surface processes and indicate the contribution of artificially induced porosity. The developed workflow consists of a combination of bulk rock measurements including helium pycnometry (He pycnometry) and mercury intrusion porosimetry (MIP) with image analysis. Image analysis was performed with high-resolution scanning electron microscopy (SEM) on broad ion beam (BIB) prepared cross sections (BIB-SEM). Different vein generations provide evidence of porosity formation at depth, as they present paleo-porosity. Towards peak metamorphic conditions (prograde path), porosity reduces to < 1 vol%, indicated by matrix porosity detected by BIB-SEM. During exhumation (retrograde path) porosity increases due to the formation of microfractures interpreted as the effect of unloading (open fractures). At the surface, porosity is further increased due to the formation of macro-fractures (fracture apertures up to 1 mm), which are interpreted as being either due to the effect of weathering processes such as freeze and thaw cycles or artificially induced by sample preparation. Additionally, porosity and pore morphology are strongly dependent on mineralogy, sample homogeneity and strain, which change dynamically in time and space.