Liquefaction is one of the most damaging functions of earthquakes in saturated sandy soil. Therefore, clearly advancing the assessment of this phenomenon is one of the key points for the geotechnical profession for sustainable development. This study presents a new equation to evaluate the potential of liquefaction (PL) in sandy soil. It accounts for two new earthquake parameters: standardized cumulative absolute velocity and closest distance from the site to the rupture surface (CAV5 and rrup) to the database. In the first step, an artificial neural network (ANN) model is developed. Additionally, a new response surface method (RSM) tool that shows the correlation between the input parameters and the target is applied to derive an equation. Then, the RSM equation and ANN model results are compared with those of the other available models to show their validity and capability. Finally, according the uncertainty in the considered parameters, sensitivity analysis is performed through Monte Carlo simulation (MCS) to show the effect of the parameters and their uncertainties on PL. The main advantage of this research is its consideration of the direct influence of the most important parameters, particularly earthquake characteristics, on liquefaction, thus making it possible to conduct parametric sensitivity analysis and show the direct impact of the parameters and their uncertainties on the PL. The results indicate that among the earthquake parameters, CAV5 has the highest effect on PL. Also, the RSM and ANN models predict PL with considerable accuracy.