A vertical soil moisture(θ) based unsaturated soil hydrodynamic numerical model was developed and used to simulate soil moisture movement at Miyun (orchard woodland), Daxing (Suburban farmland), Guantao (plain farmland) sites, which represent different underlying surface types in the Hai River Basin. The model based on one dimensional Richards’ equation (hereinafter referred to as the RE model) and the top boundary was input ground-measured rainfall and evapotranspiration data. Soil moisture during different growth stages was numerically simulated based on the fully implicit Finite Differential Method to obtain time series of soil moisture profile. Soil moisture simulated by mature HYDRUS-1D software and observed at the sites was used for crossvalidation and direct verification on the RE model simulation results, respectively. The result showed that the RE model was able to simulate the soil water dynamic changes of different underlying surfaces in Hai River Basin, with root mean square error (RMSE) compared to the observed soil moisture at three sites being  0.0313, 0.0359 and 0.0409 cm³/cm³, respectively. Compared to the simulation results of HYDRUS-1D software (whose RMSE were 0.0376, 0.0647 and 0.0467 cm³/cm³, respectively), the soil water simulated by RE model had higher precision which also showed the RE model reliability. The spatial and temporal variation of soil moisture and its impact factors were discussed. For the case of Daxing site，we analyzed the impact of plow on soil moisture modeling by optimizing the RE model parameters, under which consideration improved the simulation accuracy of the RE model was improved.