Abstract: Tides are one of the basic types of ocean water motions. Previous studies have reported that the M2 constituent exhibits seasonal variations (annual cycles) in some regions. However, based on the newly proposed method of modified two-step harmonic analysis (HA) and its application at 240 global tide gauges, we find that the M2 constituent as well as the S2 and K1 do not have significant seasonal variations at these tide gauges. The seasonal variations of the M2 constituent reported in previous studies are caused by its satellites, the H1 and H2 constituents, which are not resolved in these studies due to the short time window (one month or three months) used in HA. Because the frequency of the H1 (H2) constituent is equal to that of the M2 minus (plus) the frequency of annual cycles, the superposition of the M2, H1 and H2 constituents with constant amplitudes is equivalent to the M2 constituent with seasonally varying amplitudes. Compared with the new method, some adaptations to traditional HA aiming to capture variations in amplitudes and phase lags of constituents have some limitations, because they either neglect some satellites of the major constitutes or introduce spurious fluctuations resulting from an unreasonably large number of independent points. The nodal modulations of predominant constituents are also explored in this study. On the global scale, the nodal modulations of the M2, K1 and O1 constituents agree with the theoretical predictions, except a cold spot region with reduced nodal modulation in the Gulf of Maine and a hot spot region with enhanced nodal modulation in the South China Sea for the M2. Nodal modulation is also found for the S2 constituent (in theory, the S2 has no nodal modulation), which is 0.8% averaged at 164 tide gauges where the S2 is not too weak.