On 13 August 2010,a catastrophic debris flow with a volume of 1.17 million m 3 occurred in Xiaojiagou Ravine near Yingxiu town of Wenchuan county in Sichuan Province,China.The main source material was the landslide deposits retained in the ravine during the 2008 Wenchuan earthquake.This paper describes a two-dimensional hybrid numerical method that simulates the entire process of the debris flow from initiation to transportation and finally to deposition.The study area is discretized into a grid of square zones.A two dimensional finite difference method is then applied to simulate the rainfall-runoff and debris flow runout processes.The analysis is divided into three steps;namely,rainfall-runoff simulation,mixing water and solid materials,and debris flow runout simulation.The rainfall-runoff simulation is firstly conducted to obtain the cumulative runoff near the location of main source material and at the outlet of the first branch.The water and solid materials are then mixed to create an inflow hydrograph for the debris flow runout simulation.The occurrence time and volume of the debris flow can be estimated in this step.Finally the runout process of the debris flow is simulated.When the yield stress is high,it controls the deposition zone.When the yield stress is medium or low,both yield stress and viscosity influence the deposition zone.The flow velocity is largely influenced by the viscosity.The estimated yield stress by the equation,τy= ghsinθ,and the estimated viscosity by the equation established by Bisantino et al.(2010) provide good estimates of the area of the debris flow fan and the distribution of deposition depth. On 13 August 2010, a catastrophic debris flow with a volume of 1.17 million m 3 occurred in Xiaojiagou Ravine near Yingxiu town of Wenchuan county in Sichuan Province, China. The main source material was the landslide deposits retained in the ravine during the 2008 Wenchuan earthquake . This paper describes a two-dimensional hybrid numerical method that simulates the entire process of the debris flow from initiation to transportation and finally to deposition. The study area is discretized into a grid of square zones. A two dimensional finite difference method is then applied to simulate the rainfall-runoff and debris flow runout processes. The analysis is divided into three steps; namely, rainfall-runoff simulation, mixing water and solid materials, and debris flow runout simulation. rainfall-runoff simulation is acquired under the cumulative runoff near the location of main source material and at the outlet of the first branch. The water and solid materials are then mixed to creat e an inflow hydrograph for the debris flow runout simulation. the occurrence time and volume of the debris flow can be estimated in this step .Finally the runout process of the debris flow is simulated .When the yield stress is high, it controls the deposition zone .When the yield stress is medium or low, both yield stress and viscosity influence the deposition zone. Flow rate is largely influenced by the viscosity. Estimated yield stress by the equation, τy = ghsinθ, and the estimated viscosity by the equation established by Bisantino et al. (2010) provide good estimates of the area of ​​the debris flow fan and the distribution of deposition depth.