[en] 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 ⁸⁵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 saturated and unsaturated porous media with equilibrium sorption, FEMWASTE, was also used for comparison. The solution indicated that NESOR gives better simulation results than FEMWASTE for the test case