A new hybrid optimization algorithm was presented by integrating the gravitational search algorithm(GSA)with the sequential quadratic programming(SQP),namely GSA-SQP,for solving global optimization problems and minimization of factor of safety in slope stability analysis.The new algorithm combines the global exploration ability of the GSA to converge rapidly to a near optimum solution.In addition,it uses the accurate local exploitation ability of the SQP to accelerate the search process and find an accurate solution.A set of five well-known benchmark optimization problems was used to validate the performance of the GSA-SQP as a global optimization algorithm and facilitate comparison with the classical GSA.In addition,the effectiveness of the proposed method for slope stability analysis was investigated using three case studies of slope stability problems from the literature.The factor of safety of earth slopes was evaluated using the Morgenstern–Price method.The numerical experiments demonstrate that the hybrid algorithm converges faster to a significantly more accurate final solution for a variety of benchmark test functions and slope stability problems. A new hybrid optimization algorithm was presented by integrating the gravitational search algorithm (GSA) with the sequential quadratic programming (SQP), namely GSA-SQP, for solving global optimization problems and minimization of factor of safety in slope stability analysis. New new the global exploration ability of the GSA to converge rapidly to a near optimum solution. In addition, it uses the accurate local exploitation ability of the SQP to accelerate the search process and find an accurate solution. A set of five well-known benchmark optimization problems was used to validate the performance of the GSA-SQP as a global optimization algorithm and convenient comparison with the classical GSA. addition, the effectiveness of the proposed method for slope stability analysis was investigated using three case studies of slope stability problems from the literature The factor of safety of earth slopes was evaluated using the Morgenstern-Price method. Numerical experiments demonstrate that the hybrid algorithm converges faster to a significantly more accurate final solution for a variety of benchmark test functions and slope stability problems.