With geotechnical engineering the form of slope stability analysis most commonly used is the Bishop circular method. By assuming a circular shape, however, the method results in an overestimation of the factor of safety. This is because failure surfaces are often non-circular, given the heterogeneity of both hydrological and geotechnical slope characteristics. This problem is addressed by adapting an automated search algorithm for non-circular failure surfaces. The numerical method is based on a random search technique for the global minima of a constrained non-linear optimization problem. This is incorporated into a combined hydrology and stability model (CHASM) that allows the simulation of changes in pore-water pressures in response to individual rainfall events, and considers their role in maintaining slope stability. Initial results suggest that the implementation of the automated slip search algorithm leads to an improvement in the evaluation of slope stability in comparison to other less rigorous approaches. Copyright © 2000 John Wiley & Sons, Ltd.