{'fraction_non_alphanumeric': 0.033549277588548275, 'fraction_numerical': 0.02359803201317928, 'mean_word_length': 4.779722079258878, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Zodiacal light observations and its link with cosmic dust: a review

The zodiacal light is a night-glow mostly visible along the plane of the ecliptic. It represents the background radiation associated with solar light scattered by the tenuous flattened interplanetary cloud of dust particles surrounding the Sun and the planets. It is an interesting subject of study, as the source of the micrometeoroids falling on Earth, as a link to the activity of the small bodies of the Solar System, but also as a foreground that veils the low brightness extrasolar astronomical light sources. In this review, we summarize the zodiacal light observations that have been done from the ground and from space in brightness and polarization at various wavelength ranges. Local properties of the interplanetary dust particles in some given locations can be retrieved from the inversion of the zodiacal light integrated along the light-of-sight. We show that the current community consensus favors that the majority of the interplanetary dust particles detected at 1 au originate from the activity of comets. Our current understanding of the interplanetary dust particles properties is then discussed in the context of the recent results from the Rosetta rendezvous space mission with comet 67P/Churyumov-Gerasimenko.

micrometeoroids falling on Earth, as a link to the activity of the small bodies of the Solar System, but also as a foreground that veils the low brightness extrasolar astronomical light sources.
In this review, we summarize the zodiacal light observations that have been done from the ground and from space in brightness and polarization at various wavelength ranges. Local properties of the interplanetary dust particles in some given locations can be retrieved from the inversion of the zodiacal light integrated along the light-of-sight. We show that the current community consensus favors that the majority of the interplanetary dust particles detected at 1 au originate from the activity of comets.
Our current understanding of the interplanetary dust particles properties is then discussed in the context of the recent results from the Rosetta rendezvous

Introduction
The zodiacal light is a faint night-glow, brighter as the observer looks along the ecliptic and closer to the Sun. Its visual appearance usually takes the shape of Sun-colored cones of light visible at dusk and at dawn as illustrated in Fig. 1. In fact, the zodiacal light originates from the solar light scattered by the tenuous lenticular cloud of interplanetary dust particles surrounding the Sun and in which the planets orbit. Therefore, the study of its properties is of interest to better understand the origin and characteristics of interplanetary dust particles but also to better understand how planetary and exoplanetary systems evolve with time. The night sky brightness as seen from the ground on Earth originates from a combination of i. upper atmosphere airglow, ii. the zodiacal light scattered and emitted from the interplanetary dust particles in the solar system, iii. the integrated starlight from individual and unresolved stars, iv. the diffuse galactic light scattered from interstellar dust particles and emitted from the filamentary cirrus detected by IRAS (Hauser et al., 1984) and v. the extragalactic light and background (see e.g. Levasseur-Regourd, 1994). The first scientific work on the study of the zodiacal light was published by Cassini in the 17 th century (Cassini, 1693). Scientists progressively understood that the light was present in every celestial direction, presented a solar spectrum, and was partially polarized. They therefore assumed it to be due to sunlight scattered by tiny solid particles (e.g. Wright, 1874). Quantitative measurements were pursued during the 20 th century with a view to better understand its constituting particles and to characterize this foreground signal to improve extrasolar astronomical observations. Extensive historical notes on the study of the zodiacal light can be found in e.g. Fechtig et al. (2001); Levasseur-Regourd et al. (2001).
The cloud's brightness comes principally from scattering in the visible and by thermal emission in the infrared and so can be studied in both wavelength ranges. The most complete reviews of its properties were published almost 20 years ago in the report from the IAU Commission 21: Light of the night sky "The 1997 reference of diffuse night sky brightness"  in the "Interplanetary dust" book published in 2001 (Grün et al., 2001). Since then, a number of programs and space missions have improved the available measurements of the zodiacal light and our knowledge of the properties of the interplanetary dust particles that constitute it, like the space missions Solar Mass Ejection Imager (SMEI) in the visible (Buffington et al., 2016) and AKARI in the infrared (Pyo et al., 2010). In this work we will review the properties of the zodiacal light from observations and describe the properties of the interplanetary dust particles that can be deduced from them.

Intensity observations
The initial studies of the zodiacal light originated from a need to better understand the solid matter environment surrounding the Earth, together with the fact that this glow was a source of noise to the study of faint objects and extended sources (see e.g. Dufay, 1925;Dumont, 1965;Leinert et al., 1998).
The night sky brightness includes atmospheric glow, interplanetary dust light scattering and thermal emission, as well as galactic and extra-galactic point or diffuse sources. Long term surveys can disentangle these different signal sources to retrieve each component separately. The study of the zodiacal light necessitates year-long surveys to obtain a full assessment of its geometry and determine possible time variations of the dust cloud.
The zodiacal light intensity, I, depends on the position of the observer along the terrestrial orbit, the wavelength and the viewing direction with respect to the Sun. The dust distribution around the Sun essentially follows a mirror symmetry with respect to a surface close to the ecliptic. This symmetry surface is close to a plane, but may diverge somewhat from it and is inclined by i = 3.0 • ±0.3 from the ecliptic plane as measured by the Helios probes between Earth and Venus orbits (Leinert et al., 1980). Helio-ecliptic coordinates, illustrated in Fig. 2, are generally suitable to describe the zodiacal light with the ecliptic latitude β and the helio-ecliptic longitude (λ e − λ ), where λ e is the ecliptic longitude. The phase angle of observation is noted α. As shown in Fig. 2, one needs to keep in mind that any signal observed in a given direction corresponds to an integrated signal along the line of sight.

Figure 2: geometry of observation
In the early years of space exploration a better understanding of the dust material surrounding the Earth was sought, and systematic observations of the zodiacal light were performed. High altitude surveys were performed with narrow filters to limit airglow contamination over several years in Hawaii (Weinberg, 1964). At the Canary Islands (Dumont, 1965) devised a method to remove the night airglow using a correlation between the 557.7 nm O I line and the atmospheric continuum observed on multiple points over the celestial sphere. Other astronomers' surveys were also compiled in Leinert (1975) and Fechtig et al. (1981), and include data obtained by Weinberg (1964); Blackwell et al. (1967a); Van de Noord (1970); Sparrow & Ney (1972); Leinert et al. (1976). Fig. 3 represents the zodiacal light signal averaged over a year in intensity, I, and linear degree of polarization, P (which will be discussed in Section 2.2), on a diagram representing a quarter of the celestial sphere. The values presented correspond to the unification of the survey measurements published in tabular form in the review of Leinert et al. (1998). The outer circle represents the ecliptic. The minimum brightness is obtained at a point close to the ecliptic pole coordinates (Dumont & Sanchez, 1975a,b). Taking into account the fact that the zodiacal light has a solar spectrum, its intensity I at 550 nm is then found to be close to 2.4 × 10 −5 W m −2 sr −1 µm −1 at 30 • solar elongation in the ecliptic plane (which outshines the brightest part of the Milky Way). It reaches a minimum of 7.6 × 10 −7 W m −2 sr −1 µm −1 in the vicinity of the ecliptic pole. In the ecliptic plane, a minimum of 1.8 × 10 −6 W m −2 sr −1 µm −1 is reached at around 140 • , which is still 2.3 times higher than the ecliptic pole value. The resolution of the data is between 5 • and 15 • . The relative uncertainty in the intensity data remains below 5%. One notices the continuous decrease of I as the viewing direction moves away from the Sun, except for the direction directly opposite the Sun, named the gegenschein, where an increase in brightness is due to light backscattering and I reaches a value about 1.3 times larger than the ecliptic minimum.  (Levasseur & Blamont, 1973). The Helios 1 and 2 probes showed an increase of I with decreasing heliocentric distance within the ecliptic plane (Leinert et al., , 1982. This lead to the assumption that the dust particles population density, n, was increasing towards the Sun as n(r) ∝ r −1.3 , where r is the radial distance to the Sun. Toward the outer planets, observations were performed by Pioneer 10 and 11 (Hanner et al., 1976), establishing that the zodiacal light intensity decreases away from the Sun and becomes non detectable after 3 au (Matsumoto et al., 2018). These measurements confirmed and extended the surveys that were obtained from the ground.
For an observer located on Earth's orbit, the zodiacal light intensity will vary depending on the distance to the Sun due to the Earth's orbit eccentricity.
Once this effect is corrected for, I observed at 1 au remains consistently stable within 1.5% over a whole solar cycle . The remaining annual residual oscillations originate in a slight inclination of the zodiacal cloud's symmetry plane with respect to the ecliptic. At 1 au, this symmetry plane presents an ascending node at 95 • ± 20 • with an inclination of 1.5 • ± 0.4 • (Dumont & Levasseur-Regourd, 1978). The symmetry surface of the zodiacal cloud is possibly warped far away from the Earth, as an inclination of about 3 • was inferred from Helios measurements (Leinert et al., 1980) and a dust ring associated with the orbit of Venus detected (Leinert & Moster, 2007).
A re-analysis of Weinberg's observations from Hawaii has recently yielded a new intensity map with a 2 • resolution from which the ascending node position gives around 80 • while the inclination of the cloud's symmetry plane is around 2 • (Kwon et al., 2004). Later on, the whole sky daily photometric observations by the Solar Mass Ejection Imager (SMEI) space mission were used to generate an even more precise and resolved zodiacal light map (Buffington et al., 2016).
After image processing to remove the brightest stars, daily zodiacal maps with a resolution of 0.5 • were calculated over the 8.5 years of the satellite mission (from 2003 to 2011). An average zodiacal map is shown in Fig. 4 in the same reference frame as Fig. 3. Based on these measurements, an upper limit to the zodiacal light intensity changes of 0.3% was obtained for the time period of observations of 8.5 years. Using the same satellite data, (Buffington et al., 2009) determined that the gegenschein is always located at the anti-solar point, but varies by ≈ 10% of its intensity over time, with a portion of the variation repeating seasonally. Similarly, the latest high-resolution (5' per pixel) survey of the gegenschein using a wide-angle camera showed that the maximum scattered intensity is consistently located at the anti-solar point. Ishiguro et al. (2013) fitted a model to the gegenschein geometry and constrained a very low albedo of 0.06 for the particles located in that direction.

Figure 4: Buffington zodiacal map
Local short-term enhancements in the zodiacal light brightness are also detected and were assumed to be due to possible meteoroid streams (e.g. Levasseur & Blamont, 1973). Since then the discovery of faint cometary trails in the infrared by the Infrared Astronomical Satellite (IRAS) (Sykes & Walker, 1992) and the Spitzer telescope (Reach et al., 2007) as well as linear dust features in the optical domain (Ishiguro et al., 1999;Yang et al., 2012) have clearly demonstrated the existence of many local enhancements in dust density due to the activity and impacts amongst the Solar System small bodies population (See e.g. Fig. 4 in Levasseur-Regourd & Lasue, 2019).
Extensive descriptions of the results and methods of zodiacal light surveys can be found in the reviews previously published (see Leinert, 1975;Weinberg & Sparrow, 1978;Leinert & Grün, 1990;Leinert et al., 1998;Levasseur-Regourd et al., 2001;Mann et al., 2004;Lasue et al., 2015) 2.2. Polarimetric properties 2.2.1. Degree of linear polarization observations As light scattered by an optically thin cloud of dust particles, the zodiacal light presents a systematic linear polarization, P , defined as the ratio of the difference between the light intensity perpendicular, I ⊥ , and parallel, I , to the scattering plane and the total intensity P = I ⊥ −I I ⊥ +I . The linear polarization P can reach values as high as 20% in the integrated line of sight. P is normalized and positive if the electric field polarization direction is perpendicular to the scattering plane as shown in Figure 2. Observations demonstrated that the majority of the night sky presents a positive linear polarization signal (Weinberg, 1985). Fig. 3 bottom, represents the values of polarization annually averaged over the sky in differential ecliptic coordinates based on data from Dumont & Sanchez (1976) and Table 18 published in . The error in the measured polarization is estimated to be ≈ 2%. The linear polarization is clearly dependent to the first order on the helio-ecliptic longitude. It decreases to zero towards the Sun, in the forward scattering direction, and reaches a maximum value of ≈ 20% near 60 • of helio-ecliptic longitude and near the ecliptic pole.
The linear polarization signal then presents negative values around the gegenschein with a minimum of about -2% near 165 • of elongation similar to the backscattering signal obtained with light scattered by irregular dust particles of sizes of the same order of magnitude as the light's wavelength, λ (Levasseur-Regourd et al., 1997;Lumme et al., 1997). While the negative value of polarization in that region of the sky may be difficult to correctly assess as discussed in (Dumont & Sanchez, 1975b), the presence of a negative polarization region surrounding the gegenschein was confirmed from measurements done on the ground Weinberg & Mann (1968), from sounding rockets Wolstencroft & Rose (1967) as well as from balloons Frey et al. (1974), making it a robust detection. The inversion angle corresponds to 15 • ± 5 • , with a slope at inversion of (0.2 ± 0.1) percent per degree. Therefore the integrated linear polarization variation of the zodiacal light in the ecliptic plane follows variations similar to the ones of typical polarimetric phase curves of the Solar System small bodies surfaces (e.g., regolith on the Moon and on asteroids) or dust clouds (e.g., dust in cometary comae). This trend is typical of scattering by irregular particles, with a size greater than the wavelength ( (Belskaya et al., 2019)).
Polarimetric observations have been done at different wavelengths and towards different directions of the sky as summarized in Fig. 5. There are no strong trends in the visible range, from 450 nm to 800 nm, but near infrared observations from 1 µm to 4 µm by COBE could indicate a decreasing trend towards higher wavelengths . This could be due partly to contamination by the thermal dust emission, which starts at about 3 µm (Berriman et al., 1994). Also, because the observations summarized in Fig. 5 probe different regions of the zodiacal cloud complex, the scatterers must have different origins. As different spectral polarimetric trends are seen for different types of asteroidal surfaces (Belskaya et al., 2017) and cometary dust clouds (see e.g. Hadamcik & Levasseur-Regourd, 2003) further more precise observations may indicate potential genetic links to the dust sources. The size range of observed particles may change in the F-corona compared to the zodiacal light, causing a change of their average optical properties. Typical sizes of the particles are still expected to remain larger than 1 µm (Mann et al., 2004). The change in single particle scattering properties at small angles, however, is expected to change the radial slope of the F-coronal brightness compared to the zodiacal light. Combining overall measurements, it is expected that the brightness varies as R −2.5 at the solar equator and R −2.8 at the poles, with R the radial distance to the Sun (see Levasseur-Regourd et al., 2001, and references therein).
Color trends are not evidenced from measurements at wavelengths below 1 µm (Kimura et al., 1998). At larger wavelengths, the F-corona brightness presents a reddening. The color is expected to vary within the corona as a result of the superposition of light scattering and thermal emission effects (Mann, 1993). Possible dust rings around the Sun in the region where material sublimation may occur were anticipated due to the balance between the inward Poynting-Robertson drag and the outward radiation pressure. Observations of the 1991 solar eclipse have shown no evidence for the existence of a dust ring around the Sun (see Levasseur-Regourd et al., 2001, and references therein).
Typical classical coronal models expect the F-corona polarization to decrease and tend to zero towards the Sun (Blackwell et al., 1967b). Polarimetric observations of the Solar corona (Skomorovsky et al., 2012;Burkepile et al., 2017) and data obtained by Clementine during solar occulations by the Moon (Hahn et al., 2002) have confirmed the variations measured in previous works (see reviews by Levasseur-Regourd et al., 2001;Mann et al., 2004).

Circular polarization
The light scattered by the interplanetary dust cloud is, to the first order, linearly polarized. However, multiple scattering, partial alignment of the dust particles, or even optically active materials can introduce a circular polarization component in the scattered light (see e.g. Mishchenko et al., 2002). Measurements of the circular polarization have been attempted from space and on the ground but remain inconclusive due to the very low level of signal expected (Wolstencroft & Rose, 1967;Staude & Schmidt, 1972;Wolstencroft & Kemp, 1972). Observations do not indicate circular polarization values in excess of 0.1% for most of the sky except for a possible signal in excess of 0.5% around a helio-ecliptic longitude of 240 • (Wolstencroft & Kemp, 1972). In the case of the zodiacal cloud, one expects the major processes to explain a possible circular polarization to be either: i.particle alignment Lasue et al., 2015), or ii. optically active materials. Both processes have been suggested to explain the very low circular polarization measured for comets with recent measurements by Rosetta favouring a particle alignment (Rosenbush et al., 2007;Kolokolova et al., 2016). Further attempts to observe the zodiacal light circular polarization may be of importance as it could constrain the origin of the particles and specific dynamic properties in the solar magnetic field.

Spectral properties
In the visible part of the spectrum, the zodiacal light shows a solar spectrum.
Looking at the spectral properties of the zodiacal light integrated along the line-of-sight at larger wavelengths, a number of interesting measurements have been performed, especially with space missions such as IRAS, COBE, ISO, AKARI, WISE and Planck, to deduce the physical properties and the dynamics of the dust particles in the interplanetary dust cloud.
First, many infrared measurements have put into evidence the presence of dust trails and rings associated with the impacts of asteroids and the activity of comets (e.g. Dermott et al., 1984;Sykes, 1988;Reach, 1997) The spectral analysis of the Doppler shifted scattered solar light Fraunhofer lines gives insight into the dynamics of the dust particles constituting the zodiacal cloud. Ipatov et al. (2008) have shown the measurements to be consistent with a majority of the dust particles associated with cometary activity close to the Sun (< 1 au).
Then, a red spectral gradient of the scattered light is determined by many measurements in the range from 200 nm to 2 µm, as shown by COBE (Berriman et al., 1994), and as summarized in , and references therein).
The spectral gradient can be calculated to be S = 8.5 ± 1.0%/100 nm at 460 nm (Yang & Ishiguro, 2015). Combining this value with the albedo of A = 0.06 ± 0.01 measured at the same wavelength in the direction of the gegenschein by Ishiguro et al. (2013) shows that such particles scattering properties are most consistent with the properties of primitive small bodies such as comets or D-type asteroids (Yang & Ishiguro, 2015).
Measurements (from balloons, rockets and satellites) of the thermal emission originating from the zodiacal cloud have also been performed. The zodiacal thermal emission, which gradually takes over the solar-spectrum by about 2 µm, is actually the most prominent source of brightness in the 5 to 100 micrometer range, at least away from the galactic plane (e.g. Leinert et al., 1998).
For even larger wavelengths, the sky brightness gets dominated by the interstellar medium, the cosmic diffuse infrared background and finally the cosmic microwave background.
This thermal range can be modelled using blackbody emission curves for specific mixtures of materials (like silicates, amorphous carbon and graphite).
A change of slope in the spectrum at about 150 µm is interpreted to represent a change in the dust size distribution at a radius of about 30 µm (Fixsen & Dwek, 2002), which is consistent with the size distribution variations measured in situ (Grün et al., 1985). Comparison of the model with the data by Fixsen & Dwek (2002) indicates that replenishment of the interplanetary dust is ≈ 10 11 kg yr −1 from asteroids and comets. Dust emission rate from asteroids impacts and activity appears to be lower and thus, a significant contribution from comets would be needed to replenish the interplanetary dust cloud at 1 au.
The work by Reach et al. (2003) with the ISO satellite data also indicates a dust distribution dominated by large (> 10 µm radius), low-albedo (< 0.08), rapidly-rotating, gray particles at 1 au from the Sun. Spectral signatures of amorphous and crystalline silicates are also detected in the infrared from 9 to 11 µm indicating also the presence of small particles similar to the ones ejected by active comets, such as C/1995 O1 Hale-Bopp comet, and to proto-planetary disks emission spectra. Similar spectral signatures of crystalline silicate features were obtained with the AKARI spacecraft (Ootsubo et al., 2009).
At even larger wavelengths, in the microwave domain, Planck measurements have shown that the zodiacal cloud is a minor correction to the Cosmic Microwave Background emission (Ade et al., 2014).

Interplanetary dust local properties from local inversion
So far, we have discussed only the properties of the zodiacal light, which is an integrated signal (Fig. 1) along the line of sight of the observer (Fig. 2).
However, techniques have been tentatively developed to invert this signal in some cases, thus giving information on the local properties of the brightness and the polarization of the zodiacal light. This is what we describe in this section.

Inversion method
The dust in the zodiacal cloud cannot be assumed to have the same optical properties everywhere and in all directions. In the data integrated over the line of sight, contributions from light scattering by dust particles corresponding to different distances to the Sun, different observing phase angles and with likely different physical properties are mixed. In order to obtain specific properties of the dust particles, it is mandatory to invert the integrated intensity and polarization. Knowing that the zodiacal cloud is seen to be relatively homogeneous and in an apparent steady state, which is consistent with the mixing of injected Assuming that the zodiacal cloud has a plane of symmetry, within which the observer is located, and that the density of dust particles varies as a power-law with the heliocentric distance and is in an equilibrium state, a rigorous inversion of the zodiacal light intensity integrated along the line-of-sight can be derived from Earth-based observations as a function of the phase angle at 1 au (Dumont, 1973). Given the local values for the light intensity polarized parallel and perpendicular to the scattering plane, the associated local linear polarization can be directly calculated. This method is used to determine that at a phase angle α = 90 • at the symmetry plane, one cubic kilometer of interplanetary space scatters sunlight with an intensity that is 4 × 10 −34 times smaller than the solar brightness (Schuerman, 1979). The local linear polarization is calculated to be about 0.30 ± 0.03 (Levasseur-Regourd et al., 1990). Dumont (1983); Dumont & Levasseur-Regourd (1985) developed the method of the "nodes of lesser uncertainty" in order to invert the measurements outside the Earth's orbital path. This method assumes rotational symmetry for the symmetry plane of the zodiacal cloud, a radial power-law decrease in the dust density from the Sun with a quasi-equilibrium state of the cloud density on the time-scale of a few weeks. Knowing that the local intensity scattering function is relatively smooth and regular, it can be represented by continuous mathematical functions (e.g. polynomials, or Fourier series) with five parameters out of which one is free. By using two lines of sight with solar elongations and 180 • − , these mathematical functions present two regions where the local brightness varies over a very limited range allowing to retrieve its value with a small uncertainty in those positions. These mathematical regions correspond to : i. the socalled"radial node", located near the helio-ecliptic meridian with a variable solar distance, R, below 1 au and a constant phase angle, α = 90 • , and ii. the socalled "martian node" located at a quasi-constant solar distance of R ≈ 1.5 au but covering a phase angle range from 0 • to 42 • giving information on the backscattering function of the dust particles. Figure 6 illustrates the method and the position of the nodes. Note that the regions were named for illustration purposes but do not correspond to regions where measurements were done.
They were retrieved with limited assumptions by mathematical inversion of the line of sight integral. Some further methods of inversion have also been devised to retrieve local contributions in the plane tangential to the Earth's orbit, out of the ecliptic (Renard et al., 1995) and at other positions in the solar system (Schuerman, 1979). Improvements in the zodiacal light data collection will lead to further work on the local intensity and polarization characterization with a better spatial and angular resolution.   Table 1: Variation of the local properties of the dust in the symmetry plane. The properties are described as a function of solar distance with a power law assumption. The optical properties (linear polarization, albedo, and intensity in W m −2 sr −1 µm −1 rad −1 at 550 nm) are retrieved at α = 90 • Levasseur-Regourd et al., 2001). The density of the particles is inferred from the inverted values in intensity and polarization.
A summary of the local properties of the interplanetary dust (intensity, polarization, temperature, albedo, number density) in the symmetry plane as inverted from the data and at α = 90 • is presented in Table 1. The local values for the albedo at a given phase angle (around 90 degrees) is retrieved by combining the local values of the flux emission, the temperature and the intensity (Hanner et al., 1981).

Heliocentric variations
At the radial node, the inverted values for the local variation in intensity and polarization are presented in Fig. 6 at α = 90 • as a function of the heliocentric distance. The figure is expressed with the traditional zodiacal light brightness values given in S 10 (V ), i.e. the equivalent number of 10 th visual magnitude solar type stars per square degree. This unit was historically useful since the zodiacal light presents a solar spectrum in the visible wavelengths and so was valid over the visual domain. While SI units, i.e., W m −2 sr −1 µm −1 are easier to understand, they nevertheless need to be associated to a precise wavelength (1 S 10 (V ) is equal to 1.261 × 10 −8 W m −2 sr −1 µm −1 at 550 nm, see Levasseur-Regourd et al., 2001).
The brightness follows the power-law variation given in Table 1 and is con-sistent with the local dust density distribution under a homogeneous zodiacal cloud hypothesis as described in Leinert et al. (1977). However, taking into account the variation in albedo with respect to the heliocentric distance gives a local dust particle density value following a power-law with an exponent of −0.93±0.07 (significantly different from the -1.3 value obtained under the homogeneity assumption) (Levasseur-Regourd et al., 1991). The Poynting-Robertson effect predicts a 1/R dependence for the density of dust particles in circular orbits located within the region of production of dust particles (ranging from 0.1 au to about 10-20 au), which is consistent with the power-law coefficient thus determined (Hanner, 1980).
The degree of linear polarization, P , variation includes data obtained from a compilation by Fechtig et al. (1981), data inverted by , data from the Dumont & Sanchez (1975b) survey and data observed close to the Sun inverted by . Numerical simulations from Lasue et al. (2007) and laboratory measurements from Hadamcik et al. (2018) are also superposed to the observations. P decreases with decreasing solar distance.
The sharp decrease noticed below 0.3 au suggests that drastic changes to the dust particles properties take place closer to the Sun, in agreement with trends obtained in the F-corona (Mann & MacQueen, 1996).
The dust properties variation with the solar distance could be related to temporal evolution of the dust particles, with possible changes in their composition and/or their size distribution, as they spiral towards the Sun under the Poynting-Robertson effect and get progressively warmer. This will be discussed in Section 4.

Phase angle variations
The local brightness and polarization phase curves of the zodiacal light at the martian node (≈ 1.5 au) are plotted in Fig. 7. One can notice that around the backscattering domain (α < 20 • ), the intensity decreases by 0.015 ± 0.05 magnitude per degree. This demonstrates the existence of an opposition effect for the interplanetary dust population (Muinonen & Lumme, 1991).
The overall shape of the polarization phase curve is smooth, with a slightly negative polarization in the backscattering region, an inversion angle around (15 ± 5) • and a slope at inversion of (0.2 ± 0.1) percent per degree. This shape is typical of light interaction with irregular dust particles possibly in the shape of fluffy aggregates made of submicronic monomers (Levasseur-Regourd et al., 1997;Lumme et al., 1997). The behavior of both intensity and polarization phase curves is comparable to that of cometary dust or C-type asteroids.
The slope at inversion of the polarimetric phase curve can be related to the geometric albedo by the 'Umov empirical law' (Umow, 1905). From the inverted data, a geometric albedo value of 0.15 ± 0.08 can be calculated at 1 au (Levasseur-Regourd, 1998), which is consistent within the error bars with the value given in Table 1 and also consistent with the more recent analysis done by Ishiguro et al. (2013).
Finally, inversion of the local properties of the dust particles out of the symmetry plane of the zodiacal cloud in the plane perpendicular to the ecliptic and tangential to Earth's orbit, have also been performed and normalized to the data inverted at 1 au and for α = 90 • . The brightness decreases with the elevation, indicating a flattened geometry for the interplanetary dust cloud.
The polarization also decreases with elevation by about 10% at 1 au while the albedo increases simultaneously from 0.8 to 0.11 (Renard et al., 1995). These variations can be explained by the presence of two contributing dust populations, one coming from the asteroid belt and short period comets activity close to the ecliptic, and the other one, more isotropic, and originating from long period comets and with different properties.  (Giese et al., 1978;Haudebourg et al., 1999;Nakamura & Okamoto, 1999). Furthermore, a red polarimetric color variation of the scattered light by irregular particles is linked to a lower porosity of the dust particles. It is therefore expected that compact dust particles, possibly in the shape of elongated spheroids to present less symetries than spheres, could be significant contributors to the interplanetary dust distribution Kolokolova & Kimura, 2010).
As shown in Fig. 7, the variation in local scattered intensity and polarization gives information on the heliocentric changes in dust density and properties. To the first order, the scattered intensity can be explained by the dust density variation consistent with the Poynting-Robertson effect.
From numerical simulations in Lasue et al. (2007), that have been confirmed by laboratory simulations , the decrease in polarization sunward around 1 au can be explained to first order by a change in the optical index of the dust particles population scattering the light. In that hypothesis, the proportion of absorbing materials, such as carbonaceous compounds, to weakly absorbing ones would decrease sensibly towards the Sun (from 50-100% at 1.5 au to about 0% at 0.6 au). A plot of the numerical simulations and experimental measurements data points are shown on Fig. 7b. This explanation would be consistent with an increase in the dust population albedo sunward, which is expressed by the variation equation shown for the albedo in Tab. 1.
Other changes in the dust population properties may also occur.
Still in Fig. 7, closer to the Sun, below 0.3 au, one notices a sharp decrease in polarization as well as a change of slope in the scattered intensity. This change in behavior indicates that some properties of the dust cloud (size distribution, morphology, composition, etc.) is abruptly changing in that region. Many materials may sublimate in that region and are certainly contributors to changes in either size distribution, shape or optical indices of the dust particles (Mann et al., 2004).
As interplanetary dust particles spiral inwards and vanish in the Sun, with a typical lifetime of a couple of million years at 1 au (Burns et al., 1979), the interplanetary dust cloud requires constant replenishment. Comet nuclei activity as well as asteroidal collisions in the Main Belt can provide substantial amounts of dust particles and have long been considered to be major contributors to the inner zodiacal cloud (Whipple, 1950;Dermott et al., 1984). Other minor contributions come from the giant planets environment, the Kuiper Belt, and from interstellar dust particles (see e.g. Grün et al., 1993). Observations of the visual zodiacal intensity within 1 au by the Clementine space probe fitted with a model of dust distribution including asteroids, Jupiter Family Comets (JFC), Halley-type comets and Oort cloud comets indicates that at least 89% of the interplanetary dust located in a sphere within 1 au from the Sun originate from comets (Hahn et al., 2002). Numerical simulations of the IRAS observations in the thermal domain of the zodiacal cloud based on the dynamical properties of the dust ejected from comets and from asteroids indicates that 85-95% of the mid-infrared emission visible from 1 au is produced by JFCs, while asteroidal and long-period comets are both respectively smaller than 10% (Nesvorný et al., 2010(Nesvorný et al., , 2011. Similar conclusions are reached with the dynamical study of the evolution of aggregate dust particles injected in the inner solar system by active comets and with a size frequency distribution like the one of 67P/Churyumov-Gerasimenko, as they can reproduce the interplanetary dust properties detected around the Earth's orbit (Yang & Ishiguro, 2018). Consistent results were also obtained from the study of the Doppler shift from the Fraunhofer lines of the scattered solar light (Ipatov et al., 2008). An improved model of the IRAS and COBE thermal emission observations by Rowan-Robinson & May (2013) determined similar numbers with, in addition, less than 1% interstellar dust contribution at 1 au. Finally, lidar observations of Na and Fe fluxes in the Earth's atmosphere coupled with the cosmic spherule accretion rate measured at the South Pole also indicates that JFCs contribute 87 ± 17% of the total extraterrestrial matter mass input on the Earth (Carrillo-Sánchez et al., 2016).
In summary, based on those diverse lines of evidence, and others described in the previous sections, a consensus appears to emerge that the majority of the dust particles located within 1 au should be contributed from the activity of comets, and in particular the JFCs.

Properties of interplanetary dust and their link with comets
As we have seen above, the dust particles at the origin of the zodiacal light are expected to present specific characteristics based on the properties of the scattered and emitted zodiacal light in the visible and the thermal wavelengths.
Several lines of evidence based on different arguments make a compelling case to link the properties of the interplanetary dust particles to those of cometary dust particles.
First, the dust density distribution models and related dynamics indicate that the majority of dust particles within 1 au from the Sun probably originate in major part from the activity of short period comets. This observation is also consistent with some of the properties detected for the meteoroids population impacting the Earth. For example, almost all meteor showers have been linked with cometary parent bodies (Jenniskens, 2006). The properties of interplanetary dust particles collected in the stratosphere of the Earth by airplanes are also consistent with primitive material and represent a family of samples that are chemically unlike the extraterrestrial matter represented by the meteorites (originating from asteroids apart from possibly one or two exceptions) (Flynn et al., 2016;Koschny et al., 2019). A population of micrometeorites named as UCAMM (Ultra-Carbonaceous Antarctic MicroMeteorites) for their very high carbon content (from 48% to 85%) has also been studied and present a D/H ratio more than one order of magnitude larger than the terrestrial values strongly suggesting a cometary origin (Duprat et al., 2010).
Second, the properties of the zodiacal light are consistent with a mixture of irregular fluffy particles and compact particles formed of silicates and carbonaceous materials. These properties are consistent with the properties of interplanetary dust collected in the stratosphere (Flynn et al., 2016). They are also coherent with the properties of cometary dust particles as observed from the ground, analyzed in situ by sample return from Stardust and more recently from the in situ analysis performed by the Rosetta space mission (see e.g. Levasseur-Regourd et al., 2018). Most notably, the presence of both fluffy aggregates of grains and compact particles have been evidenced by both Stardust and Rosetta space missions (Hörz et al., 2006;Güttler et al., 2019). Similarly, both the Stardust samples and the Rosetta measurements have shown a distinct population of carbonaceous material forming a major fraction of the ejected dust with characteristics different from meteorites and possibly atmospheric IDPs (Sandford et al., 2006;Fray et al., 2016). Finally, the scattered light phase curves present a similar behavior between ejected cometary dust and interplanetary dust particles analogues .
Overall the study of the zodiacal light can give us access to properties of the cometary particles that may not be accessible by other means, since samples that arrive on Earth may be biased, and cometary space missions are limited in numbers. Continued monitoring of the zodiacal light over time will provide a better understanding of the origin and evolution of the interplanetary dust environment. Going one step further, the study of the similarities between the light scattering behavior of interplanetary dust particles, asteroids, cometary dust clouds and protoplanetary disks will help us understand planetary formation processes and the origin of the solid particles in these different environments (Levasseur-Regourd et al., 2020).

Conclusions and future observations
In this work, we attempted to give an overview of the properties of the Intensity and polarization values of the zodiacal light are also of major importance for outer Solar System observations, since they provide an estimate of the foreground noise induced by the zodiacal light, together with an optimization of the epochs of observations of faint extended astronomical objects Levasseur-Regourd & Lasue, 2019) We can expect the next generations of sensitive whole sky surveys to provide many additional constraints on the zodiacal light with the possible detection of additional local heterogeneities and global time variability. Such surveys will include ground-based telescopes such as the Large Synoptic Survey Telescope (LSST), an optical 8.4m telescope designed to survey the visible sky every week starting in 2020 (Abell et al., 2009), but also space missions such as the Messier mission, a space telescope that will scan the night sky in the 200-1000 nm range, reaching surface brightness levels of 34 and 37 mag arcsec −2 in the UV and optical domains (Valls-Gabaud & collaboration, 2016;Muslimov et al., 2017). These results will need to be compared with the results obtained by the space missions studying the small bodies of the Solar System, such as Hayabusa2 and OSIRIS-REx missions for asteroids and the newly selected ESA Comet Interceptor mission to flyby a new pristine active comet (Gater, 2019). Further data out of the Earth's orbit could be obtained closer to the Sun from the recently launched Parker Solar Probe, and by other probes tangentially to the orbits of planets or small bodies they may be orbiting. Also, future observations of the zodiacal light may be sensitive to circular polarization, which could constrain the origin of the particles and specific dynamic properties in the solar magnetic field.
Some important unsolved questions related to the zodiacal light that need to be studied in the future include: • the exact origin and physical properties of the dust particles contributing to the zodiacal cloud, • the possible time evolution of the zodiacal cloud over decades, • the circular polarization that can probe the interplanetary magnetic field structures, • the dust evolution in the vicinity of the Sun.
• how the properties of the zodiacal cloud could help us constrain the properties of exo-zodiacal dust clouds, which may behave quite differently depending on the presence and the properties of orbiting exoplanets.
Answers will progressively be found as future observations continue to focus on higher resolution photopolarimetry of the zodiacal light at all wavelengths and possibly observed out of the ecliptic. Improvements on the inversion methods based on such measurements will also likely further improve zodiacal light studies.

Acknowledgements
This work was supported by the Programme National de Planétologie (PNP) of CNRS/INSU, co-funded by CNES. We acknowledge the organizers of the Meteoroids 2019 conference for inviting the authors to present this work. We thank two anonymous referees for providing interesting discussions and improvements on the manuscript.    Skylab at the North Celestial Pole (Weinberg & Hahn, 1980); rocket Astro 7 (λe − λ = 25 • ) (Pitz et al., 1979); Helios (λe − λ = 90 • ) (Leinert et al., 1982); COBE (λe − λ = 90 • ) (Berriman et al., 1994).  Renard et al., 1995). Note that these regions correspond to places where mathematical inversion of the line-of-sight integral is possible with minimal assumptions, and not where measurements were made. See text for additional details.