Two modal decomposition techniques, including proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD), are used to identify the wake patterns past single and two crossing cylinders in 60° and 90° arrangements with gap ratio G = 4. The flow is simulated using direct numerical simulations (DNS) for Reynolds numbers Re = 100. The spatial scale of flow decreases with increasing modal frequency from the modal analysis. Two main modes are identified in the wake of the cylinders, namely spatially antisymmetric and symmetric modes. Antisymmetric and symmetric modes are related to cylinders’ vortex shedding and shedding vortices’ shift motion, respectively, whose frequencies are odd and even multiples of cylinders’ lift force frequency. In addition, a low-frequency mode concerning the shadowing effect of the downstream cylinder (DC) in 90° arrangement is found in the wake of the DC centre.