1. Introduction
The theory of sea level change includes glacial eustasy, ocean volume changes, geoidal eustasy, thermal-hyaline expansion/contraction of the water column, and changes of the rotational ellipsoid of the globe [1] - [6] . In the discussion of present and future sea level changes, the issue has generally become limited to glacial eustasy and thermal expansion [7] . Only observational facts in nature itself are good enough for a trustworthy picture, however.
Past sea level positions can be measured by watermarks, tide gauges, satellite altimetry, shore morphology, coastal stratigraphy, and historical and archaeological information. All of these factors have their limitations and margins of error, however (e.g. [8] ).
There are more than 2000 tide gauge stations [9] , but there are only a few long and continuous records. In formerly glaciated areas, the tide gauge records are dominated by glacial isostatic uplift. In other areas of the globe, the tide gauge records are usually affected by subsidence because the tide gauges are located at river-mouths where the accumulation of sediments has induced subsidence. The extraction of water from beneath cities generates sediment compaction and hence subsidence. Furthermore, many tide gauges are attached to heavy harbour constructions, which generate site-specific compaction (i.e. subsidence). Many tide gauge records are, in fact, fragmentary. Gaps in the recording may not be so serious. Gaps due to destruction (e.g. by storms and tsunamis) may generate serious discontinuities, however. This is certainly the case with relocations of the actual site of measurements. Such discontinuities invalidate straight-line trend analyses (e.g. [10] ). There are also cycles to consider; the 18.6 lunar-tidal cycle and longer “lunar-tidal super cycles”. Finally, there is the problem with the search for “global mean sea level change”, when, in fact, eustatic sea level is a regional factor to be determined separately from region to region [11] . How can the personal selection of 6 sites [12] or 184 sites [13] provide anything but meaningless average values?
Satellite altimetry (e.g. [14] ) is a new and important tool of recording sea level changes starting in late 1992. Whilst NOAA [15] gives a mean sea level rise of +2.9 ± 04 mm/yr, UC [16] gives a means sea level rise of 3.3 ± 0.4 mm/yr. Both values are strongly affected by subjective “corrections” or “manipulation” [17] , however. Remove these improper “calibrations” and the values change to +0.45 and +0.65 mm/yr, respectively [18] .
In this paper, I will focus on shore morphology and stratigraphy, the latitudinal differences in sea level changes during the last 500 years, the occurrence of a +60 - 70 cm high sea level in the equatorial regions during the Little Ice Age in the 17th century, a low sea level (peat and buildings below present sea level) in the warm period of the 18th century, and the concept of “rotational eustasy” and its interpretation with respect to physics. The concept of “rotational eustasy” has been identified and successively built up on the basis of extensive field observations in the Indian Ocean and the Pacific [19] . In the true sense of science, it presents the logical steps from observations, via interpretations, to conclusions.
A complementary review of the observational facts is presented separately [19] , and an additional a special account on biological and shore morphological criteria in [20] .
2. Material and Methods
The concept of rotational eustasy is a direct product of the accumulation of new observational facts from a number of sites in the Indian Ocean and the Pacific. In this case, I refer to the following personal field works: 1) extensive and detailed investigations by an international team of specialist of a large number of islands in the Maldives in 2000-2005, 2) a study of the Sundarban delta in Bangladesh in 2010, 3) a detailed and multidisciplinary study of the sea level changes during the last 500 years in Goa (India) in 2011, 4) a detailed documentation of the sea level changes during the last 500 years at 10 sites in the Yasawa Islands of the Fiji Islands in 2016, and 5) a study of the sea level changes in Ouvéa Island of New Caledonia in 2018. This implies a very extensive database consisting of about 100 individual observational sites and some 80 C14-dates. Each site has been presented separately in peer reviewed scientific papers. A summary is presented below (Section 3), and an extensive review is presented in [19] . In a separate paper [20] , I have addressed evidence from biological criteria (trees, corals, etc.) and shore morphology, because these data do not lie and hence seem irrefutable.
In addition to my own field documentations and sea level reconstructions, I want to refer to additional investigations of 709 atoll islands in the Pacific and Indian Ocean summarized by Duvat [21] . It refers to the changes in size and area of the individual islands over the past decade (Figure 1). It is based on independent studies by [21] - [34] . This database is interesting because it is totally objective with respect to my own. As illustrated in Figure 1 hardly any atoll has decreased in size (as would have been the case if sea level really was rising), and the vast majority have remained stable during the last decade, indicating stable sea level conditions [35] .
This is consistent with a modern understanding of coastal dynamics [35] which implies that the shore (land/sea interface) can be deformed both vertically by sea level changes (besides crustal movements and compaction) and horizontally by a number of dynamic factors (Figure 2). Consequently, atoll islands my change their size laterally despite no changes in sea level.
My personal investigations in the Indian Ocean and the Pacific concerns the reconstruction of the vertical changes in sea level in Figure 2 (i.e. eustasy vs tectonics and compaction), whilst the studies summarized in Figure 1 [21] primarily concerns the horizontal (lateral) changes (i.e. the dynamic changes). The absences of a rapid rise in sea level and any tendency of a recent acceleration have also been noted by others (e.g. [36] [37] ).
3. Observational Facts
The concept of rotational eustasy is a direct product of the accumulation of new observational facts from a number of sites in the Indian Ocean and the Pacific. In order to not include an extensive review of all the field data in this paper, I limit myself to a short review of the main findings, and present a review of the observational material separately [19] .
Studies in the Maldives revealed that the 18th century was represented by a sea level below the present level with peat below sea level, dated at 1720-1790 BP
Figure 1. Changes in size of 709 atoll islandsin the Pacific and Indian Ocean (from [35] , redrawn in upside-down view from [21] ), providing additional evidence that the equatorial eustatic component has remaind stable for the past decades.
Figure 2. Coastal dynamics deform the shore laterally (horizontally), whilst the changes in sea level or land level deforms the shrelevel vertically according to Encyclopedia of Coastal Science [38] .
( [19] [39] [40] [41] , Figs 3-6). During the 17th century, sea level was 50 - 60 cm higher than at present. This was our first record suggesting opposed sea level changes in the Indian Ocean with respect to the general climatic changes [42] .
The next place of sea level investigation was the coast of Bangladesh ( [43] ; [19] , Fig 9). The 18th century is represented by a distinct sea level low [39] [41] [43] with salt kilns in submarine position [44] .
In Goa, India, it was possible to reconstruct the sea level change during the last 500 years with a high precision ( [19] [41] , Figs 5-8). The 17th century is represented by clear morphological evidence of a +60 cm higher sea level (with a C14-date giving 1615 ±12 AD). The 18th century, on the other hand, is represented by a low sea level with buildings (dated 1700-1750 AD) now below sea level. The records are clear and straightforward: a +60 cm high in the 17th century followed by a −50 - 100 cm low in the 18th century. Archaeological data and an ancient painting back up the geological data.
Consequently, the sea level changes in the equatorial Indian Ocean record a sea level story that is opposite to the general climatic history of a cold 17th century and a warm 18th century.
In the Fiji Islands in the Pacific, we recorded a +70 cm high sea level in the 17th century and a −100 cm low sea level in the 18th century ( [19] [45] [46] , Figs 10-14); i.e. a picture quite similar to the one observed in the Indian Ocean.
The latest record comes from Ouvéa Island of New Caledonia ( [19] [35] , Figs 15-18). There are very excellent records of a +70 cm sea level in sub-recent time. This level is fresh and young (sub-recent) and recorded over the entire island. It is also seen in a sub-recent sandy shore level preceding at +70 cm. Figure 3 shows a former shore cliff with under-cut notches located 70 cm above the present high-tide level (HTL). The continuity between the present and +70 cm shore marks indicates a sub-recent age of the +70 sea level. But most impressive is a rock-cut shore platform 70 cm above the present HTL (Figure 4). It is fresh, cut into the old reef deposits and lacking any sign of weathering, indicating a young or sub-recent age ( [19] , Figs 15-16). This level must correlate to the +70 cm sea level recorded in the Fiji Islands and dated to the 17th century (just as was the case in the Maldives and in Goa, India).
In Figure 5, I summarize the observed sea level changes during the last 500 years as observed and documented in Bangladesh [43] , Goa [41] , the Maldives [39] , Fiji [46] , New Caledonia [35] and the present trend on a number of Pacific
Figure 3. Double shore notches and undercut HTL shores on Ouvéa Island (photo taken at high tide). The present high tide level (HTL) is represented by extensive under-cut notches (here used as local ±0.0 cm datum). In direct continuation of the present beach there is a dead cliff with under-cut caves representing a former HTL located 70 cm above the present one. The age is “sub-recent” in perfect agreement with the +70 cm sea level recorded in the Fiji Islands and dated at the 17th century [45] [46] .
Figure 4. A very clear and distinct rock-cut platform on Ouvéa Island levelled at +70 cm above the present HTL. The platform is fresh and un-weathered and cut into the older strongly weathered reef masses. The age must be sub-recent, fitting well with the +60 - 70 cm high sea level in the 17th century recorded and dated in Fiji, Goa and the Maldives.
Figure 5. Sea level change documented during the last 500 years in key-sites in the Indian Ocean and the Pacific [19] [41] [48] . Red stars represent an Indian Ocean tsunami event in 1733 [49] .
islands [47] . It shows 1) a high sea level during the 17th century of generally cold climate, Little Ice Age conditions and the Maunder Grand Solar Minimum, 2) a low sea level during the 18th century of generally warm climate and a Grand Solar Maximum, 3) a rise and high sea level in the early 19th century of generally cold climate, Little Ice Age conditions and the Dalton Grand Solar Minimum, and 4) a fall in sea level in the 1960-1970 when global climate changed from cooling to the present global warming. This implies changes directly opposite from what one would expect from the glacial eustatic (and thermal expansion) point of view (cf. Figure 18 in [19] ). Therefore, we must seek a new explanation.
4. Interpretation in Terms of Physics
The observations (Figure 5; [19] ) documenting a high sea level during the cold period of the 17th century, a low sea level during the warm period of the 18th century, and a rise to a moderately high sea level in the early 19th century of cold climate conditions can only be understood in terms of variations in Earth’s rate of rotation; speeding-up during Grand Solar Minima and slowing-down during Grand Solar Maxima [50] [51] [52] [53] .
This implies that water masses are dislocated from the north to the equatorial bulge at the speed-ups during Grand Solar Minima. Arctic water is pulled down all the way to central Portugal initiating Little Ice Age climate conditions in Europe and the North Atlantic [51] . At the same time the equatorial bulge increases (by about 60-70 cm in the 17th century and by about 20-30 cm in the early 19th century as shown in Figure 6).
At Grand Solar Maxima, the Gulf Stream penetrates far up into the Arctic providing warm and nice climate conditions in the whole of NW Europe [51] . The equatorial bulge, however, decreased by up to 70 - 100 cm in the 18th century, and by 10 - 20 cm at about the climate cooling in 1960-1970.
These changes in sea level are a novel finding for which I have coined the term “rotational eustasy” [5] [42] [45] . It refers to the north-south redistribution of water masses and the increase and decrease in the equatorial bulge as a function of changes in the Earth’s rate of rotation (Figure 6).
If temperature continues to rise, sea level is likely to rise by about 10 cm up to year 2100 in the northern hemisphere and remain more or less stable in the equatorial region. If, on the other hand, we will face a new Grand Solar Minimum by about 2030-2050 [51] [52] [54] , sea is likely to fall in the north and rise in the equatorial regions (Figure 6).
For the period 1850-1930, there may have been a general glacial eustatic rise in sea level as a function of the temperature rise and ice melting after the Little Ice Ages (in Figure 6). In this case, the rise would have an internal cause as proposed before [4] [55] .
The geographical distribution of changes in sea level and ocean circulation (or oscillation) at the Grand Solar Minimum (the Maunder Minimum) in the 17th century and at the Grand Solar Maximum in the 18th century are shown in Figure 7 [5] [42] [45] .
Figure 6. Latitudinal differentiation in sea level changes during the last 400 years. During the 17th century, the Sun experienced the Maunder Grand Solar Minimum, climate was cold with Little Ice Age conditions, and the Earth’s rate of rotation speeded up generating sea level fall in the north and a significant sea level rise in the equatorial region. The 18th century, on the other hand, was a period of a Grand Solar Maximum with generally warm climate conditions and the Earth’s rate of rotation slowed down generating sea level rise in the north and a significant sea level fall in the equatorial region. In the early 19th century, the Sun experienced the Dalton Grand Solar minimum, climate was cold (a new Little Ice Age), Earth’s rate of rotation speeded up and sea fell in the north and rose in the equatorial region. These changes give evidence of a predominance of rotational eustasy.
Figure 7. Sea level changes (blue ± signs), ocean circulation (arrows), temperature changes (red ± signs) in the equatorial Indian Ocean and Pacific, and in the Atlantic and Pacific during Grand Solar Maxima (a) with rotational slowing-down and Grand Solar Minima (b) with rotational speeding-up.
In the northern hemisphere, the changes in rotation between Grand Solar Maxima and Minima are revealed by the beat of the Gulf Stream in the North Atlantic and the Kuroshio Current in the North Pacific [51] . In the equatorial region, the same changes in rotation are seen in the rise and fall in sea level; i.e. rotational eustasy [5] [19] [42] [45] . The Pacific Decadal Oscillation (PDO) and the North Atlantic Oscillation (NAO) are parts of the Gulf Stream Beat (GSB) in the North Atlantic and the Kuroshio Current Beat (KCB) in the North Pacific [51] [53] .
The expansion and contraction changes of the equatorial water bulge (Figure 4) between the 17th and 18th century must be a function of a speeding-up and a slowing-down of the Earth’s rate of rotation. This is in full agreement with the observed changes in the Gulf Stream beat [51] [52] . The observed expansion/contraction changes are illustrated in Figure 8.
The changes in Earth’s rate of rotation generating rotational eustatic changes in sea level are not internally but externally driven [56] . The changes in rotation were first observed in the beat of the Gulf Stream [57] . Later, it was understood that this beat (variability) ultimately had to be driven by changes in the Solar Wind [50] ; i.e. solar variability. The changes in Solar Wind emission (and simultaneously also in solar irradiance emission) were found to be the function of planetary beat (the 60-yr cycle, the 84-yr Gleissberg cycle, etc.). Figure 9 gives an integrated system of the changes of planetary beat on the Sun, the Earth and the Earth-Moon system [42] [53] [58] [59] .
Scafetta [60] showed that the 60-yr cycle is a prime cycle in solar variability and that it exhibits excellent correlation with the planetary cycle of the combined effect of 5 Jupiter (59.30 yr) and 2 Saturn (58.90 yr) rotations around the Sun. The same applies for the Gleissberg cycle, which is the combined effect of 7 Jupiter (83.02 yr), 3 Saturn (88.35 yr) and 1 Uranus (84.3 yr) revolutions around the Sun. Neptune (164.79 yr) and Pluto (248 yr) should also be integrated in the beat systems.
The planetary cycles have a direct effect on the Sun’s position with respect to the centre of mass of our planetary system [61] [62] , and the internal motions of planetary-Sun barycentre [59] [63] . Obviously, the planetary beat cycles drive both the 60-yr cycle and the Grand Solar Maxima/Minima cycles (besides longer-term cycles). The impact on the terrestrial system is manifold as illustrated in Figure 9.
Of course there are other planetary cycles, both longer and shorter, that are likely to affect ocean circulation and sea level changes. The Gulf Stream exhibits 30 pulses in the last 13,000 years with cycles varying in length between 230 and 1000 years [48] . This paper, however, is about rotational eustasy during the last 500 years.
Figure 8. Observed changes in height of the equatorial water bulge between a +70 cm size during the Maunder Grand Solar Minimum in the 17th century and a −70 cm size during the Grand Solar Maximum in the 18th century.
Figure 9. Integrated terrestrial effects (including rotational eustasy) as driven by planetary beat on the Sun, the Earth and the Earth-Moon system. Brown boxes give processes that are involved in the changes in sea level and ocean “oscillations” (3). The 60-yr geomagnetic field cycle (2) provides evidence of a back-ground in changes in Solar Wind [50] . The changes in Earth’s shielding (1) are a well-known characteristic of Solar Wind interaction with the Earth’s geomagnetic field.
The remarkable thing with respect to rotational eustasy is that the observed changes in sea level during the last 500 years now can be traced back to an origin in the motions of the planets around the Sun and the multiple planetary beats on the Sun, the Earth and the Earth-Moon system as illustrated in Figure 9.
5. Conclusions
The concepts of glacial eustasy and thermal expansion/contraction cannot explain the sea level changes observed in the equatorial belt of the Indian Ocean and the Pacific. This calls for an alternative explanation.
The explanation presented is termed “rotational eustasy” [35] [42] [45] . It implies that the water masses are re-located between the high latitudes and the equatorial region in response to changes in Earth’s rotation. The changes in rotation originate from the planetary beat on the Sun generating variations in emission of Solar Wind which interacts with the Earth’s magnetosphere, and direct planetary beat on the Earth and the Earth-Moon system [5] [42] [45] [59] as illustrated in Figure 9.
During Grand Solar Minima, the Earth’s rate of rotation speeds up and the equatorial bulge increases generating sea level rises in the Pacific and Indian Ocean.
During Grand Solar Maxima, the Earth’s rate of rotation decreases and the equatorial bulge decreases generating low sea levels in the equatorial Pacific and Indian Ocean.
The process may also be seen as a new type of “lunar-tidal super cycles”. Besides the Grand Solar Maxima/Minima cycle, we also have the 60-yr cycle, which is a basic climatic-eustatic cycle [42] with an origin in the Jupiter-Saturn motions around the Sun [60] .
The latitudinal changes in sea level during the last 400 years are summarized in Figure 6 and Figure 7. The north-south differentiation is obvious. The driving forces can only be understood in terms of rotational eustasy (Figure 8).
The concept of rotational eustasy is, of course, also valid for periods prior to the last 500 years here analysed. I leave this for future analyses, but just note that the sea level curve of the Maldives [39] [40] ) includes a number of high amplitude changes in sea level. The changes are in the order of 0.6 to 1.2 m and they are likely to be the function of changes in rotational eustasy (because of their amplitude, rates and anti-correlation with climatic-eustatic records from northern latitudes).
The concept is, of course, also valid for the future. Therefore, it is to be expected that sea level will rise in the equatorial regions during the next Grand Solar Minimum likely to occur at about 2030-2050 [54] .
Acknowledgements
When I in 1999 was elected president of the INQUA Commission on Sea Level Changes and Coastal Evolution, I initiated an international research project on sea level changes in the Maldives. We were soon able to document a number of surprises; a 10 - 20 cm fall in sea level at around 1970, submarine peat from the 18th century and a +60 cm high sea level from the 17th century. This was the beginning of our understanding that glacial eustasy and thermal expansion were concepts not able to explain our observations. After additional documentation in Bangladesh, Goa, Fiji and New Caledonia, it is now time for this presentation of the concept of rotational eustasy. I declare no conflict of interest. I thank all colleagues involved and all local persons so generously assisting us in the field. I am indebted to Pamela Matlack-Klein for skilful linguistic review of both this paper.