Velocity fields and density fluctuations of edge turbulence are studied in I-mode [F. Ryter et al., Plasma Phys. Controlled Fusion 40, 725 (1998)] plasmas of the Alcator C-Mod [I. H. Hutchinson et al., Phys. Plasmas 1, 1511 (1994)] tokamak, which are characterized by a strong thermal transport barrier in the edge while providing little or no barrier to the transport of both bulk and impurity particles. Although previous work showed no clear geodesic-acoustic modes (GAM) on C-Mod, using a newly implemented, gas-puff-imaging based time-delay-estimate velocity inference algorithm, GAM are now shown to be ubiquitous in all I-mode discharges examined to date, with the time histories of the GAM and the I-mode specific [D. Whyte et al., Nucl. Fusion 50, 105005 (2010)] Weakly Coherent Mode (WCM, f = 100–300 kHz, Df/f approx. 0.5; and k_theta approx. 1.3 cm-1) closely following each other through the entire duration of the regime. Thus, the I-mode presents an example of a plasma state in which zero frequency zonal flows and GAM continuously coexist. Using two-field (density-velocity and radial-poloidal velocity) bispectral methods, the GAM are shown to be coupled to the WCM and to be responsible for its broad frequency structure. The effective nonlinear growth rate of the GAM is estimated, and its comparison to the collisional damping rate seems to suggest a new view on I-mode threshold physics.