Thermodynamic sorption modeling can enhance confidence in assessing and demonstrating the radionuclide sorption phenomena onto various mineral adsorbents. In this work, Ca-montmorillonite was successfully purified from Bentonil-WRK bentonite by performing the sequential physical and chemical treatments, and its geochemical properties were characterized using X-ray diffraction, Brunauer-Emmett-Teller analysis, cesium-saturation method, and controlled continuous acidbase titration. Further, batch experiments were conducted to evaluate the adsorption properties of Cs(I) and Sr(II) onto the homoionic Ca-montmorillonite under ambient conditions, and the diffuse double layer model-based inverse analysis of sorption data was performed to establish the relevant surface reaction models and obtain corresponding thermodynamic constants. Two types of surface reactions were identified as responsible for the sorption of Cs(I) and Sr(II) onto Ca-montmorillonite: cation exchange at interlayer site and complexation with edge silanol functionality. The thermodynamic sorption modeling provides acceptable representations of the experimental data, and the species distributions calculated using the resulting reaction constants accounts for the predominance of cation exchange mechanism of Cs(I) and Sr(II) under the ambient aqueous conditions. The surface complexation of cationic fission products with silanol group slightly facilitates their sorption at pH > 8.