Clay rocks such as the Opalinus Clay (OPA) in Switzerland are considered as potential host rocks for deep geological disposal of nuclear wastes. Radionuclide migration in clay rocks is dominated by molecular diffusion highly dependent on the pore network geometries. In the sandy facies of OPA (SF-OPA), this pore network is critically modified due to compositional variability and owing to diagenetic reaction products, for example, carbonate and sulfide minerals (Philipp et al., 2017). Such spatial variability is responsible for heterogeneous diffusion patterns (Kulenkampff et al., 2015). At the continuum scale, the representative elementary volume (REV) is an important parameter for understanding and modeling the reactive transport (RT) of contaminants and radionuclides. An accurate estimation of the REV for diffusion is decisive for a meaningful continuum-scale simulation of radionuclide migration in heterogeneous clay rocks like SF-OPA, thus playing an important role in the safety analysis of deep geological repositories for nuclear waste.In continuum mechanics for a heterogeneous material, the REV is the smallest volume over which a measurement can be made to produce a representative value of the macroscopic property in three-dimensional space (3D) (Bear Abstract Accurate modeling and simulation of radionuclide migration in clay rocks such as the Opalinus Clay (OPA) play a key role in the safety assessment of deep geological repositories for nuclear wastes. At the continuum scale, the representative elementary volume (REV) is a fundamental constraint to quantify the effective diffusivity, which is a key parameter in reactive transport (RT) models. Therefore, an accurate estimation of the REV is essential for a meaningful continuum-scale RT simulation in heterogeneous clay rocks. This study presents a comprehensive analysis of the heterogeneities of porosity and effective diffusivity in clay rocks by using the classical sampling theory and pore-scale simulations. First, in this study, the two-dimensional representative elementary area is correlated with the REV for porosity via a characteristic length. Next, it is shown that the REV for diffusivity is larger than the REV for porosity. Moreover, these two REVs can be correlated using Archie's law. In such a way, the REV for diffusivity can be determined by the developed correlations through analyzing two-dimensional microstructures, thus significantly reducing the computational cost. Finally, the applicability of our approach for clay rocks is validated by experimental data on the diffusion of tritiated water in the heterogeneous sandy facies of OPA. From both the experimental data and the modeling prediction, the REV for diffusivity in the sandy facies of OPA is in the order of cubic centimeters. This study provides critical insights into the diffusion in heterogeneous clay rocks toward an enhanced predictability of radionuclide migration.Plain Language Summary Contaminant migration in clay rocks is dominated by molecular diffusion due to their low permeability....