The local buckling behavior of welded square box section columns subjected to pure compression is investigated. Local buckling represents a crucial failure mode in thin-walled structures, exerting a significant impact on their overall stability and load bearing capacity. The primary objective of this research is to perform an extensive literature review considering the theoretical background of buckling phenomena and encompassing key findings and methodologies reported in previous studies. Additionally, the development and validation of a novel numerical model is presented, capable of accurately predicting the ultimate buckling capacity. Two different calculation methods are applied in the present study: (i) a numerical model using equivalent geometric imperfections to cover the residual stresses and out-of-straightness of plates, (ii) realistic geometric imperfections combined with an assumed residual stress pattern which has an experimental-based background. The objective of the numerical investigation is to investigate the accuracy of the numerical model by using different residual stress and imperfection patterns taken from the international literature. Many test results are collected from the international literature, to which the computational results are compared, and the effect of the residual stresses and geometric imperfections are analyzed. Based on the numerical analysis, the accuracy of the imperfection models is assessed and the imperfection model leading to the most accurate resistance is determined. The calculated buckling capacities are also compared to analytical design approaches, in which accuracy is also analyzed and evaluated. The current investigation proved the buckling curve developed by Schillo gives the most accurate results to the numerically calculated buckling resistance.