17.4 Theory of Ocean Tides
Tides have been so important for commerce and science for so many thousands of years that tides have entered our everyday language: time and tide wait for no one, the ebb and flow of events, a high-water mark, and turn the tide of battle.
Mariners have known for at least four thousand years that tides are related to the phase of the moon. The exact relationship, however, is hidden behind many complicating factors, and some of the greatest scientific minds of the last four centuries worked to understand, calculate, and predict tides. Galileo, Descartes, Kepler, Newton, Euler, Bernoulli, Kant, Laplace, Airy, Lord Kelvin, Jeffreys, Munk and many others contributed. Some of the first computers were developed to compute and predict tides. Ferrel built a tide-predicting machine in 1880 that was used by the U. S. Coast and Geodetic Survey to predict nineteen tidal constituents. In 1901, Harris extended the capacity to 37 constituents.
Despite all this work important questions remained: What is the amplitude and phase of the tides at any place on the ocean or along the coast? What is the speed and direction of tidal currents? What is the shape of the tides on the ocean? Where is tidal energy dissipated? Finding answers to these simple questions is difficult, and the first, accurate, global maps of deep-sea tides were only published in 1994 (LeProvost et al., 1994). The problem is hard because the tides are a self-gravitating, near-resonant, sloshing of water in a rotating, elastic, ocean basin with ridges, mountains, and submarine basins.
Predicting tides along coasts and at ports is much simpler. Data from a tide gauge plus the theory of tidal forcing gives an accurate description of tides near the tide gauge.
Note that many oceanographic books state that the tide is produced by two processes: i) the centripetal acceleration at earth's surface as earth and moon circle around a common center of mass, and ii) the gravitational attraction of mass on earth and moon. However, the derivation of the tidal potential does not involve centripatal acceleration, and the concept is not used by the astronomical or geodetic communities.
To calculate the amplitude and phase of the tide on an ocean planet, we begin by calculating the tide generating potential. This is much easier than calculating the forces. Ignoring for now Earth's rotation, the rotation of moon about Earth produces a potential VM at any point on Earth's surface
where the geometry is sketched in Figure 17.10, g is the gravitational constant, and M is moon's mass. From the triangle OPA in the figure,
Using this in (17.5) gives
r/R » 1/60, and (17.7) may be expanded in powers of r/ R using Legendre polynomials (Whittaker and Watson, 1963: §15.1):
The tidal forces are calculated from the spatial gradient of the potential, so the first term in (17.8) produces no force. The second term produces a constant force parallel to OA. This force keeps Earth in orbit about the center of mass of the Earth-moon system. The third term produces the tides, assuming the higher-order terms can be ignored. The tide-generating potential is therefore:
The tide-generating force can be decomposed into components perpendicular P and parallel H to the sea surface. The tides are produced by the horizontal component. "The vertical component is balanced by pressure on the sea bed, but the ratio of the horizontal force per unit mass to vertical gravity has to be balanced by an opposing slope of the sea surface, as well as by possible changes in current momentum" (Cartwright, 1999: 39, 45). The horizontal component is shown in Figure 17.11. It is:
The tidal potential is symmetric about the Earth-moon line, and it produces symmetric bulges.
If we allow our ocean-covered Earth to rotate, an observer in space sees the two bulges fixed relative to the Earth-moon line as Earth rotates. To an observer on Earth, the two tidal bulges seems to rotate around Earth because moon appears to move around the sky at nearly one cycle per day. Moon produces high tides every 12 hours and 25.23 minutes on the equator if the moon is above the equator. Notice that high tides are not exactly twice per day because the moon is also rotating around Earth. Of course, the moon is above the equator only twice per lunar month, and this complicates our simple picture of the tides on an ideal ocean-covered Earth. Furthermore, moon's distance from Earth R varies because moon's orbit is elliptical and because the elliptical orbit is not fixed.
Clearly, the calculation of tides is getting more complicated than we might have thought. Before continuing on, we note that the solar tidal forces are derived in a similar way. The relative importance of the sun and moon are nearly the same. Although the sun is much more massive than moon, it is much further away.
where R sun is the distance to the sun, S is the mass of the sun, R moon is the distance to the moon, and M is the mass of the moon.
Coordinates of Sun and Moon
A natural reference system for an observer on Earth is the equatorial system described at the start of Chapter 3. In this system, declinations δ of a celestial body are measured north and south of a plane which cuts the Earth's equator.
where φp is latitude at which the tidal potential is calculated, d is declination of moon or sun north of the equator, and τ1 is the hour angle of moon or sun. The hour angle is the longitude where the imaginary plane containing the sun or moon and Earth's rotation axis crosses the Equator.
The period of the solar hour angle is a solar day of 24hr 0min. The period of the lunar hour angle is a lunar day of 24hr 50.47min.
Earth's axis of rotation is inclined 23.45° with respect to the plane of Earth's orbit about the sun. This defines the ecliptic, and the suns declination varies between δ = ± 23.45° with a period of one solar year. The orientation of Earth's rotation axis precesses with respect to the stars with a period of 26,000 years. The rotation of the ecliptic plane causes δ and the vernal equinox to change slowly, and the movement called the precession of the equinoxes.
Earth's orbit about the sun is elliptical, with the sun in one focus. That point in the orbit where the distance between the sun and Earth is a minimum is called perigee. The orientation of the ellipse in the ecliptic plane changes slowly with time, causing perigee to rotate with a period of 20,900 years. Therefore Rsun varies with this period.
Moon's orbit is also elliptical, but a description of moon's orbit is much more complicated than a description of Earth's orbit. Here are the basics. The moon's orbit lies in a plane inclined at a mean angle of 5.15° relative to the plane of the ecliptic. And lunar declination varies between δ = 23.45 ± 5.15° with a period of one tropical month of 27.32 solar days. The actual inclination of moon's orbit varies between 4.97°, and 5.32°.
The shape of moon's orbit also varies. First, perigee rotates with a period of 8.85 years. The eccentricity of the orbit has a mean value of 0.0549, and it varies between 0.044 and 0.067. Second, the plane of moon's orbit rotates around Earth's axis of rotation with a period of 17.613 years. Both processes cause variations in Rmoon.
Note that I am a little imprecise in defining the position of the sun and moon. Lang (1980: §5.1.2) gives much more precise definitions.
Substituting (17.15) into (17.9) gives:
Equation (17.16) separates the period of the lunar tidal potential into three terms with periods near 14 days, 24 hours, and 12 hours. Similarly the solar potential has periods near 180 days, 24 hours, and 12 hours. Thus there are three distinct groups of tidal frequencies: twice-daily, daily, and long period, having different latitudinal factors sin2 θ, sin 2θ, and (1 - 3 cos2θ) / 2, where is the co-latitude ( 90° - φ ).
Doodson (1922) expanded (17.16) in a Fourier series using the cleverly chosen frequencies in Table 17.1. Other choices of fundamental frequencies are possible, for example the local, mean, solar time can be used instead of the local, mean, lunar time. Doodson's expansion, however, leads to an elegant decomposition of tidal constituents into groups with similar frequencies and spatial variability.
Using Doodson's expansion, each constituent of the tide has a frequency
where the integers ni are the Doodson numbers. n1 = 1, 2, 3 and n2 - n6 are between -5 and + 5. To avoid negative numbers, Doodson added five ton2 · · · 6. Each tidal constituent, sometimes called a partial tides, has a Doodson number. For example, the principal, twice-per-day, lunar tide has the number 255.555. Because the very long-term modulation of the tides by the change in sun's perigee is so small, the last Doodson number n6 is usually ignored.
If the tidal potential is expanded in Doodson's Fourier series, and if the ocean surface is in equilibrium with the tidal potential, the largest tidal constituents would have frequencies and amplitudes given in table 17.2. The expansion shows that tides with frequencies near one or two cycles per day are split into closely spaced lines with spacing separated by a cycle per month. Each of these lines is further split into lines with spacing separated by a cycle per year (Figure 17.12). Furthermore, each of these lines is split into lines with a spacing separated by a cycle per 8.8yr, and so on. Clearly, there are very many possible tidal constituents.
Why are the tidal lines in Figure 17.12 split into so many constituents? To answer the question, suppose moon's elliptical orbit was in the equatorial plane of Earth. Then d = 0. From (17.16), the tidal potential on the equator, where jp = 0, is:
If the ellipticity of the orbit is small, R = R0 ( 1 + e ) , and (17.18) is approximately
where a = ( g Mr 2 ) / ( 4 R3 ) is a constant. e varies with a period of 27.32 days, and we can write e = b cos( 2 p f2) where b is a small constant. With these simplifications, (17.19) can be written:
which has a spectrum with three lines at 2 f1 and 2 f1 ± f2. Therefore, the slow modulation of the amplitude of the tidal potential at two cycles per lunar day causes the potential to be split into three frequencies. This is the way amplitude modulated AM radio works. If we add in the the slow changes in the shape of the orbit, we get still more terms even in this very idealized case of a moon in an equatorial orbit.
If you are very observant, you will have noticed that the spectrum of the tide waves in Figure 17.12 does not look like the spectrum of ocean waves in Figure 16.6. Ocean waves have all possible frequencies, and their spectrum is continuous. Tides have precise frequencies determined by the orbit of Earth and moon, and their spectrum is not continuous. It consists of discrete lines.
Doodson's expansion included 399 constituents, of which 100 are long period, 160 are daily, 115 are twice per day, and 14 are thrice per day. Most have very small amplitudes, and only the largest are included in table 17.2. The largest tides were named by Sir George Darwin (1911) and the names are included in the table. Thus, for example, the principal, twice-per-day, lunar tide, which has Doodson number 255.555, is the M2 tide, called the M-two tide.
|Department of Oceanography, Texas A&M University
Robert H. Stewart, email@example.com
All contents copyright © 2005 Robert H. Stewart,
All rights reserved
Updated on November 24, 2008