The convection dispersion equation is coupled with non equilibrium sorption equation and a corresponding two dimension numerical model, NESOR, is developed. The calculating algorithm for sorption/desorption rate coefficient adopted in the model is presented. To validate the model, NESOR was tested using the field monitoring data of 85 Sr migration in unsaturated loess in test pit F at the CIRP's Field Test Site. Another two dimension finite element model for radionuclide transport through saturat...