| Literature DB >> 28386426 |
Gábor Horváth1, Péter Takács1, Balázs Kretzer1, Szilvia Szilasi1, Dénes Száz1, Alexandra Farkas2, András Barta3.
Abstract
If a human looks at the clear blue sky from which light with high enough degree of polarization d originates, an 8-shaped bowtie-like figure, the yellow Haidinger's brush can be perceived, the long axis of which points towards the sun. A band of high d arcs across the sky at 90° from the sun. A person can pick two points on that band, observe the yellow brushes and triangulate the position of the sun based on the orientation of the two observed brushes. This method has been suggested to have been used on the open sea by Viking navigators to determine the position of the invisible sun occluded by cloud or fog. Furthermore, Haidinger's brushes can also be used to locate the sun when it is below the horizon or occluded by objects on the horizon. To determine the position of the sun using the celestial polarization pattern, the d of the portion of the sky used must be greater than the viewer's degree of polarization threshold d* for perception of Haidinger's brushes. We studied under which sky conditions the prerequisite d > d* is satisfied. Using full-sky imaging polarimetry, we measured the d-pattern of skylight in the blue (450 nm) spectral range for 1296 different meteorological conditions with different solar elevation angles θ and per cent cloud cover ρ. From the measured d-patterns of a given sky we determined the proportion P of the sky for which d > d*. We obtained that P is the largest at low solar elevations θ ≈ 0° and under totally or nearly clear skies with cloud coverage ρ = 0%, when the sun's position is already easily determined. If the sun is below the horizon (-5° ≤ θ < 0°) during twilight, P = 76.17 ± 4.18% for [Formula: see text] under clear sky conditions. Consequently, the sky-polarimetric Viking navigation based on Haidinger's brushes is most useful after sunset and prior to sunrise, when the sun is not visible and large sky regions are bright, clear and polarized enough for perception of Haidinger's brushes.Entities:
Keywords: Haidinger's brushes; Viking navigation; human polarization vision; imaging polarimetry; sky polarization; sunstone
Year: 2017 PMID: 28386426 PMCID: PMC5367315 DOI: 10.1098/rsos.160688
Source DB: PubMed Journal: R Soc Open Sci ISSN: 2054-5703 Impact factor: 2.963
Figure 1.Principle of the sky-polarimetric Viking navigation on the basis of Haidinger's brushes. (a) Clear sky. (b) Cloudy sky. Yellow bars: local direction of polarization of skylight. Red double-headed arrow: axis of mirror symmetry of the yellow Haidinger's brush being perpendicular to the direction of polarization. The long axis of the blue Haidinger's brush is parallel to the direction of polarization. Owing to Rayleigh scattering of sunlight in the atmosphere, the celestial great circle passing through the observed sky point parallel to the yellow Haidinger's brush crosses the sun.
Figure 2.(a)(A) A clear sky photographed on 7 June 2014 at 18.36.29 (=GMT + 2 h, where GMT means Greenwich Mean Time) with solar elevation θ = 0° and per cent cloud cover ρ = 0%. (B) Same as A where the white areas are under- and/or overexposed regions as well as landmarks on the horizon, which are not taken into consideration. The white sky regions were masked out during calculations. (C) Pattern of the degree of linear polarization d of skylight measured with full-sky imaging polarimetry in the blue (450 nm) spectral range. (D) Points of the clear sky are marked with black where . (b) Same as on the left for a cloudy sky photographed on 11 January 2015 at 15.48.59 (=GMT + 1 h) with solar elevation θ = 0° and per cent cloud cover ρ = 100%. (D) The circular area is blank, because there is no celestial point with d > 23%.
Figure 3.Proportion P (%) of the clear and overcast skies in figure 2 with degrees of skylight polarization not lower than d in the blue (450 nm) spectral range. The vertical lines show the values d* = 23, 56, 87% being the three characteristic perception thresholds of Haidinger's brushes.
of P, where P is the proportion of the sky for which d > d* (= 23, 56, 87%).
measures the appropriateness of Viking navigation with Haidinger's brushes under a meteorological situation belonging to a given cell (ρ, θ). In figure 4b, each cell of a given matrix contains the relative (or normalized) standard deviation ΔP/(ΔP)max of P. Tables 1–3 and 4–6 contain the numerical data visualized by grey shades in the matrices of figure 4a and b, respectively.
Figure 4.(a) The (θ, ρ) matrices calculated for 1296 meteorological situations (characterized by sun elevation θ and per cent cloud cover ρ) in the blue (450 nm) spectral range for the perception thresholds of Haidinger's brushes , and . In every cell of a given matrix there are 12 different skies. Each cell contains the average of P, where P is the proportion of the sky for which d > d* (=23, 56, 87%). (b) As (a) for the ΔP/(ΔP)max(θ, ρ) matrices, each cell of which contains the relative (normalized) standard deviation ΔP/(ΔP)max of P.
Numerical values (%) of the (θ, ρ) matrix in the blue (450 nm) spectral range for shown in figure 4a. The 12 intervals of solar elevation θ are: −5° ≤ θ1 < 0°, 0° ≤ θ2 < 5°, 5° ≤ θ3 < 10°, 10° ≤ θ4 < 15°, 15° ≤ θ5 < 20°, 20° ≤ θ6 < 25°, 25° ≤ θ7 < 30°, 30° ≤ θ8 < 35°, 35° ≤ θ9 < 40°, 40° ≤ θ10 < 45°, 45° ≤ θ11 < 50°, 50° ≤ θ12 ≤ 55°. The nine intervals of cloud per cent ρ are: ρ0 = 0% (clear sky), 0% < ρ1 < 12.5%, 12.5% ≤ ρ2 < 25%, 25% ≤ ρ3 < 37.5%, 37.5% ≤ ρ4 < 50%, 50% ≤ ρ5 < 62.5%, 62.5% ≤ ρ6 < 75%, 75% ≤ ρ7 < 87.5%, 87.5% ≤ ρ8 ≤ 100% (nearly and totally overcast sky).
| 1.16 | 4.53 | 1.46 | 1.97 | 4.43 | 1.36 | 2.51 | 2.54 | 4.99 | 2.04 | 7.14 | 8.11 | |
| 3.37 | 4.91 | 4.10 | 1.72 | 9.95 | 10.93 | 8.09 | 10.23 | 11.73 | 9.73 | 13.54 | 4.29 | |
| 26.44 | 16.76 | 12.64 | 15.90 | 10.48 | 14.48 | 18.03 | 8.82 | 12.78 | 20.68 | 21.18 | 30.47 | |
| 45.24 | 29.65 | 9.95 | 27.76 | 30.51 | 16.91 | 23.65 | 23.92 | 19.16 | 26.48 | 26.10 | 29.93 | |
| 39.87 | 34.17 | 48.33 | 42.58 | 31.78 | 31.56 | 30.61 | 34.07 | 16.21 | 20.20 | 30.90 | 29.32 | |
| 62.32 | 59.97 | 59.23 | 56.16 | 48.26 | 42.36 | 36.07 | 22.97 | 30.15 | 34.93 | 29.14 | 24.48 | |
| 72.20 | 72.05 | 64.12 | 69.28 | 59.37 | 53.03 | 43.55 | 33.01 | 40.63 | 45.95 | 34.79 | 33.59 | |
| 76.99 | 74.76 | 72.55 | 75.80 | 69.12 | 67.98 | 59.98 | 52.76 | 49.94 | 57.47 | 51.95 | 44.57 | |
| 76.17 | 75.80 | 75.37 | 74.46 | 69.52 | 69.37 | 63.30 | 59.61 | 45.51 | 41.20 | 54.71 | 46.98 |
Numerical values (%) of the (θ, ρ) matrix in the blue (450 nm) spectral range for shown in figure 4a.
| 0.0047 | 0.0169 | 0.0051 | 0.0044 | 0.0070 | 0.0010 | 0.0053 | 0.0089 | 0.0044 | 0.0014 | 0.0098 | 0.0031 | |
| 0.0079 | 0.0070 | 0.0051 | 0.0031 | 0.0110 | 0.0128 | 0.0067 | 0.0119 | 0.0118 | 0.0052 | 0.0069 | 0.0011 | |
| 0.0160 | 0.0178 | 0.0168 | 0.0101 | 0.0111 | 0.0139 | 0.0182 | 0.0103 | 0.0105 | 0.0376 | 0.0413 | 0.0164 | |
| 0.0158 | 0.0130 | 0.0084 | 0.0103 | 0.0071 | 0.0112 | 0.0169 | 0.0197 | 0.0217 | 0.0454 | 0.0642 | 0.0370 | |
| 0.0235 | 0.0131 | 0.0160 | 0.0128 | 0.0037 | 0.0224 | 0.0230 | 0.0249 | 0.0097 | 0.0255 | 0.0640 | 0.0629 | |
| 0.0288 | 0.0288 | 0.0196 | 0.0232 | 0.0206 | 0.0141 | 0.0210 | 0.0083 | 0.0134 | 0.0579 | 0.0317 | 0.0365 | |
| 0.0354 | 0.0360 | 0.0385 | 0.0172 | 0.0710 | 0.0105 | 0.0128 | 0.0097 | 0.0029 | 0.0402 | 0.0296 | 0.0191 | |
| 0.0848 | 0.0451 | 0.0450 | 0.0189 | 0.0683 | 0.0417 | 0.0216 | 0.0162 | 0.0013 | 0.0297 | 0.0197 | 0.0328 | |
| 0.1052 | 0.1375 | 0.0653 | 0.0397 | 0.0125 | 0.0109 | 0.0268 | 0.0377 | 0.0038 | 0.0018 | 0.0551 | 0.0082 |
Numerical values (%) of the ΔP(θ, ρ) matrix in the blue (450 nm) spectral range for shown in figure 4b.
| 0.46 | 4.65 | 0.71 | 1.22 | 6.13 | 0.47 | 2.16 | 3.13 | 2.75 | 1.63 | 7.09 | 14.96 | |
| 6.84 | 7.01 | 3.71 | 1.15 | 10.70 | 9.85 | 6.28 | 8.85 | 7.92 | 9.77 | 13.27 | 3.98 | |
| 26.00 | 20.01 | 13.07 | 15.65 | 12.87 | 10.82 | 12.22 | 7.67 | 8.92 | 12.57 | 11.55 | 18.51 | |
| 25.60 | 19.83 | 18.37 | 17.07 | 20.88 | 15.65 | 13.62 | 10.43 | 14.52 | 14.91 | 10.37 | 12.09 | |
| 29.66 | 18.36 | 10.42 | 12.40 | 9.87 | 17.58 | 17.39 | 15.88 | 9.53 | 9.33 | 7.71 | 11.19 | |
| 22.81 | 19.19 | 7.60 | 9.05 | 7.44 | 12.48 | 17.65 | 15.86 | 9.85 | 20.55 | 13.27 | 10.29 | |
| 4.07 | 6.58 | 7.33 | 3.27 | 7.20 | 10.12 | 16.03 | 14.95 | 7.98 | 11.04 | 9.84 | 10.45 | |
| 2.63 | 3.08 | 2.68 | 1.25 | 6.07 | 4.866 | 9.12 | 7.49 | 5.57 | 3.18 | 1.98 | 6.90 | |
| 4.18 | 1.97 | 2.66 | 2.41 | 3.21 | 2.12 | 3.26 | 5.63 | 8.91 | 3.73 | 4.69 | 3.52 |
Numerical values (%) of the ΔP(θ, ρ) matrix in the blue (450 nm) spectral range for shown in figure 4b.
| 0.0036 | 0.0145 | 0.0040 | 0.0057 | 0.0101 | 0.0012 | 0.0129 | 0.0093 | 0.0038 | 0.0026 | 0.0069 | 0.0028 | |
| 0.0060 | 0.0050 | 0.0039 | 0.0024 | 0.0136 | 0.0212 | 0.0058 | 0.0242 | 0.0130 | 0.0064 | 0.0177 | 0.0011 | |
| 0.0065 | 0.0244 | 0.0143 | 0.0101 | 0.0156 | 0.0265 | 0.0278 | 0.0088 | 0.0121 | 0.0624 | 0.0517 | 0.0295 | |
| 0.0095 | 0.0102 | 0.0104 | 0.0039 | 0.0047 | 0.0160 | 0.0143 | 0.0198 | 0.0248 | 0.0392 | 0.1011 | 0.0390 | |
| 0.0222 | 0.0128 | 0.0045 | 0.0069 | 0.0036 | 0.0243 | 0.0320 | 0.0412 | 0.0195 | 0.0179 | 0.0403 | 0.0745 | |
| 0.0165 | 0.0108 | 0.0080 | 0.0261 | 0.0150 | 0.0152 | 0.0188 | 0.0107 | 0.0379 | 0.0518 | 0.0391 | 0.0426 | |
| 0.0207 | 0.0252 | 0.0201 | 0.0061 | 0.0630 | 0.0138 | 0.0137 | 0.0132 | 0.0040 | 0.0297 | 0.0301 | 0.0234 | |
| 0.0393 | 0.0223 | 0.0186 | 0.0055 | 0.0357 | 0.0232 | 0.0144 | 0.0287 | 0.0017 | 0.0134 | 0.0135 | 0.0265 | |
| 0.1315 | 0.2681 | 0.0293 | 0.0324 | 0.0097 | 0.0034 | 0.0194 | 0.0731 | 0.0103 | 0.0021 | 0.0371 | 0.0034 |
(θ, ρ) matrices calculated for 1296 meteorological situations (characterized by sun elevation θ and per cent cloud cover ρ) in the blue (450 nm) spectral range for the perception thresholds of Haidinger's brushes , and . In every cell of a given matrix there are 12 different skies. Each cell contains the average
of P, where P is the proportion of the sky for which d > d* (=23, 56, 87%). (b) As (a) for the ΔP/(ΔP)max(θ, ρ) matrices, each cell of which contains the relative (normalized) standard deviation ΔP/(ΔP)max of P.
(θ, ρ) matrix in the blue (450 nm) spectral range for shown in figure 4a. The 12 intervals of solar elevation θ are: −5° ≤ θ1 < 0°, 0° ≤ θ2 < 5°, 5° ≤ θ3 < 10°, 10° ≤ θ4 < 15°, 15° ≤ θ5 < 20°, 20° ≤ θ6 < 25°, 25° ≤ θ7 < 30°, 30° ≤ θ8 < 35°, 35° ≤ θ9 < 40°, 40° ≤ θ10 < 45°, 45° ≤ θ11 < 50°, 50° ≤ θ12 ≤ 55°. The nine intervals of cloud per cent ρ are: ρ0 = 0% (clear sky), 0% < ρ1 < 12.5%, 12.5% ≤ ρ2 < 25%, 25% ≤ ρ3 < 37.5%, 37.5% ≤ ρ4 < 50%, 50% ≤ ρ5 < 62.5%, 62.5% ≤ ρ6 < 75%, 75% ≤ ρ7 < 87.5%, 87.5% ≤ ρ8 ≤ 100% (nearly and totally overcast sky).
(θ, ρ) matrix in the blue (450 nm) spectral range for shown in figure 4a.
(θ, ρ) matrix in the blue (450 nm) spectral range for shown in figure 4a.
has relatively large values (
≤ 76.99%) for ρ ≤ 75%, especially at solar elevations θ ≤ 35° (figure 4a, table 1). Therefore, in many such meteorological situations this person is able to perceive Haidinger's brushes and to guess the position of the sun occluded by clouds in relatively large sky regions. For an average person with a perception threshold of Haidinger's brushes,
has the largest values (
≤ 28.14%) for completely or nearly clear skies with ρ ≤ 25% and lower sun elevations θ ≤ 15° (figure 4a, table 2). A person having low polarization sensitivity (e.g. ) can only detect Haidinger's brushes under very few meteorological conditions (figure 4a, table 3;
≤ 0.1375%, being practically zero) and would, therefore, be an unlikely choice for a Viking navigator.
Numerical values (%) of the (θ, ρ) matrix in the blue (450 nm) spectral range for shown in figure 4a.
| 0.31 | 0.39 | 0.27 | 0.29 | 0.32 | 0.21 | 0.29 | 0.30 | 0.25 | 0.20 | 0.36 | 0.27 | |
| 0.30 | 0.31 | 0.29 | 0.27 | 0.35 | 0.32 | 0.26 | 0.33 | 0.31 | 0.29 | 0.28 | 0.21 | |
| 0.89 | 0.53 | 0.34 | 0.48 | 0.29 | 0.36 | 0.35 | 0.34 | 0.36 | 0.55 | 0.59 | 0.56 | |
| 4.60 | 1.21 | 0.31 | 0.28 | 0.38 | 0.25 | 0.38 | 0.42 | 0.53 | 0.74 | 0.93 | 0.63 | |
| 8.77 | 2.04 | 1.41 | 0.78 | 0.25 | 0.49 | 0.58 | 0.66 | 0.31 | 0.57 | 0.76 | 0.81 | |
| 12.76 | 7.14 | 2.97 | 1.52 | 0.44 | 0.47 | 0.41 | 0.48 | 0.35 | 1.31 | 0.65 | 0.58 | |
| 16.72 | 16.73 | 6.46 | 4.82 | 3.22 | 1.48 | 0.40 | 0.46 | 0.41 | 0.77 | 0.56 | 0.50 | |
| 27.79 | 23.64 | 15.39 | 26.65 | 10.48 | 6.93 | 3.10 | 0.73 | 0.67 | 3.05 | 1.44 | 0.69 | |
| 26.40 | 28.14 | 26.51 | 22.85 | 11.19 | 6.99 | 2.40 | 5.08 | 0.29 | 0.21 | 0.74 | 0.30 |
are practically zero (ΔP ≤ 0.2681%; figure 4b, table 6). The larger the standard deviation ΔP at a given meteorological situation (θ–ρ cell), the smaller is the reliability of the navigation method based on Haidinger's brushes under that sky condition.
Numerical values (%) of the ΔP(θ, ρ) matrix in the blue (450 nm) spectral range for shown in figure 4b.
| 0.0597 | 0.1204 | 0.0506 | 0.0729 | 0.0700 | 0.0746 | 0.1226 | 0.0688 | 0.0554 | 0.0500 | 0.0856 | 0.0513 | |
| 0.0600 | 0.0381 | 0.0473 | 0.0376 | 0.1199 | 0.2273 | 0.1082 | 0.2035 | 0.1325 | 0.1540 | 0.1660 | 0.0392 | |
| 1.4282 | 0.4922 | 0.0736 | 0.4471 | 0.1367 | 0.3437 | 0.1912 | 0.1696 | 0.2246 | 0.5242 | 0.3197 | 0.5472 | |
| 6.0898 | 2.4652 | 0.0492 | 0.0773 | 0.3000 | 0.1133 | 0.2028 | 0.1513 | 0.4601 | 0.5565 | 0.9953 | 0.2599 | |
| 14.057 | 3.8582 | 1.8149 | 0.9965 | 0.1290 | 0.3830 | 0.3181 | 0.4286 | 0.1642 | 0.1898 | 0.2064 | 0.4319 | |
| 11.902 | 8.7061 | 4.0724 | 2.5452 | 0.1976 | 0.4637 | 0.2381 | 0.2886 | 0.3158 | 1.2297 | 0.3296 | 0.2045 | |
| 7.5923 | 9.9259 | 6.1184 | 3.1408 | 5.2772 | 2.7204 | 0.3430 | 0.4612 | 0.8336 | 0.5078 | 0.3720 | 0.1434 | |
| 9.4660 | 9.6793 | 8.6202 | 6.4879 | 11.2379 | 6.0032 | 4.7326 | 0.9265 | 0.7885 | 2.3646 | 1.4952 | 0.3427 | |
| 10.2622 | 7.7071 | 7.4711 | 8.4811 | 7.93141 | 3.7781 | 4.1875 | 6.5901 | 0.2668 | 0.0396 | 0.1969 | 0.0774 |
of P (averaged for 12 different skies possessing very similar values of ρ and θ) is a good measure of the appropriateness of a given sky for Viking navigation with Haidinger's brushes.
and its standard deviation ΔP were determined for all cells of the (θ, ρ) matrices (figure 4, tables 1–6) with the use of the patterns of the degree of polarization d measured with full-sky imaging polarimetry in the blue part of the spectrum. Although we measured the sky polarization for 1296 meteorological situations with different solar elevations θ and cloudinesses ρ in Hungary at latitude 47°15'29.83″ N, this might not affect considerably the sky polarizaton, because for a given ρ, celestial polarization depends predominantly on θ, and we restricted our measurements to solar elevations –5° ≤ θ ≤ +55° occurring at 61° latitude, where the main Viking sailing route ran. The range 0 okta ≤ ρ ≤ 8 oktas used by us covered the whole range of cloudiness that can also occur at 61° latitude.
of P is large only for the most sensitive-eyed human observers (Viking navigators) having a perception threshold of Haidinger's brushes. Viking navigators were probably those that had the highest polarization sensitivity (i.e. having the lowest threshold d*). A diet rich in fish, which is a good source of macular carotenoids, may have contributed to Vikings' general ability to see Haidinger's brushes (S. Temple 2016, personal communication).
of the sky suitable for sky-polarimetric navigation is between 26.4% and 76.2%, depending on the cloud cover (ρ < 75%).