The characteristics of unsteady wakes behind a stationary square cylinder and another upstream vibrating square cylinder have been investigated numerically with the help of a developed computational code. The effect of longitudinal as well as transverse vibrations of the upstream cylinder is studied on the coupled wake between the two cylinders, which is found to control the vortex shedding behavior behind the downstream stationary cylinder. Computations are carried out for a fixed value of Reynolds number (Re = 200) and three different values of excitation frequencies of the upstream cylinder, namely less than, equal to and greater than the natural frequency of vortex shedding corresponding to flow past a stationary square cylinder. The vortex shedding characteristics of the unsteady wakes behind the vibrating and stationary cylinders are found to differ significantly for longitudinal and transverse modes of vibration of the upstream cylinder. The wake of the downstream stationary cylinder is found to depict a synchronization behavior with the upstream cylinder vibration. The spacing between the two cylinders has been identified to be the key parameter influencing the synchronization phenomenon. The effect of cylinder spacing on the wake synchronization and the hydrodynamic forces has been examined. In addition, a comparison of the drag forces for flow past transversely vibrating square and circular cylinders for similar amplitudes and frequencies of cylinder vibration has been presented while employing the tested computational code. © 2009 The Japan Society of Fluid Mechanics and IOP Publishing Ltd.