To assess the robustness of a safety case for a deep geological repository (DGR), it is necessary to analyze a range of scenarios covering likely, less likely, and hypothetical future developments. Crystalline rock can, under ideal conditions, provide a suitable hydrogeologic barrier due to its extremely low matrix permeability. However, this host rock is often fractured, which can compromise its hydro-mechanical (HM) barrier function. We quantify how faults that are prone to reactivation during glacial events can affect radionuclide migration around a DGR in a crystalline host rock. We extend a previously developed finite element model of coupled fluid flow and radionuclide transport to numerically solve the component transport problem before and after fault reactivation. Assuming that fault reactivation is triggered by changes in mechanical boundary conditions, we derive heterogeneous permeability distributions in the reactivated faults by evaluating the Coulomb failure stress criterion of finite element solutions of a complementary hydro-mechanical problem. Specifically, we evaluate the consequences of glacial isostatic adjustment (GIA) during a glacial cycle. We find that the increased permeability in the reactivated faults accelerates the migration of radionuclides along the fault by channeling the flow, while it is reduced in the direction perpendicular to the fault. The channeling observed is also a result of heterogeneous permeability enhancement, and the flow fields differ from those of the previous model which postulated a homogeneous permeability enhancement. Although the proposed numerical workflow has been applied to the case of GIA, it is adaptable to study hydro-mechanical processes induced by seismic events or by hydrofracking in enhanced geothermal systems.
