Changes in Barents Sea ice Edge Positions in the Last 440 years: A Review of Possible Driving Forces

This is the first report of the Barents Sea Ice Edge (BIE) project. The BIE position has varied between latitude 76°N and above 82°N during the last 440 years. During the period 10,000 to 6000 years ago, Arctic climate was significantly warmer than today. We review various oceanic and atmospheric factors that may have an effect on the BIE position. The Gulf Stream beat with respect to alternations in flow intensity and N-S distribution plays a central role for the changes in climate and BIE position during the last millennium. This occurred in combination with external forcing from total solar irradiation, Earth’s shielding strength, Earth’s geomagnetic field intensity, Earth’s rotation, jet stream changes; all factors of which are ultimately driven by the planetary beat on the Sun, the Earth and the Earth-Moon system. During the last 20 years, we see signs of changes and shifts that may signal the end of the late 20th century warm period. The BIE position is likely to start advancing southward in next decade.


Introduction
When the warm Gulf Stream water in the North Atlantic Drift meets the cold polar water in the Norwegian Sea and in the Labrador Sea, it sinks forming the North Atlantic Deep Water, which is a part of the huge global thermohaline circulation system (Figure 2). This is a global system of surface water circulation and deep-water circulation. According to [34] "heat and freshwater fluxes at the ocean's surface play a key role in the forming of ocean currents, which in turn  The global thermohaline circulation system with the approximate areas of down-welling of warm surface water into cold deep water in the Arctic and Antarctic marked by orange dots (modified from [35]). International Journal of Astronomy and Astrophysics Seidov et al. [42] for the last 60 years.
Grand Solar Maxima and Minima affect the Earth's shielding capacity as a function of the Solar Wind interaction with the Earth's geomagnetic field (e.g. [43]). The shielding controls the penetration of particles responsible for the formation of 14 C in the atmosphere, the atmospheric 10 Be content and cosmic ray particles [44]. The variations in atmospheric 14 C content are well documented and provide a continual record of the variations in Solar Wind intensity and its effect on the Earth's geomagnetic field strength as recorded by the changes in Earth's shielding capacity [37] [38] [45]. Figure 4 gives the changes in 14 C production during the last millennium (redrawn from [46]). Weak shielding and high 14 C production are shown downwards with the names of the Grand Solar Minima given in blue. During these periods the North Atlantic experienced type-c (or b) circulation ( Figure 3). Strong shielding and low 14 C production (upwards) occurred when the Sun experienced Grand Solar Maxima and the North Atlantic had a type-a circulation ( Figure 3). Filtering of the 14 C production curve (red line in Figure 4) gives a main period of 210 years (grey line in Figure 4), which corresponds to the main solar cycle identified by Abreu et al. [47] for the last 9400 years.
Easterbrook [44] [48] has documented a close correlation between Grand Solar Minima (Wolf, Spörer, Maunder, Dalton, 1880Dalton, -1915Dalton, , 1945Dalton, -1977 and climatic records of severe cooling events. Mörner [37] showed that the Spörer, Maunder and Dalton Minima are very clearly identified in ocean temperature records indicating type-c circulation ( Figure 5). Even the 1930-1940 warming in the north and cooling in the south were followed by a cooling in the north and a warming in the south ( Fig. 9 in [37]) indicates a change from type-a to type-b/c circulation.  14 C production in the atmosphere (red curve) 950-1950 as a function of the Solar Wind induced changes in Earth's shielding capacity (strongly modified from [46]). Blue fields mark decreased shielding trends and proposed cooling periods of grand solar minima. Yellow fields mark increasing shielding trends and proposed warming periods of grand solar maxima. The 1970-2000 warming is missing because the records end in 1950. The thin grey line shows a filtering period of ~210 years with scale to the right.  . Ocean circulation and Gulf Stream beat during the Spörer (type-c), Maunder (type-b followed by type-c) and Dalton (type-c) Grand Solar Minima (from [37]). Data in orange refer to the time of temperature changes recorded in the database used by Mörner [37]. Figure 4 suggests that the Medieval Warm Optimum may have lasted from 1050 to 1250, which is significantly longer than originally proposed by Lamb (e.g. [49]).
The period AD 950 -1050 is an interesting and important period. It is well covered in Figure 4 with a warm peak around 970 and a cold peak around 1050 (the Oort Grand Solar Minimum in Figure 4). In Sweden, the late 10 th century was warm and the early 11 th century cool and exceptionally dry with sand drift and dying peat bogs over in northern Europe. Figure 6 illustrates the situation in Greenland [50] [51]. In the period AD 900 -1000, quite warm climatic conditions prevailed and there was a free sailing route from Iceland to Greenland. Already a century later icebergs and sea ice prevented a direct route to Greenland, International Journal of Astronomy and Astrophysics and the sailing route was displaced to the south. This was just in the time of the Oort Grand Solar Minimum ( Figure 4). By AD 1200 -1400 the sailing route was even more displaced to the south, and the Viking settlements on Greenland were left to become extinct. Obviously, this was due to the cooling events during the Wolf and Spörer Grand Solar Minima. Figure 6. The changes in sailing route to Greenland as a function of changes in climate and ocean circulation (modified from [50] [51]).
Isotopic data from the GISP2 ice core suggest multi-year intervals of low temperature especially in the early and mid-14 th century, presumably representing the coldest period in central Greenland for the last 700 years [52], among other things exposing the remaining Norse settlement in Greenland to disaster [53].
The growth of open system pingos in Disko Island (central W Greenland) also points towards multi-year cold intervals during the Little Ice Age [54] [55].
Analysis of the isotopic composition of ice samples from a Little Ice Age glacier-derived rock glacier in Disko Island [56] indicate a mean annual air temperature about 2˚C -4˚C below present values, presumably reflecting net ice accumulation during a cold interval of the Little Ice Age, tentatively dated to AD 1600 -1775.
The drift-ice records of Bond et al. [57] [58] and "deglacial oscillations in ocean circulation" [59] seem directly to add on to the Gulf Stream beat picture of Mörner [37]. Cores from both sides of the North Atlantic (off Newfoundland and off Ireland) record cyclic peaks in ice-rafted material from drift-ice (IRD), indicating cyclic alternations between two modes: 1) low IRD, warmer SST and a circulation as today (i.e. the type-a circulation of Figure 3), and 2) ice-rafting events with peaks in IRD, colder SST and a significant southward advection of colder and fresher water (i.e. the type-c circulation of Figure 3). As a conse- quence, the thermohaline circulation was reduced [57] [60]. This is consistent with the model of Mörner [51], where the surface water circulation affects the deep-water circulation, not vice versa as illustrated in Figure 7. Furthermore, this southward shift of colder, ice-bearing surface water far into the sub-polar North Atlantic had an ocean surface circulation origin, which was coupled with the atmospheric circulation over Greenland [57]. Consequently, we can safely assume that the observed cooling events recorded far south along the European coasts (Figure 3(c)) also extended over the central and western parts of the North Atlantic (as recorded by the drift-ice events). The correlation by Bond et al. [58] between periods of increased drift-ice and periods of reduced Solar activity (as recorded by cosmogenic nuclides which, in fact, refers to Solar Wind activity [38], not irradiance as claimed by Bard et al. [61] and Sejrup et al. [62]. This implies that the Gulf Stream pulses must represent Solar Wind effects on the Earth's total rate of rotation as proposed by Mörner [36] [37] [43] [63] [64]. The presence of an external rotational pulse was recently endorsed by [65]. Figure 7. Effects and feedback interaction of changes in Earth's rate of rotation (from [51]). Ocean surface circulation plays the dominant role in redistribution ocean-stored energy (as recorded by paleoclimate) and water masses (as recorded by sea level). The external forcing refers to the planetary beat on the Sun, the Earth and the Earth-Moon system [36] [43] [63]. Internal forcing comes from the redistribution of mass and the changes in Earth's radius due to changes in sea level [66] [67] [68].
Kerr [69] wrote about "a slow, multidecadal climate pulse that beats the Atlantic Ocean" (his AMO) and noted that already Bjerknes [70] wrote about "a surge of warm water up the Gulf Stream" in order to explain the early 20-century's warm climate over Northwest Europe. Similarly, the North Atlantic changes recorded by Levitus [71] had to be related to changes in the strength of the Gulf Stream. Furthermore, the major shifts in climate and sea level recorded over the last 150 years in Northwest Europe correspond with the major shifts in LOD (length of the day) [67] [68], indicating an interaction between rotation and Gulf Stream beat. In the Florida Straits, the water masses transported by the Florida Current (28)(29)(30)(31)(32) varied significantly through the last millennium with a general lowering of about 3 Sv during the "Little Ice Age" [72]. International Journal of Astronomy and Astrophysics During LGM the land ice had its maximum expansion and the polar front was located far down in the North Atlantic forcing the Gulf Stream along its southern branch (Figure 8(A)). In the Pacific, the Kuroshio Current was forced to lower latitudes as a function of a considerable speeding-up of Earth's rate of rotation due to the sea level lowering and decrease in Earth's radius [36] [68]; i.e. internal forcing ( Figure 7). On the Faroe shelf, grooves and furrows from grounded icebergs have been recorded at several 100 m depth indicating drift from the NE to the SW of polar water full of floating icebergs, some of which must have been of gigantic size to scratch the seafloor as deep down as 300 -500 m [73] (green arrow in Figure 8
Ocean circulation changes as presented in Figure 9 imply that the observed data may record planetary cycles (e.g. [1]), solar cycles (e.g. [92]), lunar-solar tidal cycles (e.g. [90] [95]) and lunar cycles [86] [96]. Figure 9. Integrated effects of planetary beat on solar variability in luminosity (e.g. electromagnetic radiation) and, via the Solar Wind, on a number of fundamental terrestrial processes, and, via direct effects on the Earth-Moon system, on gravity, rotation, wind, ocean circulation, and sea level changes and oceanic oscillation systems [43] [63] [64]. International Journal of Astronomy and Astrophysics  [97]) indicating the redistribution of ocean surface water between the Arctic and the Equatorial regions as proposed by Mörner [37] [38]. The alternating trend may be a signal for the transition between type-a and type-b/c circulation. (B) Argo temperature changes (˚C) 2004-2019 [97] [98]. Whilst the Arctic Oceans were cooling (down to 1300 m), the Equatorial Oceans were worming (in their uppermost 50 -100 m).
A recent manifestation of the Gulf Stream beat (i.e. the alternations between type-a and type-b/c circulation in Figure 3) and the concept of rotational eustasy [94] seems to be recorded in the ocean temperature variations between the Arctic from 2011 and the equatorial regions from 2012 as illustrated in Figure   10  The Argo temperature changes with depth from 2004 to 2019 [97] demonstrate that the Arctic Oceans were cooling considerably in the upper 1300 m at the same time as the Equatorial Ocean were warming in the upper 50 -100 m ( Figure 10(B)). This is consistent with a decrease in the northern branch of the Gulf Stream (i.e. a change from type-a towards a type-b or type-c circulation; Figure 3), and a decrease in the Arctic down-welling and thermohaline circulation ( Figure 2). The Arctic cooling recorded is likely-if continued (cf. Figure   10(A)) -to initiate a southward advance of the BIE position in the near future.
The temperature changes shown in Figure 10(A) and Figure 10(B) are linked to a considerable freshening of the North Atlantic after 2012 as excellently shown by Holliday et al. [99]. Their analysis is based on the Argo measurements, too.
According to them, the freshening was both rapid and extensive, and reached a low salinity level not recorded before in the last 120 years.
The concept of Gulf Stream beat [37] and the observations of latitudinal changes in temperature in the North Atlantic ( Figure 10(A) and Figure 10(B) at about 2011 has recently got an important support in observed changes in northward ocean heat transport and cooling (about 2˚C) and freshening in the Eastern Subpolar Gyre [100], which may be taken to signify the onset of a general change from type-a to type-c circulation.
Another important implication of Figure 9 scheme is that the solar impact on planet Earth primarily goes via heliomagnetic Solar Wind emission and not via total solar irradiance (TSI). When changes in the Earth's geomagnetic shielding capacity are involved, we can be sure that the main forcing function comes from the Solar Wind interaction with the magnetosphere and the Earth's rate of rotation ( Figure 9). Therefore, observed variations in 14 C and 10 Be are the manifestation of changes in the shielding capacity (Solar Wind) not in TSI (as proposed by, for example, Bard et al. [61], Scafetta [1], and Veretenenko and Ogurtsov [101]).

Temperature and Ice Conditions during the Holocene
The Last Glaciation ended at about 11,800 BP, which is the boundary between the Pleistocene and the Holocene and also the boundary between the Younger Dryas and the Preboreal. The warming following the last glaciation at the end of the Younger Dryas period 11.8 Ka ago was abrupt. In central Greenland, temperatures increased by 7˚C or more within 10 -30 years [103] [104] [105]. Most of the changes in wind-blown materials, snow and some other climate indicators registered by ice cores, including a doubling of snow accumulation, played out within a few years [106].
In the Northern Hemisphere, summer insolation peaked 10 -9 Ka ago when the last remnants of the large ice sheets in the northern hemisphere retreated rapidly. At that time the incoming solar radiation was approximately 8% greater than at present [107]. This may have been especially significant at high latitudes International Journal of Astronomy and Astrophysics  [126], and (c) according to Humlum [130] and Humlum et al. [87].
in summer when daylight persists for 24 hours [108]. Since then, the Northern Hemisphere summers have gradually seen decreasing incoming solar radiation. The peak insolation is reflected in paleoecological data from the Arctic, where the Holocene thermal maximum, resulting in northward movement of the tree-line and shrub tundra replacing herbaceous tundra north of the tree line, occurred close to the Pleistocene-Holocene transition, broadly in the interval between 13 -7.5 Ka BP depending on latitude [109].
The Holocene changes in glacier extension and temperature in the Svalbard region have been studied by several authors (e.g. [87] [110]- [128], all of which indicate an early Holocene with climatic conditions significantly warmer than today, and a late Holocene of cooling conditions as illustrated in Figure 11. Figure 11 gives the general changes in temperature in the Svalbard region during the last 12,000 years. Curve (a) gives the changes in water temperature at a depth of 100 m in the Fram Strait according to Werner et al. [129]. Curve (b) gives air temperature change in Svalbard according to Mangerud and Svendsen [126]. Curve (c) gives the major changes in ice conditions with arrows indicating the increase to a maximum at about 2300 cal. yrs BP (cf. above; Humlum et al. [87]), the major sea ice expansions in the periods 1660 -1720 and 1800 -1818 and retreat at around 1769 [130] (see Figure 12). The main trend is clear; a rapid rise in temperature in early Holocene followed by a temperature well above the present from 10,500 to 5500 cal. yrs BP (with a period of low temperature 8500 -8200 cal. yrs. BP) and values below the present from about 3000 to 100 cal. yrs BP. Curves (a) and (b) agree well in general trend, but differ in details. The Late Holocene is best covered by the data of curve (c), too few data for a detailed analysis of the last 2500 years, however. International Journal of Astronomy and Astrophysics In much of the Nordic Seas the mid Holocene period 6 -4 Ka BP encompasses a transition to reduced Atlantic Water influence and lower sea surface temperatures. The late Holocene period from ca. 4 Ka BP to the present is characterized by sea-surface cooling caused by decreased Atlantic Water influence in the surface waters of the Nordic Seas, including coastal Svalbard (e.g. [116]). Atlantic Water inflows to the Arctic Ocean may also have been affected by bathymetric changes elsewhere in the Arctic (i.e., Nares Strait, Bering Strait and other eastern Canadian Arctic Channels). Of these, the opening of Bering Strait about 10,000 cal. yrs BP might be the most important because it introduces relatively fresh Pacific waters into the Arctic Basin. This may have strengthened the meridional Atlantic Water circulation through the Fram Strait, at the expense of Barents Sea inflow [79]. Jennings et al. [81], working on the East Greenland Shelf, found evidence for late Holocene cooling from 4.7 Ka BP, associated with southward expansion of the Arctic sea ice and increased polar water influence. A general shift towards cooling and return to the present tundra environments at many Eurasian arctic and subarctic sites apparently started around 4.5 Ka BP [109] [131].
Additional evidence for Late Holocene glacier growth on Svalbard was presented by Snyder et al. [121]. In Svalbard permafrost near sea level was reforming shortly before 3 Ka BP [87]. Furrer [112] [113] found evidence in NW Spitsbergen for a late Holocene glacier advance, reaching maximum size around 2.3 Ka BP and again during the Little Ice Age. Röthe et al. [123] studied the Karlbreen glacier in northwestern Svalbard. They documented a small to absent glacier from 9500 to 4000 cal. yrs. BP, a significant expansion at 3600 BP, a glaciation peak at 2400 BP, a significant peak 1300 -1800 BP, a peak 700 -100 BP and LIA peaks in the period 100 -500 cal. yrs BP. Snyder et al. [121], working in western Spitsbergen, recovered laminated sediments in a core in the proglacial lake Linnévatnet (W Spitsbergen) to reconstruct the glacial history of a presently ice-free cirque. Terrestrial macrofossils from the core were radiocarbon dated to provide a limit on the onset of glaciation, and indicate that a nearby short glacial valley (a cirque) was not glaciated until in the Little Ice Age. The absence of glaciation during earlier Neoglacial periods suggests local climatic conditions "favouring glacier expansion were either of greater magnitude or of longer duration during the 'Little Ice Age' than at any other time during the Holocene". Humlum et al. [87]  Island of about 5˚C [135]. International Journal of Astronomy and Astrophysics In the late 16 th century, climate conditions were favourable and Dutch sailors attempted to find the northeast-passage to the Pacific. During one of the expeditions, Willem Barents in 1596 discovered Svalbard, which obviously then must have been accessible, indicating warmer climatic conditions (cf. peak at around 1600 in Figure 4). His discovery route is shown on a sea map from 1599. The discovery was challenged by the English, who claimed Svalbard was discovered by Hugh Willougby already in 1573. In our BIE series we have data back to 1579, almost 20 years before the discovery of Svalbard by Barents.
The whaling industry flourished in the 18 th century, and in the year 1769 there are firm evidence that the spring maximum ice edge was north of Svalbard, and even north of the 1966 position as indicated in Figure 12, and further discussed in Part 2 [136]. The 18 th century was, on the whole, a period of warmer climatic conditions in Northern Europe (cf. Figure 4).    Figure 12 map. High-resolution sediment cores from the Fram Strait [129] provide too imprecise data to record this, however. This is also the case with a record of the last 2000 years [138].

Observations Referring to the Last 150 Years
Only a limited number of Arctic meteorological stations were in operation be-  The Svalbard-northern Russia region was exposed to relatively warm autumns (SON), suggesting increased cyclonic activity, leading to retreat of sea ice in the Kara Sea. This may be explained by the onset of a prolonged weakened state in the continentally based Siberian High in this period, as suggested by Meeker and Mayewski [82]. Hanssen-Bauer and Førland [140], however, suggested that the early 20 th century temperature increase in the Norwegian Arctic cannot be explained as a result of changes in advection alone. It may therefore be suggested that increasing cyclonic activity along the northern coasts of Russia and Siberia by wind stress may have contributed to opening of the sea ice cover, releasing heat from warm water below to the atmosphere, and thereby providing an additional source of heat. In order to illustrate the potential effect of this on air-temperatures, the average heat loss over the ice covered Arctic basin is about 6 W/m 2 , while the heat loss to the atmosphere from the ice-free north eastern Greenland Sea averages is no less than 200 W/m 2 [141] [142].
Whatever the reasons for this warming may be, several very mild winters were recorded in W Greenland, most notably in 1929, where no land-fast sea ice formed at all and hunters were able to use kayaks throughout the whole winter [56]. At The most important new species in the Disko Island region was the Greenland cod, locally known as the big-headed cod, which presumably is identical to the Icelandic cod [56]. Most likely, larvae of the Icelandic cod spread to SE Greenland waters in the early 1920s, driven by an enhanced Irminger Current, and for some decades were able to reproduce and survive in SW Greenland waters [143].
Whatever the precise background, this species represented a very important resource in Greenlandic economy for about 15 -20 years, after which a decline set in, ending with a number of disastrous fishing seasons.
The changing ice conditions after 1920 also made the first navigation of the NE Passage without wintering possible in 1932 for the Russian trawler "Sibiryakow" [144]. Ten years later, in August 1942, the reduced sea ice extension prompted the German naval high command to order the heavy cruiser "Admiral Sheer" into the Kara Sea, to intercept suspected allied convoys from US and Canada with supplies to the Red Army, making use of the extraordinary open water conditions along the Russian and Siberian coasts at that time [145]. During this operation "Wunderland", famed "Sibiryakow" was regrettably sunk by "Admiral Sheer".
The background for the German concern was the fact that a German raider, the "Komet", without great difficulties passed through the Northern Sea Route to the Pacific Ocean in summer 1940 [146] [147]. However, the Allies never attempted to make use of the NE Passage during the war.
Special interest has often been attached to the Svalbard region, as this part of the Arctic apparently displays a high climatic variability, as was recognized by both Ahlmann [148] and Lamb [132], and later by IPCC [149]. Significant climatic variations in the Svalbard region during the 20 th century are well documented by meteorological data since 1911 [139] [150] [151]. As an example, a marked warming around 1920 within 5 years changed the mean annual air temperature at sea level from about −9.5˚C to −4˚C. This is one of the most pronounced increases in air temperature documented anywhere during the instrumental period. Later, from 1957 to 1968 the air temperature sharply dropped about 4˚C, with a subsequent more gradual increase towards the end of the 20 th century ( Figure 13). The high climatic variability of the Svalbard region probably International Journal of Astronomy and Astrophysics derives from three mechanisms. Firstly, the islands are located directly in the main transport pathway for air masses into the Arctic Basin [78] [152]. Secondly, Svalbard is located near the confluence of air masses and ocean currents with very different temperature characteristics. Finally, the variability is further enhanced by rapid variations in the sea-ice extent, coupled with both atmospheric and oceanic circulation [80] [153] [154] [155] [156].
Moving on to the period 1940-1965, in contrast to 1915-1940, this represents a period of cooling in many parts of the Arctic ( Figure 13). Again, the temperature change was by far most pronounced during the winter and weaker during other seasons. Especially Alaska, NW Canada, Russia and parts of Siberia were affected by this development. The North Atlantic sector with Scandinavia at that time had relatively cold autumns, suggesting an increase of northerly winds.
North Atlantic sea ice also increased in this period, and a greatly increased flow of the cold East Greenland Current brought polar water far south (cf. Figure 3(c)).
In several years (1965 and especially 1968 and 1969) more Arctic sea ice came to the coasts of Iceland than seen for fifty years [137]. In April-May 1968 and1969 Iceland was half surrounded by sea ice, which has not occurred since 1888. At the same time, the Greenland cod migrated to Iceland waters, and for a few years (1967)(1968)(1969)(1970)(1971) offset the declining stocks there. The 1960s also saw the abandonment of attempts at grain growing in Iceland, which were resumed during the previous relatively warm period ) after a lapse of some hundreds of years [137]. At the beginning of the 21 st century, grain (barley) was only grown at one farm near Skoga in southernmost Iceland, and only for making hay (H. Björnsson, pers.com. 2003).

The Last 50 Years
The final ca. 30 years of the 20 th century has been a period of renewed warming in many parts of the Arctic, lending observational support to the common notion of 'Global Warming' [93] [149] [157] [158] [159] [160]. Again, this warming has been most pronounced in the winter season and less so in other seasons. The late 20 th century warming has affected especially western Canada, Russia and Siberia.
This was accompanied by retreating sea ice in both the arctic North Atlantic sector, and in the Kara Sea, after reaching a maximum extension in the early 1970s (e.g. [161] [162]). A study by Winsor [163] concluded that the mean sea ice thick- An apparent paradox is that warming registered in the Atlantic water layer below the sea ice in the Arctic Ocean became noticeable in the early 1990s, whereas the major ice thinning occurred prior to 1990 [163], suggesting the dynamics of sea ice extension in the Arctic Ocean to be more complex than usually thought. Polyakov et al. [164] proposed that the Arctic Ocean has experienced an "atlantification" (i.e. increased role of Atlantic inflow) during the last decades. This seems contradicted by Klyashtorin et al. [165] [166]). This AW warming has slowed slightly since 2008. Ocean ventilation may have had pronounced effects on the ice edge position (a question to be returned to in Part-2 [136]).
Whilst the Arctic region, in general, records a warming during the last 25 years, the Eurasian region records a cooling [159]. According to Roy [93], this is the function of the weaker sunspot cycles in the last 40 years. The influence of the Sun on the Arctic winter climate goes via changes in Arctic Oscillation (AO).
During winter solar minima, there is a warming over most of the Arctic (in December-February), especially over the Kara and Barents Seas, but a cooling over Siberia and Scandinavia [93].

Life and Resources in the Arctic
The ice cover in the European Arctic Seas affects the biologic productivity and the marine life in many ways and on all trophic levels. This includes the marine lightscape [167], the total primary production [168] [169] [170], the timing and quality of the primary and secondary production [171] [172] [173] [174], and the total production of commercial and non-commercial fish, seals and whale stocks [175] [176]. In this paper, we present data on the sea ice cover in the European Arctic from 1579 to 2018 (below). After a warm period that characterized the Viking time (see Figure 6) the winter temperature decreased rapidly from about 1030 and stayed relatively low in the Svalbard area until 1890 [177]. Despite that the period 1550 to 1890 in general was a period with very low winter temperature in Svalbard area, there has been periods with higher temperatures and a large variation in the summer (August) ice cover over the period 1580 to 2016 [177] [178]. The periods with little sea ice cover are associated with high biological production, fishing activity and hunting for whales and seals.
From Iceland, we have fish landings statistics covering the period 1600-1862 [179]. It shows "fluctuations of 50-60 years" with deep lows in the 1690s, 1750s and 1810s. International Journal of Astronomy and Astrophysics

BIE in the Whaling Period
Due to the pioneering work of Torgny Vinje [180] on assessing the location of the summer ice edge between Svalbard and Frans Josef Land (20˚E -45˚E), based on log books of whalers, warships, early explorers, and scientist (e.g. [181]), and later airplanes and satellites, we have a record of the sea ice cover in this region, dating back to the 1580s. Falk-Petersen et al. [178]  The area north of Svalbard was called the Whalers Bay by the Dutch whalers, who carried out intense whaling here for more than 100 years. The peak of the whaling period, between 1660 and 1730, coincided with a period where the summer ice edge was north of 80˚N for extended periods. Interesting, the short period with the expansion of the ice edge between 1730 to 1750 also showed a reduced catch of the bowhead whale as demonstrated in Figure 14 showing the reconstructed aggregate catch [183]. Comparing Figure 5, Figure 12 and Figure   14, we find that the maximum catch of more than 3000 bowheads happened in 1701 (i.e. at the end of the Maunder Minimum) when there was a type-c ocean circulation with cold water all the way south to Portugal ( Figure 5). In historical time, it was the period with most spring ice in the North Atlantic and the sea temperature along the Norwegian coast was so low that the cod did not reproduce. The Whalers Bay in the Fram Strait is a high productive hotspot today, and we suggest that this also was the case in the period with the historical whaling from 1670 to 1790, and that the high primary and secondary production sustaining the large stocks of whales and seals in this area [178] [183] [184] [185]. Bowheads are slow swimmers that feed by filtering zooplankton out of engulfed water masses with their 4.5 -5 m long baleen plates. The main prey is the energy rich Calanus copepods such as Calanus glacialis and Calanus hyperboreus, but also other lipid-rich zooplankters [186] [187]. Weslawski et al. [184] estimated that bowhead whales consumed as much as 4 million tons of plankton per season from June to September in the Svalbard area at the peak of their abundance.

BIE in the Post Whaling Period
After the intense Dutch bowhead whaling period 1680 to 1790, the sea ice expanded rapidly, and large parts of the Nordic Sea and the Barents Sea were covered by ice during spring and summer until 1910. Then the ice retreated rapidly until the 1940s, before the ice expanded again in the 1960s and 70s. In the long period with extensive ice cover after 1790, there was a short period with less ice around the 1850s and in the decade that followed. In 1873 the cod appeared on the banks of Spitzbergen, and a fishery was established the year after. The cod appeared in dense aggregations where the cod were feasting up on sea butterflies (Limacina helicina) [189] [190]. The fishery in this period was done by hand-line from small rowboats. The intense cod fishery off Svalbard (Spitsbergen) lasted from 1874 to1882 [189]. Then there was a rapid retreat towards north by the ice edge between 1910 and 1940 ( Figure 15) and already in 1931 H. U. Sverdrup [181] with the submarine Nautilus did marine research as far north as 82˚N in the Fram Strait. In the 1930s there was an intense cod fishery off Svalbard with yearly landings of several million tons [191]. The Northeast Arctic cod stock in the Nordic Seas and the Barents Sea was strong during the mid-1950s, with a stock of about 4 -5 million tons. However, during the cold 1960s and 1970s, with an expansion of the sea ice in the Barents Sea (see Figure 15), the stock was International Journal of Astronomy and Astrophysics dramatically reduced to 1 million tons by 1964 [190] [192]. The cold period lasted to approximately 1985 in the Nordic Seas [169] [193], with persistent large ice cover in the European Arctic. This period falls together with a collapse both in the herring and in cod stocks in the Norwegian Sea and Barents Sea [190] [192]).
The expansion or contraction of distributional range of cod will often coincide with long-term temperature changes [194]. Northward expansion of the cod stock with increasing temperature and reduced ice cover has previously been documented during the warm periods of the late 1800 and during the 1920s and 1930s.
During the warming event in the 1920s and 1930s and associated decreasing ice cover, the cod changed its distribution from southern Greenland up to Disko Island [195][196][197]. The cod was present as far north as the Disko Bay until the 1970s when they disappeared following a rapid decline in seawater temperature, with a slight recovery after the 1990s [197].
After 1990 there has been a rapid increase in temperature and decrease in ice cover in the European Arctic, followed by a large increase in the commercial fish stock in the Nordic Sea and the Barents Sea. In the Whalers Bay a mesopelagic layer of fish and zooplankton was recorded with an increased presence of mammals north of Svalbard [176] [185] [198] [199].  The mechanism and phase lock of these events will be discussed in more detail in Part 2 [136]. Since the 179 yrs period is a stable period, the José cycle may confirm an ice minimum in 1580 and predict a future ice minimum in 2120.

The Atlantic Ocean Oscillations
There are a number of ocean oscillations systems like PDO (Pacific Decadal Oscillation), AMO (Atlantic Multidecadal Oscillation) and NAO (North Atlantic Oscillation), which sets the character of regional changes in oceanic environment. They are all included in the "stadium wave wheel" of Wyatt and Curry [202], indicating coupling and feedback interaction (Section 9.5 and Appendix 1). Easterbrook [44] [203] [65]. They [65] modelled SST based on LOD and INAO back to 1650 and found that the cold sea was dominating in the beginning of 18 th century, as shown in Figure 12, even if there were warm land temperature periods in Europe. International Journal of Astronomy and Astrophysics Another anomaly was the period 1900-1930 with an increasing positive NAO and a cold SST [205]. Lansner and Pedersen [206] showed that in the warm period 1920-1950, it was only stations with little ocean influence, which showed warming, while the warming at the end of the century took place both at ocean affected and shielded stations. The result is a slight cooling for inland stations 1900-2010 measuring −0.03˚C/century, while ocean affected stations showed a warming 0.78˚C/century. This is interpreted as a sign of a slow ocean warming since the cold LIA.
The AMO index ( Figure 16) is defined as the pattern of sea surface temperature variability when the linear trend is removed [207]. It has a 65 -70 years' quasi-periodic variation [208]. The updated AMO series in Figure 16 has a period of 67 years. The BIE position variation is roughly in phase with the AMO index. During the period 1900-1925 the SST was cold, AMO negative, and NAO posive, which shows that AMO is dominating over NAO [205]. Both Alheit et al. [209] and Drinkwater et al. [194] discussed the AMO variations and recorded its effects on fish faunas and North Atlantic temperature and ice-edge variations.
Mette et al. [156]] record a close correlation between AMO and changes in ocean temperature in northern Norway. Figure 16. The AMO oscillation in temperature variation from 1856 to 2018 (11 yr running mean) ( [207]; from http://www.esrl.noaa.gov/psd/data/timeseries/AMO/). [210] and Klyashtorin et al. [165], the stocks of herring and cod in the Barents Sea fluctuate with a 60-year periodicity (the cods lagging by 8 -10 years), which correlates with the changes in Arctic temperature, the ocean water temperature and the ice cover conditions in the Barents Sea, all exhibiting a ~60-year periodicity. Mörner [43] noted that the changes in ice cover reflect the "delivery of warm Atlantic water to the region" and that "the main source of heat delivered to the Arctic basin is warm water inflow from the North Atlantic Stream" as recorded by Nikolaev and Alexeev [211].

According to Klyashtorin and Lyubushin
It agrees with the observed changes in Earth's rotation [36] and Gulf Stream beat [37]. This is illustrated in Figure 17.

Alternative External Driving Forces
The Arctic sea ice grows in the winter, when it is cold and dark, reaching a maximum area in March/April. It then retreats to a summer minimum in late August or September. A recent re-calibration of Arctic sea ice extent datasets 1901-2015 ( Fig. 10 in [212]) shows that the Arctic summer sea ice cover was growing until about 1920, then retreating slightly until 1960, when there suddenly became more ice, which retreated slowly until 1990, then faster. At the same time the annual average surface air temperatures in the Arctic changed in opposite directions, indicating that air temperature variations may be the first cause for ice cover variations.
For the Whalers Bay area, the winter sea ice loss has since 1979 been larger than for the Arctic in general. The winter ice loss trend has been close to 10% per decade. The reason may be 0.3˚C warmer water in the Atlantic water pathway. Wind pattern may control the year-to-year changes, but influx of warmer water controls the decadal variations according to Onarheim et al. [213].
A new estimate of Northern Hemisphere and Arctic weather stations based air temperature variation, well supported by sea surface temperature trends, glacier length variations, and three rings proxy temperatures [214], suggest that surface  In order to investigate the proposed solar-Arctic climate connection on multidecadal to centennial timescales, longer time series have to be analysed. It is important to keep track of systematic patterns as: seasonal and annual cycles, local and regional air pressure systems: topography, landscape, and storage and exchange of heat/energy through atmospheric and oceanic circulation and glaciation; and delayed actions [217]. An example of a coupled interaction is the 13 years' delay of the North Pacific Multidecadal Oscillation (NPMO) to the Atlantic Multidecadal Oscillation (AMO) [218].
In addition to direct solar irradiation, the Arctic warming is related to the atmospheric heat transport from the tropics to the Arctic, which is modulated by TSI variations with little or no delay [219]. This results in a warming of the mid and high latitudes while the tropics do not warm much. One example is the recent temperature "pause" (or "hiatus"), which is dominated by a slight cooling in the equatorial regions while the high latitudes are warming [220]. In the following we will discuss various driving mechanisms covering millennial, centennial, multi-decadal and shorter time scales.

Millennial, Centennial and Multi-Decadal Solar and Planetary Cycles
The interaction between Earth, the Sun and other planets in the solar system results in variations in the Earth's orbit on many time scales. Secular changes in eccentricity, axial tilt and precession result in long climate cycles, also called Milankovitch cycles [221] [222], which in generally are held to be the main reason for the Quaternary Ice Ages. The rapid changes of cooling and warming periods observed are in phase with 10 Be variations, however, which makes it more likely that solar activity is the main reason for the alternations between Ice Ages and Interglacial, between stadials and interstadials, and short-term changes in climate [48].
Analysis of solar activity in a 9400-year long record of solar modulation potential φ determined from cosmogenic radionuclides such as 10 Be and 14 C shows distinct periodicities of 88, 104, 150, 208 and 506 yrs [47]. They find excellent agreement between the long-term cycles and periods in planetary torque, suggesting a physical connection (see Figure 4). It should be noted, however, that they did not detect a 60-yrs cycle in their analysis.  [223] show that in all series exists a strong cycle with duration 55 -60 years in the 18 th and 19 th centuries. It is less pronounced during the end of the 19 th and the beginning of the 20 th centuries. In all series except Rz there is also strong and stable quasi 100 -120 years and ~200 years cycles. In the Rz series a strong ~95 years period was present.
Also quasi 40 -45 and 30 years cycles are detected, but they are unstable during long time scales and have low amplitudes. An interesting difference is that the simulated aa-index shows a period 180 years while the Rh shows 202 years. Cliver et al. [224] found a strong correlation between decadal aa and global temperature values (r = 0.90) and conclude that 50% -100% of the temperature increase since Maunder Minimum is due to solar variability.
The solar dynamo theory [225] explains that in the beginning of a sunspot cycle, the dynamo process transforms the Sun's polar field into an internal toroidal field that increases during the maximum phase of the solar cycle, then affecting Rz, when completing the sunspot cycle the toroidal field transforms into a poloidal affecting the aa-index. However, analysing wavelet power spectra Duhau and Chen [226] find that Rz and aa undergo cyclical oscillations whose amplitudes vanished in 1923. After that the ratio between Rz and aa increases from 5.3 to 12.4 and the phase difference decreases from 0.12 to 0.05 radians.
They find that the secular Rz modulation steadily rise until the next Gleissberg International Journal of Astronomy and Astrophysics cycle maximum in 1971. The temporal change in the coupling function between Rz and aa is synchronous with a loss in the relationship between temperature and sunspot number in the period 1920-1955 [227]. For the secular changes in Rz 1844 -2000 Duhau and Chen [226] determined the following periods: 85.5, 121.1, and 171.2 years.
Finally, it may be significant to note that Blanter et al. [228]  In order to forecast future variations, it is necessary to identify stationary oscillations. This was done for the first time for TSI variations and group sunspot series by Yndestad and Solheim [92], who searched for stationary components using the wavelet autocorrelation method. They found stationary components of the solar variability controlled by the 12 years Jupiter period and the 84 years Uranus period and its 3/2, 5/3 and 9/2 harmonic components giving periods 125, 210, and 373 years. The TSI record of Hoyt and Schatten [229] showed only the stable periods of Jupiter, Uranus and Neptune.
According to Kuznetsova [230], the changes of the terrestrial climate during the last millennium are driven by the interaction of Jupiter, Saturn, Venus, the Sun and the Moon.
Le Mouël et al. [231], analysed the Ri series by singular spectrum analysis A global temperature harmonic analysis [236] based on paleotemperature reconstructions and instrumental records for the period AD 1 to 2015 gave three significant periods: ≈1000, ≈460 and ≈190 years in addition to a less significant period ≈60 years. A 190-year period was also found modulating the length of the sunspot cycles ([237] [238]).
Harmonic analysis of the HadCRUT4 global temperature (1850-2012) [239] showed that a linear trend of 0.5˚C per century and periods 9. During solar cycles with negative polarity of the Sun's northern polar field, the Galactic Cosmic Ray (GCR) variation has a peaked form. In the other phase it has a plateau. This is an effect of the differences in cosmic ray drift in the positive and negative phases of the magnetic cycle. Integrated GCR counts are higher in plateau cycles compared with peaked cycles [242]. This may be the reason for the amplification of the 22 years component in the global temperature curve.
The difference between 11 and 22 years climate response, is also seen in the latitudinal difference in the rhythm of growth in pine trees, (as shown in [239]; Fig.   6), where the 20 years period dominates north of 65 degrees latitude, while the 10 years period dominates at lower latitudes. This may be a result of differences in atmospheric circulation or effects of cosmic rays of lower energies reaching deeper at higher latitudes. For Svalbard at 78˚N, an analysis of temperature series by Humlum et al. [88] shows that periods 17 and 26 years are much stronger than those at 10 -12 years.
Kasatkina et al. [243] showed that dendrochronological data from the Kola Peninsula and Finnish Lapland recorded the effect of low-latitude volcanic activity with "significant decrease in polar tree-ring growth over 7 years after eruption", and "the main cycles of solar activity (11, 20 -25 and ~100 years)".
The 66 yr period can be explained as a harmonic solar (3 * 22 = 66 yr) or lunar (3.5 × 18.6 = 65.1 yr) period. A maximum entropy power spectra analysis by Scafetta [244] showed that 9 different global and regional temperature series Periods mentioned above may well be present in global time series, but the sensitivity of the local (high latitude) climate is different from low latitude sensitivity, and the combination of thermal-salinity-cryospheric interactions, also involves lunar and orbital forcing, acting different at seasons. A framework for the solar-lunar-orbital interaction for high latitudes is presented by Cionco et al. [246] and demonstrated in Figure 19. The latitudinal temperature gradient (LTG) drives the (wind driven) ocean circulation between the tropics and the polar regions. Analyses of recent changes in LTG suggest that the Earth's climate system is unusually sensitive to changes in LTG [247]. LTG varies according to seasons. In the NH summer the variation in obliquity is dominating with 41 Ka orbital and 18.6-yr lunar cycle. In NH winter the TSI is dominating with 11yr cycle [90] [246]. The LTG is clearly modulated by the TSI variations [219] and represent a channel both for solar and lunar forcing of Arctic climate. Figure 19. Mean daily irradiance at 65˚N for the summer solstice based on a variable TSI and variation in the Earth's orbit due to interaction from the Moon and other planets (modified from [246]). The Spörer, Maunder and Dalton Minima are easily observed, like the 1920-2000 high and a coming minimum 2020-2085.

Lunar Cycles
Treloar [95] identified that lunar-solar forces produce multiple tidal ocean ef- Correlations are also found between the lunar nodal tide effect and temperature and salinity in hydrographic time series for the Faroe-Shetland channel and the Kola Section, suggesting that the temperature and salinity variations are forced by harmonics of the nodal tide [96] [250]. Dominating periods were the 18.6 yr lunar nodal period and its 9.3 year harmonics. Davis and Brewer [90] stated: "empirical evidence suggests that extraterrestrial forcing influences the Earth's climate" (a notion that we certainly agree with), and they presented data in the effort of "unifying orbital, solar and lunar forcing". It is interesting to notice that a major reorganisation in atmospheric circulation took place about 1420. During the preceding Medieval Warm Period, winter storminess across the North Atlantic was at a minimum. This is confirmed in Na + deposits in the GISP2 borehole in Greenland [252]. After 1420 until about 1925 cells of winter Icelandic low pressure have been strongly anchored within northern hemisphere circulation, which has resulted in increased frequency and severity of North Atlantic winter storms. This corresponds to a period of reduced June 21 solar Insolation at 65˚N (Figure 19). A reduction in winter storminess is documented for Grimsey and Stykkisholmur weather stations in Iceland after 1920 [253].

Changes in the Earth's Shielding Capacity
The Earth's magnetic field is modulated by the solar wind, which varies in phase with the solar activity or sunspot cycles, with an average period of 11.06 years.
At maximum solar activity the shielding capacity of the Earth's magnetic field against cosmic rays from exploding stars is maximum, and this reduces the number of cosmic rays that hit the Earth. Since cosmic rays may be related to formation of clouds [254], this has a climate effect. We can therefore expect less clouds and higher temperatures at solar maxima. The variations in the solar wind amplify the TSI-effect on the Earth's climate, and longer periods with reduced solar wind intensity are in general cold periods on Earth (cf. Figure 4).
A linear relation between sunspot numbers and the rotation of the Earth, or the length of the day (LOD), shows that increased solar activity or strong solar wind decelerates the Earth's rotation (Fig. 2 in [37]). This affects the oceanic and atmospheric circulation systems (cf. Figure 9), both of which have a significant effect of terrestrial weather and climate. The relations among solar activity, irradiance, solar wind, variations in the Earth's atmospheric shielding and variation in LOD are visualized in Figure 20.
The changes in Earth's shielding capacity are recorded in the changes in atmospheric 14 C production. Figure 4 gives the changes from AD 950 to 1950. The production is modulated by a 210-yr periodic fluctuation, which corresponds to the 208-yr period observed by Abreu et al. [47] during the last 9400 years. The    The figure shows that the solar wind has decreased more than 30% from the peak in 2004. Since 1990 the neutron counts (GCR) has increased and the solar wind speed decreased indicating that Earth's shielding capacity has decreased considerably since 1990 ( Figure 22). Bucha [75] demonstrated that increased corpuscular radiation (as recorded by the neutron counts) implied enhanced penetration of electrons into the middle atmosphere and stratosphere, resulting in increased temperature and atmospheric pressure in the stratosphere and troposphere in the auroral oval. This leads to a change from a meridional to a zonal type of atmospheric circulation over the Arctic. Visa versa, low corpuscular magnetic activity leads to low temperatures in the stratosphere. Further effects on the Solar Wind interaction with the terrestrial system are given in Figure 9

Internal Coupling and Feedback Processes
We have previously discussed possible external driving forces for the main Gulf Stream beat as well as for the BIE variations. In this section, we discuss internal coupling and feedback processes because they provide insights into the dynamics on the regional dimension.

The ENSO Events
The equatorial current lags behind the rotation of the solid Earth, which lead to

The Gulf Stream Pulsation and Earth's Rate of Rotation
Mörner [255] noted that the Gulf Stream exhibited a pulsation of 16 cycles within the Holocene, recorded in high-resolution deep-sea cores, sea level changes in the Kattegatt-North Sea region, and temperature changes in Northwest Europe.
He interpreted this in terms of a feedback coupling between water masses redistributed from low latitudes to middle-to-high latitudes, and changes in Earth's rotation in order to compensate for this. We now know, that the primary driving force is to be found in external planetary interaction with the Sun and the Earth The Gulf Stream beat system proposed by Mörner [37] has a close relation to the alternations between NAO+ and NAO− stages recorded by Wanner et al. [256]. Their main illustration (Fig. 9 in [256]) has been endorsed and reprinted by Pinto and Raible [257] and Boström [258]. Whilst the Gulf Stream beat theory has a driving energy input from external forces (Figure 7 and Figure 9), the NAO theory of Wanner et al. [256] is based on internal interaction and feedback process. It seems significant that  concluded that he had presented facts, indicating, "that the phases of two of the Earth's major climate systems, the North Atlantic Oscillation (NAO) and the Pacific Decadal Oscillation (PDO), are related to changes in the Earth's rotation rate". As shown by Mörner [36] [37] the changes in rotation originate from planetary-solar beat as illustrated in Figure 9, i.e. from external sources (Figure 7).

The Ocean Oscillation Interaction
In his account on "the Lunar nodal cycle influence on the Barents Sea", Yndestad [259] [260] presented both an internal and an external view of the climate feedback loops of the Arctic; i.e. the coupling and interaction between temperature, water, ice extent and light reflection in the Arctic region of the Atlantic sector.
He noted that his first lunar view model was good to illustrate the fundamental dynamics operating with respect to feedback loops in an autonomous dynamic system. This model represents the internal interaction and feedback processes.
He noted that this model was not enough, but called for external oscillating forces.
In his external forced oscillators model, he included solar/orbital influx (from the Milankovitch variables), Moon orbital forces (tides and rotation) and CO2 effects.
Today, we know that the Sun varies in quite another way, and that the planetary beat effects the Sun, the Earth and the Earth-Moon system ( Figure 9).
Yndestad's thesis [259] [260] is deeply relevant, because it at the same time as it shows the fundamental interaction and coupling between different internal parameters, it also calls for external forcing functions (today in improved understanding, however).

The Rotating Wine Glass Mechanism
The polar position oscillates with periods 6, 18, and 74 yrs cycles, all related to the 18.6 lunar nodal cycle.
Yndestad [259]  Yndestad [259] presents the detection of lunar nodal oscillation or harmonics in all Arctic climate time series analysed (Barents Sea ice extent, Kola section ocean temperatures, Hammerfest tidal variations, Greenland Sea ice extent, Røst winter temperature, and NAO winter index). He concludes that these stationary cycles may enable us to predict the future Arctic climate. International Journal of Astronomy and Astrophysics 9.5. The "Stadium-Wave" Mechanism Wyatt [261] and Wyatt and Curry [202] have presented a quite remarkable "Stadium-wave" mechanism of a special and temporal wave of 4 segments propagating across the Northern Hemisphere activating a network of synchronized climate variables. They identified four temporal groups (I -IV) changing between negative and positive phases. It offers a synchronization and feedback coupling between ocean, ice and wind factors. In Appendix 1, it is shown that the "stadium wave" mechanism of Wyatt and Curry [202] can be combined with the planetary-solar modulation, Gulf Stream beat and the external rotation forces here discussed.
In this paper, we acknowledge that the mechanism seems to imply a logical manifestation of "terrestrial changes in climate & environment" in Figure 9, but calls for external driving forces (as indicated in Figure 7 and Figure 9).

Discussion
The yearly cycle of freezing and melting of the Arctic Ocean ice cover depends on the variation of heat influx during the seasons. This is modulated by many different internal and/or external factors, as the direct solar radiation, the cloud cover modulated by cosmic ray particles and other particles (pollution or volcanic), transport of warm water by ocean currents and advection of air by the wind directed by pressure system, again modulated by the solar wind. Our working hypothesis is that the planetary gravity beat is somehow preserved as a weak signal through a chain of processes schematically shown in Figure 9 and in some cases amplified in local or regional synchronized systems: especially Earth rotation [36] [262] and the polar wobble [86]. Stream, which both are controlled by the solar wind. The solar wind also acts on the Earth's speed of rotation, which again forces the ocean currents (cf. Figure 7).
The summer Artic sea ice was nearly gone during the Holocene temperature maximum 10,500 to 5500 cal. yrs BP (Figure 11), but reappeared during the Late Holocene cooling period, which culminated with maximum ice in the period 1660-1720 as demonstrated in Figure 12. Thereafter a rapid warming took place in the 18 th century (in 1769 the ice edge was even north of Svalbard) before another cold period occurred around 1800-1818 (in agreement with stages 1 -3 in Figure 21

The Maunder and Dalton Minima
The Gulf Stream circulation system is strongly driven by the main alternations between Grand Solar Maxima and Minima (Figures 1-6, Figure 8, Figure 12, Figure 19). During the Maunder and Dalton Minima, the inflow of warm Atlantic water to the Arctic was strongly reduced, and cold Arctic water reached all the way south to central Portugal ( Figure 5). In the GIN Sea region, severe climate conditions and periods of extensive ice expansion occurred ( Figure 12).
Both periods are well recorded in the Solar wind intensity record of Figure 21 ( Another problem is that the high level of solar wind through the 19 th century did not move the BIE north. This did not happen until the start of the solar Grand Maximum in the 1920s, and was accompanied by a rapid temperature increase at many Arctic stations ( Figure 13). One possible explanation will be explored in Paper 2 [136]: The Sun was stressed after the producing short sunspot-cycles resulting in the 18 th century maximum, and it relaxed by producing International Journal of Astronomy and Astrophysics long cycles in the 19 th century. Only at the start of cycle 16 (1900) it was "refreshed" and able to produce short cycles [238]. This pushed the BIE north in the 20 th century, still modulated by the planetary cycles.
In the BIE position estimates ( Figure 14) we find minimum ice 180 years apart, corresponding to the José period of alignment of the planets in the solar system.
It may also seem significant that Yndestad and Solheim [92] were able to document the planetary beat of Jupiter, Saturn, Uranus and Neptune in TSI and sunspot datasets in addition to the recorded 210-yr (de Vries/Suess) cycle, which is the predominant cycle in the 14 C production record of the last 9400 years [47].

The 60-Yr Cycle
The most prominent feature in the BIE record ( Figure 15) are the two deep minima 1630-1660 and 1800-1900 and gradual move north after 1900. A series of maxima 180 years apart suggest a modulation following the José cycle of orbital synchronization. Since this is a sum of three 60-year cycles we may also find the 60-year cycle in the data (to be searched for in Paper 2 [136]).
Scafetta (e.g. [244] [266]) and Solheim [239] showed that the 60-yr cycle had an origin in the combined 60-yr beat of Jupiter and Saturn. Consequently, we must place planetary beat-in this case the combined beat by Jupiter and Saturn-at the top of the forcing (Figure 9). This beat simultaneously affects the Sun, the Earth and the Earth-Moon system. The solar variability generates changes in luminosity (irradiance) and Solar Wind, the last one of which controls a number of fundamental Earth's processes as illustrated in Figure 9.
The planetary beat effects on Earth's rotation and the Earth-Moon system imply that there must also be a generation of tidal effects. It therefore seems sig-  [276], and around 2040 [165]. The 60-yr cycle ( Figure 19) and the lunar cycles must interact with the Gleissberg (~84 years) and the de Vries/Suess (208 years) cycles as assumed in the concept of a New Grand Solar Minimum peaking some 20 -30 years in the future (above). The effect of the positive trend of a millennium long cycle [235] should also be considered. It seems significant, however, that the combined Gleisssberg-de Vries cycles provides a cyclic trend which agrees well with the observed changes in 14 C prediction and the Gulf Stream beat between type-a and type-c circulation patterns ( [38] [39]; Appendix 2).

Conclusions
We have presented a new 440 long record of the changes in BIE position ( Figure   15). It includes an observed ice maximum in the mid 17 th century, extreme north- During the last 200 years, the BIE record seems predominantly forced by a long-period maybe of millennium scale steady northward move, and shorter cycles, and the 60-yrs cycle, ultimately driven by the combined Jupiter and Saturn beat on the Sun, the Earth and the Earth-Moon system (Figure 9). This planetary beat seems also to be the external force needed to drive the "stadium wave" mechanism of Wyatt and Curry (202) as further discussed in Appendix 1.
We note signs suggesting the ending of the present warm period; e.g.  Figure 10(B)) and the North Atlantic is recording a freshening [99].
As for the future changes in climate and environment of the Arctic, there are two scenarios: 1) a linear, maybe even accelerating, warming leading to ice-free conditions in 2061-2088 according to Onarheim and Årthun [160], or 2) cyclic changes, in response to planetary beat leading to a New Grand Solar Minimum at about 2030-2050 with colder climatic conditions. Our data are in favour of the second scenario. The climatic evolution during the next couple of decades will provide the final test.
Finally, we note with pleasure that the bowhead whale (Balaeana mysticetus), once feared to have become extinct (Figure 14), has now reappeared. International Journal of Astronomy and Astrophysics

Appendix 1 (by Nils-Axel Mörner)
The "stadium wave wheel" mechanism in need of external forcing Wyatt [1] and Wyatt and Curry [2] have presented a quite remarkable "Stadium-wave" mechanism of a special and temporal wave of 4 segments propagating across the Northern Hemisphere activating a network of synchronized climate variables. They identified four temporal groups (I -IV) changing between negative and positive phases. It offers means of synchronization and feedback coupling between ocean, ice and wind factors, but it lacks driving input energy forces. Therefore, I analysed the system in view of our paper on the changes in Arctic sea ice edge positions during the last 440 years with respect to planetary-solar modulation and Gulf Stream beat.
Wyatt [1] proposed that the variables of the Northern Hemisphere surface tem-   In Figure 24, we compare the temporal zonation given by Wyatt and Curry I therefore propose that the stadium wave mechanism by Wyatt and Curry [2] of a number of inter-connected internal factors of coupling and feedback, in fact, is ultimately driven by external forces going back to a planetary beat ( Figure 9 in the main paper). There are strong arguments for seeking a primary forcing of AO in geomagnetic changes [6], and of AMO, NAO and PDO in external rotation [7] [8] [9]. International Journal of Astronomy and Astrophysics

Appendix 2 (by Nils-Axel Mörner)
The Gulf Stream Beat Theory: Origin and Perspectives The close correlation observed among the changes in Gulf Stream beat [1] [2], the changes in a combined curve of the Gleissberg and de Vries cycles (Fig. 5 in [2]) and the observed changes in atmospheric 14 C production ( [3] [4]; Figure 4), which is a measure of the strength of the Earth's geomagnetic shielding capacity, provide a chain of facts indicating that solar variability driven by planetary beat controls the Earth's geomagnetic field strength driving the latitudinal beat of the Gulf Stream. Additional data [5] indicate that there is a close correlation between Gulf Stream beat and Earth's rate of rotation. The Solar Wind changes observed ( Figure 21; [6]) are indicative of the solar-terrestrial interaction [7] [8].
This implies that we see the following linkage between solar variability as seen in the alternation between Grand Solar Maxima and Minima), Earth's rotation and the beat of the distribution of water masses by the Gulf Stream ( Figure 25). Over the last 800 years the Gulf Stream has switch between circulation patterns of type-a (at Grand Solar Maxima) and type-c (at Grand Solar Minima) in  ). Observed periods of type-a and type-c Gulf Stream circulation patterns are marked with dots and label 2 (type-a) and 4 (type-c) according to Mörner [1]. Vertical scale gives temperature in centigrade.
The changes in Gulf Stream beat ( Figure 26) also correlates well with recorded changes in 14 C production, indicating a correlation between changes in the Earth's shielding strength and the ocean circulation changes ( Figure 28). Figure 28. Comparison between Gulf Stream circulation patterns according to Figure 27 (above) and a curve of 14 C production (below) by Bard et al. [9] termed "solar irradiance" though it, in fact, is a curve of geomagnetic shielding primarily controlled by the Solar Wind effects [1] [2] [7] [8]. S, M, D and F stands for Spörer, Maunder, Dalton and Future Grand Solar Minima. The corresponding alternation between type-a (orange wedges) and type-c (blue wedges) North Atlantic circulation is given in the top graph (cf. Figure 26). Figure 30 gives an update version of the geomagnetic shielding during the last 450 years and Figure 4 covers the last 1000 years.
In 1996, I suddenly realised that the Earth's rate of rotation ( Figure 29) was strongly controlled by external changes in the Earth's geomagnetic shielding strength due to interaction with the Solar Wind [5]. The external forcing has later been confirmed by Mazzarella [10], and Mazzarella and Scafetta [11].    14 C production since AD 1500 (a longer record is given Figure 4).
Whilst the cycles from 1450 up to 1790 are in the order of 80 -90 years, the cycles after 1790 are in the order of 60 years. This is in agreement with the LOD record ( Figure 29) and the Solar Wind record (Figure 21). It also agrees with the Gulf Stream changes between type-a and type-c circulation patterns.
In conclusion, this means that we see the following chain of cause-and-effect behind the observed cyclic pattern in changes of North Atlantic circulation and its stepwise origin in causal changes in rotation, shielding, solar wind intensity, solar variability and planetary beat as shown in Figure 31 (cf. Figure 9): Figure 31. Chain of cause-and-effects from Gulf Stream Beat back to its ultimate origin in Planetary Beat (cf. Figure 9) [14] [15].