Crystalline rocks such as granites have been investigated as potential host rocks for the geological disposal of radioactive waste in many countries (e.g., EUR, 2005; NEA, 2012; NASEM, 2015). Radionuclide (RN) transport in fractured crystalline rocks can be conceptualized using a dual-porosity model where RNs are transported by advective water flow through a fracture and are retarded by diffusion and sorption into the surrounding rock matrix (Neretnieks, 1980). RNs can then diffuse from the fracture surfaces into the pore network of the fracture fillings and rock matrix and sorb onto pore walls. RN sorption on fracture surfaces and diffusion into the adjacent rock matrix are a key process influencing the safety of the geological disposal. To develop a realistic model and reliable parameters for long-term safety assessments, it is necessary to understand and quantify RN diffusion and sorption processes in the heterogeneous rock matrix and fracture systems.The different types and scales of heterogeneities that must also be considered for RN transport in fractured crystalline rocks include (a) heterogeneous distribution of minerals and pores in the rock matrix, (b) heterogeneity in mineral and pore distribution near fracture surfaces (e.g., fracture fillings and altered zones near fracture surfaces), and (c) heterogeneous flow distribution in the complex channel structures along fracture openings. As discussed in detail in Tachi et al. (2018), the first heterogeneity (a) in a rock matrix affects sorption and diffusion at a relatively smaller scale. The domain pore spaces between mineral grains and intragranular secondary pores form water-filled, heterogeneous networks with varying diffusion-related geometric factors, such as tortuosity and constrictivity, and accessible reactive surfaces related to sorption (e.