Landslides, which are triggered by heavy rainfall or seismic activity, pose serious threats in mountainous regions, necessitating the implementation of structural mitigations, early warning systems, and hazard maps. Traditionally, hazard maps have been empirically derived based on topographical parameters; however, recent advancements have seen the integration of numerical models with digital elevation models (DEMs) to calculate hazard maps. These numerical models are grounded in a comprehensive understanding of landslide dynamics. Large-scale landslides often inundate a larger area than predicted by existing models, which is attributed to the impact of suspended fine sediment reducing bottom friction. To address the inadequacies of existing models in reflecting these effects, we developed a new numerical model for largescale landslides that considers the suspension of fine sediment. First, we investigated the relationship between the quantity of suspended fine sediment and the kinematic conditions of a landslide through flume experiments, resulting in a regression equation for estimating the quantity of suspended fine sediment. This equation was incorporated into a two-dimensional depth-averaged model, forming the basis for the new large-scale landslide model. When this model is integrated with a DEM, it accurately calculates the inundated area associated with a landslide. We tested the model using data for a largescale landslide event in Japan, where pre-and post-landslide DEMs were available for model validation. The results showed that our model successfully replicated the actual inundation area, demonstrating its potential utility in generating hazard maps based on numerical simulations.