1 Introduction

Solar eclipses have fascinated humankind since the beginning of recorded history. As valuable resource for historians they allow certain historical events to be dated precisely thereby improving the accuracy of historical chronology. For example a solar eclipse of June 15 in 763 BC mentioned in an Assyrian eponym list (Rawlinson 1867) is an important timestamp for the chronology of the Ancient Orient. In pre-scientific times, before chemistry separated from alchemy, solar eclipses (aka “chemical weddings”) were thought to be special insofar that certain (al-)chemical reactions would only take place during such events (Huang 2014). In the beginning of the 20th century some anecdotal evidence (Kolisko 1927) suggested a possible correlation between a solar eclipse and the capillary forces of an aqueous gold colloid solution; in 1975 investigations by Schütz and Schütz (1975) suggested a correlation between a solar eclipse and the mixing dynamics of water and glycerol. It should be mentioned that the solar eclipse was not visible at the location where the experiments took place, though. According to the knowledge of the authors none of the above mentioned experiments satisfy the standards of modern science and should therefore not be considered as scientific literature. However, they can serve as inspiration to explore actual effects that occur during a solar eclipse next to the obvious changes in illumination, temperature and weather conditions. Apart from changes to water physics, various experiments have been performed in order to detect anomalous gravitational behavior during solar eclipses. Some observations gave positive results and some failed to detect any noticeable effect. Allais reported in 1954 and in 1959 that the plane and/or amplitude of oscillation of a pendulum shifted during an eclipse (Allais 1975, 1999); however, an investigation by Slichter et al. (1965) failed to report any gravitational effect. In Saxl and Allen (1970) again confirmed the observations of Allais (1970), whereas in 1990 Kuusela confirmed the conclusion of Slichter reporting that no effect could be observed (Kuusela 1991). In Van Flandern and Yang (2003) provided a possible explanation for these effects. Rapid air movement for the bulk of the atmosphere above normal cloud levels are said to create a rise of atmospheric pressure causing the anomalies. The topic of gravitational effects during solar eclipses therefore seems to be explained as barometric variation effect.

1.1 Effects on Earth’s Ozone Layer

There are several studies claiming that a solar eclipse has an impact on Earth’s ozone layer (Bezverkhny et al. 1956; Fournier 1956; Kawabata1936; Khrigian et al. 1961; Jerlov et al. 1954; Pariisky and Chen-Chao 1958; Steblova 1961; Stranz 1961; Svensson 1968; Chakrabarty et al. 1997; Zerefos et al. 1999). These claims are based upon measurements of the amount of UV-B radiation on the surface during the eclipse, and they vary between a 40% decrease (Jerlov et al. 1954) to a 20–30% increase (Khrigian et al. 1961; Steblova 1961) and even an unlikely increase of 300–400% (Bezverkhny et al. 1956). However, comparable to the gravitational study situation, there are papers rejecting any relation between the solar eclipse and the ozone layer. Initially in Bojkov (1968) and subsequently in Blumthaler et al. (2006) showed convincingly that the differences observed can be explained by limb darkening (Bojkov 1968; Blumthaler et al. 2006), an effect that deals with the difference in spectral composition of radiation coming from either the limb or the center of the solar disk. Therefore, any secondary effects of a solar eclipse based upon changes in the ozone layer should be treated with caution. Nevertheless, some of these investigation should be mentioned here, for example the study by Chakrabarty et al. (1997) of the fluctuation in ozone column over Ahmedabad during the solar eclipse of 24 October 1995 , and the study by Mims and Mims about the same topic during the solar eclipse of July 11th, 1991 (Mims and Mims 1993). Whereas the interpretation of their results might be disputable, their measurements show effects with maxima 10 min and 25 min after 3rd and 2nd contact, respectively, which is a similar time window compared to the findings in this work.

1.2 Environmental Studies

Solar eclipses can have effects on terrestrial winds, temperature and weather in general. For example during the 11 August 1999 total solar eclipse transient a temperature decreases of up to 3 °C and a mean regional wind speed decrease of 0.7 m s−1 during the maximum eclipse hour with a mean anticlockwise wind direction change of 17° were observed (Gray and Harrison 2012). Also the electric parameters of the atmosphere change: the ambient electric field dropped by up to 65% during the eclipse, and potential gradient showed epochs of enhancements during and after the eclipse until after Sunset (Anil Kumar et al. 2013). The authors conclude that during the course of an eclipse, radiative cooling enables the formation and development of water droplets resulting in scavenging or coalescence of free ions/free space charges in the atmosphere causing the decreases in electric currents and field strengths with the reduction in galactic cosmic rays flux as possible additional effect. Moreover it has been shown that the pH value of seawater during solar eclipse drops by 20% of the difference between ordinary water and seawater (Santhosh et al. 2014); and this drop is said to be due to the disputed effect of an eclipse on the ozone layer (see 1.1). Also effects on the marine microbiome have been reported (Backus et al. 1963). For example, solar eclipse induced effects on the marine environment have been shown to be driven mainly by the change of underwater irradiance (Economou et al. 2006). In another study it was shown that locking of the Sun ray during an eclipse at sea seems to promote the progeny of microbial predators in the race of better acclimatization and survival in the natural and changing environmental conditions (Shriyan et al. 2011).

2 Theoretical Aspects

2.1 Surface Tension

Surface tension is the elastic tendency of a fluid surface to create the least surface area possible. The field dates back to the days of Galileo (Galilei 1612); the first papers were written by Young (1805); a good overview can be found in the review by Good (1992). Briefly, surface tension is caused by cohesive forces among liquid molecules. In the bulk each molecule is pulled equally in every direction by neighboring liquid molecules, resulting in a net force of zero, whereas at the surface there is a net force pointing inwards. This force creates internal pressure and is responsible for liquid surfaces to contract to the smallest area possible. One method to measure this force is called the “Wilhelmy plate method” (Wilhelmy 1863; Adamson and Gast (1967). Here a thin plate is submerged into the liquid perpendicular to the interface and carefully pulled out until the retracting force reaches a maximum. This position is then kept constant, and the force exerted on it is measured with a microbalance using the Wilhelmy equation,

$$\gamma = \frac{F}{lcos\theta }$$
(1)

where γ is the surface tension, F the force, l the wetted perimeter and θ the contact angle between liquid phase and the plate. This method has been used successfully for several decades (Neumann et al. 1979; Jordan and Lane 1964; Lane and Jordan 1970; Princen 1970; Penn and Miller 1980; Bayramli et al. 1981; Seebergh and Berg 1992; Hayes and Ralston 1992, 1993a, b; Sauer and Kampert 1998). The apparatus used in this study, a Kibron AlphaPi Plus automatic tensiometer, uses a small needle instead of a plate and assumes complete wetting (θ = 0), which simplifies the measurement so that the downwards force measured by the instrument’s built-in micro balance becomes directly proportional to the surface tension.

2.2 Astronomical Facts About Solar Eclipses

A solar eclipse is a perceived darkening of the Sun that occurs when the Moon passes between the Sun and Earth whereby fully or partially blocking the Sun’s light. Such a situation can only occur at new Moon, when the Sun and the Moon are in conjunction as seen from the Earth in an alignment astronomically referred to as syzygy. In a total eclipse, the disc of the Sun is fully covered by the disc of the Moon. In partial and annular eclipses, only part of the Sun is obscured. The event is therefore characterized by the optical coverage C, the ratio of the area of the solar disc occulted by the lunar disc compared to the total area of the solar disc, C = ASun,obscured/ASun,total. A complete derivation of this parameter for the physics relevant for the MSW effect is given elsewhere (Narayan et al. 1999). Other important parameters for an eclipse are the ratio of the apparent radii of the Moon and the Sun’s discs seen from Earth, a = rm/rs and the impact parameter p defined as the apparent distance of closest approach between the center of the Sun’s disc and the center of the Moon’s disc. An eclipse is total if a ≥ 1 and it is annular (as in the present case) if a < 1. For an annular solar eclipse the Moon has to be in or close to its apogee. The progress of an eclipse can be divided into four stages:

  1. (I)

    First contact: the Moon’s disc starts to impinge on the Sun’s disc (C > 0)

  2. (II)

    Second contact: the Moon’s disc has moved completely within the Sun’s disc (C = 100 for a total eclipse, and C = 94.33%, C = 97.31% and C = 35.52 for the annular eclipses in Réunion and Puerto Chacabuco and the partial eclipse in Rothera, respectively)

  3. (III)

    Third contact: the Moon’s disc starts to move out of the Sun’s disc

  4. (IV)

    Fourth contact: the Moon’s disc leaves the Sun’s disc entirely

These stages are marked in the following figures using the same Latin numbers as above. For the annular eclipse studied in this work the coverage C can be calculated as follows (Narayan et al. 1999):

$$C\left( d \right) = f\left( u \right) + a^{2} \left[ {1 - f\frac{u}{a}} \right]$$
(2)

where

$$f\left( u \right) = \frac{1}{\pi }\left( {\arcsin u - u\sqrt {1 - u^{2} } } \right),$$
(3)

and d is the apparent separation between the centers of the Sun and the Moon (in units of solar radii),

$$d\left( t \right) = \left\{ {p\left[ {1 - \left( {\frac{2t}{T}} \right)^{2} } \right] + \left( {1 + a^{2} } \right)\left( {\frac{2t}{T}} \right)^{2} } \right\}^{{\frac{1}{2}}}$$
(4)

with t being the time (set to 0 for the point of maximum coverage) and T the time between first and fourth contact.

For the calculation of the coverage in this work a computer program was used (Win-Eclipse and copyright 2002). The paths of the lunar shadow for the three eclipses are depicted in Fig. 1., relevant data about the eclipses is summarized in Table 1.

Fig. 1
figure 1

Paths (red lines) of the lunar shadow for the solar eclipses in Réunion (a), Puerto Chacabuco (b), Oregon (c), and Rothera (d). Image created with WinEclipse software (Win-Eclipse and copyright 2002)

Table 1 Data on the eclipses investigated in this campaign

3 Experimental

3.1 Réunion

Experiments took place at the Observatoire Les Makes, 18, rue G. Bizet, Les Makes, 97421 La Rivière, Ile de La Réunion. An HP laptop running Windows 7 was used to control all devices and store the data. The water used in these experiments was “EDENA l’eau pure du cirque de mafate”, a normal potable water available on the island. The ion content is given in Table 2. This water was chosen over distilled or purified water due to its stability concerning ion strength and CO2 content and associated parameters such as surface tension or refractive index, thereby minimizing a measurement drift induced by chemical reactions.

Table 2 Ion content of the sample water EDENA eau de source (2016)

Temperature and humidity inside the lab were measured using a TSP01 2-channel compact USB temperature and humidity logger with an accuracy of 0.05 °C and 0.1% relative humidity, respectively (Thorlabs GmbH, Dachau, Germany). The measurement rate was one per second.

Refractive index measurements were conducted using a Rudolph J47-WR S5 refractometer in the continuous measure mode connected to the laptop via a virtual COM port. The J47 has a built-in temperature correction using the latest ICUMSA tables (ICUMSA Methods Book and Bartens 2016). The accuracy of the measurements was 0.0001 nD, the measurement rate was one per 5 s. Before the actual measurements the instrument was calibrated using Milli-Q water (one-point calibration). For the measurements, several mL of the sample water were put onto the J47 sapphire window, then the lid was closed to prevent interference from outside light sources and to minimize evaporation.

Background magnetic fields were measured using a 3-axis gaussmeter (Alphalab Vector/Magnitude Gauss Meter Model VGM). The model VGM is designed for rapid identification and measurement of areas that are magnetized. It is a DC gaussmeter with X, Y and Z sensors at the end of a probe with a range from 0.00 to 799.99 gauss in 0.01 gauss steps, noise level (jitter) of 0.02 gauss for X, Y and Z and magnitude (compare to the Earth field of about 0.5 gauss); and a 1% accuracy according to the manufacturer. All data (X, Y, Z and magnitude) were sampled at a 2 Hz rate.

Surface tension measurements were carried out using an AquaPi Plus tensiometer (Kibron Inc., Helsinki, Finland) with an accuracy of 0.1mN/m. Surface tension was measured using the system’s Wilhelmy method (Wilhelmy 1863; Adamson and Gast 1967). In this method the needle is first submerged into the sample (by moving the sample) and then slowly pulled out until the point where pulling force does not increase any further. At this position the sample is held in place and the force is measured continuously at a rate of 10 values per second. Since the software delivered with the instrument was not designed for long-term measurements software was developed using said Wilhemy method which would store all measured data directly to the hard disk. Before the measurement the system was calibrated using two points (water and air) with 72.8 and 0.0 mN/m, respectively. After initial stabilization of the system (variation of the readout < 5%) average and standard deviation values of 200 single measurements were recorded. Whenever the standard deviation of that measurement was larger than 0.05 mN/m the measurement period was extended until such a precision was reached. Because the Kibron AquaPi Plus sample chamber is not gastight after some minutes of measurement the value started to drop due to evaporation of the sample. This process led to a detachment of the needle from the surface after 6–7 h. This time period was nevertheless long enough to measure continuously during the whole solar eclipse, from first contact to fourth.

The AquaPi + tensiometer shows very stable operation when measuring for the normal duration of surface tension measurements which is a few minutes at most. When measuring over hours, a linear drift can be observed, mostly downwards with different inclinations. The same drift is visible when just the weight of the needle is recorded, allowing the conclusion that the reason for this behavior is an electric drift of the balances of the tensiometer. For that reason, all surface tension measurements have been baseline corrected before testing their correlation to the eclipse.

3.2 Puerto Chacabuco

Experiments took place at a hotel room in the Loberias del Sur hotel, Jose Miguel Carrera 50 XI region, Puerto Chacabuco, Chile, since no scientific institute was available in the totality zone. An HP laptop running Windows 7 was used to control all devices and store the data. The water used in these experiments was “PUYEHUE Patagonia mineral water”, a normal potable water available at the hotel. The ion content is given in Table 3. Temperature and humidity inside the hotel room were measured using the same equipment as in Réunion. Refractive index was not measured in Puerto Chacabuco. Surface tension measurements and DC magnetic field measurements were carried out as in Réunion. The air pressure was measured with a barometric pressure USB data logger type B1100-2, Gulf Coast Data Concepts, LLC, 1 Hz sample rate, ± 250 Pa absolute accuracy.

Table 3 Ion content of the sample water in Puerto Chacabuco (2017)

3.3 Oregon

Experiments took place in the open field in Oregon 97874, USA (exact coordinates latitude: 44.684277, longitude: − 120.033875), The open field was chosen for reasons of maximum duration and also to obtain one result without the possible interference of building roofs and walls. A Dell XPS 15 9550 laptop running Windows 10 was used to control all devices and store the data. The water used in these experiments was “Arrowhead Mountain Spring Water”, a normal potable water available at most supermarkets. The ion content is given in Table 4. Temperature and humidity at the measurement site were recorded using the same equipment as in Réunion. Refractive index was again not measured. Surface tension measurements were carried out as in Réunion. The air pressure and DC magnetic fields were not measured during this eclipse.

Table 4 Ion content of the sample water in Oregon (2017)

3.4 Rothera

Experiments took place at the Bonner labs at the Rothera research station in Antarctica. An HP laptop running Windows 7 was used to control all devices and store the data. The water used in these experiments was “cachantum ORIGEN NATURAL” mineral water, a normal Chilean potable water. The only ion content mentioned on the bottle is a sodium concentration of 20 mg/L. Temperature and humidity inside the lab were measured using the same equipment as in Réunion. Refractive index was not measured at Rothera. Surface tension measurements were carried out as in Réunion, but this time also with a second tensiometer (Aqua Pi + , Kibron, newer model with micro USB interface) in parallel. The second tensiometer, although bearing the same name, had a higher sensitivity and an increased digital resolution, but was also more unstable and had higher relative errors. The air pressure was measured with a barometric pressure USB data logger type B1100-2, Gulf Coast Data Concepts, LLC, 1 Hz sample rate, ± 250 Pa absolute accuracy. In addition, the background magnetic field (AC) was measured using an AC Milligauss Meter Model UHS2 in a frequency range of 13 Hz to 75 kHz. The meter measures the true 3-axis magnitude of the AC field. The background electric field was measured using a Trifield Natural EM Meter NEM which can detect small changes (3 V/m) in the natural DC fields and some radio transients.

3.5 Remark on Instrument Precision

After the end of the measurement campaign the authors discovered an issue with the data transfer routines for the surface tension meters and the magnetometers which had resulted in the occasional addition of artificial noise to the measurement data, visible in the large error bars of some data points, but fortunately not so frequent or large that would obscure the effects reported. The authors could correct this issue to ensure that planned future measurements on this phenomenon will be more precise than the ones presented.

4 Results

In all locations a positive cross-correlation between the eclipse and the baseline-corrected surface tension was found. These results are summarized in Table 5.

Table 5 Cross-correlation of the baseline-corrected surface tension measurements with the eclipse at the respective locations during the solar eclipse experiments

In all locations except for Oregon reference measurements were done at the same place at the same time but one day before or after the eclipse. As well as surface tension, environmental parameters (temperature, humidity and air pressure) were recorded. Depending on availability and equipment functionality not all parameters could be measured at all locations. Whenever there are no error bars visible the standard deviation is smaller than the diameter of the gray dot. The linear baseline drift was taken into account in all graphs by subtraction of said baseline. The baseline itself was determined from the initial behavior of the surface tension before the eclipse or, in case of the Rothera measurements, where the behavior itself was much more unstable than anywhere else, as a linear fit over the whole measurement. After starting the measurement the tensiometer needs a few minutes for stabilization (see e.g. steep rise Puerto Chacabuco data without eclipse, Fig. 5); these data points were not considered in the baseline calculation. The (spread of the) scales for all parameters are the same in the following graphs in order to facilitate comparison, except for Oregon where the temperature drop increased the surface tension rise almost by a factor of 2 (Figs. 2, 3, 4).

Fig. 2
figure 2

Surface tension, temperature and humidity during the annular eclipse in Réunion. The scales for humidity and surface tension are on the left, the scales for eclipse and temperature are on the right side, respectively

Fig. 3
figure 3

Surface tension, temperature and humidity one day after the annular eclipse in Réunion. The scales for humidity and surface tension are on the left, the scale temperature is on the right side, respectively

Fig. 4
figure 4

Surface tension, temperature, humidity and air pressure during the annular eclipse in Puerto Chacabuco. The scales for humidity, surface tension and air pressure are on the left, the scales for eclipse and temperature are on the right side, respectively

Fig. 5
figure 5

Surface tension, temperature, weight of an open water bottle (purple crosses) and water evaporation rate (thick yellow line) during the annular eclipse in Puerto Chacabuco. The scales for surface tension, weight and evaporation rate are on the left, the scales for eclipse and temperature are on the right side, respectively

In Puerto Chacabuco also the evaporation rate of water, the weight of a filled (and sealed) and an empty water bottle were monitored during the eclipse. Whereas the weight of both the filled and the empty bottle did not change during the eclipse, the evaporation rate did, as is shown in Fig. 5.

The evaporation rate was calculated by differentiation of the weight data every 10 min. It is noteworthy that the water evaporation rate peak, visible also as steepest increase in the water weight curve, coincides with the surface tension peak. This result is quantitatively reflected in the cross correlation analysis (see Table 6, Figs. 6, 7, 8, 9).

Table 6 Cross-correlation of the baseline-corrected surface tension and evaporation rate measurements with the eclipse in Puerto Chacabuco
Fig. 6
figure 6

Surface tension, temperature and humidity one day after the annular eclipse in Puerto Chacabuco. The scales for humidity and surface tension are on the left, the scale temperature is on the right side, respectively

Fig. 7
figure 7

Surface tension, temperature and humidity during the total eclipse in Oregon. The scales for humidity, surface tension and air pressure are on the left, the scales for eclipse and temperature are on the right side, respectively. The gray crosses show the calculated change of surface tension due to the temperature drop

Fig. 8
figure 8

Surface tension, temperature, humidity and air pressure during the partial eclipse in Rothera. The scales for humidity and surface tension are on the left, the scales for air pressure, eclipse and temperature are on the right side, respectively. The gray crosses are the data from the second, more sensitive tensiometer

Fig. 9
figure 9

Surface tension, temperature, humidity and air pressure one day before the partial eclipse in Rothera. The scales for humidity and surface tension are on the left, the scales for air pressure and temperature are on the right side, respectively. The gray crosses are the data from the second, more sensitive tensiometer

Refractive index, electric field and background radiation measurements did not show any correlation to the solar eclipse (data not shown). Additional correlations found in the data analysis are summarized in Table 7.

Table 7 Additional correlations (|P| > 0.5) of measured parameters with the eclipses

The results from all other measurements described in the experimental section did not reveal a correlation with the solar eclipses.

5 Discussion

Measuring and understanding a hitherto unknown effect can be a very difficult and delicate task. The force that causes the surface tension to change in correlation to the solar eclipse might not be constant in the absence of the eclipse, so getting good reference data is a challenge. Not knowing its nature makes it impossible to shield it, and the best option is to keep all parameters constant except for the eclipse. That, however, is also only possible within certain margins. Nevertheless, the results of the measurements presented show a general common behavior: the rise of the surface tension in correlation with the eclipse event, albeit between 20 and 85 min later.

5.1 Steps and Linear Retreat of the Surface Tension Measurement

A closer look at the surface tension data reveals some downward steps mostly accompanied by short periods of instability (see e.g. Figure 1, shortly after the maximum of the eclipse). These instabilities caused the instrument to continue measuring until a standard deviation of 0.05 was reached, since during these moments the values were fluctuating so much that a reliable readout (σ < 0.05) was not possible. Such events are visible in the data as “blank periods” followed by large, sudden jumps. They can be explained by a step-wise film retreat from the needle, a known phenomenon for droplet evaporation (McHale et al. 2005) and is due to irregularities and associated energy barriers on wetted surface. In the present case some of these irregularities are probably created during the experiment by microscale calcium carbonate crystallization on the needle or by vibrations. Especially on Réunion, an island with an active volcano, it is likely that little tremors caused the irregular behavior. For the overall increase of surface tension observed in this work, however, we conclude that these steps are irrelevant and represent an independent phenomenon unrelated to the solar eclipse since they occurred both during and after the eclipse and also during both increase and decrease of the surface tension during the eclipse.

5.2 Terrestrial Parameters Capable of Influencing the Surface Tension of Water

5.2.1 Humidity

Pérez-Díaz et al. (2012) have shown that the effective surface tension on a water droplet decreases as humidity in the air increases. Their experimental data is in agreement with their model based on Newton’s laws leading to

$$\gamma_{e} = \gamma + \frac{{a\left( T \right)\left\langle {V_{n} } \right\rangle }}{2K}\left( {1 - h} \right),$$
(5)

where γe is the effective surface tension, γ is the surface tension, a(T) is a function of the temperature T and h is the unitary humidity valued from 0 to 1. In this model, it is assumed that evaporating molecules of water are ejected from liquid to gas with a mean normal component of the speed of ‘‘ejection’’ < vn> greater than zero. As a result, depending on the temperature, the surface tension of water decreases with increasing humidity (negative correlation). At 20 °C the value ranges from ~ 75.5 mN/m at 20% humidity linearly down to 73 mN/m at 100% humidity. Thus a humidity induced change in surface tension would be visible as a negative correlation. Table 8 shows the correlation of baseline-corrected surface tension and humidity for the four locations.

Table 8 Correlation of the baseline-corrected surface tension measurements with humidity at the respective locations during the solar eclipse experiments

If humidity were (also) responsible for the observed effects, humidity drops of 20% negatively correlated to the surface tension changes should have been observed. That was not the case for any of the experiments performed. On the contrary, a moderate positive correlation was found for the measurements at Réunion, which probably weakened the observed effect.

5.2.2 Temperature

Next to humidity, from Fig. 1 of Pérez-Díaz et al. (2012) the temperature dependence of surface tension can be extracted using linear regression. In the present experiments, the humidity was different for each location, ranging between close to 50% in Réunion to about 20% in Rothera. The corresponding temperature coefficients for this humidity range are shown in Table 9.

Table 9 Temperature coefficients according to (Pérez-Díaz et al. 2012) obtained by linear regression from measurements at 5, 10, 15 and 20 °C

In principle, as with humidity, there is a negative correlation between temperature and surface tension. The correlations found in the present experiments are given in Table 10.

Table 10 Correlation of the baseline-corrected surface tension measurements with temperature at the respective locations during the solar eclipse experiments

If temperature was responsible for the observed effects, a negative correlation to the surface tension changes should have been observed. That was not the case for any of the experiments performed. On the contrary, again a moderate positive correlation was found for the measurements at Réunion, which probably weakened the observed effect.

In Oregon, where the experiments took place in the field, both temperature and humidity in that experiment changed drastically during the eclipse. There is also correlation of the eclipse with temperature (P = − 0.623). In order to estimate how much of the effect is due to the temperature change, the theoretical change of the surface tension due to the temperature drop at the measured humidity was calculated and shown as grey x symbols. The temperature coefficient for each data point was obtained from a linear regression from the data in Table 9 and the humidity data. About one-third of the surface tension rise can indeed be explained by the temperature drop.

Since humidity and surface tension are negatively correlated and the humidity rose during the experiments in Oregon, any correction for a humidity change would only slightly increase the observed surface tension rise, but cannot be the cause of it.

5.2.3 Air Pressure

Van Flandern and Yang (2003) have shown that rapid air mass movement for the bulk of the atmosphere above normal cloud levels increased the atmospheric pressure by at least 0.6% during an eclipse (from 979 to 985 g/cm2) and returned to pre-eclipse level afterwards. This change was induced by rapid air mass movement for the bulk of the atmosphere above normal cloud levels and was said to be sufficient to explain gravitational anomalies reported previously (Van Flandern and Yang 2003). This change is much bigger than any air pressure change observed in the experiments in this work. Massoudi and King (Massoudi and King 1974) have investigated the effect of air pressure on the surface tension of water. According to their measurements, the interfacial tension as function of pressure for water can be approximated by a 3rd grade polynomial. Assuming the same relation can be used as good approximation for the potable water used in this study and an atmospheric composition of 78.08% N2, 21.88% O2 and 0.04% CO2, a pressure drop of 9.8 atm would be necessary to explain the increase of surface tension observed in this work. Although air pressure was not monitored during all eclipses it is safe to say that such a drop did not occur, and that any atmospheric pressure increase as reported by Van Flandern and Yang (2003) would have caused a reduction of the surface tension in the µN/m range, which is beyond the precision of the instrument used. A more direct approach is the calculation of the volume change due to air pressure change. With a compressibility of ~ 5 × 10−10 Pa it is straightforward to calculate that the volume of water would decrease by 0.25 ppm due to the 5 mbar pressure increase. The sample beaker can be approximated by a cylinder of 1 cm height, which would be reduced by 2.5 nm due to the increase pressure. This distance is equivalent to a layer of about 10 water molecules and far beyond the precision of the instrument used in this study.

5.3 Cosmic Parameters Capable of Influencing the Surface Tension of Water

5.3.1 Gravity and Tides

The tidal acceleration at the Earth’s surface along the Moon-Earth axis is about 1.1 × 10−7 G (Ministry of Defence (NAVY) 1987). According to a NASA report (Weislogel et al. 1998) the surface tension of water is dependent on gravity, and a linear interpolation of the values given in this report allows an estimation of the effect during a solar eclipse, which is an increase of ~ 1.1 × 10−7 mN/m. This change is 6 orders of magnitude lower than the effect measured in this work.

Another possible cause would be the change of capillary rise due to said change of gravity. The height of a liquid column is given by Batchelor et al. (1967)

$$h = \frac{2\gamma cos\theta }{\rho Gr}$$
(6)

where γ is the surface tension, θ is the contact angle, ρ is the density of liquid, G is the local acceleration due to gravity, and r is the radius of tube. A change of gravity by 1.1 × 10−7 G results in a change of h by ~ 10−13 m independent from the other parameters. Or, a small change of, for example, 0.8 mN/m would require an increase of gravity by 0.11 G. The first effect is negligible, and the second is unrealistic and carries the wrong sign.

Finally, it may be mentioned that the water level in a small beaker cannot be directly compared to an ocean concerning tidal range. The huge tidal range in the oceans is possible because large bodies of water move from one location (where there is less attraction from the Moon) onto another (where there is lunar attraction)—water is not expanding or shrinking, its compressibility is extremely low (~ 5 × 10−10 Pa−1). In a beaker the circumstances are different. Displacement of water is not possible, and the volume expansion is negligible, and tilting of the surface would also be in the nm range. In addition this change will take place over a period of 12 h in a sinusoidal manner, so the change during the period of an eclipse, when it will be at the max of the sinusoidal variation will be even more negligible.

5.3.2 Solar Wind

The continuous flow of particles emanating from the Sun and permeating the solar system is often referred to as solar wind. It consists mostly of electrons, protons and helium nuclei (alpha particles). Its existence was first suggested by Arthur Eddington (2006) and Biermann (1951) suggested it to be reason for a comet’s tail pointing always away from the Sun independent from its heading. In 1959 solar wind was measured directly for the first time by the Soviet satellite Luna 1 (2007). The solar wind consists of about 1.3 × 1036 particles per second (Kallenrode 2004) resulting in a total annual solar mass loss of ~ (2–3) × 10−14 solar masses (Carroll et al. 1995) equal to about a billion kilograms per second. The velocity distribution of the solar wind shows two maxima and is therefore divided into the slow solar wind and the fast solar wind, with their differences extending beyond their speeds. In near-Earth space the slow solar wind has a velocity of 300–500 km s−1 and a composition comparable to the corona. The fast solar wind has a typical velocity of 750 km s−1 and its composition is comparable to the Sun’s photosphere (Geiss et al. 1995). The slow solar wind is twice as dense and more variable in nature than the fast solar wind (Kallenrode 2004; Suess 1999).

The temporal offset of the surface tension effect suggests a dependence on a cause that travels from Moon to Earth significantly slower than the speed of light, and on the first glance solar wind seems to be a promising candidate. The observed time differences yield velocities between 193 and 386 km s−1, which is at the lower end of the slow solar wind component. However, predicting the occurrence of a solar wind shadow precisely is not a simple task: When a syzygy of Earth, Moon and Sun results in an eclipse the solar wind is passing by the Moon, which is creating a shadow. However, the slow solar wind component will take about 16 min to reach Earth. Whereas the solar wind shadow that was being cast at the time of eclipse moves towards the Earth, the Earth is sliding sideways as it orbits the Sun. After 16 min, when the shadow reaches the position that the Earth was in at the time of the eclipse, the Earth has been travelling for 16 min on its orbit around the Sun at ~ 30 km s−1 so will have moved about 28830 km or 2.3 Earth diameters. In a first approximation, as a result the Moon shadow would completely miss the Earth. However, solar wind consists of charged particles which are strongly deflected by the Earth’s and Sun’s magnetic fields and attracted (and repelled, depending on the charge) by the Earth’s electric field. Without knowing the exact magnitude and geometry of these fields at the time of the eclipse it is impossible to predict whether a solar wind shadow will reach Earth or not.

The impact of the solar wind on Earth is measured by NASA satellites and commonly known as the so-called K-index. It quantifies disturbances in the horizontal component of Earth’s magnetic field with an integer in the range 0–9 with 1 being a calm and 5 or more indicating a magnetic storm, and is derived from the maximum fluctuations of horizontal components observed on a magnetometer during a three-hour interval. The index was introduced already in Bartels et al. (1939). Table 11 summarizes the K-indices for the time periods of the solar eclipses studied in this work.

Table 11 K-indices and surface tension rise during solar eclipse events 2015–2018. The surface tension rise in Oregon due to the temperature decrease is subtracted from the total amount

There is a seeming correlation between the K-index and the magnitude of the effect, corroborating the solar wind hypothesis, taking into account that the eclipse at Rothera was only partial. In addition, a solar eclipse in a region with high solar wind density such as Rothera should produce a stronger effect and a more strongly fluctuating base line. Again, taking into account that the eclipse in Rothera was only partial, both of these effects were found. Keeping in mind that a correlation alone can never prove a causality relation, the data at hand makes the solar wind hypothesis the most likely one of all hypotheses mentioned here. This statement, however, requires a model as how the solar wind physically changes the surface tension, which is presented in the next paragraph.

5.3.2.1 Terrestrial Muons Disrupt the H-bond Network At the Water/Air Interface

Muons are abundantly produced by the interaction of high energy cosmic rays with Earth’s upper atmosphere according to

$$CR + A \to \pi + K \to \mu + \nu$$
(7)

where CR are the cosmic rays, A is Earth’s atmosphere, π are pions, K kaons, µ muons and ν neutrinos. Kaons decay again into muons and neutrinos (Marinho et al. 2014). The cosmic ray flux has a broad energy spectrum extending from 103 up to 1020 eV (Yao et al. 2006), most of which, however, is perfectly shielded by the geomagnetic field (Abraham et al. 2010). Muons, however, are unaffected by this shielding since they are created inside Earth’s atmosphere. Therefore they are the dominant component of the radiation flux at sea level for energies above 100 MeV (Alpen 1998). Interestingly, it has been shown that the muon flux on Earth’s surface is dependent on the K-index (Brancus et al. 2009), pointing at a correlation between the intensities of the solar wind and the cosmic radiation. Without going into detail about how exactly these two parameters are related, this finding (Brancus et al. 2009) allows to conclude a correlation between the muon density on Earth and the solar wind intensity measured by the K-index. Furthermore, it has been shown that when other cosmic rays shower in the Earth’s atmosphere and create muons, there is a relative deficit of muons from the direction of the Moon (Boersma et al. 2009). In this hypothesis we suggest a similar “Moon shadow” event in conjunction with the solar wind and its relation to cosmic rays (Brancus et al. 2009).

The main reaction caused by muon interactions is ionization, giving rise to several electrons and low-energy photons that will propagate and interact further (Marinho et al. 2014). The energy necessary for ionizing a water molecule, along their paths through water is about 13 eV, thus many orders of magnitudes lower than the typical energies deposited by a muon, which are in the keV range. Muons leave a trace of secondary electrons, ionized water molecules, radicals and the products thereof along their penetration path, thereby creating additional charge carriers in the liquid. It has recently been shown that water with positive excess charge does exist (Fuchs et al. 2016) and that this charge will accumulate at the surface, lowering the surface tension (Fuchs et al. 2019). Moreover, it has been shown (Chen et al. 2016) that the ionisation of a single water molecule can result in collective effects which can be linked to macroscopically detectable distortions like the surface tension. These findings support the hypothesis that the absence of muons is the source for the surface tension rise shortly after a solar eclipse. If the solar wind and the cosmic radiation are blocked by the Moon during a solar eclipse event, the muon production in the upper atmosphere should be reduced as well. Fewer muons then impact on the water/air interface, which leads to fewer disruptions and a higher surface tension—and a lower evaporation rate as observed during the Puerto Chacabuco eclipse.

Considering that the average muon flux is about 10,000 muons/m2 per minute (Grieder 2001), it is straightforward to calculate the number of molecules in the beaker, approximated as a perfect cylinder with 1 cm radius and height, with which these muons interact assuming total muon penetration and a density of 1 kg/dm3), which is 1.1 × 108 molecules/minute. We could now assume that time lag of the surface tension rise is partly due to velocity of the solar wind and partly due to the dynamics of the muons in solution. The average solar wind velocity close to Earth is roughly 500 km/s, which could explain about 13 min of the delay. If we average the time lags of the sites with similar K-index (all except Puerto Chacabuco), the average lag time is (31 ± 4)min (standard error), which would mean that it takes 18 min for muon influence to wear of, or to start again. This time constant is comparable to the time constant for discharging protonically charged water (22 min) (Fuchs et al. 2019). That effect would increase the number of ionized molecules causing the effect to 2 × 109, which is equivalent to a concentration in the order of 1 pico mol/L. Surface tension changes due to free charges have been found in the nano mol range (Fuchs et al. 2016). Thus, in order to explain the surface tension rise with this effect, we have to assume that the impact of a muon on a water molecule will influence (ionize) at least a thousand water molecules in the neighbourhood, and that the charged species created by this process will not immediately recombine, but survive in the liquid for a matter of minutes. The latter has been shown recently (Fuchs et al. 2016, 2019); and given the energy of a muon is in the GeV range, a mass ionisation event seems possible. Theories that support the idea of the coherent behaviour of thousands of water molecules that would allow such processes do exist (Del Giudice et al. 1985), and the authors are planning to devise a fully developed quantum field theory explanation describing in detail the interaction between a muon and a coherence domain. In addition, the authors plan to further test this hypothesis by measuring the muon intensity directly during an eclipse event.

For the sake of completeness, additional less likely (but nevertheless possible) hypotheses explaining the result are presented in the appendix.

5.4 Summary

We have shown that the blocking of the solar wind and cosmic radiation with the subsequent absence of atmospheric muons during an eclipse event is a plausible working hypothesis for the observed changes in surface tension and evaporation rate during solar eclipses. Alternative hypotheses including a a novel neutrino/water interaction are presented in the appendix. None of the obvious environmental parameters (temperature, humidity and pressure) can explain the effects observed.

It is noteworthy that many solar eclipse effects reported in literature, such as the hypothesized change in the ozone layer (Bezverkhny et al. 1956; Fournier 1956; Kawabata 1936; Khrigian et al. 1961; Jerlov et al. 1954; Pariisky and Chen-Chao 1958; Steblova 1961; Stranz 1961; Svensson 1968; Chakrabarty et al. 1997; Zerefos et al. 1999), or the pH of the sea (Santhosh 2014), happened at approximately the same time offset as the effects in our study, some tens of minutes after the eclipse. An influence much slower than the speed of light seems plausible at first glance, pointing at solar wind and cosmic radiation blocking as possible cause for these effects as well. Although the involved magnetic and electric fields make an actual simulation of the effect very difficult, it remains to be the most likely hypothesis.

6 Conclusions

Studies of surface tension of water during four solar eclipse events have resulted in an increase between 0.6 and 2.4 mN/m between 20 and 85 min after the maximum eclipse, respectively, independent of environmental parameters like temperature, air pressure and humidity. The magnitude and times can be correlated to the K-index which reflects the impact of solar activity on Earth. During one of the eclipses water evaporation experiments were conducted and showed a temperature independent reduction of water evaporation at the same time as the surface tension increase was observed. A disruption of the water/air interface by muons which is reduced during a solar eclipse due to blocking of both cosmic radiation and solar wind is presented as working hypothesis. As alternative hypothesis, we propose a novel neutrino/water interaction based on chiral preferences of water interactions. More studies and theoretical approaches are called for to further investigate and explain the phenomenon observed.