The safe disposal of high-level radioactive waste (HLW), including the discharged spent nuclear fuel (SNF) and contaminated by-products produced from relevant chemical treatments, has become a serious pending problem for numerous countries that operate the nuclear power plants. The deep geological disposal (DGD) has thus far been considered the most proven and viable solution for isolation of the HLW and preventing any significant release of radionuclides into the biosphere. The DGD system consists of the multiple engineered and natural barrier components. Among them, the montmorillonite-based buffer and tunnel backfills are designed to perform the two major geochemical functions: 1) preventing the ingress of groundwater and any chemicals that compromise the safety of waste canister and 2) retarding the migration of released radionuclides by providing sufficient chemisorption sites. Therefore, it is essential to investigate the sorption mechanism of radionuclides onto montmorillonite and develop a thermodynamic reaction model in advance in order to accurately predict the long-term performance of engineered barriers and to reduce the uncertainties in the safety assessment of a deep geological repository (DGR) ultimately; thus far, sorption of chemical species onto mineral adsorbents has been widely described based on the concept of sorption-desorption distribution coefficient (Kd), the value of which is intrinsically conditional, and active scientific efforts have been made to develop robust thermodynamic sorption models which offer the potential to improve confidence in demonstration of radionuclide migration under a wide range of geochemical conditions. The natural montmorillonites are generally classified into Na-type or Ca-type according to its exchangeable cation, and the Ca-montmorillonite containing clays are being considered as candidate materials for the engineered barriers of DGR in several countries; they generally have advantages of higher thermal conductivity and lower price than the Na-montmorillonite based clays, but their sorption capacities are still comparable. In this framework, we aimed to investigate the chemical interactions of Ca-montmorillonite with selenite [Se(IV)], which is a major oxyanionic species in terms of HLW disposal, and develop a reliable thermodynamic sorption model (TSM). The present work summarizes the characterization of Ca-montmorillonite separated from the newly adopted reference bentonite (Bentonil-WRK) by means of XRD, BET, FTIR, CEC measurement, and acid-base titration. Further, its sorption behaviors with aqueous selenite species under aqueous conditions of S/L = 5 g/L, I = 0.01-0.1 m CaCl2, pH = 4.5-8.5, pCO2 = 10-3.5 atm, and T = 25°C were examined, and the resulting thermodynamic data are discussed as well.