Grid search method for locating the critical failure surface is extended by deriving additional analytical expressions for slip center grid (x(min), z(min); x(max), z(max))(,) where global minimum of safety factor occurs, including the prediction of minimum and maximum values for safety factor (Fs(min,max)) and for slip circle radius (R-min,R-max). Derived models are proposed in a form of nonlinear functions of geometrical parameters (slope height H, depth to bedrock d and slope angle beta) and soil factors (bulk density gamma, cohesion c, angle of internal friction phi and pore water pressure coefficient tau(u)). Research was performed using Box-Behnken experimental design, for which the input data were provided by Spencer limit equilibrium analyzes of different slopes with circular slip surface. Reasonable predictive power of the proposed models was verified both by internal and external validation, latter of which included the analyzes of slopes with random geometrical and soil properties. Regarding the impact of input parameters, beta has the strongest influence on response values (Fs, R, x, z), with the predominant linear and quadratic effect. As for the influence of remaining factors, c and phi also have strong impact on Fs, while H and phi have significant influence on slip circle radius and the location of slip center grid. However, due to existence of two-factor interactions, it is shown that the effect of beta on Fs, x, z and R is highly dependent on the values of c, phi, tau(u), H and d/H, including the significant effect of tau(u) x phi, c x H and c x gamma. When compared to traditional grid search method, proposed approach could be used to locate the circular slip surface with global minimum of safety factor, without the need for additional slope stability analyzes. (C) 2015 Elsevier Inc. All rights reserved.