Large-scale landslides hazard analysis is based, among other methods, on the numerical evaluation of the safety factor using physical and mechanical parameters measured in the field. The safety factor is determined as the ratio between the forces acting in the favor of sliding phenomenon and the resistance of the earthen massive. The main purpose of the paper is to create two bi-dimensional models of landslides hazard assessment based on the estimation of the stability factor, under different degrees of soil saturation and considering different scenarios of seismic acceleration. Slope stability analysis allows the simultaneity of two exceptional loads such as saturation and earthquake. The first model follows the methodology provided by the national law. In this case, the safety factor is determined by applying limit equilibrium methods using specific geotechnical engineering software. The second model aims to create a deterministic model for safety factor assessment implemented in a GIS system. This one is based on the infinite slope model. The secondary objectives of the paper are: providing theoretical principles, attaining the comparative analysis between the methodologies mentioned above, identifying the critical points of the created models and the reciprocal validation of the results. The study area is a hilly area located in the NorthWestern part of the Iasi City, Romania. The geotechnical parameters were obtained from the laboratory tests carried out on samples taken from 22 boreholes. The geomorphological parameters resulted from the high quality digital elevation model with 1m resolution. The final maps representing the spatial distribution of the safety factor values are reclassified using a common scale. Similarity analysis of the results indicates a good mutual validation.