The large-scale landslides that have occurred around the world, mainly due to extreme rainfall events, serve as a warning signal for the need to understand the triggering mechanisms of these mass movements. However, the influence of fundamental variables, such as soil thickness and initial moisture conditions, in the study of watershed stability has not been fully clarified. In this dissertation, numerical flow simulations and slope stability analysis are carried out for the hydrographic basins of the Quitite and Papagaio rivers. In this study, two analytical soil thickness models are applied to saturated-unsaturated three-dimensional finite element flow simulations. The initial moisture condition is analyzed using precipitation data for a spin-up period. Subsequently, the safety factor (SF) is analysed, using the infinite slope method, after the imposition of an extreme rainfall event. From these results, it is possible to verify that the ranges of soil thicknesses from 0 to 3 m are the ones that most resemble the mapped scars, suggesting a more realistic approach for considering the average basin thickness. In addition, it is verified that, in extreme events, the use of spin-up time for initial variable humidity conditions may not have much influence on the safety factor, since the intensity and duration of the rain will be the triggering agents.