Technical note: On seasonal variability of the M2 tide
Seasonal variability of the M2 ocean tide can be detected at many ports, perhaps most. Examination of the cluster of tidal constituents residing within the M2 tidal group can shed light on the physical mechanisms underlying seasonality. In the broadest terms these are astronomical, frictional–advective interactions, and climate processes; some induce annual modulations and some semiannual, in amplitude, phase, or both. This note reviews how this occurs and gives an example from each broad category. Phase conventions and their relationship with causal mechanisms, as well as nomenclature, are also addressed.
It has long been noticed (Darwin, 1907; Corkan, 1934) that ocean tide constituents at some ports may experience significant seasonal variations. Especially noteworthy are large modulations discovered in some polar regions (e.g., Godin, 1986; Rotermund et al., 2021; Bij de Vaate et al., 2021). Significant modulations can also occur at lower latitudes, both regionally (Kang et al., 2002) and especially locally (e.g., Foreman et al., 1995). In fact, nearly all coastal tides show at least a small seasonal modulation (Pugh and Woodworth, 2014).
Climate-driven processes capable of inducing seasonal changes in barotropic tides are myriad: variability in ocean stratification (Müller, 2012), variability in ice cover (Prinsenberg, 1988), seasonal runoff and changes in river discharge (Guo et al., 2015), and tide–surge interactions from predominantly wintertime storms (Pugh and Woodworth, 2014). Similar processes and others can induce seasonal perturbations in baroclinic tides detectable in surface measurements (Ray and Zaron, 2011; Zhao, 2021). The extent the astronomical tidal potential plays in seasonality – usually small, but potentially important where tide amplitudes are large – is often overlooked. Some confusion on the issue was recently clarified by Du and Yu (2021).
To help unravel observations, it is useful to revisit the spectral characteristics behind seasonal variability. Understanding how certain spectral lines in observed sea level arise often points to possible physical mechanisms at work. The purpose of this note is to review this topic, focusing solely on the principal semidiurnal tide M2. Much of what follows is hardly new; the purpose is to review and clarify, including even the nomenclature used. Whether seasonal variability is annual or semiannual, and whether it occurs in amplitude or phase or both, are obviously important aspects of variability; seeing one type of modulation rather than another can narrow the list of causes.
Recall the technical definitions (e.g., Munk and Cartwright, 1966) of “tidal group” – a cluster of spectral lines with the same first two Doodson numbers – and “tidal constituent” – a cluster with the same first three Doodson numbers. Tidal groups are separated in frequency by about one cycle per month, constituents by about one cycle per year. So when one speaks of seasonal variability of the M2 tide, one must be speaking, in a sense, of the M2 tidal group. It is the variability seen, for example, in a series of monthly estimates of M2.
The main constituents within the M2 group are listed in Table 1 in order of frequency. Included are tides generated by the astronomical potential (Cartwright and Edden, 1973; Hartmann and Wenzel, 1995), compound tides generated by shallow-water processes, and two annual sidelines (MA2 and MB2) commonly employed to effect an annual modulation from the wide range of possible climate processes. Climate processes are broadband, but they do commonly display a large spectral peak at once per year, and the two sidelines attempt to account for that. An example of these spectral lines in a real sea-level spectrum is shown in Fig. 1, computed from the tide-gauge measurements at Saint-Malo, France (and used further below).
Since all tides in the table are lunar, they all have 18.6-year nodal sidelines. These will be ignored in the present context, but they are obviously important for tidal prediction.
τ – mean lunar time, h – mean longitude of the Sun, p – mean longitude of lunar perigee, p′ – mean longitude of solar perigee, currently about 283∘.
The constituent arguments given in the table are in the form of the standard mean longitudes that Doodson (1921) found so useful for ordering the whole tidal catalog. Simple expressions, linear in time for present-day tides, are readily available to evaluate the longitudes, and thus the tidal arguments, for any particular time (e.g., Pugh and Woodworth, 2014, p. 68).
All the tides of Table 1 differ in frequency from the central M2 tide by either one or two cycles per year (cpy). They thus act to modulate the amplitude and/or phase of M2 at those frequencies – the source, at least in spectral terms, of the seasonality of M2. Most of the tidal arguments differ from M2 by either ±h or ±2h, so the modulation is tied to the Sun's declination, zero at the equinox; the speed of h is one cycle per tropical year. The two astronomical tides, α2 and β2, differ from M2 by , where is the Sun's mean anomaly, zero at perihelion, so the modulation is tied to the Sun's distance; the speed of ℓ′ is one cycle per anomalistic year. The astronomical tide Γ2 also depends on the position of the moon's perigee p. In an average sense, over many years, the contribution from Γ2 to M2 seasonality mostly cancels out, because p (period 8.85 years) varies from one year to the next.
The frequencies of Table 1 are given with sufficient precision to show the tiny differences stemming from tropical versus anomalistic years, which are almost identical because the slowly moving p′ requires 209 centuries to complete one revolution. So the annual pairs MA2, MB2 and α2,β2 take practically the same frequencies, and either pair may be used in a tidal analysis. However, because their phases differ, the pairs cannot be interchanged; analysis and any subsequent prediction must maintain consistency.
Each class of constituent is now addressed in more detail.
2.1 Astronomical tides
The four gravitational sidelines are all very small, only about 0.3 % of the primary, or smaller. With arguments depending on the mean longitude of the Sun, they evidently arise from the Sun's third-body perturbations of the lunar orbit. In particular, α2,β2 arise from what in lunar theory is called the “annual equation”, which refers to an expansion or compaction of the moon's orbit depending on whether the Sun is at perihelion or aphelion, respectively. With a change in the orbital radius there is a corresponding change in the moon's angular velocity and thus longitude. The variation in longitude is given approximately by (Brouwer and Clemence, 1961, p. 329). The radial variation is 49cos ℓ′ km.
Because the ocean cannot support extremely high-Q resonances, the ocean's response (admittance) to gravitational forcing must be nearly constant across the small frequency band of the M2 group. In that case, the gravitational components of the four sidelines can all be inferred from measurements of the central M2 primary. Their amplitudes must be in the same proportion as the relative amplitudes of Table 1, and their phases must be nearly identical to the phase of M2. It may be possible that this rule is violated in locations where large M2 currents act through nonlinear dissipation to suppress the sidelines; this has been observed to happen for the M2 nodal sideline (Ku et al., 1985).
Aside from this one exception, the induced seasonal modulations in M2 from the astronomical sidelines are easily worked out. At a location with mean M2 amplitude A and phase lag G, the combined elevation of the group is
where ci are the relative amplitudes from Table 1, and the Γ2 constituent has been dropped for the reasons noted above. In analogy with the usual approach for handling nodal modulations, this expression may be written as a single modulated wave:
where functions f,u vary throughout the year as periodic functions of h (or ℓ′). Expanding the trigonometric functions and gathering like terms in cos G and sin G lead to
for the amplitude modulation and
for the phase modulation. The amplitude modulation is insignificant; because the main coefficients in Eq. (2) cancel, the second term from δ2 is actually larger than the first term but still very small. For the phase modulation in Eq. (4), the first term peaks at ∘, which is early April, and the smaller second term shifts that peak to mid-April. So, the observed phase lag of M2 (u subtracts from G) takes its minimum value in mid-April and its maximum in mid-September. It is interesting to note that the first term in Eq. (4) is exactly twice the Moon's longitude modulation of , alluded to above as arising from the annual equation in the Sun's perturbation of the lunar orbit (it is twice, because M2 is semidiurnal).
At most tide gauges, the gravitational components of any seasonal variability in the observed M2 will thus comprise only a minor part. Observed variability will likely differ significantly from what might be inferred from the M2 admittance or as modeled by Eqs. (3)–(4). Nonetheless, when only the astronomical modulation is important and the M2 amplitude is moderately large, a phase modulation of 0.7∘ is easily seen; an example is given in Sect. 3.1.
2.2 Compound tides
Aside from the two compound tides listed in Table 1, there are possibly others that fall within the M2 group (Simon, 2013). The most important is OP2, exactly coinciding with MSK2, and KO2, exactly coinciding with M2. The OP2 and MSK2 constituents likely arise from different nonlinear aspects of the hydrodynamics (Parker, 1991), although both are generated in shallow water.
The KO2 tide is important only when M2 is anomalously small; it can potentially induce an unusual nodal modulation (relative to the standard M2 modulation), but it has no bearing on seasonality, with one exception: any seasonal variations in the primary K1 or O1 tides would induce seasonality in the compound KO2, but cases where this is significant relative to M2 must be rare. The few additional compound tides within the group noted by Simon (2013) are from interactions involving the diurnal S1, which is normally small and will be ignored. Thus, when nonlinear shallow-water processes are acting sufficiently to produce a noticeable MSK2 (or OP2), the main effect of this will be equivalent to semiannual modulations in M2 amplitude and/or phase.
2.3 Climate-induced modulations
The astronomical sidelines in Table 1 sit at discrete known frequencies. Compound tides are similar, although the number of them can sharply increase in shallow water. In contrast, climate-driven modulations are broadband. Although an annual cycle typically dominates, which is the justification for MA2, MB2, one might expect higher harmonics in some cases and also possibly a general smearing of observed spectral lines across the group (e.g., Munk et al., 1965).
Müller et al. (2014), building on an earlier study by Foreman et al. (1995), discussed the case of Victoria (Canada). There was a clear annual modulation in both amplitude and phase of M2 but with an apparent second harmonic. Multi-year tidal analysis of the Victoria data (not shown) does show energy at the MSK2 (or OP2) frequency. Whether this is due to the compound tide(s) or to a true second harmonic below MA2 is not clear. During tidal analysis it is sometimes possible to separate two constituents of identical frequency by exploiting their different nodal modulations (Parker et al., 1999), but in this case the amplitudes appear too weak to allow it.
The question of frictional compound tides versus higher harmonics of MA2 and MB2 is only one instance of the more general, and difficult, problem of deciphering the causes of climate-induced modulations, even when manifested by MA2 or MB2 acting alone (e.g., Pugh and Vassie, 1994).
The constituent names in Table 1 follow both historical and current international conventions – for the latter, see the table maintained by the Tide, Water Level and Current Working Group of the International Hydrographic Organization (https://iho.int, last access: 17 April 2022). The exception is β2, which does not (as yet) appear in the working group's table. Godin's 1988 book has the same usage as the IHO; his earlier 1972 book left the smaller lines unlabeled (Godin, 1972, 1988). Regarding β2, however, it is commonly employed in the Earth tide community (e.g., Hartmann and Wenzel, 1995; Calvo et al., 2016; Ducarme and Schueller, 2018), and β2 falls into the obvious pattern of using the first four letters of the Greek alphabet for the four astronomical constituents.
The two climate constituents apparently began as MA2 and Ma2 (Corkan, 1934). That lost favor, probably because a speaker cannot distinguish the two and also because early computers employed only upper-case letters. MB2 has been used for decades, at least in British work1. Amin (1976) still used Ma2, but by 1983 he too had switched to MB2.
Not included in Table 1 is the pair, H1, H2, for the two annual sidelines, which is starting to appear fairly often in the literature, probably because it is used in a popular open-source software package. The use of H1 for a semidiurnal wave is certainly an oddity, as the tide community has used an integer subscript to denote tidal species since at least the end of the nineteenth century (Darwin, 1883; Harris, 1895; Cartwright, 1999). It is not clear where the symbol originated. An early (possibly first) use was by Pugh and Vassie (1976). Before that, in a discussion of shallow-water tides, Godin (1972, Table 2.14) used the label “(Horn)” for both sidelines, and he cited a discussion by Horn (1960), although Horn himself left the lines unlabeled. Perhaps “(Horn)” has morphed into H1, H2, or the labels merely reflect the argument differences of ±h. An important point, however, is that both Godin and Horn, as well as Pugh and Vassie, were referring to the two climate lines – none of their arguments involved p′ – whereas the H1, H2 in present-day use are substituting for the gravitational tides α2,β2, as their arguments do involve p′. In any event, the use of a wrong subscript should be discouraged.
In this section an example of M2 seasonality arising from each of the three categories – gravitation, frictional interaction, and climate processes – is presented. For each tide gauge analyzed, a tidal solution based on a single inversion of many years of data was computed, sufficient to obtain reliable estimates of all constituents in Table 1. Based on an appropriate set of the estimated constituents (different for each category), the implied modulation of M2 over one year was computed by complex demodulation. This is then compared with results of a second tidal inversion (or rather a set of inversions) in which estimates of M2 were obtained for every month of the multi-year time series, from which monthly means were then computed. The monthly calculations accounted for the conventional 18.6-year nodal modulation of M2, but no other modulations. (Standard errors in these monthly means were estimated from the standard deviations for each month, scaled by for n years of data.)
The goal is to confirm that seasonality of M2, as delineated by monthly mean estimates of amplitude and phase, can – at least in these cases – be accurately reproduced by the modulations from a particular set of spectral lines. Which lines are in play differ depending on the category of causation.
3.1 Astronomical modulations
It is actually not easy to find good examples of seasonal variability stemming solely from the purely astronomical constituents of Table 1. Any potential case must display fairly constant admittances across the group of gravitational constituents. Yet at most tide gauges one sees perturbations in the admittance or one sees significant amplitudes in the compound tides.
The time series at Port Orford, Oregon, is one of the better examples. Tidal constants estimated from 26 years of data2 (1994–2021) are listed in Table 2. The magnitudes of tidal admittances are all consistent within error limits, all phases are close, and the compound MSK2 is very small. Combining the harmonic constants of the three astronomical constituents implies a seasonal modulation of M2 given by the solid lines of Fig. 2. These are in good agreement with the monthly mean estimates, aside perhaps for the January amplitude.
The theoretical seasonal modulations, based on Eqs. (3)–(4), are shown as the dashed lines. The solid and dashed lines agree well in phase and less well in amplitude at first glance. However, note that the amplitude vertical axis spans only 1 cm, so in fact the amplitude agreement is also quite good, with all data implying very little amplitude modulation. The small differences in amplitude curves occur because the estimated admittances in Table 2 are not identical, simply due to unavoidable estimation error.
The analysis at Port Orford confirms that the astronomically induced seasonal modulation of M2 results in almost no amplitude modulation and a phase modulation of about 0.7∘, with a minimum in April and a maximum in September.
3.2 Frictional/advective modulations
An example of modulations dominated by one or both of the compound tides in the M2 group is Saint-Malo, France, whose spectrum was shown in Fig. 1. The two annual constituents (α2,β2 or MA2, MB2) are smaller, as the spectrum reveals, but are still too large to ignore. The total modulation, based on four constituents, is shown in Fig. 3. The amplitude clearly reveals the presence of a semiannual modulation, with a slow rise during the beginning of the year and then a rapid decay between July and October. The phase is dominated by the semiannual effect, with an annual contribution responsible for the September phase lag exceeding the earlier peak in March.
3.3 Annual climate modulations
Tide gauges with annual variations in M2, and thus with sidelines dominated by MA2 and/or MB2, are easy to find. The case of Victoria was already noted (Müller et al., 2014), partly for having a second harmonic in its modulation. However, there are many tide gauges where only the annual modulation presents itself. A fair number can be found along the coast of Japan, even though M2 itself is not especially large there. The example chosen here is in fact one of the largest modulations discovered anywhere: Chittagong, along the coast of Bangladesh. At that location the annual sea-level term, Sa, is also anomalously large, presumably reflecting large discharge from the Ganges. Several other nearby tide gauges, with somewhat smaller M2 modulations, were studied by Tazkia et al. (2017).
Harmonic analysis of 11 years of hourly data (2008–2018) at Chittagong reveals astonishing large amplitudes for MA2 and MB2 of 15.5 and 10.1 cm, respectively, with the M2 constituent at 171.6 cm. The resulting seasonal modulation is shown in Fig. 4. The monthly mean amplitudes range from 146 cm in February to a high of 195 cm in August. In comparison, the phase modulation is not very large, about 4∘. The monthly mean phases appear slightly erratic and fit the demodulated curve only moderately well.
The large modulation in amplitude at this location closely mimics the large oscillation in annual sea level, which is minimum in February and maximum in July, with a mean range of 71 cm. Tazkia et al. (2017) developed a tide model for the region that explores this interdependence.
One indication of the variety of processes responsible for seasonal variability of the M2 tide is the variety of different mechanisms generating spectral lines within the M2 group: astronomical motions of the moon, frictional and other nonlinear interactions between tidal waves, and climate processes. The astronomical contribution is predictable given good long-term mean estimates of the M2 constants; it is mostly a ±0.37∘ modulation in phase, precisely double the solar perturbation in the moon's longitude arising from the “annual equation” of lunar theory. On the other hand, when a tide gauge is found to be affected by substantial seasonality in M2, it usually arises from one or both of the constituents MA2, MB2. The variety of climate processes responsible for those two constituents – annual changes in stratification, sea level, ice cover, etc. – is where the real complication lies when attempting to understand seasonal variability.
The Port Orford and Chittagong tide-gauge data are available from the University of Hawaii Sea Level Center (https://uhslc.soest.hawaii.edu, University of Hawaii Sea Level Center, 2022). The Saint-Malo tide-gauge data are available from Le SHOM (https://doi.org/10.17183/REFMAR, SHOM, 2022), the French national hydrographic service.
The supplement related to this article is available online at: https://doi.org/10.5194/os-18-1073-2022-supplement.
The author has declared that there are no competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
It is a pleasure to thank Philip Woodworth and David Pugh for discussions. This work was supported by the U.S. National Aeronautics and Space Administration through the Sentinel-6 and Sea Level Change projects.
This paper was edited by Joanne Williams and reviewed by Qian Yu and one anonymous referee.
Amin, M.: The fine resolution of tidal harmonics, Geophys. J. Roy. Astr. Soc., 44, 293–310, 1976. a
Bij de Vaate, I., Vasulkar, A. N., Slobbe, D. C., and Verlaan, M.: The influence of Arctic landfast ice on seasonal modulation of the M2 tide, J. Geophys. Res.-Oceans, 126, e2020JC016630, https://doi.org/10.1029/2020JC016630, 2021. a
Brouwer, D. and Clemence, G. M.: Methods of Celestial Mechanics, Academic Press, New York, 1961. a
Calvo, M., Rosat, S., and Hinderer, J.: Tidal spectroscopy from a long record of superconducting gravimeters at Strasbourg, in: International Symposium on Earth Sciences for Future Generations, edited by: Freymueller, J. and Sanchez, L., 131–136, Springer, https://doi.org/10.1007/1345_2016_223, 2016. a
Cartwright, D. E.: Tides: A Scientific History, Cambridge Univ. Press, Cambridge, 1999. a
Cartwright, D. E. and Edden, A. C.: Corrected tables of tidal harmonics, Geophys. J. Roy. Astr. Soc., 33, 253–264, 1973. a
Darwin, G. H.: The harmonic analysis of tidal observations, British Association Report for 1883, 49–118, reprinted in Darwin's Scientic Papers, 1907, Vol. 1, 1–70, Cambridge Univ. Press, 1883. a
Darwin, G. H.: On the Antarctic tidal observations of the “Discovery”, in: Scientific Papers. Volume 1: Oceanic Tides and Lunar Disturbance of Gravity, 372–388, Cambridge Univ. Press, 1907. a
Doodson, A. T.: The harmonic development of the tide-generating potential, Proc. Royal Soc., 100, 305–329, 1921. a
Du, Z. and Yu, Q.: Comment on “Seasonal and nodal variations of predominant tidal constituents in the global ocean” by Cao et al., Cont. Shelf Res., 227, 104524, https://doi.org/10.1016/j.csr.2021.104524, 2021. a
Ducarme, B. and Schueller, K.: Canonical wave grouping as a key to optimal tidal analysis, Marees Terr. Bull. Inf., 150, 12131–12244, http://maregraph-renater.upf.pf/bim/BIM/bim150.pdf (last access: 15 July 2022), 2018. a
Foreman, M. G. G., Walters, R. A., Henry, R. F., Keller, C. P., and Dolling, A. G.: A tidal model for eastern Juan de Fuca Strait and the southern Strait of Georgia, J. Geophys. Res., 100, 721–740, 1995. a, b
Godin, G.: Modification by an ice cover of the tide in James Bay and Hudson Bay, Arctic, 39, 65–67, 1986. a
Godin, G.: Tides, CICESE, Ensenada Baja California, Mexico, 1988. a
Guo, L., van der Wegen, M., Jay, D. A., Matte, P., Wang, Z. B., Roelvink, D., and He, Q.: River-tide dynamics: exploration of nonstationary and nonlinear tidal behavior in the Yangtze River estuary, J. Geophys. Res.-Oceans, 120, 3499–3521, https://doi.org/10.1002/2014JC010491, 2015. a
Harris, R. A.: Manual of tides. Part III, in: Report of the Superintendent, U.S. Coast and Geodetic Survey, Part II, 125–262, U.S. Gov't. Printing Office, Washington, 1895. a
Horn, W.: Some recent approaches to tidal problems, Int. Hydrogr. Rev., 37, 65–88, 1960. a
Kang, S. K., Foreman, M. G. G., Lie, H.-J., Lee, J.-H., Cherniawsky, J., and Yum, K.-D.: Two-layer tidal modeling of the Yellow and East China Seas with application to seasonal variability of the M2 tide, J. Geophys. Res., 107, 3020, https://doi.org/10.1029/2001JC000838, 2002. a
Ku, L.-F., Greenberg, D. A., Garrett, C. J. R., and Dobson, F. W.: Nodal modulation of the lunar semidiurnal tide in the Bay of Fundy and Gulf of Maine, Nature, 230, 69–71, 1985. a
Müller, M.: The influence of changing stratification conditions on barotropic tidal transport and its implications for seasonal and secular changes of tides, Cont. Shelf Res., 47, 107–118, https://doi.org/10.1016/j.csr.2012.07.003, 2012. a
Munk, W. H. and Cartwright, D. E.: Tidal spectroscopy and prediction, Philos. T. R. Soc. A, 259, 533–581, 1966. a
Munk, W. H., Zetler, B., and Groves, G. W.: Tidal cusps, Geophys. J. Roy. Astr. Soc., 10, 211–219, 1965. a
Parker, B. B.: The relative importance of the various nonlinear mechanisms in a wide range of tidal interactions (review), in: Tidal Hydrodynamics, edited by: Parker, B. B., 237–268, John Wiley, New York, 1991. a
Parker, B. B., Davies, A. M., and Xing, J.: Tidal height and current prediction, in: Coastal Ocean Prediction, edited by: Mooers, C. N. K., 277–327, American Geophysical Union, Washington, 1999. a
Prinsenberg, S. J.: Damping and phase advance of the tide in western Hudson Bay by the annual ice cover, J. Phys. Oceanogr., 18, 1744–1751, 1988. a
Pugh, D. T. and Vassie, J. M.: Tide and surge propagation off-shore in the Dowsing region of the North Sea, Dt. Hydrogr. Z., 29, 163–213, 1976. a
Pugh, D. T. and Vassie, J. M.: Seasonal modulations of the principal semidiurnal lunar tide, in: Mixing and Transport in the Environment: A Memorial Volume for Catherine Allen, edited by: Beven, K. J., Chatwin, P. C., and Millbank, J. H., chap. 13, 247–267, John Wiley, 1994. a
Rotermund, L. M., Williams, W. J., Klymak, J. M., Wu, Y., Scharien, R. K., and Haas, C.: The effect of sea ice on tidal propagation in the Kitikmeot Sea, Canandian Arctic archipelago, J. Geophys. Res.-Oceans, 126, e2020JC016786, https://doi.org/10.1029/2020JC016786, 2021. a
Tazkia, A. R., Krien, Y., Durand, F., Testut, L., Papa, F., and Bertin, X.: Seasonal modulation of M2 tide in the northern Bay of Bengal, Cont. Shelf Res., 137, 154–162, https://doi.org/10.1016/j.csr.2016.12.008, 2017. a, b
David Pugh has a copy of a November 1977 memorandum from David Cartwright, then director of the Bidston Observatory, proposing use of the modified symbols MA2 and MB2, in “deference to Corkan's original work”. He suggested similar notation for other constituents affected by seasonal variability, such as NA2, NB2; these are now included in the IHO tables, as are MA4, MB4.
Data at Port Orford are available from 1978 to present, but the data before 1994 yield tidal estimates too erratic to use. Data after 1993 appear to be of good quality.