Sensitivity and robustness of sky-polarimetric Viking navigation: Sailing success is most sensitive to night sailing, navigation periodicity and sailing date, but robust against weather conditions.

Although Viking sailors did not have a magnetic compass, they could successfully navigate with a sun-compass under a sunny sky. Under cloudy/foggy conditions, they might have applied the sky-polarimetric Viking navigation (SPVN), the high success of which has been demonstrated with computer simulations using the following input data: sky polarization patterns measured with full-sky imaging polarimetry, and error functions of the navigation steps measured in psychophysical laboratory and planetarium experiments. As a continuation of the earlier studies, in this work we investigate the sensitivity of the success of SPVN to the following relevant sailing, meteorological and navigational parameters: sunstone type, sailing date, navigation periodicity, night sailing, dominance of strongly, medium or weakly cloudy skies, and changeability of cloudiness. Randomly varying these parameters in the simulation of Viking voyages along the latitude 60° 21' 55'' N from Norway to Greenland, we determined those parameters which had strong and weak influences on the success of SPVN. The following intrinsic parameters of the simulation were also randomly changed: sailing speed, visibility distance of Greenland's southeast coastline and start time of diurnal sailing. Our results show that the sailing success is sensitive to the night sailing, navigation periodicity and sailing date, while it is robust against the sunstone type, dominance of strongly, medium or weakly cloudy skies, and changing cloudiness.

Introduction

Vikings navigated successfully on the North Atlantic Ocean for three hundred years without a magnetic compass [1-5]. Under sunny weather conditions, they might have navigated by a sun-compass and when the sun was occluded by clouds/fog or was below but near the horizon, they used this compass and the sky polarization detected by sunstone crystals [6-8]. Although all aspects of this hypothetical sky-polarimetric Viking navigation (SPVN) were purely speculative, many scientists [9-13] accepted and cited it as if it were a fact. Some sceptic researchers, however, expressed contra-arguments [14] disputing why this navigation method could not have functioned under certain meteorological conditions. Between the two groups of believers and sceptics, a third group was formed by researchers who tried to reveal which components/aspects of this hypothesis may be valid or unfounded. Ropars et al. [15] studied a calcite crystal that might have been a precise depolarizing sunstone. Le Floch et al. [16] analysed the sixteenth century Alderney calcite crystal that could be an efficient reference optical Viking compass.

In the last decades the most thorough and systematic studies of SPVN have been performed by a Hungarian group cooperating with German, Swedish and Swiss researchers: They measured the polarization patterns of the sky with full-sky imaging polarimetry, and investigated the atmospheric optical and meteorological prerequisites of SPVN [17-26]. The error functions of the steps of SPVN have been measured on numerous test persons in psychophysical laboratory and planetarium experiments [27-31]. Száz et al. [31] determined the accuracy of SPVN as functions of the solar elevation and cloudiness. Bernáth et al. [32] proposed an alternative interpretation of the Uunartoq (Viking) sundial artefact and suggested that it might have been an instrument to determine the latitude and local noon. Bernáth et al. [33] also interpreted the Uunartoq fragment as a twilight board and demonstrated in a field experiment how the Viking sun-compass could have been used with sunstones before and after sunset. Using celestial polarization patterns and psychophysically measured errors, Száz and Horváth [34] revealed with computer simulations the chance that Vikings could reach Greenland.

After these studies the logical question of the sensitivity or robustness of SPVN has arisen. Using an improved version of the software of Száz and Horváth [34], in this work we investigate the sensitivity/robustness of the success of SPVN to the following relevant sailing, meteorological and navigational parameters: sunstone crystal type, sailing date, navigation periodicity, sailing or staying at night, dominance of strongly, medium or weakly cloudy skies, and cloudiness changeability. Randomly varying these parameters in the simulations of Viking voyages, we determined those parameters which have strong and weak influences on the success of SPVN. The sailing speed and start time of diurnal sailing were also randomly changed, and the visibility distance of Greenland’s southeast coastline depended on the current cloudiness situation.

In this work, the terms ‘navigation’ and ‘sailing’ are used in the following contexts: Navigation means the four-step process of SPVN, during which a Viking navigator determines the geographical north and then the ship’s intended moving direction. Sailing means the journey of a Viking ship, the navigator of which performs numerous navigation processes.

Methods

Computational methods

We simulated 1 000 000 Viking voyages from Bergen (Norway) to Hvarf (Greenland). The route of the Viking voyage between Norway and south Greenland (Hvarf) ran along the 60° 21’ 55’’ northern latitude and started from Hernam (nowadays the Norwegian Bergen) [2, 3]. One of the possible explanations of this latitude could be that the prevailing ocean currents might have been taken into consideration by Vikings. On the other hand, this straight route was the shortest between Norway and Hvarf. If the departure had not been from Bergen, but much further north or south in Norway, then the duration of a voyage would have been longer, and other shapes of the gnomonic lines should have been used as the straight (for spring equinoxes) or hyperbolic (for summer solstices) lines found on the sun-compass artefacts. Our computations had the following five internal parameters: solar elevation, north error, sailing speed, start time of diurnal sailing, and visibility distance of Greenland’s southern coastline. The sensitivity/robustness of the successful SPVN was determined for the following six variables: sailing date, sunstone type, night sailing, cloudiness dominance, cloudiness changeability, and navigation periodicity. Before each run of a simulated voyage, the values of the following parameters were chosen randomly and independently of each other:

sailing date: spring equinox or summer solstice

sunstone type: calcite, cordierite or tourmaline

night sailing: yes or no

The spring equinox (21 March) and the summer solstice (21 June) meant two important dates of the Viking sailing season: the former and latter was approximately the start and the middle of this season. In our simulations we chose only the spring equinox and the summer solstice as sailing dates, because the two artefacts of the Viking sun-compass (discovered in the Vatnahverfi and Uunartoq Fjords of Greenland) had only two clearly discernible gnomonic lines (trajectories of the tip of the shadow cast by the vertical gnomon on the horizontal disc of the sun-compass in sunshine) [3]: a straight line for navigation at spring equinoxes and a hyperbolic line for summer solstices. For other sailing dates the gnomonic lines are different hyperbolic curves.

The values of the following continuous parameters were also chosen randomly from a uniform distribution:

cloudiness dominance: -1 ≤ mdominance ≤ +1

cloudiness changeability: 0 ≤ σchangeability ≤ 4

navigation periodicity: 0.5 hour ≤ Δt ≤ 6 hours

The temporal step of each simulated voyage was ε = 1 minute, while the cloudiness ρ was changed every 60 minutes. During an ε period, the ship advanced Δ = ε, where is the velocity vector of the ship in the current step. We performed the simulations with a custom-developed Python script. The variables during the simulations were the following:

Sailing date: Spring equinox or summer solstice

Voyages were simulated at spring equinox (21 March) and summer solstice (21 June). Before each run we randomly selected the sailing date.

Sunstone type

In the simulations, we used calcite, cordierite and tourmaline sunstones, for which the error functions of the steps of SPVN have been published earlier [28, 30]. Before each run, we randomly selected a crystal, which was used during the full voyage.

Solar elevation

Always the actual solar elevation θ was used in every point of time of the sailing. Solar elevation data were taken from the following web site: https://www.sunearthtools.com/dp/tools/pos_sun.php.

Navigation periodicity

During a given time period τ, the simulated Viking ship moved in a constant direction along a straight line with a constant sailing speed w until the navigator determined the new sailing direction. In the simulation, the time segment τ was chosen randomly with a uniform distribution from the interval Δt − Δt/6 ≤ τ ≤ Δt + Δt/6, where Δt is the navigation periodicity (= the elapsed time since the last navigation) [34]. Before each run, Δt was randomly selected from the following range: 0.5 hour ≤ Δt ≤ 6 hours.

Dominance and changeability of cloudiness

Száz and Horváth [34] assumed an equal distribution/chance of good (with cloudiness 0 ≤ ρ < 3 oktas) and bad (5 oktas < ρ ≤ 8 oktas) weathers during their simulated Viking voyages. This assumption is considered here as a medium (normal) weather with cloudiness dominance mdominance = 0 meaning a medium cloudy sky. We tried to simulate the temporal change of cloudiness as realistic as possible. The cloudiness ρ ranged from 0 to 8 oktas (0 okta: no clouds, 8 oktas: completely overcast sky), where okta means 1/8 area of the sky hemisphere. The first ρ-value at the start of a simulated voyage was chosen from a Gaussian distribution with a median of 4 and a deviation of 2. Each succeeding ρ-value was calculated by adding a discrete change (increment) in okta to the previous ρ-value according to a second Gaussian distribution with median mdominance and deviation σchangeability, where mdominance denotes the dominance of clouds on the sky (-1 ≤ mdominance ≤ -0.75 for weakly cloudy skies, -0.125 ≤ mdominance ≤ +0.125 for medium cloudy skies, +0.75 ≤ mdominance ≤ +1 for strongly cloudy skies) and σchangeability denotes the cloudiness changeability (the higher the σchangeability, the higher the changeability of oktas in the next hour), so that ρ had to remain in the range of 0 and 8 oktas: for ρ < 0 okta or ρ > 8 oktas, the cloudiness was set to 0 and 8 oktas, respectively. The cloudiness ρ was recalculated hourly.

North error

The angle of deviation from the geographical north, called north error, was determined as follows: For a given solar elevation angle θ we calculated the actual cloudiness ρ as described above, then we used the 1080 data files from the θ-ρ matrix determined in [31]. In this matrix, every θ-ρ pair contains 12 different sky situations with a known north error distribution calculated with the use of the errors of the four steps of SPVN measured in psychophysical experiments. One of these 12 sky situations was selected randomly according to uniform distribution, then from the corresponding error distribution we used a randomly chosen error value as the actual north error.

Sailing speed

We characterized the ship’s sailing speed w (which was constant 11 km/h in [34]) with three parameters: (i) maximum speed wmax = 21 km/h, (ii) average speed wave = 11 km/h, and (iii) speed’s standard deviation Δwave = 2 km/h. For each change of the north error the simulation generated a new speed value.

Night sailing

Considering the night sailing, in the simulations of sailing routes there were two options:

Between sunset and sunrise, Vikings lowered their sails and stopped their voyage as in [34]. In this case, at night the ship’s position was constant. With a non-lowered sail, the ship would have randomly drifted due to the random winds, which drift would have made difficult the safe daytime continuation of the voyage. With a lowered sail, the ship did not drift at night. The influence of water currents and wind was neglected, because the average water speed in the Atlantic Ocean is 0.54–0.72 km/h [35]. On the other hand, the average wind speed is more than ten times larger (50–60 km/h, [35]).

At night, the ship continued its route in the last direction determined by the navigator at sunset.

Before the start of a simulated voyage, we randomly turned on or off the night sailing option.

Start time of diurnal sailing

The start time of diurnal sailing was selected randomly on every sailing day after sunrise within the range 0 and 1 hour from a uniform distribution. Note that in [34] the navigator started to navigate always at sunrise with zero sun elevation.

Visibility distance of Greenland’s southeast coastline (end of simulation)

A given simulation was stopped, when the ship reached the northeast coastline of North-America (unsuccessful sailing route) or the visibility distance d of the southeast coastline of Greenland’s south tip (successful sailing route), depending on the actual cloudiness value ρ (measured in okta) as follows: d(ρ) = -16.02875·ρ + 128.23 km. This formula is derived in [34]. Assuming that in overcast and foggy weather with ρ = 8 the visibility distance d(ρ = 8) = 0 km, the d(ρ) function was calculated with a linear interpolation between d(ρ = 8) = 0 km and d(ρ = 0) = 128.23 km. The linear function d(ρ) modelled how the cloudiness may affect the visibility distance of the coastline, where ρ changed stochastically. The map of the North Atlantic region was generated by a software written by us. The contours of continents and islands were manually digitalized from the map available as open-source data from http://www.gnuplotting.org/plotting-the-world-revisited/ (the raw data points of the contours can be freely downloaded in text format from: http://www.gnuplotting.org/data/world_10m.txt). These open-source data can be freely used without permission/licence. The coastline of Greenland (Fig 1) as the direct goal of voyages became limited to a southeast section, instead of the full east borderline in [34]. When a ship reached shore on the northern coastline of Greenland where there was no Viking settlement, the voyage was considered unsuccessful.

Fig 1

Examples for simulated sailing routes using SPVN from Norway to Greenland.

Successful (green) and unsuccessful (red) routes of 1000 Viking voyages from Bergen (Norway) to Hvarf (Greenland) at summer solstice using cordierite sunstone with navigation periodicity Δt = 1 hour (A, B), Δt = 3 hours (C, D) and Δt = 6 hours (E, F) without (A, C, E) and with (B, D, F) night sailing. The values of the sailing success rate s are: (A) 99.8%, (B) 99.8%, (C) 19.0%, (D) 100.0%, (E) 0.0%, (F) 24.3%. The blue curve is the borderline of visibility from which the southeast mountains of Greenland can already be seen from a Viking ship. The simulation of a voyage stops when the navigator sees the southeast coastline where the visibility distance is determined by the current cloudiness value ρ. Some simulated sailing trajectories pass through Iceland and/or North Scotland. In these cases, it was assumed that the Vikings continued their voyage towards Greenland. The maps are generated by our software after a manual selection of the contours of continents and islands from the open-source data originating from http://www.gnuplotting.org/plotting-the-world-revisited/.

Statistical analyses

We used logistic regression to obtain the significant variables determining the probability π of successful SPVN as functions of the variables sailing date, sunstone type, cloudiness dominance mdominance, cloudiness changeability σchangeability, navigation periodicity Δt and night sailing. Száz and Horváth [34] showed that the sailing success does not change monotonously with increasing Δt. Therefore, for the logistic regression we divided the continuous variable Δt into 100 equal intervals. In the 1st logistic model we supposed interferences between mdominance and σchangeability, as well as between sunstone type, Δt and night sailing. We split the simulation results to train (900 000 simulations) and test (100 000 simulations) datasets. Using the train dataset, we built the 1st logistic model and determined the significant variables by applying ANOVA test for this model. The test dataset was used to calculate the accuracy of the model. On the basis of the results of the 1st model, we applied a 2nd logistic regression (model) with only the significant variables to determine the probability of successful sailing under different conditions. Statistical analyses were performed with the R statistics package 3.6.3. [36].

Results

Some visualized Viking voyages

Fig 1 visualizes some simulated Viking sailing routes from Bergen (Norway) to Hvarf (Greenland) at summer solstice using cordierite with navigation periodicity Δt = 1, 3 and 6 hours without and with night sailing. The sailing success rate s is larger for Δt = 1 hour (without night sailing: s = 99.8%, with night sailing: s = 99.8%) than for Δt = 6 hours (without night sailing: s = 0.0%, with night sailing: s = 24.3%), independently of the night sailing. However, for Δt = 3 and 6 hours, with night sailing s (sΔ = 100.0%, s Δ = 24.3%) is larger than without night sailing (sΔ = 19.0%, sΔ = 0.0%). Fig 1 demonstrates that night sailing has a significant influence on the sailing success, that is explained in the Discussion.

Robust parameters having minimal effect on the sailing success

Fig 2A shows the success rate of simulated voyages at summer solstice (21 June) as a function of the navigation periodicity Δt for calcite, cordierite and tourmaline sunstones without night sailing. The sunstone type has only a minor effect on the sailing success rate s. The effect of sunstone type on s is similarly negligible for different sailing dates, independently of night sailing. Száz and Horváth [34] obtained that the accuracy of sunstone adjustment in the first step of SPVN slightly depends on the crystal type. Here, we demonstrate that the sunstone type has a negligible influence on the success rate s of voyages using four-step SPVN. Thus, SPVN is robust against the sunstone type.

Fig 2

Success rate s of simulated voyages using SPVN at summer solstice (21 June) as a function of the navigation periodicity Δt without night sailing.

(A) s(Δt) for calcite, cordierite and tourmaline sunstones. (B) s(Δt) with calcite sunstone for weakly cloudy (-1 ≤ mdominance ≤ -0.75), medium cloudy (-0.125 ≤ mdominance ≤ +0.125) and strongly cloudy (+0.75 ≤ mdominance ≤ +1) skies. (C) s(Δt) with calcite sunstone for weak/slow (1 ≤ σchangeability ≤ 1.5), medium (2.25 ≤ σchangeability ≤ 2.75) and strong/rapid (3.5 ≤ σchangeability ≤ 4) cloudiness changeabilities of the cloudiness ρ between simulation steps. The sunstone type, cloudiness dominance mdominance and cloudiness changeability σchangeability have only a minor influence on s, thus SPVN is robust against these parameters.

Cloudiness dominance

Fig 2B displays the success rate s of voyages using calcite sunstone at summer solstice (21 June) without night sailing as a function of the navigation periodicity Δt for weakly cloudy (-1 ≤ mdominance ≤ -0.75), medium cloudy (-0.125 ≤ mdominance ≤ +0.125) and strongly cloudy (+0.75 ≤ mdominance ≤ +1) skies. The effect of mdominance variable is negligible on s. The influence of cloudiness dominance on s is similarly minimal for different sailing dates, independently of night sailing. Hence, SPVN is robust against the cloudiness dominance mdominance.

Cloudiness changeability

Fig 2C illustrates the success rate s of voyages using calcite sunstone at summer solstice (21 June) without night sailing as a function of the navigation periodicity Δt for weak/slow (1 ≤ σchangeability ≤ 1.5), medium (2.25 ≤ σchangeability ≤ 2.75) and strong/rapid (3.5 ≤ σchangeability ≤ 4) change possibilities σchangeability of the cloudiness ρ between simulation steps. σchangeability has only a minor effect on s. The effect of sunstone on s is similarly negligible for different sailing dates, independently of night sailing. Consequently, SPVN is robust against the cloudiness changeability σchangeability.

Sensitive parameters significantly influencing the sailing success

Sailing date

Fig 3 illustrates the success rate s of simulated voyages using SPVN with calcite sunstone at spring equinox (21 March) and summer solstice (21 June) as a function of the navigation periodicity Δt with and without night sailing. Our simulations show a strong dependence of s on the sailing date: The sailing successes at spring equinox and summer solstice differ significantly from each other, independently of the night sailing. At summer solstice, the length of the day suitable for SPVN is longer than at equinox, thus the advancement of the ship is definitely longer for the former date. If, however, we assume continuous advancement of the ship at night based on the last set direction during the last navigation before sunset, the advancement time is equal for both dates. The significant differences between the sailing successes at spring equinox and summer solstice mean that the success s of SPVN is sensitive to the sailing date.

Fig 3

Success rate s of simulated voyages using SPVN with calcite sunstones at spring equinox (21 March, A-B) and summer solstice (21 June, C-D) as a function of the navigation periodicity Δt without (A, C) and with (B, D) night sailing.

Since the effect of sunstone type, cloudiness dominance mdominance and cloudiness changeability σchangeability was negligible on s, we did not sort the data for these variables when created this figure.

In their simulations, Száz and Horváth [34] assumed that the Viking ships did not sail during the night. In our present simulations, we tested the effect of night sailing on the sailing success. The success rate s strongly depended on night sailing. As Száz et al. [31] showed, if the navigation is in the afternoon, then the determined direction always points toward north. When night sailing was allowed, the Vikings sailed toward north all night (because even with Δt = 6 hours, the last navigation of the day always happened in the afternoon). If other parameters were set so that the ship deviated toward south (e.g. in summer solstice with Δt = 3 and 6 hours), night sailing could help to increase the success rate s (see Fig 1C versus 1D and 1E versus 1F). However, for deviating too much toward north (i.e. anything above 65.37° latitude at Greenland’s coastline), the sailing was regarded as unsuccessful. According to Fig 3, the sailing success rate s versus Δt is more or less chaotic, depending strongly on the initial conditions of the simulation (sailing date and night sailing). Thus, SPVN is sensitive to the night sailing.

Fig 3 demonstrates that the sailing success rate s strongly depends on the navigation periodicity Δt. Thus, SPVN is sensitive to Δt.

Statistics

The accuracy of the 1st logistic model is 0.8874. According to Table 1, in this model the variables sailing date, crystal type, cloudiness dominance mdominance, navigation periodicity Δt and night sailing are highly significant (p < 0.05), while the cloudiness changeability σchangeability is not significant. The interacting variable pairs are also significant. The high significance of most of the variables is not surprising because of the high number of simulation (900 000 runs in the train dataset).

Table 1

Results of the ANOVA test for the 1st logistic model.

variables	df	deviance	resid. df	resid. dev	p	significance
NULL			999999	1348182
crystal type	2	264	999997	1347918	< 2.2e-16	***
m dominance	1	74	999996	1347844	< 2.2e-16	***
σ changeability	1	0	999995	1347843	0.496711
Δt	99	141068	999896	1206775	< 2.2e-16	***
night sailing	1	45475	999895	1161300	< 2.2e-16	***
sailing date	1	161	999894	1161139	< 2.2e-16	***
mdominance & σchangeability	1	8	999893	1161131	0.004813	**
Δt & night sailing	99	248383	999794	912748	< 2.2e-16	***
Δt & sailing date	99	280489	999695	632260	< 2.2e-16	***
night sailing & sailing date	1	100488	999694	531772	< 2.2e-16	***
Δt & night sailing & sailing date	99	18739	999595	513033	< 2.2e-16	***

df: degree of freedom, resid. df.: residual degree of freedom. The interacting variables are connected by &.

For the 2nd logistic model we used only the significant variables determined by the 1st logistic model: the sailing date, crystal type, cloudiness dominance m, navigation periodicity Δt and night sailing. We supposed interactions between the following variables: sailing date, crystal type and navigation periodicity Δt. Table 2 contains the results of the ANOVA test applied for the 2nd logistic model. All variables and variable interactions show highly significant effects. The accuracy of the 2nd model is 0.8873, practically the same as that of the 1st model.

Table 2

Results of the ANOVA test for the 2nd logistic model.

variables	df	deviance	resid. df	resid. dev	p	significance
NULL			999999	1348182
crystal type	2	264	999997	1347918	< 2.2e-16	***
m dominance	1	74	999996	1347844	< 2.2e-16	***
Δt	99	141067	999897	1206777	< 2.2e-16	***
night sailing	1	45475	999896	1161302	< 2.2e-16	***
sailing date	1	161	999895	1161141	< 2.2e-16	***
Δt & night sailing	99	248378	999796	912763	< 2.2e-16	***
Δt & sailing date	99	280470	999697	632293	< 2.2e-16	***
night sailing & sailing date	1	100481	999696	531811	< 2.2e-16	***
Δt & night sailing & sailing date	99	18736	999597	513075	< 2.2e-16	***

df: degree of freedom, resid. df.: residual degree of freedom. The interacting variables are connected by &.

Figs 4–6 show the predicted probability π of successful SPVN as functions of the significant variables. It is clearly seen that the crystal type and cloudiness dominance mdominance do not affect π. The shorter the navigation periodicity, the higher the π, but this effect is not monotonous as seen in Figs 4–6. However, if the voyage happens at spring equinox with night sailing, the longest navigation periodicity Δt = 6 hours leads to higher π.

Fig 4

Probability π of successful sky-polarimetric Viking navigation (SPVN) as functions of the navigation periodicity Δt and cloudiness dominance mdominance for calcite sunstone predicted by the 2nd logistic model.

It is clearly seen that π is independent of mdominance.

Fig 6

Probability π of successful SPVN as functions of the navigation periodicity Δt and cloudiness dominance mdominance for tourmaline sunstone predicted by the 2nd logistic model.

It is clearly seen that π is independent of mdominance.

Probability π of successful SPVN as functions of the navigation periodicity Δt and cloudiness dominance mdominance for cordierite sunstone predicted by the 2nd logistic model.

Discussion

This work is a logical sequel of our previous investigation [34] in the topic of sky-polarimetric Viking navigation (SPVN). We provide a further elaboration of the potential success rate of sailing by SPVN using dichroic/double-refracting sunstone crystals to the southeast shore of Greenland from Norway. Our new findings are obtained via a computer simulation which is refined compared to its original version [34]. This improved model has subtler features, and therefore, is capable to take into consideration more precisely also the weather that may influence the success of SPVN. One of the main findings of this paper is that SPVN can be very robust against weather conditions, which could be a fundamental factor in verifying the operability and veracity of this navigational method.

Both the earlier [34] and the present simulations of North Atlantic voyages using SPVN share the common feature that the same

sunstone crystals (cordierite, tourmaline, calcite),

sailing dates (spring equinox, summer solstice) and

north error datasets

were used, furthermore

the small quasi-random drift of the ship with lowered sails (when the voyage was stopped at night) due to water currents and wind was neglected.

In order to make more reliable the model of Viking voyages, the simulations of Száz and Horváth [34] were improved in our present work as follows:

The discrete values of the navigation periodicity Δt (ranging from 1 to 6 hours with an increment of 1 hour) were changed to continuous values chosen randomly between 0.5 and 6 hours from a uniform distribution.

The following new variables were introduced:

night sailing,

cloudiness dominance mdominance and

cloudiness changeability σchangeability

The sailing speed w of the ship changed around an average value wave = 11 km/h with a Gaussian normal distribution (wmax = 21 km/h, Δwave = ±2 km/h), instead of using the earlier constant value (11 km/h).

Instead of using the earlier constant visibility distance d of the south coastline of Greenland, d of Greenland’s southeast coastline changed randomly between 0 and the maximum value of 128.23 km, which is the theoretical limit of visibility considering the height of the ship. d = 0 km and d = 128.23 km mean that the weather is extremely cloudy/foggy (ρ = 8 okta) and absolutely clear (ρ = 0 okta), respectively.

The start time of diurnal sailing changed randomly between 0 and 1 hour after sunrise, instead of the earlier constant zero value (start at sunrise).

The coastline of Greenland as the direct goal of voyages became limited to a southeast section, opposite to the earlier south borderline (Fig 1).

The aim of this work was to determine the sensitivity/robustness of the success of sky-polarimetric Viking navigation (SPVN) to variables sunstone type, sailing date, navigation periodicity Δt, night sailing, cloudiness dominance mdominance and cloudiness changeability σchangeability. Száz and Horváth [34] have already shown that SPVN is robust (almost insensitive) against the change of sunstone type, but depends strongly on the sailing date and navigation periodicity. Using a more sophisticated, improved simulation of SPVN, here we corroborated these earlier findings and demonstrated that SPVN is also very sensitive to night sailing, but is robust against the changing weather conditions (dominance and changeability of cloudiness).

If we look at the north error dataset published in [31] and used also in this work, we can see that there are periods within a day, when the sign of this error changes from negative to positive. Negative and positive north error means that the ship will steer more towards south and north, respectively. The success of a voyage is basically depending on whether or not the ship sails longer periods with southern or northern components of the moving direction. Thus, the ratio between these two direction increments will give the net average direction. For example, if a simulated navigator measures the intended sailing direction when the simulation uses a randomly set negative north error just before a period with a positive north error, then the ship will move towards south, and it will only move towards north again after a new measurement of the sailing direction. This situation can create such simulation days, when the ship advances more towards south than north (Fig 7). This is the reason why the navigation periodicity, sailing date and night sailing have a profound effect on the sailing success (Fig 7).

Fig 7

North error accumulated on a day of sailing versus time in a day at solstice (A, B, C) and equinox (D) using calcite sunstone.

The black curve on each figure represents the north error averaged for cloudiness as a function of time. Each rectangular column represents a sailing period where the area of the rectangle is the accumulated north error for a measurement cycle. The area under the curve AUC is the averaged north error a simulated Viking route suffers on that day (= sum of the area of rectangles).

In our simulations, variable mdominance drives the increment/decrement of cloudiness ρ (okta) in skewing the change of ρ-values towards more or less negative or positive changes. In order that this scheme can work without automatically converging to and then getting stuck at zero cloudiness (ρ = 0 okta) or total cloudiness (ρ = 8 okta), variable σchangeability is needed (with an appropriately greater range of values than that of mdominance) to allow variation to the preset cloudiness dominance. As examples, Fig 8 shows the cloudiness distributions during simulated voyages using cloudiness dominances mdominance = -1, 0 or +1 (meaning dominantly weakly, medium or strongly cloudy skies) and cloudiness changeabilities σchangeability = 1, 2 or 4 (meaning weakly/slowly, medium or strongly/rapidly changing cloudiness).

Fig 8

Cloudiness distributions of Viking sailing routes.

Cloudiness distributions during simulated voyages using calcite sunstones at summer solstice with navigation periodicity Δt = 1 hour for cloudiness dominances mdominance = -1, 0 or 1 (meaning dominantly weakly, medium or strongly cloudy skies) and for cloudiness changeabilities σchangeability = 1, 2 or 4 (meaning weakly/slowly, medium or strongly/rapidly changing cloudiness).

It is striking that the success rate s of simulated sailing routes are practically independent of the dominance mdominance and changeability σchangeability of cloudiness (Figs 2B and 4–6). This can partly be explained with the cumulated navigation (north) error due to cloudiness having only a low scatter near the end point of the sailing route so that it practically does not change the success rate s. Note that s depends mainly on the (a)symmetry of the particular north error cumulated in the morning and the afternoon of every day. Using the same data series as we in this work, Száz et al. [31] have already studied the dependence of the north error on the cloudiness ρ, taking into consideration the solar elevation and the sailing date (spring equinox or summer solstice). They also found (see Fig 4 of [31]) that both the net (cumulated) morning and the net afternoon north errors are practically independent of ρ, but depend on the sailing date (spring equinox or summer solstice), the solar elevation and the time of day (morning or afternoon). Hence, since the north error does not vary significantly with changing ρ, ρ has no great effect on the morning-afternoon asymmetry, and therefore it influences only slightly the success rate s of SPVN.

The night sailing results in a big difference compared to the situation without night sailing, since at night the ship advances towards the last direction measured immediately before sunset. This could create situations when the ship steers towards a bad direction (deviating from the intended sailing direction), however, it can also create situations when the direction errors made during the day are actually compensated. Fig 1B, 1D and 1F show that the simulated voyages with night sailing headed more northward compared to those without night sailing (Fig 1A, 1C and 1E). The reason for this is that in afternoon the direction error heads tendentiously northward [31], and with night sailing the simulated ship continued its way to this northward direction all night.

The graphs in Figs 2 and 3A–3C show smaller- or higher-amplitude oscillation of the success rate s versus the navigation periodicity Δt for certain periods with lower and higher local extrema. The main reason for this phenomenon is the alternating sign of the strongly oscillating north error accumulated on a day (Fig 7): At higher navigation periodicities than 1 hour, the absolute accumulated north error can greatly vary, depending on the actual point of time of navigation. Thus, for a certain periodicity Δt the compensation for the accumulated north error can be good, while for a slightly different Δt-value it can break down quickly by accumulating north error skewed into one of northern or southern directions.

In the future, it could be studied how a progressive restriction of the coastline section defined as successful destination influences the success rate s of Viking navigation. It is expected that s gradually decreases as the desired destination becomes more and more localized. Thus, the chance of a precise localization of the destination is small in all probability.

Conclusion

We found that among the investigated meteorological parameters and navigation/sailing variables, the success of SPVN is robust against the sunstone type, as well as the cloudiness dominance and changeability. Contrary to this, SPVN is sensitive to the sailing date, navigation periodicity and night sailing, which can optimally be chosen/selected by the sailors/navigators. Remarkably, the accuracy of this navigation method is practically not affected by the dominance and changeability of cloudiness, though one could expect that weather plays the largest role in the sailing success.

29 Oct 2021

16 Nov 2021

Our detailed response to the comments of the two reviewers is uploaded as a separate file (+VikingNavigation_PLoS-One_response.docx).

Submitted filename: +VikingNavigation_PLoS-One_response.docx

Click here for additional data file.

5 Jan 2022

Sensitivity and robustness of sky-polarimetric Viking navigation Sailing success is most sensitive to night sailing, navigation periodicity and sailing date, but robust against weather conditions

PONE-D-21-15021R1

Dear Dr. Horvath,

We’re pleased to inform you that your manuscript has been judged scientifically suitable for publication and will be formally accepted for publication once it meets all outstanding technical requirements.

Within one week, you’ll receive an e-mail detailing the required amendments. When these have been addressed, you’ll receive a formal acceptance letter and your manuscript will be scheduled for publication.

An invoice for payment will follow shortly after the formal acceptance. To ensure an efficient process, please log into Editorial Manager at http://www.editorialmanager.com/pone/, click the 'Update My Information' link at the top of the page, and double check that your user information is up-to-date. If you have any billing related questions, please contact our Author Billing department directly at authorbilling@plos.org.

If your institution or institutions have a press office, please notify them about your upcoming paper to help maximize its impact. If they’ll be preparing press materials, please inform our press team as soon as possible -- no later than 48 hours after receiving the formal acceptance. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information, please contact onepress@plos.org.

Kind regards,

Simone Lolli

Academic Editor

PLOS ONE

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. If the authors have adequately addressed your comments raised in a previous round of review and you feel that this manuscript is now acceptable for publication, you may indicate that here to bypass the “Comments to the Author” section, enter your conflict of interest statement in the “Confidential to Editor” section, and submit your "Accept" recommendation.

Reviewer #1: All comments have been addressed

Reviewer #2: All comments have been addressed

**********

2. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #1: Yes

Reviewer #2: Yes

**********

3. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #1: Yes

Reviewer #2: Yes

**********

4. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.

Reviewer #1: Yes

Reviewer #2: Yes

**********

5. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #1: Yes

Reviewer #2: Yes

**********

6. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #1: Congratulations - this has been a valuable contribution: an interesting paper, well planned and executed..

Reviewer #2: The authors have clearly addressed all outstanding questions and minor issues in a satisfying way, the paper in its now improved form is most definitely suitable for publication.

**********

7. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: Yes: Victor Benno MEYER-ROCHOW

Reviewer #2: No

10 Jan 2022

PONE-D-21-15021R1

Sensitivity and robustness of sky-polarimetric Viking navigation: Sailing success is most sensitive to night sailing, navigation periodicity and sailing date, but robust against weather conditions

Dear Dr. Horvath:

I'm pleased to inform you that your manuscript has been deemed suitable for publication in PLOS ONE. Congratulations! Your manuscript is now with our production department.

If your institution or institutions have a press office, please let them know about your upcoming paper now to help maximize its impact. If they'll be preparing press materials, please inform our press team within the next 48 hours. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information please contact onepress@plos.org.

If we can help with anything else, please email us at plosone@plos.org.

Thank you for submitting your work to PLOS ONE and supporting open access.

Kind regards,

PLOS ONE Editorial Office Staff

on behalf of

Dr. Simone Lolli

Academic Editor

PLOS ONE