Laminar to turbulent transition induced by spanwise periodic arrays of cylindrical roughness elements over a NACA 0012 airfoil is investigated by hotwire anemometry and infrared thermography. The roughness elements are placed in the flow under adverse pressure gradient. Three configurations are investigated, namely an isolated roughness element, a spanwise array of roughness elements, and a pair of arrays in stagger. The streamwise and spanwise interactions between roughness wakes are addressed, focusing on the evolution of mean flow features and mechanisms for the subsequent process of laminar-turbulent transition. The spanwise interaction between roughness elements involves the connections and merging of neighboring low-speed regions (MLS) in the wake, which affects the spanwise distribution and amplitude of the velocity streaks. The maximum effect on promoting transition is observed when two neighboring low-speed regions overlap with each other in the near wake (within 6 times roughness height). The addition of a second roughness array promotes transition when the spanwise spacing is larger than two times the roughness diameter. Spectral analysis of the streamwise velocity fluctuations reveals that the number of roughness elements within the spanwise array affects the number of MLSs and the dominant instability mechanism. For an odd number of MLSs, the Kelvin–Helmholtz instability dominates the growth of velocity fluctuations around the three-dimensional shear layers. For an even number of MLSs, both Kelvin–Helmholtz and asymmetric instabilities appear in the wake. In this case, the dominant mode that leads to transition depends on the spanwise spacing between roughness elements.