Vegetation plays a critical role in regulating the catchment water balance and enhancing soil stability through root reinforcement. The dynamic nature of vegetation, particularly its seasonal change, significantly affects the magnitude of this influence. However, quantifying the long-term impacts of dynamic vegetation on both flood and landslide occurrences at the catchment scale remains challenging due to the complexity of root structures and the varying dimensions of landslides. In this study, we improved the coupled hydrological-geotechnical model iHydroSlide3D v1.0 by incorporating key vegetation components, such as Leaf Area Index (LAI), root characteristics, and their seasonal dynamics. The improved model was validated using historical observations and applied to a 100-years simulation driven by a weather generator. Three computational scenarios were employed to assess the influence of vegetation on key hydrological and slope-stability variables. Results show that vegetation reduces soil moisture and runoff during low to moderate rainfall events but has a limited impact during larger rainfall events. Additionally, slope stability is found to be more influenced by root reinforcement than soil water uptake. The dynamic nature of vegetation plays a decisive role in modulating its effects on hydrological processes and soil stability, depending on the growth or decay trend of vegetation. This modeling framework offers a robust tool for assessing long-term flood and landslide risks in vegetated catchments.

Physically-based hydrological-geotechnical modeling at large scales is difficult, especially due to the time-consuming nature of flow routing and 3D soil stability models. Although parallelization techniques are commonly used for each model individually, there is currently no concurrent parallelization strategy for both. This study proposed an open-source, Parallelized, and modular modeling software for regional Hydrologic processes and Landslides simulation and prediction (PHyL v1.0). It offers parallel computation in both hydrological and 3D slope stability modules, cross-scale modeling ability via a soil moisture downscaling method, and advanced input/output (I/O) and post-processing visualization. Additionally, PHyL v1.0 is flexible and extensible, making it compatible with all mainstream operating systems. We applied PHyL v1.0 in the Yuehe River Basin, where the computational efficiencies, parallel performance, parameter sensitivity analysis, and predictive capabilities were evaluated. The PHyL v1.0 is therefore appropriately used as an advanced software for high-resolution and complex simulations of regional floods and landslides.