Fluctuations of the southern sea ice edge and the Marginal Ice Zone in the Arctic Ocean in the North Atlantic occupy a key function for the understanding of climate changes. The aim of this paper is to investigate a possible connection between the orbits of the planets and Earth’s climate as documented in our revised 440 years long BIE position record. This is an old hypothesis, based on the complex synchronization structure of the solar system, which since Pythagoras of Samos (ca. 570 - 495 BC) is known as the music of the spheres . As science evolved with Kepler’s elliptical orbits and Newton’s mechanics, and in particular with the discovery by Christian Huygens in the 17th Century, that entrainment or synchronization between coupled oscillators requires very little energy exchange if enough time is allowed, the hypothesis of a connection is strengthened. Since the solar system is about 5 billion years old, there is plenty of time to synchronize.
Synchronization is observed in a large number of numerical relations in the solar system . Since short periods may be washed out by stochastic noise in the climate system it is important to find long time series where the long periods may be preserved. In this paper we discuss what we can learn from a 440 years long series of estimated positions of Barents Sea Ice Edge (BIE).
In this series of papers, we present analysis of a revised data set, covering 440 years, of the estimated position of the late August ice edge position in the Barents Sea Ice Edge (BIE). The data set is unique in the sense that it covers the two last grand solar minima, the Maunder Minimum (1640-1720) and the Dalton Minimum (1790-1820) , the first one of which includes the coldest period of the last millennium in the Northern hemisphere (1690-1700) according to .
The Barents Sea Ice Edge (BIE) project was launched in 2015. The present paper, termed Paper 1, provides an introductory review of active variables and potential forcing functions. The second report, termed Paper 2, provides a more detailed analysis of correlations with other data sets and possible drivers of the BIE position variation.
The present paper begins with a general outline (Sections 2 and 3) of the geological and paleoclimatical background. The Gulf Stream plays a central role and is discussed in Section 4 with respect to its temporal and spatial variations during the last millennium. Section 5 discusses multiple integrated effects. Temperature and ice conditions during the Holocene are presented in Section 6. The life and resources in the Arctic including a graph of the BIE changes during the last 440 years are given in Section 7. Section 8 is devoted to a discussion of possible external driving forces. Section 9 covers internal coupling and feedback processes. The paper ends with a discussion in Section 10, and the conclusions drawn in Section 11.
2. North Atlantic Ocean
The opening of the North Atlantic Ocean began in the Early Eocene (e.g.  ). Mörner  gave a general account of the major Cenozoic events in the North Atlantic with respect to plate tectonics, crustal-lithospheric movements, paleoceanography and paleoclimatology. As stressed by Hopkins , the GIN Sea (i.e. the Greenland, Iceland and Norwegian Seas) “plays a unique and vital role in the circulation of the world’s ocean, in coupling the Polar Sea (Arctic Ocean) to the Atlantic Ocean”.
The rate of modern sea floor spreading is in the order of a couple of centimetres per year . The Faroe-Iceland Ridge underwent a stepwise long-term subsidence ending in uplift in the Mid-Pleistocene . The Fennoscandian Shield was suddenly uplifted at about the Oligocene/Miocene boundary around 23 Ma (Ma = million years), initiating “The Baltic River” debouching in the Netherlands   and slump deposits west of Svalbard  ending in a considerable subsidence of the Baltic region in the last million years. The Barents Shelf experienced uplift pulses at about 38 Ma, 22.5 Ma and 5 - 4 Ma, ending in a subsidence in the Mid-Pleistocene. Svalbard was uplifted at 22.5 Ma.
In the Paleocene to Eocene, Svalbard had a temperate climate  . The temperate flora on Svalbard was originally assigned a Miocene age , which posed some paleoclimatological problems.
At about 3 - 4 Ma (million year), the Panama Isthmus closed   , which initiated the formation of the Gulf Stream bringing warm, high-salinity equatorial water to high latitudes in the North Atlantic (cf. below) and the Arctic cross-polar flow. By this, regional evaporation increased, especially from the ocean between Greenland, Svalbard and Norway, leading to increasing precipitation over the adjoining land areas, partly in the form of snow, resulting in significant glacier growth.
At about 3 Ma, there was a short period of quite temperate climatic conditions with a boreal forest in the Cape Copenhagen Formation of northeast Greenland. This period represents a global warm peak and global sea level high stand spanning the Kaena Reversed Geomagnetic Polarity zone within the Gauss Normal Polarity zone  .
At the Gauss/Matuyama boundary at about 2.6 Ma   or rather at during Marine Isotope Stage (MIS) G10 at around 2.8 Ma    the alternation between continental glaciations and interglacials commenced. In Europe the sequence of continental glaciations started with the Pretiglian Glaciation at around the Gauss/Matuyama boundary and in the early Matuyama    . Jakobsson et al.  presented a useful overview of the glacial history of the Arctic Ocean during the last 150,000 years. The probable impact from previous interglacials on our understanding of present-day climate conditions has been assessed by Kandiano et al. . The deglaciation of the Barents Sea region is covered by   .
4. The Gulf Stream
The Gulf Stream sets the character of the North Atlantic (e.g.   ). From the Gulf of Mexico, the Gulf Stream follows the American east coast up to New Foundland where it turns east and splits into two main branches (Figure 1); the northern branch (known as the North Atlantic Drift) and the southern branch (known as the Canary Current). The northern branch (NB) reaches all the way up into the Barents Sea and the Arctic Ocean, making northwestern Europe mild and inhabitable. The southern branch (SB) affects the Gibraltar area, northwest Africa and the Canary Islands.
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  “heat and freshwater fluxes at the ocean’s surface play a key role in the forming of ocean currents, which in turn
Figure 1. The Gulf Stream and its division into a northern (NB) and southern (SB) branch.
Figure 2. 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 ).
have major effect on climate” and he finishes by concluding: “understanding its past and future behaviour is crucial to our understanding of climate change”. This is, of course, exactly the reason why we have undertaken the present study of the ice edge behaviour and its forcing functions.
Figure 3. The beat (pulsation) of the Gulf Stream as a function of changes in the Earth’s rate of rotation   redistributing ocean water masses and ocean stored heat along the northern (a) and southern (c) branches of the Gulf Stream. The type-a circulation is typical for Grand Solar Maxima, and the type-c (and type-b) is typical for Grand Solar Minima.
The Gulf Stream Beat
The Gulf Stream splits at about 40˚N Lat. and 30˚W Long. into a northern branch (NB) and a southern branch (SB) as illustrated in Figure 1. The distribution of oceanic water along those two branches is not constant with time but subjected to a “beat” or pulsation in its intensity and especially in the alternation of water distribution along the two branches    as illustrated in Figure 3. When the northern branch increases in intensity and warms Scandinavia and the Arctic, the southern branch decreases and cools the Gibraltar-Canary area (Figure 3(a)). When southern branch increases and warms the Gibraltar-Canary region, the northern branch decreases and Scandinavia and the Arctic cools (Figure 3(b) & Figure 3(c)). The forcing of the alternations in beat along the northern and southern branches comes from changes in Earth’s rate of rotation  ; deceleration in the first case (Figure 3(a)) and acceleration in the second case (Figure 3(c)).
The type-a circulation is typical for Grand Solar Maxima and is recorded at the known warm peaks at about AD ~970, ~1250, ~1320, ~1590 - 1600, ~1720 - 1790, and after 1900. The type-b and type-c circulations are typical for Grand Solar Minima and is well recorded and documented at the “little ice ages” around AD ~1050, ~1300, ~1450, ~1690, ~1810, and there are strong reasons to expect that it will re-occur at about AD 2040   .
The alternations between type-a and type-c circulation have wide environmental effects in the North Atlantic. Significantly, Yasuhara et al.  find a correlation between cooling events and periods of reduced Atlantic Meridional Overturning Circulation (AMOC) during the last 10,000 years, Thibodeau et al.  between the Little Ice Age and a slowdown of the North Atlantic circulation, and Seidov et al.  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. ). The shielding controls the penetration of particles responsible for the formation of 14C in the atmosphere, the atmospheric 10Be content and cosmic ray particles . The variations in atmospheric 14C 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   .
Figure 4 gives the changes in 14C production during the last millennium (redrawn from ). Weak shielding and high 14C 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 14C production (upwards) occurred when the Sun experienced Grand Solar Maxima and the North Atlantic had a type-a circulation (Figure 3). Filtering of the 14C 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.  for the last 9400 years.
Easterbrook   has documented a close correlation between Grand Solar Minima (Wolf, Spörer, Maunder, Dalton, 1880-1915, 1945-1977) and climatic records of severe cooling events. Mörner  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 ) indicates a change from type-a to type-b/c circulation.
Figure 4. Changes in 14C 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 ). 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.
Figure 5. 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 ). Data in orange refer to the time of temperature changes recorded in the database used by Mörner .
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. ).
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 10th century was warm and the early 11th century cool and exceptionally dry with sand drift and dying peat bogs over in northern Europe. Figure 6 illustrates the situation in Greenland  . 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, 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  ).
Isotopic data from the GISP2 ice core suggest multi-year intervals of low temperature especially in the early and mid-14th century, presumably representing the coldest period in central Greenland for the last 700 years , among other things exposing the remaining Norse settlement in Greenland to disaster . The growth of open system pingos in Disko Island (central W Greenland) also points towards multi-year cold intervals during the Little Ice Age  . Analysis of the isotopic composition of ice samples from a Little Ice Age glacier-derived rock glacier in Disko Island  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.   and “deglacial oscillations in ocean circulation”  seem directly to add on to the Gulf Stream beat picture of Mörner . 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 consequence, the thermohaline circulation was reduced  . This is consistent with the model of Mörner , 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 . 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.  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 , not irradiance as claimed by Bard et al.  and Sejrup et al. . 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     . The presence of an external rotational pulse was recently endorsed by .
Figure 7. Effects and feedback interaction of changes in Earth’s rate of rotation (from ). 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   . Internal forcing comes from the redistribution of mass and the changes in Earth’s radius due to changes in sea level   .
Kerr  wrote about “a slow, multidecadal climate pulse that beats the Atlantic Ocean” (his AMO) and noted that already Bjerknes  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  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)  , indicating an interaction between rotation and Gulf Stream beat. In the Florida Straits, the water masses transported by the Florida Current (28 - 32 Sv) varied significantly through the last millennium with a general lowering of about 3 Sv during the “Little Ice Age” .
Figure 8. Changes in polar front position (red lines) and Gulf Stream distribution (hatched flow arrows) during the main stadials and interstadials of the last 25 Ka (thousand years): (A) at LGM, (B) at 18 Ka, (C) at 14 - 16 Ka and (D) at 12 - 13 Ka (modified from ).
In Figure 8, we compare the polar front positions and Gulf Stream distribution at the major stadial/interstadial changes of the last 25 Ka (from ); viz. at LGM (Last Glacial Maximum, 22 - 24 Ka), the 18 Ka glacial maximum, the Bölling/Alleröd (BÖ/AL) interstadials at about 16 - 14 cal. Ka BP and the Younger Dryas (YD) stadial at around 13 - 12 cal. Ka BP (Ka BP = thousand years before present).
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  ; 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  (green arrow in Figure 8(A)).
At the Pomeranian Stadial around 18 Ka, the polar front was just south of Iceland. The Gulf Stream seems to have penetrated up to the edge of this front (Figure 8(B)).
During the Bölling-Alleröd Interstadials, the Gulf Stream penetrated all the way up into the Barents Sea (Figure 8(C)). A boreo-arctic mollusc fauna entered the North Sea and the Kattegatt. At the Younger Dryas Stadial the climatic conditions returned to arctic and glacial conditions. The polar front and the Gulf Stream were displaced to the south (Figure 8(D)). This was also the case with the Kuroshio Current in the Pacific   indicating increasing rate of Earth’s rotation. The drastic changes in ocean water distribution, climate and rotation in association with the BÖ/AL to YD periods pose complicated forcing  .
5. Multiple and Integrated Effects
The changes in ice edge positions within the Svalbard and Greenland regions during the last 400 - 500 years—the topic of this project—must be the function of multiple variables (e.g. Bucha  ; Kellogg ; Dickson et al. ; Forman et al. ; Humlum  ; Jennings et al. ; Meeker and Mayewski ; Vinje and Goosse ; Soon ; Divine and Dick ; Yndestad ; Humlum et al.  ; Sejrup et al.  ; Davis and Brewer ; Solheim et al. ; Yndestad and Solheim ; Roy ; Årthun et al. ). Because the Gulf Stream beat records a close correlation with the Grand Solar Maxima/Minima alternations (above), there are reasons to consider the multiple effects of the planetary beat on the Sun, the Earth and the Earth-Moon system as illustrated in Figure 9 (from ). The observed changes in Gulf Stream beat, North Atlantic Oscillation (NAO), Pacific Multidecadal Oscillation (PDO), other oceanic “oscillations” system and the new concept of rotational eustasy seem to lead their origin in processes within Figure 9 system .
Ocean circulation changes as presented in Figure 9 imply that the observed data may record planetary cycles (e.g. ), solar cycles (e.g. ), lunar-solar tidal cycles (e.g.  ) and lunar cycles  .
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   .
Figure 10. (A) Global ocean temperature 0 - 1900 m according the Argo measurements (reworked from ) indicating the redistribution of ocean surface water between the Arctic and the Equatorial regions as proposed by Mörner  . 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  . 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  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(A).
The Argo temperature changes with depth from 2004 to 2019  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. . 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  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 , 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 14C and 10Be are the manifestation of changes in the shielding capacity (Solar Wind) not in TSI (as proposed by, for example, Bard et al. , Scafetta , and Veretenenko and Ogurtsov ). The so-called 60-yr cycle  must primarily be a function of heliomagnetic Solar Wind emission  .
6. 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   . 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 .
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 . This may have been especially significant at high latitudes
Figure 11. Temperature changes in the Svalbard region during the last 12,000 years; (a) according to Werner et al. , (b) according to Mangerud & Svendsen , and (c) according to Humlum  and Humlum et al. .
in summer when daylight persists for 24 hours . 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 .
The Holocene changes in glacier extension and temperature in the Svalbard region have been studied by several authors (e.g.   - , 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. . Curve (b) gives air temperature change in Svalbard according to Mangerud and Svendsen . 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. ), the major sea ice expansions in the periods 1660 - 1720 and 1800 - 1818 and retreat at around 1769  (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.
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. ). 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 . Jennings et al. , 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  . Additional evidence for Late Holocene glacier growth on Svalbard was presented by Snyder et al. .
In Svalbard permafrost near sea level was reforming shortly before 3 Ka BP . Furrer   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.  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. , 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.  described in situ subglacial vegetation beneath the glacier Longyearbreen. Dating of the plants recovered demonstrated a glacier advance of about 2 km during the last 1100 years.
Professor Hubert Lamb stands out as one of the key persons in the understanding of climate change in historical times (e.g.     ). Besides his excellent analyses on the changes in temperature, wind and ocean circulation, he also gave vivid examples on the cooling during the Little Ice Ages in the form of the sudden appearance of Eskimos in kayaks at the coast of Scotland and the drop of ocean surface temperature between Iceland and the Faroe Island of about 5˚C .
In the late 16th 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 18th 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 . The 18th century was, on the whole, a period of warmer climatic conditions in Northern Europe (cf. Figure 4).
Figure 12 by Humlum  provides an expressive paleoenvironmental map of the Northeast Atlantic region during the 1660 - 1720 period with ice positions compared with 1769, 1899 and 1966. The period of cool climatic conditions corresponds to the Maunder Minimum shown in Figure 4.
Figure 12. Palaeoenvironmental map of the Northeast Atlantic region (from ) during the Little Ice Age in the 17th century with the spring maximum ice edge position in 1660-1720 and comparative positions in 1769, 1866 and 1966.
Figure 12 map is based on various written sources, including Lamb      . Red arrows indicate prevailing wind directions. Symbols for animals and humans with black cross indicate reported periods of absence/extinctions. Symbols for ice crystals indicate periods with reported periods of glacier growth. Yellow symbols (Europe) indicate periods with reported coastal flooding and/or strong coastal sand drifting.
The ice edge displacements documented in Figure 12 imply a manifestation of the beat between the cold East Greenland Current and the warm West Spitsbergen Current in pace with the Gulf Stream beat (Section 4); viz. an intensification of the East Greenland Current (and decrease of the West Svalbard Current) and a northward displacement of the ice edge during Grand Solar Maxima, and an intensification of the North Atlantic Current and the West Spitsbergen Current (and decrease of the East Greenland Current) and a southward displacement of the ice edge during Grand Solar Minima. Both the Maunder Minimum cold phase and the Dalton Minimum cold phase are well recorded in Figure 12 map. High-resolution sediment cores from the Fram Strait  provide too imprecise data to record this, however. This is also the case with a record of the last 2000 years .
6.1. Observations Referring to the Last 150 Years
Only a limited number of Arctic meteorological stations were in operation before 1910, and in few places have observations been carried out before the late 19th century. In Fennoscandia, the oldest systematic climatic observations north of 65˚N were in Tornio (Finland) during a period from 1737 to 1749, while regular weather stations were initiated around 1850. The first station in northern Russia was at Arkhangelsk in 1813. Early meteorological observations were initiated in Greenland in the period 1866-1895. Meteorological observations in the Arctic Ocean began at the voyage headed by F. Nansen on-board the “Fram” in 1894-1896 and were again taken up during the forced drift in the ice of other ships in the 1930s. The first real meteorological station in the High Arctic was Green Harbour (Svalbard) in 1911. A new phase of the Arctic data coverage started in the mid-1930s with North Pole Ice stations that continuously operated in the Arctic Ocean up to 1991. Economic reasons later caused a disappointing reduction in the number of existing Arctic meteorological stations during the 1990s, especially in northern Russia and Canada. The available meteorological records, however, clearly conveys the impression of high Arctic temperature variability; both on an annual basis and between individual regions in the Arctic.
Long Arctic surface air temperature records (Figure 13) illustrate conditions during the last cold interval of the Little Ice Age, which was terminated by a marked increase in mean annual air temperature around 1920. This rise was mainly derived from rising winter temperatures, while the temperature increased less in other seasons, especially summer. The termination of the Little Ice Age was followed by a warm period, recorded at many Arctic meteorological stations, continuing to around 1940-1950, where after a widespread cooling commenced. In some parts of the Arctic this cooling culminated around 1970-75 and was followed by a corresponding temperature recovery towards the end of the 20th century. In some regions, such as Greenland and Iceland, cold conditions prevailed almost to the end of the 20th century, while in other regions (e.g. Svalbard) temperatures have been rising without interruption since the cold period around 1970 .
Figure 13. Long Arctic surface annual air temperature series (redrawn from ). The thin blue line represents the mean annual air temperature, and the thick blue line is the running 5-year average. Data source: NASA Goddard Institute for Space Studies (GISS) and http://www.rimfrost.no/.
The period 1915-1940 (Figure 13) represents a period of general warming in the Arctic especially influencing northern Alaska, northern Canada, Greenland and the northern parts of Russia and Siberia. The warming was by far most pronounced during the winter season (DJF) and least during the summer (JJA). 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 . Hanssen-Bauer and Førland , however, suggested that the early 20th 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/m2, while the heat loss to the atmosphere from the ice-free north eastern Greenland Sea averages is no less than 200 W/m2  .
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 . At the same time new fish species appeared in the central W Greenland region, causing a significant change in Inuit hunting habits towards open sea fishing. 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 . 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 . 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” . 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 . 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  . 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  and Lamb , and later by IPCC . Significant climatic variations in the Svalbard region during the 20th century are well documented by meteorological data since 1911   . 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 20th century (Figure 13). The high climatic variability of the Svalbard region probably derives from three mechanisms. Firstly, the islands are located directly in the main transport pathway for air masses into the Arctic Basin  . 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     .
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 . In April-May 1968 and 1969 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-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 (1920-1950) after a lapse of some hundreds of years . At the beginning of the 21st 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).
6.2. The Last 50 Years
The final ca. 30 years of the 20th century has been a period of renewed warming in many parts of the Arctic, lending observational support to the common notion of ‘Global Warming’      . Again, this warming has been most pronounced in the winter season and less so in other seasons. The late 20th 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.  ). A study by Winsor  concluded that the mean sea ice thickness in the Arctic Ocean remained constant from 1986 to 1997. Satellite mapping showed the late summer 2012 sea ice extension in the Arctic to be the smallest since initiation of satellite mapping in 1978. However, the sea ice expanded again in 2013 and 2014, and despite the season 2018-2019 showed the next smallest ice-covered area, there is no trend in September ice-cover 2007-2019.
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 , suggesting the dynamics of sea ice extension in the Arctic Ocean to be more complex than usually thought.
Polyakov et al.  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.  who documented an increased Atlantic inflow through the Faroe-Shetland Strait from 1910 to 1970, a fall to low levels in 1980 and continued low levels at least up to 2000. Humlum  recorded a drop in the North Atlantic Current temperature from 2004 to 2019 with a peak-low in 2016. This can hardly be combined with an “atlantification” in the last two decades (nor is the general decreasing trend in warm water to the Arctic since 2011; Figure 10(A)).
According to Polyakov et al.  “the role of Atlantic water heat in sea-ice reduction is not negligible”. Observations from the 1990s and 2000s documented two warm, pulse-like Atlantic water (AW) temperature anomalies on the order of 1˚C (relative to the 1970s), entering the Arctic through Fram Strait and occupying large areas of the Arctic Ocean. The strength of the 2000s warming peaked in 2007-2008, with no analogy since the 1950s (cf. ). 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 ).
Whilst the Arctic region, in general, records a warming during the last 25 years, the Eurasian region records a cooling . According to Roy , 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 .
7. 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 , the total primary production   , the timing and quality of the primary and secondary production    , and the total production of commercial and non-commercial fish, seals and whale stocks  . 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 . 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  . 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 . It shows “fluctuations of 50-60 years” with deep lows in the 1690s, 1750s and 1810s.
Figure 14. Reconstructed aggregate harvest of bowhead whales from the Greenland-Spitzbergen stock (from ).
7.1. BIE in the Whaling Period
Due to the pioneering work of Torgny Vinje  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. ), 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.  updated this time series to the present and related it to the primary and secondary production in the Whalers Bay north of Svalbard. They hypothesized that in periods with little ice the total bio-production in the European Arctic increased due to two main factors: a) the increase of the total area of open water available for primary production during spring and summer in the Barents Sea and Nordic Seas, and b) the winter upwelling of nutrient rich Atlantic water along the shelf break in the area north of Svalbard and Frans Josef Land fuels production hotspots in Whalers Bay as soon as the Sun returns. The high primary production attracts the lipid-rich Arctic grazing zooplanktons such as Calanus glacialis, Calanus finmarchicus and Calanus hyperboreus, the key species in the high latitude lipid-based ecosystems  .
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 . 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    . 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  . Weslawski et al.  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.
The documentation on the “French Naval Operations in the Spitsbergen during Louis XIVs Reign” further illustrates the productivity of this area in the Dutch whaling period. The commercial value of whaling activity was high and the commercial output of interest for several European nations . In 1693, the French king sent a fleet of warships north of Svalbard to take the Dutch whaling fleet as prizes. On the 6th August 1693, the two French frigates Agil and Favori encountered and took a fleet of more than 40 ships, all flying the Dutch standard, under fire at the mouth of the Hinlopen strait (La Baye aux Ours—today Sorgbukta). The French frigates sized 13 ships and sunk a large part of the Dutch ships. The remaining Dutch whaling fleet escaped by sailing northeast in the ice-free area in the Hinlopen Strait. Despite the losses, the Dutch ship owners still were able to make a profit. The Dutch whaling stopped in some periods due to naval conflicts in Europe, for instance 1665-1667, 1691 and 1710 .
7.2. 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)  . The fishery in this period was done by hand-line from small rowboats. The intense cod fishery off Svalbard (Spitsbergen) lasted from 1874 to1882 . 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  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 . 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 dramatically reduced to 1 million tons by 1964  . The cold period lasted to approximately 1985 in the Nordic Seas  , 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  ).
The expansion or contraction of distributional range of cod will often coincide with long-term temperature changes . 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-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 .
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    .
Figure 15. The estimated position of the BIE 1579-2018 (black dots) and an 11-point running mean (red curve). Some important events (marked 1-11): 1) 1611-51 approx. 250 whales/yr, 2) 1662-71 increase to 2000/yr, 3) 1689-1696 reduced whaling, 4) 1701 max catch 3180 whales, 5) 1723-1735 reduced whaling, 6) 1750-1800 declining catch, 7) 1831 commercial whaling ends, 8) 1879-1882 rich cod fisheries in Svalbard waters, 9) Eastern Arctic bowhead assumed extinct, 10) 1955-1960 strong Barents Sea cod stock, 11) 1990-cod recovery, the bowhead is back.
Figure 15 shows the revised BIE position estimates for the period 1579-2018 (black dots) and an 11-point running mean. Some important biological events are marked with numbers (Figure 15, points 1 - 11). The red curve documents the presence of all the major solar maxima/minima known for the last 440 years (Figure 4); viz. little ice in the late 16th century solar maximum, ice expanding from 1600 to 1650 corresponding to the first half of the Maunder Minimum. The retreat of BIE in the second half will be discussed in Part 2 . Then we observe less ice during the 18th century solar maximum, expanding ice from 1780 to 1800 referring to the onset of the Dalton Minimum, less ice in a short period around 1850, and rapid retreat of the ice edge 1910-1940. Finally, we observe more ice 1940-1970, and a retreat 1970-2010. Since 2010, except for the record low ice in 2012, the summer ice cover has not shown any significant decrease (http://nsidc.org/arcticseaicenews/). We may also notice the rapid expansion of the ice in the years 1614-1616, 1664-67, 1690-95, 1780-85, 1864-69, and 1960-62, which may be due to special weather events. Apart from the cycles, the dominating part of the curve is the gradual movement of the BIE northwards for more than 200 years (1800-2018).
Another interesting point is that the ice minima in 1760 and 1940 are 180 years apart. This may be related to the José cycle of 179 years , which is the fundamental orbital repetition period of the Jovian planets and is also synchronized with the sidereal orbital periods of Venus, Earth, Mars and Jupiter as well as the Schwabe (11.1 yrs), Hale (22.2 yrs) and Gleissberg (90 yrs) periods . The mechanism and phase lock of these events will be discussed in more detail in Part 2 . 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.
7.3. 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 , indicating coupling and feedback interaction (Section 9.5 and Appendix 1). Easterbrook   noted that glaciers on Mt. Baker in Washington State showed a regular pattern of advances and retreats, which matched the PDO variations. This has a bearing on the effects of changes in the North Atlantic (Gulf Stream, NAO, AMO) with respect to the changes in ice-edge position on the Barents Sea region.
The NAO index is defined as the normalized pressure difference between subtropical high and Icelandic low pressure. Its connection to the Arctic Ocean oscillation (AO) is well described by D’Aleo and Easterbrook ( , Fig. 11.16) indicating alternations between stages of zonal jet stream and warm positive NAO conditions, and meridional jet streams and cold NAO conditions in the North Atlantic. The NAO-index was increasing from 1895 and had maxima between 1908-1926 and 1985-1997 and negative values 1955-75. It is about 70 years between its two maxima. The Earth’s rotation (LOD) is negatively correlated with the integral function of NAO (INAO), and also with global sea surface temperature (SST) . They  modelled SST based on LOD and INAO back to 1650 and found that the cold sea was dominating in the beginning of 18th century, as shown in Figure 12, even if there were warm land temperature periods in Europe. Another anomaly was the period 1900-1930 with an increasing positive NAO and a cold SST . Lansner and Pedersen  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 . It has a 65 - 70 years’ quasi-periodic variation . 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 . Both Alheit et al.  and Drinkwater et al.  discussed the AMO variations and recorded its effects on fish faunas and North Atlantic temperature and ice-edge variations. Mette et al.  ] 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) ( ; from http://www.esrl.noaa.gov/psd/data/timeseries/AMO/).
According to Klyashtorin and Lyubushin  and Klyashtorin et al. , 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  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 . It agrees with the observed changes in Earth’s rotation  and Gulf Stream beat . This is illustrated in Figure 17.
This implies that the BIE position variations may be synchronized with the movements of planets in the solar system and follow the grand solar maxima/minima alternations and the 60-yrs cycle during the last 440 years. This suggests a primary planetary—solar forcing (further discussed below, and in Part 2 ).
Figure 17. From planetary beat via a number of North Atlantic variables to the changes in fish faunas as observed by Klyashtorin et al.  (from Mörner ).
8. 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 ) 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. .
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 , suggest that surface temperatures in the Arctic increased during the period 1915-1940, decreased 1940-1970, and then increased from 1975, in good agreement with the proposed “Total Solar Irradiance” (TSI) according to Hoyt and Schatten  and updated by Scafetta and Willson , using the post 1980 ACCRIM satellite measurements for calibration (Figure 18).
Figure 18. Comparison between Total Solar Irradiation (TSI) according to the Scafetta and Willson  calibration (red curve) and the composite temperature curve for Arctic stations (a) and a rural composite for the Northern Hemisphere (b) . The Arctic temperature curve is based on data from 77 Global Historical Climatology Network stations above the Arctic Circle. The NH-composite is based on rural stations in the US, China, Ireland and Arctic. (Reproduced with permission from the authors W. Soon, R. Connolly & M. Connolly).
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 . An example of a coupled interaction is the 13 years’ delay of the North Pacific Multidecadal Oscillation (NPMO) to the Atlantic Multidecadal Oscillation (AMO) .
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 . 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 .
In the following we will discuss various driving mechanisms covering millennial, centennial, multi-decadal and shorter time scales.
8.1. 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  , 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 10Be 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 .
Analysis of solar activity in a 9400-year long record of solar modulation potential φ determined from cosmogenic radionuclides such as 10Be and 14C shows distinct periodicities of 88, 104, 150, 208 and 506 yrs . 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.
An analysis of the Zürich mean annual sunspot number Rz (1749-2009 AD), the group sunspot number Rh, (1610-1995 AD), the Extended time series of Solar activity Indices (ESAI) (1090-2003 AD), and the simulated extended geomagnetic index aa (1619-2003) by Komitov et al.  show that in all series exists a strong cycle with duration 55 - 60 years in the 18th and 19th centuries. It is less pronounced during the end of the 19th and the beginning of the 20th 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.  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  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  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 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 . For the secular changes in Rz 1844 - 2000 Duhau and Chen  determined the following periods: 85.5, 121.1, and 171.2 years.
Finally, it may be significant to note that Blanter et al.  gave 4 possible explanations to the low values of the Solar Cycle 23: it may mark: 1) the beginning of a New Grand Solar Minimum, 2) the beginning of a new Gleissberg cycle, 3) the end of a great maximum episode, or 4) a random variation.
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 , 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  showed only the stable periods of Jupiter, Uranus and Neptune.
According to Kuznetsova , 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. , analysed the Ri series by singular spectrum analysis (SSA) and found 4 components: a multi-secular trend with a slight acceleration around 1900, two Hale periods of 11.0 and 10.0 years, the fist increasing in amplitude, the second with secular amplitude variation with maximum around 1740 and 2010, and a clean Gleissberg period component which varied from 85 to 94 years with a decreasing amplitude—signalling a longer period beat. The multi-secular trend confirms the modern grand solar maximum from 1920  .
Analysis of the last 4000 years of the reconstructed central Greenland surface temperature, derived from the bi-decadal 𝑑18O record from the Greenland Ice Sheet Project II (GISP2)   revealed periods ≈3500, ≈1130, ≈790, ≈580 and ≈370 years . A 3-period model based on non-linear optimization of the longest periods resulted in a good fit that clearly shows the timing of historical warm periods ( , their Fig. 8).
A global temperature harmonic analysis  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 (  ).
Harmonic analysis of the HadCRUT4 global temperature (1850-2012)  showed that a linear trend of 0.5˚C per century and periods 9.1, 21.6 and 66 years. A 9.3 years period is a harmonic of the 18.6 years lunar tide period. This period is also identified as a standing wave in the Atlantic Ocean by Mazimov and Smirnov  . The 21.6 years period is close to the Sun’s magnetic period, the Hale period of 22 years.
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 . 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 ; 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.  shows that periods 17 and 26 years are much stronger than those at 10 - 12 years.
Kasatkina et al.  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  showed that 9 different global and regional temperature series contained the same 10 spectral cycles from the Sun and the Planets, in addition to a 9.1 years period from the Moon.
Scafetta  demonstrated multi-frequency spectral coherence between planetary and global surface temperature oscillations by advanced techniques and statistical significance tests, with large coherence peaks of 20 and 60 years quasi-periods. The latter has maxima around 1880, 1940 and 2000.
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.  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 . 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  . The LTG is clearly modulated by the TSI variations  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 ). The Spörer, Maunder and Dalton Minima are easily observed, like the 1920-2000 high and a coming minimum 2020-2085.
8.2. Lunar Cycles
Treloar  identified that lunar-solar forces produce multiple tidal ocean effects with periods of 86.795, 46.22, 13.53, 18.02 and 20.295 years. Lunar Nodal oscillation components of ~56 and ~74 years are observed in the North Sea, Kattegatt and Baltic tide gauge data  .
An analysis of time-series for the North polar position, extent of Arctic ice, sea level, sea temperature, air temperature, and the NAO winter index shows that a harmonic spectrum from the 18.6-year lunar nodal cycle is present in all time-series  . The cycles in the harmonic spectrum have stationary periods but not stationary amplitudes and phases. A sub-harmonic cycle of 74 years may introduce a phase reversal of the 18.6-year cycle.
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  . Dominating periods were the 18.6 yr lunar nodal period and its 9.3 year harmonics.
8.3. Solar and Lunar Cycles and Earth’s Orbit
Davis and Brewer  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”.
Figure 19 shows the summer insolation at 65˚N during the last 2000 years . It gives the combined effect of variable TSI and orbital variations due to lunar and planetary cycles, including solar variability. The rapid decline 1350-1400 coincides with the start of the Little Ice Age (LIA), which had its maximum cold period 1687-1703, according to the longest instrumental temperature record from Central England . Recovery from this cold period was first visible in the Arctic from about 1920 as demonstrated in Figure 12 and Figure 15. The summer insolation at 65˚N increased in three steps after 1700 with the recent 1920-2000 maximum being the highest. After 2005 the June 21 TSI at 65˚N is decreasing with a minimum in 2045 and a recovery at a lower level about 2085.
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 . 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 .
8.4. 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 , 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 ). 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 14C 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.  during the last 9400 years. The cycle had its maximum in 1910 and a minimum in 2005. This signals the end of the warming period.
The shielding capacity is a function of the solar wind, which brings magnetic plasma, which strengthens the Earth’s magnetic field. Weaker solar wind means less shielding. 30 years of output from a data-constrained magneto-hydrodynamic model of the solar corona has made it possible to model solar wind variations over the last four centuries  as shown in Figure 21 for the equatorial component that hits the Earth. The global solar wind variations show main evolution in four zones: 1) the Maunder Minimum of lower solar wind speed values, 2) the 18th century of generally higher values, 3) the Dalton Minimum of lower values, and 4) the last 180 years of surprisingly similar values.
Figure 20. Variations in solar activity affect the Earth’s rotation, cloud formation and solar irradiance). The isotopes 14C and 10Be may be used as indicators for past solar activity .
Figure 21. Estimated helio-equatorial mean solar wind speed 1617-2013 at the distance of one astronomical unit (AU) from the Sun, the red curve is an 11-yr running mean. Four zones are identified: 1) the Maunder Minimum, 2) the 18th century high, 3) the Dalton Minimum, and 4) the last 180 years with a fairly constant level (modified from ).
Figure 22 shows solar wind measured by spacecraft in geocentric orbits, made available through the omniweb.gsfc.nasa.gov website, and observations at the Oulu neutron count detector. Neutron counts are a good measure of cosmic galactic rays (GCR) particles and are clearly anti-correlated with the solar wind. 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).
Figure 22. (a): Solar wind measured at the Earth’s orbit (6 solar rotations running mean; https://spdf.gsfc.nasa.gov/pub/data/omni/low_res_omni/extended/omni_27_av.dat), (b): Oulu neutron counts (6 months average; http://cosmicrays.oulu.fi). Both graphs show a change in mean trend at about 1992 indicating the onset of decreased shielding capacity.
Bucha  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 (cf.   . Kuznetsova  identified a 198 yrs magnetic field cycle and a 54 yrs solar wind cycle in terrestrial temperature data.
A closer study of the geomagnetic activity over the past century of the solar wind flow speed Vsw, the interplanetary magnetic field strength B, and the open solar flux Fs, indicate that the grand maximum started in 1920, peaked around 1955 and 1986 and then declined. Extrapolating the linear decline trends indicate that the grand maximum of the three parameters will end in 2013, 2014 and 2027, respectively .
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.
9.1. The ENSO Events
The equatorial current lags behind the rotation of the solid Earth, which lead to an accumulation of warm water-masses outside East Asia rising local sea level by about 30 cm. At El Niño events, 0.4 ms of the rate of rotation (LOD) is transferred from the solid Earth to the hydrosphere, which is manifested in an eastward migration of the water-masses from East Asia towards the west coasts of the two American continents. When the water hits the west American coasts, sea level rises by about 30 cm, and 0.4 ms in rotation (LOD) is transferred back to the solid Earth. This provides a significant documentation of the interaction and feedback coupling between the redistribution of mass (water) and changes in rotation with the interchange of angular momentum between the solid Earth and the hydrosphere  .
9.2. The Gulf Stream Pulsation and Earth’s Rate of Rotation
Mörner  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 (Figure 9;   ).
The Gulf Stream beat system proposed by Mörner  has a close relation to the alternations between NAO+ and NAO− stages recorded by Wanner et al. . Their main illustration (Fig. 9 in ) has been endorsed and reprinted by Pinto and Raible  and Boström . 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.  is based on internal interaction and feedback process. It seems significant that Wilson (2011) 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   the changes in rotation originate from planetary-solar beat as illustrated in Figure 9, i.e. from external sources (Figure 7).
9.3. The Ocean Oscillation Interaction
In his account on “the Lunar nodal cycle influence on the Barents Sea”, Yndestad   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   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).
9.4. 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  explains that in the cold Arctic Ocean, water accumulates in circulating layers, the Arctic Surface Water, the Atlantic Water, the Deep Water, and the Bottom Water, all of which have different densities and circulation patterns. The Atlantic Water residence time is about 25 years, while the Deep Water residence time is about 75 years, which is the same period as the polar motion. His model is that the gravity force from the 18.6 years lunar nodal cycle, influences the polar position and most likely the circulating water in the Arctic Ocean, the Arctic ice extent, and Arctic atmospheric condition. A continuous source of energy from the stationary polar position cycle is transformed into a spectrum of oscillating water circulation periods, where some harmonic cycles are preserved by resonance. This is the base for “The wine glass theory”, that the rotating Arctic Ocean behaves like liquid in a rotating wine glass, where the Deep Water oscillation has a resonance of about 74.5 years. The denser Bottom Water has a resonance of 4 × 74.5 = 290 years, while the less dense Atlantic water has a resonance 74.5/3 = 25 years.
Yndestad  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.
9.5. The “Stadium-Wave” Mechanism
Wyatt  and Wyatt and Curry  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  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).
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   and the polar wobble .
The climate of the North Western coast of Europe and the position of the ice edge in the North Atlantic Ocean are controlled by the Northern branch of the Atlantic current (Gulf Stream) and by advection modulated by the Northern Jet 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 18th 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 solar wind reconstruction). From then on, we observe a successive northward shift of the ice edge (Figure 12).
The explorations of the Arctic that started late in the 16th century led to the discovery of Svalbard and rich stocks of seals and bowhead whales (Balaeana mysticetus). By inspection of logbooks from expeditions and hunters and later from airplanes and satellites, it has been possible to estimate the BIE August position since 1579. This series is not complete, but it makes it possible to investigate secular trends, which are important to understand the evolution of the Arctic ice cover and Earth’s climate.
When the BIE expanded during the cold 19th century, and the stock had been hunted down, commercial whale hunting stopped. The bowhead disappeared (Figure 14). It came back in the 1990s when the BIE position retreated north of Svalbard (Figure 15).
The position of the Sea Ice Edge in the Barents Sea (BIE) is a result of many dynamic factors acting on different time scales. Signatures of periodic and quasi-periodic variations can give indications on relations preserved during a chain of processes as indicated in Figure 9 and Figure 18. Periods with much ice correspond to deep minima in the solar activity, suggesting a causal connection between solar activity and the BIE positions.
10.1. 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 (zones 1 and 3). The solar maximum of the 18th century, is recorded in strongly intensified inflow of warm Atlantic water to the Arctic, high solar irradiance at 65˚N (Figure 19), high values of solar wind (Figure 21, zone 3), and generally warm climate conditions and strongly reduced ice extension in the GIN Sea region with the April maximum ice edge of year 1769 located even north of Svalbard (Figure 12).
A particular problem with this interpretation is that the maximum catch of bowhead whales took place in the period 1685-1703 (Figure 14) when the coldest climate in the millennium was recorded in Europe, and the global sea surface temperature was 1.0˚C - 1.5˚C cooler than in the 1950-1980 period . From Figure 15 and Figure 21 it is clear that the BIE moved north at least 50 years before the recent high solar wind level was reached. Gulf Stream circulation type-a dominated the 18th century ( , Fig. 7a in ). The water north of Svalbard and in the Fram Strait was apparently warm enough for rich marine life. This seems to call for a local source of heat or special wind pattern, maybe recorded in the weak positive NAO-index in 1705-1740 .
Another problem is that the high level of solar wind through the 19th 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 : The Sun was stressed after the producing short sunspot-cycles resulting in the 18th century maximum, and it relaxed by producing long cycles in the 19th century. Only at the start of cycle 16 (1900) it was “refreshed” and able to produce short cycles . This pushed the BIE north in the 20th 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  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 14C production record of the last 9400 years .
10.2. 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 ).
A 60-year quasi-periodic cycle is detected in global temperature (Figure 19;     ), in TSI , in LOD   , in geomagnetism  , in Jet Stream changes between zonal and meriodional circulation   , in sea level changes     and ocean oscillations systems like PDO, AMO and NAO (Figure 18; ).
Scafetta (e.g.  ) and Solheim  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 significant that Yndestad was able to identify the 18.6 and 74.4 tidal cycles in different Norwegian Sea parameters    , and Hansen  was able to document the 55.8 yrs and 74.4 yrs tidal cycles in the North Sea and Kattegatt records.
10.3. The Near Future
Cyclic forcing functions imply that we would have means of predict changes in the near future. The deterministic model analyses of sunspots and TSI recorded by Yndestad and Solheim  suggest that we will have a new Dalton-type minimum by about 2025-2050 and 2040-2065, respectively. This fits well with other estimates of a New Grand Solar Minimum; e.g. at about 2030-2050  , starting in 2014 and reaching its minimum in 2055-2060 ), in 2030-2040  , and around 2040 .
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  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 14C prediction and the Gulf Stream beat between type-a and type-c circulation patterns (  ; Appendix 2).
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 17th century, extreme northern positions in the 18th century and expansion to low latitudes in the early 19th century. Since the end of the Dalton minimum about 1820, BIE is moving north with cycles imposed on a linear trend; slow in the beginning, faster after 1920. We correlated this with the changes in Gulf Stream beat with its major changes in inflow of Atlantic water to the Arctic; low inflow, cold climatic condition and ice expansion during the Maunder Minimum, strong inflow, warm climatic conditions and ice reduction during the Solar Maximum of the 18th century, and cold climatic condition and ice expansion during the Dalton Minimum. This is consistent with a solar forcing where the changes in Earth’s shielding strength play a central role. Much ice persisted during the 19th century. Significant warming and ice melting are recorded during the 20th and early 21st centuries (except for a minor step back in 1940-1960).
This implies changes in solar wind, and an ultimate driving force in planetary beat (Figure 9). The planetary beat includes the interacting effects of Jupiter, Saturn, Uranus and Neptune on the 60-yrs, 84-yrs and 208-yrs cycles in addition to an integrated effect of almost two centuries of persistent high solar wind (Figure 21) or a long cycle perhaps of millennium length, driven by the combined outer planets, not to forget lunar cycles imposed by the polar wobble. Peaks (or ice minima) 180 years apart may be a result of planetary synchronization (the José cycle).
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. 1) the long-term solar cycles, which have set the character of the terrestrial changes in climate during the last millennium (Figure 4), are in the stage of changing back to a new Grand Solar minimum , 2) the solar cycles record a progressive decrease in sunspot number of cycles 23 and 24, 3) the atmospheric neutron count and solar wind speed have changed in the last decade (Figure 22) towards weaker shielding conditions, and 4) the Arctic Oceans are cooling in the last decade (Figure 10(A) and Figure 10(B)) and the North Atlantic is recording a freshening .
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 , 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.
This is our first report of the Barents Sea Ice Edge (BIE) project. It is a typical collaboration product. JES has acted as editor and contributed large parts of the text and figures, NAM has contributed large parts of the text and figures, OH is behind most of Section 6, and SFP is responsible for most of Section 7. None of the authors have any special funding for this project to declare. We declare no conflict of interest. Our second report will follow soon. It is entitled: Changes in Barents Sea Ice Edge Positions during the last 440 years: Sun, Moon or Planets  . Two supplementary notes (by NAM) are added to the paper as Appendix 1 and Appendix 2.
Appendix 1 (by Nils-Axel Mörner)
The “stadium wave wheel” mechanism in need of external forcing
Wyatt  and Wyatt and Curry  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  proposed that the variables of the Northern Hemisphere surface temperature (NHT), the Atlantic Multidecadal Oscillation (AMO), the North Atlantic Oscillation (NAO), the El Nino-Southern Oscillation index (NINO3.4), the North Pacific Oscillation (NPO), the Pacific Decadal Oscillation (PDO), the Aleutian Low Pressure Index (ALPI), and the index of Atmospheric-Mass-Transfer anomalies (AT) could be combined into a logical system of a wave propagating around the Northern Hemisphere in a manner resembling a stadium wave of applauding athletes at sport events. This initiated the “Stadium Wave” Wheel theory of Wyatt and Curry . Wyatt  “did not assign an exact period to the signal; rather, the signal is described as secularly varying, with an apparent 64-year period during the 20th century” (quot. ). Therefore, I think, we are allowed to elaborate with the timing of their four (I-IV) “temporal groups” as done in Figure 23 and Figure 24.
Figure 23. The “stadium wave wheel”  here divided into two halves; a warm half beginning 1910 and 1976 and a cold half beginning 1878, 1945 and 2006 (for explanations of acronyms ).
Figure 24. The warm segments I-IV (orange) and cold segments I-IV (blue) in the “stadium wave wheel”  compared to their probable positions with respect to the global temperature curve   and its division in periods of 30 years rises and falls indicating a 60-yrs external driving force .
In Figure 23, I have simply divided the “stadium wave wheel” of Wyatt and Curry  into two halves; a warm half beginning 1910 and 1976 and a cold half beginning 1878, 1945 and 2006 as given the main global temperature curve in Figure 24.
In Figure 24, we compare the temporal zonation given by Wyatt and Curry  and the main global temperature record of the last 150 years   . The similarity is clear. There is an alternation between warm and cold periods of a duration of about 30 years; i.e. a full cycle of 60 years. It predicts a cooling from 2006 to 2036. The dominant 60-yr cycle is indicative of an external force of planetary-solar origin.
I therefore propose that the stadium wave mechanism by Wyatt and Curry  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 , and of AMO, NAO and PDO in external rotation   .
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  , the changes in a combined curve of the Gleissberg and de Vries cycles (Fig. 5 in ) and the observed changes in atmospheric 14C production ( ; 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  indicate that there is a close correlation between Gulf Stream beat and Earth’s rate of rotation. The Solar Wind changes observed (Figure 21; ) are indicative of the solar-terrestrial interaction  .
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).
Figure 25. Relations among solar cycles, Earth’s rate of rotation and the observed changes in the ocean circulation in the North Atlantic  .
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 Figure 3, in the following way (Figure 26 from ):
Figure 26. Changes main Gulf Stream circulation patterns during the last 800 years between type-a (orange wedges upwards) and type-c (blue wedges downwards) (from ).
The changes in Gulf Stream circulation pattern (Figure 26) fit very well with a curve combining the Gleissberg 84 yrs cycle and the De Vries 240 yrs cycle (from ) indicating correlation between solar variability and Gulf Stream beat (Figure 27).
Figure 27. Cyclic phasing of the combined “Gleisberg” and De Vries Cycles over the last 600 years giving a new Solar Minimum at about 2040-2050 (from ). 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 . Vertical scale gives temperature in centigrade.
The changes in Gulf Stream beat (Figure 26) also correlates well with recorded changes in 14C 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 14C production (below) by Bard et al.  termed “solar irradiance” though it, in fact, is a curve of geomagnetic shielding primarily controlled by the Solar Wind effects    . 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 . The external forcing has later been confirmed by Mazzarella , and Mazzarella and Scafetta . Figure 29 record resemble the Solar Wind record (Figure 21) in showing long cycles during the 17th and 18th centuries and a dominance of a 60 yrs cycle for the last 200 years.
Figure 29. The changes in LOD (down implies acceleration and up deceleration) and corresponding changes in temperature and sea level in the Northwest European region (from ). The correlations seem good and lend support to the rotational/circulation theory proposed  . The changes in LOD as well as in temperature and sea level follow a cyclic pattern: centennial for the first 200 years and of an about 60 yrs cyclicity for the last 200 years. The PDO cycle shows a similar pattern with a cyclicity of about ~56 yr (1921-1977 and 1942-1998).
Figure 30. Delta 14C variations , with interpretations in terms of Grand Solar Maxima/Minima and terrestrial warm/cold phases (also type-a/type-c Atlantic circulation).
Figure 30 (from ) provides a very sharp and clear picture of the changes in atmospheric 14C 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)  .
*Barents Sea Ice Edge (BIE).
 Mörner, N.-A. (2010) Solar Minima, Earth’s Rotation and Little Ice Ages in the Past and in the Future: The North Atlantic-European Case. Global and Planetary Change, 72, 282-293.
 Sekonya, K.G. (2017) Development of a 6 Mv Tandem Accelerator Mass Spectrometry Facility and Its Applications. PhD Thesis at the Faculty of Science, University of the Witwatersrand, Johannesburg, 197 p.
 Dergachev, V.A. and Vasiliev, S.S. (2019) Long-Term Changes in the Concentration of Radiocarbon and the Nature of the Hallstatt Cycle. Journal of Atmospheric and Solar-Terrestrial Physics, 182, 10-24.
 Mörner, N.-A. (2012) Planetary Beat, Solar Wind and Terrestrial Climate. In: Escaropa, C.D. and Berós Cruz, A.F., Eds., Solar Wind: Emission, Technologies and Impacts, Nova Sci. Publ., New York, 47-66.
 Mazzarella, A. (2008) Solar Forcing of Changes in Atmospheric Circulation, Earth’s Rate of Rotation and Climate. The Open Atmospheric Science Journal, 2, 181-184.
 Mörner, N.-A. (1988) Terrestrial Variations within Given Energy, Mass and Momentum Budgets; Paleoclimate, Sea Level, Paleomagnetism, Differential Rotation and Geodynamics. In: Stephenson, F.R. and Wolfendale, A.W., Eds., Secular Solar and Geomagnetic Variations in the Last 10,000 Years, Kluwer Academic Publishers, Dordrecht, 455-478.
 Mörner, N.-A. (1989) Changes in the Earth’s Rate of Rotation on an El Niño to Century Basis. In: Lowes, F.J., et al., Eds., Geomagnetism and Paleomanetism, Kluwer Academic Publishers, Dordrecht, 45-53.
 Scafetta, N. (2010) Empirical Evidence for a Celestial Origin of the Climate Oscillations and Its Implications. Journal of Atmospheric and Solar-Terrestrial Physics, 72, 951-970.
 Yndestad, H. and Solheim, J.-E. (2017) The Influence of Solar System Oscillation on the Variability of the Total Solar Irradiance. New Astronomy, 51, 135-152.