Soil-bedrock interface is difficult to determine and remains essentially unknown in most Brazilian slopes. In this thesis, we present an analytic model for the spatial prediction of regolith depth built on the bottomup control on fresh bedrock topography hypothesis and high-resolution topographic data. Most of the parameters of the model represent physical entities that can be measured directly in the laboratory or field. The model includes a term which simulates the loss of regolith due to stochastic mass movements and another term that mimic the bedrock-valley morphology. We reconcile our model with field observations from boreholes using a light dynamic penetrometer at Tijuca massif, Rio de Janeiro. We use Bayesian inference, with Markov chain Monte Carlo simulation to summarize the posterior distribution of the parameters, which led to model parameters that best honor our field data as well as the stratigraphic predictive uncertainty. To test the results of the Bayesian inference in slope stability, we develop a software to integrate unsaturated flow simulations, which provide the pressure head distributions and a numerical limit analysis code, that generates the factor of safety (FS), both in three dimensions. We propagate the stratigraphic uncertainty through the developed program to quantify the FS variability and the probability of failure of a natural unsaturated hillslope in the study region. Finally, we emphasize the importance of bedrock topography in slope stability analysis.