Oceanic Tidal Mixing as a Contributor to Milankovitch-scale Climate Change

We propose that changes in the magnitude of oceanic tidal mixing on long time scales is an important, but previously unrecognized, contributor to global climate change. It is well known that Earth's orbital and rotational state changes significantly on 10⁴-10⁵ year time scales, and that this influences the spatial and temporal pattern of incident radiation. It is widely supposed that climatic variations on these same time scales are, in large part, a response of the ocean-atmosphere-cryosphere system to this radiative forcing. Our proposal is that variations in the luni-solar tidal potential, induced by these same orbital and rotational variations, influences oceanic mixing and thus modulates meridional heat transport, by amounts which are competitive with the radiative forcing. There are some obvious differences between tidal potential and insolation. First is that the Sun and Moon both contribute to tides, whereas the radiation is entirely of solar origin. Second is that the Earth is transparent to gravity but opaque to radiation. Clipping associated with this opacity makes the radiation pattern temporal spectrum rather more complex than the tidal spectrum. A third point is that solar radiation directly delivers energy to Earth's surface whereas tidal mixing will only expedite lateral transport of heat in association with oceanic thermo-haline circulation. The diurnal average insolation pattern is best parameterized via a Fourier series in time of year and Legendre polynomials in sine of latitude. Our present focus will be on the annual average terms. The Legendre degree n=0 term describes the global average insolation, and is nearly constant. The degree n=1 term describes differences between northern and southern hemispheres, and the annual mean is zero. The degree n=2 term is the main contributor to the equator to pole variations, and varies with obliquity and orbital eccentricity, with the obliquity variation dominating. The lowest order decomposition of the tidal potential recognizes 3 constituents: semi-diurnal, diurnal, and long period, with solar and lunar contributions to each. Our present focus will be on long term variations in the mean square amplitude of the semi-diurnal constituent, with averaging over all the short period variations. For the solar tide that includes the day and year. For the lunar tide it includes the day, month, year, and the apsidal (8.85 year) and nodal (18.6 year) periods. We present calculations of the variations in radiative and tidal forcing for the past 3 million years. The two signals are quite similar. Both vary by ~1% of their respective mean values, are dominated by obliquity variations, and exhibit only secondary influence from orbital eccentricity.