A novel selective mode decomposition, proper orthogonal decomposition and dynamic mode decomposition methods are used to analyse large-eddy simulation data of the flow field about a NACA0012 airfoil at low Reynolds numbers of $5\times 10^4$ and $9\times 10^4$ , and at near-stall conditions. The objective of the analysis is to investigate the structure of the laminar separation bubble (LSB) and its associated low-frequency flow oscillation (LFO). It is shown that the flow field can be decomposed into three dominant flow modes: two low-frequency modes (LFO-Mode-1 and LFO-Mode-2) that govern an interplay of a triad of vortices and sustain the LFO phenomenon, and a high-frequency oscillating (HFO) mode featuring travelling Kelvin–Helmholtz waves along the wake of the airfoil. The structure and dynamics of the LSB depend on the energy content of these three dominant flow modes. At angles of attack lower than the stall angle of attack and above the angle of a full stall, the flow is dominated by the HFO mode. At angles of attack above the stall angle of attack the LFO-Mode-2 overtakes the HFO mode, triggers instability in the LSB and initiates the LFO phenomenon. Previous studies peg the structure, stability and bursting conditions of the separation bubble to local flow parameters. However, the amplitude of these local flow parameters is dependent on the energy content of the three dominant flow modes. Thus, the present work proposes a more robust bursting criterion that is based on global eigenmodes.