Heat from radioactive decay of nuclear waste disposal in deep geological repositories may thermally load the host rocks containing fractures and faults, potentially causing fault reactivation and induced seismicity. To investigate the effect of heating rate on stability of thermally induced fault slip, we conducted three experiments of thermally induced fault reactivation on a granite sample containing an inclined sawcut fracture under conditions of varying heating rates spanning two orders of magnitude. We find that the spatiotemporal evolution of temperature and the resultant thermal stress control the timing and size of induced stick‐slip events. A fast heating rate induces earlier but smaller stick‐slip events. In contrast, a slow heating rate promotes a uniform distribution of temperature and thermal stress, facilitating larger stick‐slip events that occur later and exhibit higher apparent frictional strength. Our findings highlight the importance of controlling heating rates in deep geological repositories to mitigate seismic risks.
The radiological impact of radionuclide transport via groundwater pathways at the Wolsong Low- and Intermediate-Level Waste (LILW) Disposal Center was estimated by considering site-specific characteristics, including hydrogeology, geochemistry, and land use. Human intrusion scenarios, such as groundwater well development, were analyzed to evaluate potential pumping volumes and radionuclide migration pathways. Particular attention was given to the hydrological and geochemical aspects of radionuclide transport, with a focus on local aquifer heterogeneity, flow dynamics, and interactions with engineered barriers and surrounding rock formations that delay radionuclide migration through sorption and other retention mechanisms. Sorption coefficients (Kd), calibrated using site-specific geochemical data, were incorporated to ensure realistic modeling of radionuclide behavior. A hierarchical approach integrating scenario screening, particle tracking techniques, and mass transfer modeling was employed. Numerical simulations using FEFLOW ver. 7.3 and GoldSim ver. 14.0 software provided insights into near-field and far-field transport phenomena under well pumping conditions. The results revealed distinct spatial flux behaviors, where carbon-14 (14C) dominated near-field flux due to its high inventory, while technetium-99 (99Tc) emerged as the primary dose contributor in the far-field flux, owing to its anionic nature and limited sorption capacity. Additionally, under high-pH conditions near concrete barriers, cellulose degradation into isosaccharinic acid was identified, enhancing radionuclide mobility through complex formation. These findings underscore the importance of site-specific sorption and speciation parameters in safety assessment and highlight the need for accurate geochemical modeling to optimize waste placement and ensure long-term disposal safety. The outcomes provide valuable insights for optimizing waste placement and contribute to the development of evidence-based safety strategies for long-term performance assessment.
Long-term safety assessment of deep geological repositories for spent nuclear fuel requires explicit evaluation of thermo-mechanical (TM) processes induced by decay heat and their influence on fractured host rock. A safety-relevant, though low-probability, scenario concerns shear reactivation of fractures intersecting deposition holes, which could compromise canister integrity if displacement exceeds design limits. This study presents a three-dimensional discrete element modelling approach to analyze the thermal evolution of the Forsmark repository (Sweden) and the associated long-term response of a discrete fracture network (DFN) during the post-closure phase. The model explicitly represents repository panel, deterministic deformation zones, and a stochastically generated fracture network embedded in a bonded particle assembly representing the rock for Particle Flow Code (PFC) numerical simulations. Time-dependent heat release from spent nuclear fuel canisters is implemented using a physically based decay power function. A deposition panel-scale heat-loading formulation accounts for deposition-hole and tunnel spacing. Two emplacement scenarios are analyzed: (a) a simultaneous all-panel heating scenario, used as a conservative bounding case, and (b) a sequential panel heating scenario representing staged emplacement and closure. The simulations show that temperature and thermally induced stress evolution are sensitive to the emplacement and closure sequence. Sequential heating produces a more gradual thermal build-up and lower peak temperatures than simultaneous heating, indicating that thermal and stress perturbations in the host rock can be influenced not only through repository design, but also by operational strategy. Thermally induced fracture shear displacement displays a systematic temporal response. Fractures located within the deposition panel footprint develop shear displacement rapidly during the early post-closure period, reaching peak values at approximately 200 years, followed by gradual relaxation as temperatures decline. The average peak shear displacement on fractures is on the order of 2–3 mm, while fractures outside the panel footprint show smaller early-time displacements and a more prolonged long-term response. All simulated shear displacements remain more than one order of magnitude below the commonly cited canister damage threshold for Forsmark of approximately 50 mm, even for the conservative simultaneous heating case. These results indicate that thermally induced fracture shear is unlikely to cause direct mechanical damage to canisters. At the same time, the persistence of residual shear displacement after heating implies permanent fracture dilation, which may influence long-term hydraulic properties and indirectly affect processes such as groundwater flow and canister corrosion. The modelling framework and results presented here were conducted for review purposes independently from the Swedish safety case, and provide a mechanistic basis for evaluating thermally induced fracture deformation in crystalline rock repositories and contribute to bounding the role of thermo-mechanical processes in the safety assessment of spent nuclear fuel disposal at Forsmark.
Identification of seismically active fault zones and the definition of sufficiently large respect distances from these faults which enable avoiding the damaged rock zone surrounding the ruptured ground commonly are amongst the first steps to take in the geoscientific evaluation of sites suitable for nuclear waste disposal. In this work we present a GIS-based approach, using the earthquake-epicentre locations from the instrumental earthquake record of South-Korea to identify potentially active fault zones in the country, and compare different strategies for fault zone buffer creation as originally developed for site search in the high seismicity country Japan, and the low-to-moderate seismicity countries Germany and Sweden. In order to characterize the hazard potential of the Korean fault zones, we moreover conducted slip tendency analysis, here for the first time covering the fault zones of the entire Korean Peninsula. For our analyses we used the geo-spatial information from a new version of the Geological map of South-Korea, containing the outlines of 11 rock units, which we simplified to distinguish between 4 different rock types (granites, metamorphic rocks, sedimentary rocks and igneous rocks) and the surface traces of 1,528 fault zones and 6,654 lineaments identified through years of field work and data processing, a rich geo-dataset which we will publish along with this manuscript. Our approach for identification of active fault zones was developed without prior knowledge of already known seismically active fault zones, and as a proof of concept the results later were compared to a map containing already identified active fault zones. The comparison revealed that our approach identified 16 of the 21 known seismically active faults and added 472 previously unknown potentially active faults. The 5 seismically active fault zones which were not identified by our approach are located in the NE- and SW-sectors of the Korean Peninsula, which haven’t seen much recent seismic activity, and thus are not sufficiently well covered by the seismic record. The strike directions of fault zones identified as active are in good agreement with the orientation of the current stress field of the peninsula and slip tendency analysis provided first insights into subsurface geometry such as the dip angles of both active and inactive fault zones. The results of our work are of major importance for the early-stage seismic hazard assessment that has to be conducted in support of the nuclear waste disposal siting in South-Korea. Moreover, the GIS-based methods for identification of active fault zones and buffering of respect areas around fault zone traces presented here, are applicable also elsewhere.
Identification of seismically active fault zones and the definition of sufficiently large respect distances from these faults which enable avoiding the damaged rock zone surrounding the ruptured ground commonly are amongst the first steps to take in the geoscientific evaluation of sites suitable for nuclear waste disposal. In this work we present a GIS-based approach, using the earthquake-epicentre locations from the instrumental earthquake record of South-Korea to identify potentially active fault zones in the country, and compare different strategies for fault zone buffer creation as originally developed for site search in the high seismicity country Japan, and the low-to-moderate seismicity countries Germany and Sweden. In order to characterize the hazard potential of the Korean fault zones, we moreover conducted slip tendency analysis, here for the first time covering the fault zones of the entire Korean Peninsula. For our analyses we used the geo-spatial information from a new version of the Geological map of South-Korea, containing the outlines of 11 rock units, which we simplified to distinguish between 4 different rock types (granites, metamorphic rocks, sedimentary rocks and igneous rocks) and the surface traces of 1,528 fault zones and 6,654 lineaments identified through years of field work and data processing, a rich geo-dataset which we will publish along with this manuscript. Our approach for identification of active fault zones was developed without prior knowledge of already known seismically active fault zones, and as a proof of concept the results later were compared to a map containing already identified active fault zones. The comparison revealed that our approach identified 16 of the 21 known seismically active faults and added 472 previously unknown potentially active faults. The 5 seismically active fault zones which were not identified by our approach are located in the NE- and SW-sectors of the Korean Peninsula, which haven’t seen much recent seismic activity, and thus are not sufficiently well covered by the seismic record. The strike directions of fault zones identified as active are in good agreement with the orientation of the current stress field of the peninsula and slip tendency analysis provided first insights into subsurface geometry such as the dip angles of both active and inactive fault zones. The results of our work are of major importance for the early-stage seismic hazard assessment that has to be conducted in support of the nuclear waste disposal siting in South-Korea. Moreover, the GIS-based methods for identification of active fault zones and buffering of respect areas around fault zone traces presented here, are applicable also elsewhere.
Identification of seismically active fault zones and the definition of sufficiently large respect distances from these faults which enable avoiding the damaged rock zone surrounding the ruptured ground commonly are amongst the first steps to take in the geoscientific evaluation of sites suitable for nuclear waste disposal. In this work we present a GIS-based approach, using the earthquake-epicentre locations from the instrumental earthquake record of South-Korea to identify potentially active fault zones in the country, and compare different strategies for fault zone buffer creation as originally developed for site search in the high seismicity country Japan, and the low-to-moderate seismicity countries Germany and Sweden. In order to characterize the hazard potential of the Korean fault zones, we moreover conducted slip tendency analysis, here for the first time covering the fault zones of the entire Korean Peninsula. For our analyses we used the geo-spatial information from a new version of the Geological map of South-Korea, containing the outlines of 11 rock units, which we simplified to distinguish between 4 different rock types (granites, metamorphic rocks, sedimentary rocks and igneous rocks) and the surface traces of 1,528 fault zones and 6,654 lineaments identified through years of field work and data processing, a rich geo-dataset which we will publish along with this manuscript. Our approach for identification of active fault zones was developed without prior knowledge of already known seismically active fault zones, and as a proof of concept the results later were compared to a map containing already identified active fault zones. The comparison revealed that our approach identified 16 of the 21 known seismically active faults and added 472 previously unknown potentially active faults. The 5 seismically active fault zones which were not identified by our approach are located in the NE- and SW-sectors of the Korean Peninsula, which haven’t seen much recent seismic activity, and thus are not sufficiently well covered by the seismic record. The strike directions of fault zones identified as active are in good agreement with the orientation of the current stress field of the peninsula and slip tendency analysis provided first insights into subsurface geometry such as the dip angles of both active and inactive fault zones. The results of our work are of major importance for the early-stage seismic hazard assessment that has to be conducted in support of the nuclear waste disposal siting in South-Korea. Moreover, the GIS-based methods for identification of active fault zones and buffering of respect areas around fault zone traces presented here, are applicable also elsewhere.