In instances of damage to engineered barriers containing nuclear waste material, surrounding bedrock is a natural barrier that retards radionuclide movement by way of adsorption and delay due to groundwater flow through highly tortuous fractured rock pathways. At the Gyeongju nuclear waste disposal site, groundwater mainly flows through granitic and sedimentary rock fractures. Therefore, to understand the nuclide migration path, it is necessary to understand discrete fracture networks based on heterogeneous fracture orientations, densities, and size characteristics. In this study, detailed heterogeneous fracture distribution, including the density and orientation of the fractures, was considered for a region that has undergone long periods of change from various geological activities at and around the Gyeongju site. A site-scale discrete fracture network (DFN) model was constructed taking into account: (i) regional fracture heterogeneity constrained by a multiple linear regression analysis of fracture intensity on faults and electrical resistivity; and (ii) the connectivity of conductive fractures having fracture hydraulic parameters, using transient flow simulation. Geometric and hydraulic heterogeneity of the DFN was upscaled into equivalent porous media for flow and transport simulation for a large-scale model.

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 (K_d), 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 (¹⁴C) dominated near-field flux due to its high inventory, while technetium-99 (⁹⁹Tc) 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.

This study examines the hydrogeological response to seismic activity in the southeastern region of Korea, focusing on groundwater-level changes induced by the 2016 Gyeongju earthquake (ML 5.8). Using a 2D hydro-mechanical coupled bonded particle model, we simulated the dynamic fault rupture process to analyze stress redistribution and its impact on pore pressure and groundwater levels. The results indicate that compressional areas correlate strongly with pore pressure increases and groundwater-level rises, while extensional areas show decreases in both. Observations from the groundwater monitoring well at Gyeongju Cheonbuk, located approximately 10 km north of the earthquake’s epicenter, align well with the model's predictions, providing validation for the simulation. These findings highlight the capability of hydro-mechanical models to capture fault-induced hydrological responses and offer valuable insights into the interplay between seismic activity and groundwater systems.

This study examines the hydrogeological response to the 12 September 2016 Gyeong-Ju earthquake (ML 5.8) in the southeastern region of the Korean Peninsula. Using 2D hydro-mechanical coupled bonded particle modeling, we simulated the dynamic fault rupture process to analyze stress redistribution and its impact on pore pressure and groundwater levels (GWLs). The results indicated that compressional areas correlated strongly with pore pressure increases and GWL rises, while extensional areas showed decreases in both. Observations from the groundwater monitoring Well 5 at Gyeong-Ju San-Nae and Well 8 at Gyeong-Ju Cheon-Buk, located approximately 15 km from the earthquake’s epicenter, aligned well with the model’s predictions and interpretation, providing validation for the simulation. These findings highlight the capability of hydro-mechanical models to capture fault-induced hydrological responses and offer valuable insights into the interplay between seismic activity and groundwater systems.