The effects of obstacle-height variability on spatial characteristics of turbulent organized structures were investigated with the use of a large-eddy simulation technique for airflows over roughness obstacles. Two-types simulation cases were considered: one is uniform-height case in which uniform-height obstacles are aligned in streamwise direction, the other is height-variability case with staggered higher-height obstacles. Streaky structures were observed above the roughness sublayer (RSL) regardless of obstacle-height variability. When obstacles are uniform, flow fields within the RSL contain low- and high-speed regions along the streamwise streets. When obstacle heights vary, airflow within the RSL collides with the front-facing surfaces of taller obstacles. The statistical features of low- and high-speed structures were examined using the spatial correlations of flow fields centering on strong ejection and sweep, respectively. The ejection- and sweep-center spatial correlations extend forward and backward in the streamwise direction, respectively. Length scales were obtained from the ejection-center and sweep-center spatial correlations. The streamwise lengths vary significantly below the canopy height when obstacles are uniform. In contrast, the streamwise length scales remain nearly constant when obstacle heights vary. The horizontal aspect ratios below the canopy heights indicate that turbulent organized structures over obstacles with variable heights are more isotropic than those over uniform obstacles. The inclination angles of the organized structures were also deduced using the spatial correlations. The angles of sweep-center structures are steeper than those of the ejection-center structures. The angles of the ejection-center structures at the RSL heights become larger with obstacle-height variability.