In this paper, a mode identification technique in the context of spline finite strip method (SFSM) is presented to compute the contribution of primary (global, distortional and local) and secondary (shear/transverse extension) buckling modes. The base vectors corresponding to individual buckling modes are developed based on the principles of generalized beam theory. The buckling mode shape in SFSM is approximated as a linear combination of these orthonormal base vectors to evaluate the participation of individual buckling mode. The proposed mode identification technique is able to successfully quantify the participation of different buckling modes and the mode participation is comparable with mode identification using finite strip method (FSM) and generalized beam theory (GBT). Illustrative examples are presented to calculate the participation of individual modes in cold-formed steel sections under different loading and boundary conditions. Also the specific application of mode identification in SFSM is demonstrated. © 2017 Elsevier Ltd