Global gyrokinetic Vlasov simulations for trapped ion modes are performed by solving a Vlasov equation averaged over the cyclotron and bounce motions of trapped ions. The distribution function, for trapped ions, is then calculated in a two-dimensional phase space, parametrized by the longitudinal action (energy) and the magnetic moment in presence of magnetic shear. The physical mechanism of the saturation processes between streamerlike structures and zonal flows in relation to the suppression of turbulent transport is discussed. The magnetic shear is identified to play a key role in the dominant streamer-induced transport regime, which exhibits a Bohm-like scaling. The interaction of streamerlike structures with plasma turbulence is shown to produce the inverse cascade that condenses onto long-wavelength trapped ion structures, on the basis on wave triad interactions.