The search for the critical slip surface on a slope is an optimization problem whereby the factor of safety is minimized over a set of parameters which define the shape of the slip surface. In limit equilibrium slope stability analysis, traditional methods for searching for the critical slip surface include grid search and auto-refine search. More recently, metaheuristic optimization methods such as Particle Swarm and Cuckoo Search, among other variations, have been used to search for critical slip surfaces. These simulate natural processes that search the solution space for a minimum solution for various optimization problems encountered in a vast range of disciplines. Typically, the parameters of spheres or ellipsoids which cut the ground topography are varied to create different slip surfaces. The parameters of cutting planes and wedges can also be varied to create multi-planar slip surfaces using the same metaheuristic techniques. However, critical slip surfaces are not always spherical, ellipsoidal, or planar in nature. This paper introduces a novel method which employs the use of three-dimensional spline surfaces in a metaheuristic search to find the critical slip surface in a slope. By varying the parameters which define the location, size and curvature of the spline, the critical slip surface can be found. The proposed formulation of parameters is shown to perform better than the parameters which define the preceding shapes due to the superior flexibility of a spline surface in its curvature.