The unprecedented liquefaction-related land damage during earthquakes has highlighted the need to develop a model that better interprets the liquefaction land damage vulnerability (LLDV) when determining whether liquefaction is likely to cause damage at the ground’s surface. This paper presents the development of a novel comprehensive framework based on select case history records of cone penetration tests using a Bayesian belief network (BBN) methodology to assess seismic soil liquefaction and liquefaction land damage potentials in one model. The BBN-based LLDV model is developed by integrating multi-related factors of seismic soil liquefaction and its induced hazards using a machine learning (ML) algorithm-K2 and domain knowledge (DK) data fusion methodology. Compared with the C4.5 decision tree-J48 model, naive Bayesian (NB) classifier, and BBN-K2 ML prediction methods in terms of overall accuracy and the Cohen’s kappa coefficient, the proposed BBN K2 and DK model has a better performance and provides a substitutive novel LLDV framework for characterizing the vulnerability of land to liquefaction-induced damage. The proposed model not only predicts quantitatively the seismic soil liquefaction potential and its ground damage potential probability but can also identify the main reasons and fault-finding state combinations, and the results are likely to assist in decisions on seismic risk mitigation measures for sustainable development. The proposed model is simple to perform in practice and provides a step toward a more sophisticated liquefaction risk assessment modeling. This study also interprets the BBN model sensitivity analysis and most probable explanation of seismic soil liquefied sites based on an engineering point of view.