Large Eddy Simulations have been performed for fully developed turbulent flow within and above explicitly resolved simple cube arrays. The model (LES-CITY) was validated in comparison with laboratory experiments. The systematic influence of cubic density on turbulent flow characteristics was investigated through numerical experiments in a wide range of cubic area density (0 to 0.44).<BR>The following results were obtained. 1) The intermittency of flow: the temporally and spatially averaged flow structure confirms the existence of conventionally described canyon flow regime. However, the intermittency of the canyon flow for any cubic densities is large and the stream patterns are never persistent. 2) Turbulent Organized Structures (TOS): TOS are characterized by longitudinally elongated low speed streaks and corresponding streamwise vortices. The streaks in sparse and dense canopy flows are likely to be aligned to the street-line and to the roof-line, respectively.