The flow behind nine different arrays of cylinders is experimentally investigated via Particle Image Velocimetry (PIV) at a Reynolds number of Re ∼100 based on the diameter of the cylinders. Each array consists of a column of four cylinders in front and three in the rear. The horizontal distance between the two columns and the vertical distance between the cylinders within each column are varied for H/D=[2,4,8] and V/D=[2,4,6], resulting in nine different arrays denoted as mVnH, where m corresponds to V/D and n stands for H/D. The PIV measurements are conducted for 15 s at 200 Hz frequency, corresponding to 39 to 360 vortex shedding events for the wakes in this study. Then, proper orthogonal decomposition is applied to the velocity fields to analyze the flow dynamics. All arrays show unsteady flow, and based on their flow structures, they are classified in to three main categories of single bluff body (SBB), transitional (TR), and co-shedding (CS) flow. SBB characteristics can be seen for 2V2H and 2V4H arrays, but the latter has more steady vortex shedding as the H/D increases from 2 to 4. Then, 2V8H and 4V2H have an asymmetric flow with several vortex streets and act as an intermediary stage in the shift from SBB to CS flow structure when the distances are increased. The highest total kinetic energy values and widest probability density functions of the velocity components are observed for this group. The five remaining arrays in the CS group have symmetric flow, with three or five vortex streets present behind. However, based on the distances, the frequency and phase synchronization of the vortex streets change considerably, which might have an important effect on, for example, the heat transfer or the structural load of the cylinders.