Bedrock topography is known to affect subsurface water flow and thus the spatial distribution of pore water pressure, which is a key factor for determining slope stability. Therefore, the aim of this study is to investigate the effect of bedrock topography on the timing and location of landslide initiation using 2D and 3D simulations with a hydromechanical model and the Local Factor of Safety (LFS) method. A set of synthetic modeling experiments was performed where water flow and slope stability were simulated for 2D and 3D slopes with layers of variable thickness and hydraulic parameters. In particular, the spatial and temporal development of water content, pore water pressure, and the resulting LFS were analyzed. The results showed that the consideration of variable bedrock topography can have a significant effect on slope stability and that this effect is highly dependent on the intensity of the event rainfall. In addition, it was found that the consideration of 3D water flow may either increase or decrease the predicted stability depending on how bedrock topography affected the redistribution of infiltrated water. pore water pressure [1,2,21] and reduce matric suction [23]. Accordingly, the suction stress and the effective stress, which is the sum of the total stress and the suction stress [24] are reduced. Such a reduction in the effective stress reduces the shear strength of the soil, which may cause failure in the hillslope [12,25]. Therefore, a dynamic assessment of coupled hydromechanical processes in variably saturated hillslope is valuable to assess slope stability.A wide range of studies has focused on developing models to predict the timing and distribution of rainfall-induced landslides [26][27][28]. Many studies have used simplified representations of water flow and determined slope stability using simple limit-equilibrium methods (e.g., infinite-slope stability method) despite well-known limitations [29] such as the required a priori knowledge of the shape of the failure plane and the overestimation of slope stability [30]. Most studies have also focused on 2D slope stability modeling even though some studies have shown that 2D and 3D slope stability assessments can differ considerably [31]. Recently, the local factor of the safety (LFS) concept was proposed to assess slope stability [30]. This Coulomb stress-field based method describes the stability status of cohesive variably saturated soils at each point within hillslopes and does not require a priori assumptions with respect to the location and shape of the failure plane. By taking full advantage of modern numerical solutions for variably saturated flow and stress distribution in hillslopes, the LFS method can be applied on unstructured meshes with high accuracy. Lu et al. [30] showed that the results of stability assessment using the LFS method are in agreement with results of conventional methods while, at the same time, providing further insights into the initiation and evolution of the potential failure surface. Therefore, the LFS method ...