Physically-based models have been used to assess landslide susceptibility, hazard, and risk in many regions worldwide. They have also been regarded as valuable tools for landslide prediction and the development or improvement of landslide early warning systems. They are usually validated to demonstrate their predictive capacity, but they are not deeply studied regularly to understand the sensitivity of the input variables and the behavior of the models under many different rainfall scenarios. In this research paper, we studied two distributed physically-based models for shallow landslides: SLIP and Iverson. For this, the first-order second-moment (FOSM) method was used to calculate the contribution of random input variables (soil strength, unit weight, and permeability parameters) to the variance of the factor of safety. Different intensity and duration rainfall events were simulated to assess the response of the models to those rainfall conditions in terms of the factor of safety and failure probability. The results showed that the shear strength (cohesion and friction angle, in order of significance) parameters have the largest contribution to the variance in both models, but they vary depending on geological, geotechnical, and topographic conditions. The Iverson and SLIP models respond in different ways to the variation of rainfall conditions: for shorter durations (e.g. ≤ 8 h), increasing the intensity caused more unstable areas in the SLIP model, while for longer durations the unstable areas were considerably higher for the Iverson model. Understanding those behaviors can be useful for practical and appropriate implementation of the models in landslide assessment projects.