In this study, we experimentally investigated the flow field characteristics in the complex geometry of a tightly packed 61-rod bundle with wire spacers by employing matched-index-of-refraction (MIR) and stereoscopic particle image velocimetry (SPIV) techniques. SPIV measurements were performed on a cross-flow plane at Reynolds numbers of 500, 2500, and 6300, and the measurement area covered the interior and edge sub-channels that were located near the enclosure wall. The SPIV results showed that peaks of the velocity magnitude occurred in the edge sub-channel, indicating the wall effects on the flow symmetry in the cross-flow plane. Mean velocities at the central areas of the edge and interior sub-channels showed flat profiles for Re=500 and local peaks around 0.15<y/Drod<0.3, due to the wire’s presence and near-wall location of the edge sub-channel. The rms fluctuating profiles velocities showed that local maximal peaks linearly shifted from y/Drod=-0.44 to −0.36 and −0.28 as Re increased. Taylor’s hypothesis was applied to the SPIV velocity vector fields for reconstructing quasi-instantaneous three-dimensional (3D) vortical structures in the wire-wrapped rod bundle. The visualization of 3D vorticity iso-surfaces showed that vortical structures were generated from shear layers of neighboring rods and traveled in the streamwise direction. Finally, proper orthogonal decomposition (POD) analysis was applied to SPIV velocity snapshots taken at the edge and interior sub-channels to extract the most statistically dominant flow structures. The POD velocity decomposition revealed coherent structures with sizes and shapes comparable to instantaneous vortical structures identified by iso-surfaces. Energy fractions of the first POD modes were found to reduce from 46% and 31% to 26% and 18% when Reynolds numbers increased from Re=500 to 6300, respectively. On the other hand, kinetic energy levels of low-order POD modes (modes 2, 3, 4, etc.) increased when the Reynolds number increased.