Safety assessments of deep geological repositories for high level radioactive waste require robust quantitative and qualitative characterization of future repository environments over time horizons up to one million years. In Germany, this period spans roughly ten glacial and interglacial cycles, during which ice sheet advance, permafrost formation, and thawing can repeatedly alter the thermal and hydraulic regime of the overburden from near surface to intermediate depths. To assess radionuclide transport during glaciation, parameters for process based numerical predictions are required and must be obtained from laboratory measurements under controlled conditions. An improved rock volumetric deformation testing system was used to measure freezing induced deformation and hydraulic response. System and coolant effects were calibrated using stainless steel. Freeze thaw tests on consolidated and unconsolidated materials yielded temperature dependent deformation behavior, water migration characteristics during freezing, and permeability evolution with freeze thaw cycle number. These experimentally derived relations were used to constrain a coupled thermo hydraulic numerical framework implemented via a FEFLOW plugin. The model resolves liquid water and ice phase change in saturated porous media, incorporates nonlinear temperature dependent hydraulic conductivity and effective thermal properties, and represents lateral groundwater flow driven by hydraulic gradients. Permeability levels and their freeze thaw induced evolution were parameterized using laboratory results, and simulated stagewise redistribution of water during cooling was evaluated against measured inlet and outlet flux signatures. The calibrated framework quantified how advective heat transport associated with lateral flow controls the timing, persistence, and spatial heterogeneity of the freezing front.