Abstract
In this paper, for the dust environment in complicated atmospheric on Mars, the Martian atmosphere is vertically stratified based on the vertical distribution characteristic of Martian dust particles. Combining the number concentration and particle size distribution of each layer of dust, the transmission characteristics of polarized laser in the Martian dust are studied by using the polarized Monte Carlo method. The results show that among the six selected wavelengths, the transmittance of the wavelength of 7.46 μm in Martian dust is the largest, while the laser intensity attenuation of the wavelength of 0.49 μm in Martian dust is the largest. The horizontally polarized laser of the wavelength of 0.55 μm (effective variance b = 0.1 μm) is taken as an example; the depolarization degree of the polarized laser under stratified conditions increases from the minimum value of 3.5440157e − 4 at 0 km to the maximum value of 0.0092563 at 24 km, with an increase of 0.00890189843; while the depolarization degree under nonstratified conditions increases from the minimum value of 1.2678582e − 4 at 0 km to the maximum value of 0.0045636 at 24 km, with an increase of only 0.004436842, the degree of depolarization under the condition of stratified is greater than that of nonstratified. On the whole, for different effective variance b, the variation trend of the polarization parameter with height for different polarized laser is basically identical.
1. Introduction
As the second-closest terrestrial planet to the Earth, Mars has always been the focus of deep-space exploration in the world. Mars has a complicated atmosphere and lacks water and vegetation on its surface. Almost every Martian year (equivalent to 687 days on Earth), there will be an extraordinary dust storm sweeping the globe on Mars, with a high speed of 180 m/s and a long duration. The absorption and scattering effect of suspended dust particles in the Martian dust storm will cause laser attenuation, which will seriously affect the stability and reliability of wireless optical communication on Mars. In the near-surface atmosphere on Mars, the concentration and particle size of dust particles are unevenly distributed in vertical height, so it is necessary to perform stratification when calculating the transmission characteristics of laser. In this research field, Wu et al. [1] in 2004 used the four-flux method and Monte Carlo method to study the multiple scattering and transmission attenuation characteristics of laser in a stratified dust atmosphere, which provided a theoretical basis for the numerical calculation of laser transmission in a dust environment. Compared with the traditional laser, the polarized laser is widely used in the field of laser communication and detection due to its unique characteristics [2, 3]. In terms of polarized laser research, Plass and Kattawar and Plass et al. [4, 5] from 1972 to 1973 used the Monte Carlo method to calculate and analyze the variation in polarization degree and polarization direction after multiple scattering of polarized laser in the clouds, which laid a foundation for the later detection of aerosol parameters in the atmosphere. In 2012, Yang [6] studied the depolarization effect of polarized laser in atmospheric transmission, which furnished a basis for the application of circular polarization coding in communication systems. In 2015, Fernández et al. [7] used polarized Monte Carlo to analyze the full polarization characteristic of the Rayleigh and Compton effects in multiple scattering. In 2019, Ma et al. [8] studied the transmission characteristics of the signal laser in free-space optical communication systems, which provided theoretical support for obtaining higher mixing efficiency by using a polarization control algorithm. In 2020, Yang et al. [9] used the polarized Monte Carlo method to simulate the multiple scattering of polarized laser in the Martian atmosphere, which offered theoretical support for analyzing the effect of wavelength, particle size distribution, wind speed, height, and particle number concentration on the transmission characteristics of the polarized laser. In 2021, Alemanno et al. [10] used spectral polarization technology to measure the optical constant of Mars simulated matter, which supplied a method for the planetary science community to derive the optical constant of mineral and material (natural or synthetic).
Based on the Martian atmospheric environment, this paper uses the polarized Monte Carlo method to stratify the atmosphere of 0–24 km on Mars and calculates the transmission characteristics of polarized laser at the vertical height. Then the transmittance and depolarization degree of polarized laser transmitted in the Martian dust environment are analyzed to provide theoretical support for the realization of reliable and stable wireless optical communication on Mars in the future.
2. Laser Transmission Theory
According to the complex refractive index of Martian dust particles [11, 12], as shown in Table 1, Mie theory [13] is used to calculate the extinction efficiency factor , scattering efficiency factor , and absorption efficiency factor of Martian dust particles as follows:
Here, and represent the scattering coefficient, and the spectral scattering characteristic of Martian dust particles can be calculated by equation (1).
Based on the spectral scattering characteristic of Martian dust particles, combined with the laser transmission theory, the multiple scattering of polarized laser in Martian dust environment is imitated by the polarized Monte Carlo simulation method [14]. The specific steps are as follows [9]:(1)Set the initial state of the photon, including position coordinate, movement direction, and so on.(2)Movement of the photon. The photon moves from a point to the next collision point . The relationship between the two points is as follows: Here, is the cosine of the direction angle; is the moving step length, ; is a random number evenly distributed between 0 and 1; and is the average attenuation coefficient of the particles.(3)Selection of scattering direction. Suppose the Stokes vector of the incident wave is , then its phase function is as follows: Here, and are the two elements of the scattering matrix , is the scattering angle, and is the rotation angle between the incident meridian surface and the scattering surface. Then use the “acceptance-rejection sampling” method to sample the random variable , , and in equation (2), and the corresponding interval is , , and . If , then the selection is successful, and go to the next step; otherwise, select again.(4)Photon scattering tracking and polarization processing. After the scattering angle and azimuth are successfully selected, the Stokes parameter will be converted in three steps:(i)From the incident meridian surface to the scattering surface, the incident Stokes parameter is multiplied by the rotation matrix to obtain the rotated Stokes parameter ; the calculation of is as follows:(ii)Transformation on the scattering surface. The Stokes parameter is transferred from one point to another on the scattering surface, and the conversion matrix is Mueller matrix . At the same time, the conversion of the direction angle cosine is performed, and the new direction cosine is .(iii)Return from the scattering surface to the new meridian surface to reach the final goal. The rotation matrix at this time is , which is calculated as follows: If , equation (5) takes the plus sign. If , equation (5) takes the minus sign. The final Stokes parameter is shown in the following equation [15]: Here, represents the Stokes parameter before scattering.(5)Termination condition and final state. The tracking of the photon is terminated when the photon escapes the interface or the weight is less than the set value. The final state of the total polarized laser is as follows:
Here, is the number of photons. The quantitative value describing the degree of polarization of polarized laser is expressed by the degree of polarization as follows [16]:
Here, represents the total laser intensity, represents the difference value of laser intensity in and directions, represents the completely linearly polarized laser component with an angle of to the scattering plane, represents the right- or left-hand polarized state component of the circularly polarized laser, and represents the proportion of the laser intensity of completely linearly polarized laser in the total laser intensity. Use the degree of depolarization to describe the variation in the degree of polarization of linearly polarized laser, as follows:
Here, represents the laser intensity in the horizontal component and represents the laser intensity in the vertical component.
Because the concentration and particle size of dust particles in the near-surface atmosphere on Mars are unevenly distributed in the vertical height, it is necessary to stratify the lower atmosphere of Mars when calculating the transmission characteristics of the laser. For the polarized Monte Carlo processing of stratified Martian dust atmosphere, it is necessary to consider not only the reflection and transmission within each stratified surface but also the reflection and transmission relationship between adjacent layers. In the case of multiple scattering, it is assumed that photon moves from i layer to ith + 1 layer. At this time, it needs to be handled according to the following steps [1]: if the photon is transmitted to the ith + 1 layer, it is moved first to the boundary of the ith + 1 layer along the scattering direction; the transmission probability and weight function are calculated; and its transmission direction is decided, and then the photon travels the remaining distance along the new scattering direction. If the photon is not transmitted out and remains in the ith layer, the parameter of this layer is used to continually track the position, scattering direction and transmission path of the next scattering; If the photon is transmitted to another layer, the above steps are repeated until the photon is absorbed in one layer or transmitted out of the dust.
2.1. Vertical Distribution of Martian Dust Particles
According to the launched Mars probes Viking 1 and Viking 2, the Martian atmosphere below 300 km is divided into troposphere, mesosphere, thermosphere, and exosphere from bottom to top. Because the concentration of Martian dust particles in the Martian atmosphere above 24 km is very small and almost negligible, this paper mainly studies that the Martian lower atmospheric communication environment under the height of 24 km has an influence on the transmission characteristics of the laser.
The topography of Mars has different aerodynamics in different regions. Dust storms of different sizes will be generated when the wind speed is different. In 2001, the Mars Global Surveyor probe revealed the entire process of the global dust storm on Mars. The results showed that the dust storm had already bypassed Mars in just 15 days, and its expansion speed was 16 m/s. Chassefière et al. [12] studied the influence of height variation on the distribution of dust particles in the Martian atmosphere and found that different size distribution, effective radius, and variance have different particle numbers and concentrations at different heights. The distribution of dust particles in three different models is measured, and the relationship of its distribution with height, effective radius, and variance is shown in Table 2 [17, 18].
2.2. Attenuation Characteristic in the Vertical Direction of the Martian Atmosphere
If the number concentration of Martian dust particles is taken as a constant, the calculated transmission attenuation error is too large because the effective radius and particle number concentration of Martian dust particles vary with height [19]. When there is no dust storm, the lower atmosphere of Mars is divided into 24 layers according to Table 1. The 1st layer has the lowest height and the biggest particle concentration; the particle number concentration is 1.31 n/cm3, 2.26 n/cm3, and 5.56 n/cm3 (the effective variance is b = 0.10 μm, b = 0.25 μm, and b = 0.40 μm, respectively); the 24th layer has the highest height, but its particle number concentration is the lowest; the particle number concentration is 0.12 n/cm3, 0.33 n/cm3, and 1.20 n/cm3 (the effective variance is b = 0.10 μm, b = 0.25 μm, and b = 0.40 μm, respectively). The polarized Monte Carlo method is used to calculate and analyze the laser transmission under two conditions (stratified and nonstratified). The particle number concentration is fixed to 0.39 n/cm3 0.78 n/cm3, and 2.16 n/cm3 (the effective variance is b = 0.10 μm, b = 0.25 μm, and b = 0.40 μm, respectively) under the condition of nonstratified, and the effective radius is taken, respectively, as 1.6 μm, 1.64 μm, and 1.69 μm. The calculated results are shown in Table 2, which is the transmission transmittance corresponding to several typical wavelengths under the condition of stratified and nonstratified.
It can be seen from Table 3 that when the wavelength is 7.46 μm, the calculated transmittance is larger than that of other wavelengths because the real part of the complex refractive index of this wavelength is the smallest among several typical wavelengths, and its transmittance reaches the maximum value of 9.7939100e − 001 when the effective variance b is 0.1 μm. On the contrary, the real part of the complex refractive index of the wavelength of 10.6 μm is larger than that of other wavelengths, so the calculated transmittance at this wavelength is the smallest, and its transmittance reaches the minimum value of 8.3379000e − 002 when the effective variance b is 0.4 μm. The transmittance of different wavelengths is an obvious difference, and the calculated result is also different between the condition of stratified and nonstratified, but the calculated result under the condition of stratified is more accurate than that of nonstratified, which is more in line with the actual dust environment on Mars.
3. Transmission Characteristics of Different Polarized Laser in Martian Dust Environment
Then the transmission characteristics of the horizontal and vertical components of different polarized laser are studied in the Martian dust. The distributed parameter is reff = 2.34 μm and , and the particle number concentration is 1.8 × 109/m3 at the height of 18 km. The calculation results are shown in Figures 1–5.

(a)

(b)

(a)

(b)

(a)

(b)

(a)

(b)

(a)

(b)
Figures 1(a) and 1(b), respectively, show the variation curves of the horizontal and vertical components with the transmission distance when horizontally polarized laser [1, 1, 0, 0] transmitted in the Martian dust. Figures 2(a) and 2(b), respectively, show that the variation curves of the horizontal and vertical components with the transmission distance when vertically polarized laser [1, −1, 0, 0] transmitted in the Martian dust; it almost symmetrically coincides with the variation curves of horizontally polarized light [1, 1, 0, 0] in the Martian dust. As can be seen from Figure 2(b), with the increase of transmission distance, the laser intensity in the vertical component increases first and reaches a peak at about 1 km and then decreases continually. In summary, when horizontally and vertically polarized laser transmitted in the Martian dust, as the transmission distance increases, the vertical and horizontal laser intensity will increase first to a peak and then decrease. The intensity reduction in the horizontal component and the intensity increase in the vertical component at the wavelength of 7.46 μm are significantly smaller than those of other wavelengths [9].
Figures 3(a) and 3(b), respectively, show the variation curves of the horizontal and vertical components with the transmission distance when completely linearly polarized laser [1, 0, 1, 0] (with an angle of 45° to the scattering plane) transmitted in the Martian dust. Figures 4(a) and 4(b), respectively, show the variation curves of the horizontal and vertical components with the transmission distance when completely linearly polarized laser [1, 0, −1, 0] (with an angle of −45° to the scattering plane.) transmitted in the Martian dust.
It can be seen from Figure 3 that both horizontal and vertical components of completely linearly polarized laser [1, 0, 1, 0] (with an angle of 45° to the scattering plane) gradually decrease with the increase of the transmission distance. Among the six typical wavelengths, the laser intensity attenuation in the horizontal and vertical components at the wavelength of 7.46 μm is the lowest. From 0.4 km to 4 km, the laser intensity of the horizontal component drops from 0.1081 to 0.0576, and its laser intensity attenuation is only 0.0505. On the contrary, the laser intensity in the horizontal component drops from 0.33157 to 0.04651, and the laser intensity attenuation of the wavelength of 0.49 μm is 0.2851, which also shows that polarized laser at this wavelength is most unsuitable for transmission in Martian dust environment. As can be seen from Figure 4, it is basically consistent with the variation trend of Figure 3 because completely linearly polarized laser [1, 0, 1,0] (with an angle of 45° to the scattering plane) is symmetrical with completely linearly polarized laser [1, 0, −1, 0] (with an angle of −45° to the scattering plane), so their transmission characteristics are symmetrically coincident.
Figures 5(a) and 5(b), respectively, show the variation curves of the horizontal and vertical components with the transmission distance when circularly polarized laser [1, 0, 0, 1] (in a completely right-hand circularly polarized state) transmitted in the Martian dust. Figures 6(a) and 6(b), respectively, show the variation curves of the horizontal and vertical components with the transmission distance when circularly polarized laser [1, 0, 0, −1] (in a completely left-hand circularly polarized state) transmitted in the Martian dust. Because the completely right-hand circularly polarized laser [1, 0, 0, 1] is symmetrical with the completely left-hand circularly polarized laser [1, 0, 0, −1], their transmission characteristics are also symmetrically coincident.

(a)

(b)
Considering Figures 1–6, when horizontally polarized laser and vertically polarized laser transmitted in the Martian dust, the laser intensity variation in the horizontal and vertical components is different from that of other polarized laser; they will increase first to a peak and then decrease with the increase of transmission distance. But, when circularly polarized laser and completely linearly polarized laser transmitted in the Martian dust, their laser intensity in the horizontal and vertical components will gradually decrease with the increase of transmission distance. This is one of the reasons why horizontally and vertically polarized laser is often used to transmit in the Martian dust. In addition, among the selected wavelengths, no matter that polarized laser transmitted in the Martian dust, the laser intensity attenuation in the horizontal and vertical components of the wavelength of 7.46 μm is the smallest.
4. Transmission Characteristics of Polarized Laser in the Vertical Direction of the Martian Atmosphere
Based on the content of Section 3, the transmission characteristics of polarized laser are studied in the vertical direction of Mars. Firstly, the horizontally polarized laser [1, 1, 0, 0] at the wavelength of 0.55 μm is taken as an example; the value of different polarization parameters is calculated under the condition of stratified and nonstratified; and the calculation results are shown in Table 4. In the comparison table of polarization parameters.
It can be seen from the calculation results in Table 4 that the reflection parameter IR1 and degree of polarization under the condition of stratified are significantly less than those of nonstratified and the transmission parameter IT1 and degree of depolarization dep1 under the condition of stratified are greater than those of nonstratified. It shows that the depolarization phenomenon of polarized laser is more obvious in the stratified Martian dust environment.
Secondly, the variation of polarization parameter with height is calculated and analyzed, and the horizontally polarized laser at the wavelength of 0.55 μm is taken as an example and compared the calculation results in two cases (stratified and nonstratified). The calculation results are shown in Figures 7–9.

(a)

(b)

(a)

(b)

(a)

(b)
Figures 7(a), 8(a), and 9(a) are the variation curves of the transmission parameters with the height of the Martian atmosphere. In these figures, IT1 and QT1 represent the transmission parameters under the condition of stratified, and IT and QT represent the transmission parameters under the condition of nonstratified. It can be seen from Figure 7(a) that the transmission parameter is monotonically decreasing with the increase of height; it also can be seen that the degree of depolarization is monotonically increasing as the increase of height in Figures 7(b) and 9(b), and the degree of depolarization under the condition of stratified is obviously greater than that of nonstratified. On the whole, for different effective variance b, the variation trend of the polarization parameters with height for different lase is basically identical. And the transmittance and depolarization degree of polarized laser under the condition of stratified when it transmitted in the Martian dust on Mars are obviously greater than those of nonstratified, which indicated that the polarization parameters calculated under the condition of stratified more accurately reflect the transmission characteristics of polarized laser in the Martian dust environment.
5. Conclusion
Based on the vertical distribution characteristic of Martian dust particles, the polarized Monte Carlo method combined with the number concentration and particle size distribution of each layer of dust was used to study the transmission characteristics of polarized laser in the vertical direction of the Martian atmosphere. The transmission transmittance of several typical polarized laser under the condition of stratified and nonstratified are calculated, and the variation regulation of horizontal and vertical components of horizontally and vertically polarized laser with transmission distance is given when they are transmitted in Martian dust. When the wavelength is 0.55 μm, for different effective variance b, the depolarization degree under two conditions (stratified and nonstratified) enlarges with the increase of height identically. The depolarization degree of polarized laser under the condition of stratified is significantly greater than that of nonstratified, so it shows that the depolarization phenomenon of polarized laser is more obvious in stratified Martian dust environment. In addition, the real part of the complex refractive index of the wavelength of 7.46 μm is the smallest among the six wavelengths, so its transmittance is the biggest when it is transmitted in the Martian dust, while the laser intensity attenuation in the horizontal and vertical components of the wavelength of 0.49 μm is the largest, so the trend of intensity attenuation in the vertical and horizontal component is the fastest when it transmitted in the Martian dust.
Data Availability
Previously reported theoretical and measurement data were used to support this study and are available at [DOI: 10.1364/opex.13.004420 and DOI: 10.1029/94JE03363]. These prior studies are cited at relevant places within the text as references [1–19].
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this article.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (61405157) and the Key Industry Innovation Chain Project of Shaanxi Province (2017ZDCXL-GY-06-01).