1. Introduction
The atmosphere is an important environmental element for human beings, and it is an indispensable natural resource. The Earth’s radiation balance is determined by the combination of aerosols, clouds, atmospheric gases, and surface reflections in the atmosphere [
1,
2]. When the solar radiation passes through the atmosphere, it is scattered by the atmospheric particles, which will cause the polarization of light. In general, the primary scattering of atmospheric particles causes the sky polarization to produce a positive value, while the multiple scattering causes the sky polarization to produce a negative value. The intersection of positive and negative polarizations in the sky is the atmospheric neutral point (Neutral Point), which is the point of zero polarization in the sky [
3]. Due to the scattering absorption of incident sunlight by air molecules and aerosol particles, the sky has a relatively stable polarization mode at some points in the day and at a certain position, which is the sky polarization field [
4,
5]. In 1809, for the first time, Arago discovered the polarization phenomenon of skylight and found that there is a point in the sky where the degree of polarization is zero, that is, the atmospheric neutral point. In 1870, Strutt proposed the Rayleigh scattering theory, which scientifically explained the polarization phenomenon of sky-scattered light and can more accurately describe the polarization state distribution of scattered light in clear sky.
The human eye cannot directly perceive the polarization information of light, but we can create some fast and high-precision polarization measuring instruments [
6], which is one of the hot topics for innovation. Since polarization is one of the basic properties of light, which contains characteristic information of the object being tested, polarization detection is a basic measurement method that cannot be ignored in the field of optical measurement. G. Horváth et al. proved that the skylight polarization field under cloudy, foggy, and dusty weather is similar to that in sunny weather [
7,
8]. Nathan J. Pust and Joseph A. Shaw [
9] proved that the degree of polarization of cloudy weather is significantly lower than that of sunny weather due to the influence of multiple scattering, but the distribution of polarization in azimuth is almost unchanged. A fully automatic imaging, whole-sky, polarized-light-testing instrument is designed and manufactured. The actual distribution of sky polarization mode in cloudy weather is simulated, and the simulation results are consistent with the test results [
7]. With the continuous development of polarization instruments, polarization measurement is playing an increasingly important role in many fields [
10,
11].
When solar radiation enters the Earth–atmosphere system, it is affected by the reflection, scattering, and transmission of atmospheric particles, underlying surface, etc., which cause the polarization of the incoming neutral sunlight [
12]. Relevant research found that polarized light is broadly utilized. For example, creatures such as bats use visual polarized light for accurate environmental sensing and navigation [
13], voyagers use polarized light to identify directions in the Atlantic Ocean [
14], weakly polarized reflected light of the moon can reflect vegetation traces [
15], surface-polarized reflections can improve the accuracy of biochemical content inversion [
16], and ground-based polarization observations can obtain effective information on solar flares [
17]. These cases raise fundamental questions. Is there a skylight polarization field in nature? Can this field penetrate the atmosphere and support surface observations? Is it possible to accurately monitor the interaction between the atmosphere, ocean transpiration, and the surface ecosystem? Can this field be combined with the magnetic field and the gravity field to reach a better understanding of the natural phenomena and laws of the Earth? Solving these problems could lead people to a deeper and wider insight into what skylight polarization can do. For example, the skylight polarization field might be helpful for navigation under complicated natural conditions where signals from global positioning satellites are missing. Moreover, applied to atmospheric or public health-related issues, the polarization field is meaningful to the separation of different kinds of atmospheric particles, or even the μm-size virus existing in the aerosol particles, which polarizes light differently from other particles.
In this study, we explored the laws, characteristics, and global distribution of the skylight polarization field map and its seasonal fluctuation, we explored the stability and continuous diffusivity of the light intensity distributions through coupling land surface and atmosphere polarization effects, we studied the daily periodicity of the skylight polarization field through all-day measurements with polarization instruments, and we also compared the characteristic parameters of the skylight polarization field with the existing gravity field and geomagnetic field. Furthermore, this work provides the possibility to study the three basic vector fields under the influence of the three elements of the vast universe (universal gravitation, Earth rotation, and solar radiation) on the Earth, as well as the scientific understanding of the influence of the three vector fields on all the objects on the Earth.
2. Materials and Methods
The experiment was carried out at the Huairou Solar Observing Station National Astronomical Observatories CAS, Huairou District, Beijing, China. The observing instruments are a Solar Observing Telescope and Nikon D200 digital camera with a fisheye lens and an optical polarizer (
Figure 1). We used the digital compass and leveling instruments to find the north direction and adjust the tripod platform to the horizontal position. We then installed the digital camera on the tripod, installed the fisheye lens on the digital camera, and installed the optical polarizer on the top of the fisheye lens. By using the scale around the polarizer, the rotation angle of the polarizer can be obtained, which can be marked as 0 degrees, 60 degrees, and 120 degrees. For each acquisition, the polarizer was rotated at three different degrees, respectively, and three all-sky images with different polarizing angles at that time were obtained. Every 10 min or 30 min, we collected experimental data of the whole sky and observed the distribution of sky-polarized light.
Rayleigh scattering and Mie scattering occur in the atmosphere under clear-sky conditions, and the two are also the origin of the atmospheric polarization effects. Thus, they become the theoretical and application basis of the skylight polarization field. The degree of linear polarization (DOLP) is the most important parameter for characterizing polarized light and establishing a polarization model. The DOLP of atmosphere during Rayleigh scattering can be directly simulated by the scattering angle (the angle between the directions of incidence for sunlight and scattering).
In measuring the atmospheric polarization, a Stokes vector, S, is usually used to describe the polarization state of the polarized beam. The Stokes vector,
, can be expressed as another form,
, where
I is the total intensity of light,
Q and
U are linearly polarized light in two orthogonal directions, and
V is circularly polarized light. In atmospheric polarization measurement, circular polarized light,
V, is usually ignored because linear polarized light is the most common type of polarization in nature. The Stokes vector represents the polarization state of the light, while the Mueller matrix represents the process of the polarization device changing the Stokes vector of the incident light. If the Stokes vector of the incident light is
S and the Mueller matrix of the linear polarization device is
T, then the Stokes vector of the outgoing light,
, can be obtained by linear transformation,
:
The Mueller matrix of the ideal optical detection system is as follows:
where
α is the angle between the preferred transmission plane and the reference plane of the linear polarizer.
In the new Stokes vector,
S′, the first row is used to represent the intensity of the outgoing light passing through the optical system. If
I′ is expressed as a function of
α, the intensity of the outgoing light is as follows:
According to Formulas (2) and (3), if the light intensity values at three different
α positions are known, then the parameters of the Stokes vector, DOLP, and AOLP can be calculated.
In atmosphere polarized observation, the direction of space-borne observations is from top to bottom, whereas that of ground observations is from bottom to top. However, the DOLP, at a fixed location, ideally does not vary with the observing direction. At a specific moment, the solar altitude and azimuth angles are different at different positions on the Earth. The solar altitude angle is complementary to the zenith angle. In a horizontal coordinate system, the solar altitude angle is the angle between the solar incidence and the normal placement of the observation point, which can be calculated. From this, the solar altitude angles for various positions on Earth can be obtained at the same time so as to obtain the DOLP globally at the nadir directions.
Solar radiation is non-polarized before it enters the atmosphere. In the process of sunlight incidence on the Earth surface, through the scattering interaction with the atmospheric components (gas, aerosol particles, water droplets, and ice crystals), the skylight forms a sky polarization field with the sun as the center. The whole skylight polarization field consists of two elements: the field axis and the force line. The point with the lowest degree of polarization formed by the solar incidence point is the polarization neutral point, which is defined as the axis; the point with this as the center forms a trend of diffusion to the surrounding space (a concentric circle), which is defined as the force line. Thus, the skylight polarization field is formed. Under the condition of a non-equilibrium atmosphere, the force line of the skylight polarization field is distorted, which can be used to measure the distribution of various components in the atmosphere and the properties of the underlying surface. Using the Rayleigh model to represent the scattering of molecules in the atmosphere, we can describe and simulate the state of polarized light and the distribution characteristics of the global skylight polarization field.
The degree of polarization of the atmosphere in Rayleigh scattering can be simulated by the scattering angle
θ (the angle between the direction of incident light and the direction of scattered light) [
18]:
When a polarization observation is carried out, the satellite observation mode is from top to bottom, and the scattering angle , where is the solar zenith angle. When the ground observation mode is from bottom to top, the scattering angle . For a vertical observation, the degree of the polarization is the same whether from top to bottom or from bottom to top.
At the same time, the solar altitude angle and solar azimuth angle are different at all points in the world. The relationship between the solar altitude angle and the zenith angle satisfies
. In the horizontal coordinate system, the sun height angle
means the angle between the direct sunlight and the plane where the observation point is located:
where
ϕ is the latitude of the observer’s location,
ω is the hour angle, and
δ is the declination angle.
The time angle
ω is the angular distance from the equator of the celestial meridional circle of the observation point to the time circle of the sun. It can be calculated as follows:
where
T is the true solar time, and the declination
δ is the angle between the Earth’s equator and the connecting line between the sun and the Earth’s center:
where
n is the day of year.
The true solar time can be calculated as:
where
is the longitude of the place to be measured, “+” is applicable to the eastern hemisphere, “−” is applicable to the western hemisphere, and
is the longitude of the standard time zone.
The corrected time difference is:
where
.
Therefore, the solar altitude angle of all points in the world at the same time can be obtained, and the degree of polarization can be obtained:
The Rayleigh scattering occurs when the diameter of particles in the atmosphere is much smaller than the wavelength. For visible light, the Rayleigh scattering phenomenon is very significant. A horizontal coordinate system is established to simulate the distribution of the DOLP for the upper half of the sky. The simulated result of the DOLP is shown in
Figure 2a for Rayleigh scattering. Mie scattering occurs when the diameter of particles in the atmosphere is comparable to the wavelength of radiation. The simulated result of the DOLP is shown in
Figure 2b for Mie scattering. The atmospheric model we used is a second simulation of a satellite signal in the solar spectrum-vector (6SV) [
19]. The type of aerosol we used is ocean-type aerosol, and the AOD (Aerosol optical Depth) is 0.5. The polarization effect of the real atmosphere is caused by the combined effects of Rayleigh scattering, Mie scattering, and multiple scattering.
Figure 2c shows the measured DOLP under clear weather conditions. The white line is the meridian of the sun. The white solid circle is the polarization neutral point.
5. Conclusions
From the image of the polarization angle in a clear sky, it can be observed that the concentration point of the polarization angle is the position of the atmospheric polarization’s neutral point. The position and shape of the polarization angle’s distribution are different with different solar height angles. It can be found that with the change in the position of the sun, the convergence point of the polarization angle of the whole sky revolves around the zenith. When the solar altitude angle is high, we can only observe one convergence point of the sky polarization angle, and we can only see one convergence point for most of the day. Only when the solar altitude angle is low, such as the image at about 16:00 p.m., can we observe another convergence point.
In cloudy weather, the clouds cut off the complete polarization ring; however, the change trend of polarization can still be seen, and the neutral point area of the atmosphere in the sky is still obvious. On sunny days, the value of sky polarization is large, and the circular distribution is relatively complete; in cloudy conditions, the value of sky polarization is small, and the circular distribution is irregular.
The experimental results of the two weather conditions are in accordance with the theoretical law, and the experimental data are available. However, clouds interfere with the sky’s polarization pattern. In order to ensure the accuracy of navigation, the sky polarization mode should be measured in sunny and cloudless conditions.
- (1)
A simulation of the global skylight polarization field distribution was achieved based on the theory of the electromagnetic wave and atmospheric scattering model and global multi-point verification with the sky polarization observation sequence;
- (2)
Through experimental research on the effects of surface polarization and ground–air coupling polarization, it is demonstrated that the polarization field is centered on the incidence of the sun (that is, the axis of the polarization field, which is called the neutral region of atmospheric polarization), and the degree of polarization increases when it diffuses outward. This proves the coverage of the field to the whole sky and the continuous diffusivity of the field’s axis and force lines. The instrument constructed by polarization bionic navigation proves that the intensity of this diffusive polarization force line is stable, objective, and can be accurately measured. Its cycle time is a single day;
- (3)
At any point on the Earth, the skylight polarization field has specific value and direction, which can be measured in various forms. The physical quantities describing the magnitude and direction of polarization field are the degree and angle of polarization. The polarization field is a function of time and space. The polarization field is affected by the solar altitude angle in a day. The skylight polarization field in the same place changes over a short period (a single day) and a long period (over a year).