A self-organizing model was developed for simulating rill erosion process on slopes with particular attention to the role of gravitational erosion.For a complete simulation circle,processes such as precipitation,infiltration,runoff,scouring,gravitational erosion and elevation variation were fully considered.Precipitation time(or runoff time) was regarded as iteration benchmark in the model.To specify the contribution of gravitational erosion to the process of rill formation and development,a gravitational erosion module was inserted into the model.Gravitational erosion in rill development was regarded as a Gaussian random process.A model was calibrated by our experimental data,and further validated satisfactorily with 22 runs of experimental results from different investigators. Systematic comparison was made between sediment yields with and without consideration of gravitational erosion module.It was demonstrated that the model could reasonably simulate the rill erosion process under a variety of slope gradients,rainfall intensities and soil conditions upon the gravitational erosion being considered.However,the role of gravitational erosion on sediment yields in rill systems varies significantly under different conditions,although it is of the utmost importance in steeper slopes.The process of gravitational erosion in rill development was studied by a newly-defined parameter a>,which is defined as the volume ratio of gravitational erosion over hydraulic-related erosion.The gravitational contribution to the total erosion could be over 50%for the rill systems with higher rainfall intensity and steeper slopes. A self-organizing model was developed for simulating rill erosion process on slopes with particular attention to the role of gravitational erosion. For a complete simulation circle, processes such as precipitation, infiltration, runoff, scouring, gravitational erosion and elevation variation were fully considered. Precipitation time (or runoff time) was as as iteration benchmark in the model. To specify the contribution of gravitational erosion to the process of rill formation and development, a gravitational erosion module was inserted into the model. Gravitational erosion in rill development wastioned as a Gaussian random process. A model was calibrated by our experimental data, and further validated satisfactorily with 22 runs of experimental results from different investigators. Systematic comparison made between sediment yields with and without consideration of gravitational erosion module. could reasonably simulate the rill erosion process under a variety of slope gradients, rainfall intensities and soil conditions upon the gravitational erosion be considered. However, the role of gravitational erosion on sediment yields in rill systems is significantly under different conditions, although it is of utmost importance in steeper slopes. The process of gravitational erosion in rill development was studied by a newly-defined parameter a>, which is defined as the volume ratio of gravitational erosion over hydraulic-related erosion. gravitational contribution to the total erosion could be over 50% for the rill systems with higher rainfall intensity and steeper slopes.