Abstract Spatiotemporal variations in momentum flux spectra of convective gravity waves (CGWs) at the source level (cloud top), including nonlinear forcing effects, are examined based on calculations using an offline version of CGW parameterization and global reanalysis data for a period of 32 years (1979–2010). The cloud-top momentum flux (CTMF) is not solely proportional to the convective heating rate but is affected by the wave-filtering and resonance factor and background stability and temperature underlying the convection. Consequently, the primary peak of CTMF is in the winter hemisphere midlatitudes, associated with storm tracks, where a secondary peak of convective heating exists, whereas the secondary peak of CTMF appears in the summer hemisphere tropics and intertropical convergence zone (ITCZ), where the primary peak of convective heating exists. The magnitude of CTMF fluctuates largely with 1-yr and 1-day periods in major CTMF regions. At low latitudes and Pacific storm-track regions, a 6-month period is also significant, and the decadal cycle appears in the southern Andes. The equatorial eastern Pacific region exhibits a substantial interannual to decadal scale of variabilities. The correlation between convective heating and the CTMF is relatively lower in the equatorial region than in other regions. The CTMF in 10°N–10°S during the period of the pre-Concordiasi campaign approximately follows a lognormal distribution but with a slight underestimation in the tail of the probability density function. In Part II, the momentum flux and drag of CGW in the stratosphere will be examined.