Based on slope analysis software and finite element software,simulation analysis method and probability statistics were applied by establishing ageneralization model of Guandigou watershed in Suide Countyof Northern Shaanxi Province that is a classical slope-gully system.Concerning on the mechanical stabilization and damage by stress in the slope-gully system,the stability,land sliding probability and land sliding amount of the slope-gully system with gradual dam land rise were enhanced,weakened and reduced,respectively.Variation of the dam land rise satisfies binomial,binomial or linear relation.The stress field and displacement field in the slope-gully system were also analyzed by FEM(finite element method)strength reduction theory.When the stability factor was the lowest,the land sliding probability was the greatest,and theslope-gully system was on the edge of sliding erosion,the maximum stress and displacement fields obtainedwere as follows:(1)The maximum displacement of Xdirection:the region about 15mperpendicular to vivoslope in the middle and lower segment of slope-gully system;(2)The maximum displacement of Ydirection:the sector region of about 10mextension from the top of mound slope to both sides;(3)The maximum region of tensile stress:strip region of about 9mvertical extension from the top of mound slope to the edgeline of gully slope,through mound slope.The result may offer valuable reference for allocating the engineering measures of soil and water conservation and implementing the biological measures,and provide a certainreliability basis for evaluating the stability of slope-gully system.