Jeremie Lasue
email: jlasue@irap.omp.eu

Anny-Chantal Levasseur-Regourd

Jean-Baptiste Renard

Zodiacal light observations and its link with cosmic dust: a review

Keywords: Zodiacal light, brightness, polarization, thermal emission, interplanetary dust, cometary dust

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

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 [START_REF] Hauser | IRAS observations of the diffuse infrared background[END_REF] and v. the extragalactic light and background (see e.g. [START_REF] Levasseur-Regourd | Natural Background Radiation-the Light from the Night Sky[END_REF]. The first scientific work on the study of the zodiacal light was published by Cassini in the 17 th century [START_REF] Cassini | Découverte de la lumière celeste qui paroist dans le zodiaque[END_REF]. 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. [START_REF] Wright | On the Polarization of the zodiacal Light[END_REF]. 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. [START_REF] Fechtig | Historical perspectives[END_REF]; [START_REF] Levasseur-Regourd | Optical and thermal properties of interplanetary dust[END_REF]. 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" [START_REF] Leinert | The 1997 reference of diffuse night sky brightness[END_REF] and in the "Interplanetary dust" book published in 2001 [START_REF] Grün | Interplanetary dust[END_REF]. 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 [START_REF] Buffington | Measurements and an empirical model of the Zodiacal brightness as observed by the Solar Mass Ejection Imager (SMEI)[END_REF] and AKARI in the infrared [START_REF] Pyo | Brightness map of the zodiacal emission from the AKARI IRC All-Sky Survey[END_REF]. 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.

Integrated observations of the zodiacal light along the line of sight

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. [START_REF] Dufay | La polarisation de la lumière zodiacale[END_REF][START_REF] Dumont | Séparation des composantes atmosphériques et interplanétaire et stellaire du ciel nocturne à 5000 Å. Applications à la photométrie de la lumière zodiacale et du Gegenschein[END_REF][START_REF] Leinert | The 1997 reference of diffuse night sky brightness[END_REF].

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 [START_REF] Leinert | The plane of symmetry of interplanetary dust in the inner solar system[END_REF]. 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 [START_REF] Weinberg | The zodiacal light at 5300 Å[END_REF]. At the Canary Islands [START_REF] Dumont | Séparation des composantes atmosphériques et interplanétaire et stellaire du ciel nocturne à 5000 Å. Applications à la photométrie de la lumière zodiacale et du Gegenschein[END_REF] 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 [START_REF] Leinert | Zodiacal light-a measure of the interplanetary environment[END_REF] and [START_REF] Fechtig | Interplanetary dust and zodiacal light[END_REF], and include data obtained by [START_REF] Weinberg | The zodiacal light at 5300 Å[END_REF]; Blackwell et al. (1967a); [START_REF] Van De Noord | Observations of the zodiacal light with a balloonborne telescope[END_REF]; [START_REF] Sparrow | Observations of the Zodiacal Light from the Ecliptic to the Poles[END_REF]; [START_REF] Leinert | Interpretation of a rocket photometry of the inner zodiacal light[END_REF]. 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 [START_REF] Leinert | The 1997 reference of diffuse night sky brightness[END_REF]. 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. [START_REF] Levasseur | Satellite observations of intensity variations of the zodiacal light[END_REF]. The Helios 1 and 2 probes showed an increase of I with decreasing heliocentric distance within the ecliptic plane [START_REF] Leinert | The zodiacal light from 1.0 to 0.3 AU as observed by the Helios space probes[END_REF][START_REF] Leinert | Helios zodiacal light measurements-a tabulated summary[END_REF]. 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 [START_REF] Hanner | 1.1. 4 Pioneer 10 Observations of Zodiacal Light Brightness near the Ecliptic: Changes with Heliocentric Distance[END_REF], establishing that the zodiacal light intensity decreases away from the Sun and becomes non detectable after 3 au [START_REF] Matsumoto | Zodiacal light beyond earth orbit observed with Pioneer 10[END_REF]. 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 (Levasseur-Regourd, 1996). 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 • [START_REF] Dumont | Zodiacal light photopolarimetry. IV-Annual variations of brightness and the symmetry plane of the zodiacal cloud: Absence of solar-cycle variations[END_REF]. 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 [START_REF] Leinert | The plane of symmetry of interplanetary dust in the inner solar system[END_REF] and a dust ring associated with the orbit of Venus detected [START_REF] Leinert | Evidence for dust accumulation just outside the orbit of Venus[END_REF].

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 • [START_REF] Kwon | An observational model of the zodiacal light brightness distribution[END_REF]. 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 [START_REF] Buffington | Measurements and an empirical model of the Zodiacal brightness as observed by the Solar Mass Ejection Imager (SMEI)[END_REF].

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, [START_REF] Buffington | Measurements of the gegenschein brightness from the Solar Mass Ejection Imager (SMEI)[END_REF] 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. [START_REF] Ishiguro | High-resolution imaging of the gegenschein and the geometric albedo of interplanetary dust[END_REF] 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. [START_REF] Levasseur | Satellite observations of intensity variations of the zodiacal light[END_REF]. Since then the discovery of faint cometary trails in the infrared by the Infrared Astronomical Satellite (IRAS) [START_REF] Sykes | Cometary dust trails: I[END_REF] and the Spitzer telescope [START_REF] Reach | A survey of debris trails from short-period comets[END_REF] as well as linear dust features in the optical domain [START_REF] Ishiguro | First detection of visible zodiacal dust bands from ground-based observations[END_REF][START_REF] Yang | High-Resolution Map of Zodiacal Dust Bands by WIZARD Measurements. In Asteroids, Comets[END_REF] 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 [START_REF] Levasseur-Regourd & Lasue | Figure 7: Variation of the local near-ecliptic brightness (top) and linear degree of polarization, P[END_REF].

Extensive descriptions of the results and methods of zodiacal light surveys can be found in the reviews previously published (see [START_REF] Leinert | Zodiacal light-a measure of the interplanetary environment[END_REF][START_REF] Weinberg | Zodiacal light as an indicator of interplanetary dust[END_REF][START_REF] Leinert | Interplanetary dust[END_REF][START_REF] Leinert | The 1997 reference of diffuse night sky brightness[END_REF][START_REF] Levasseur-Regourd | Optical and thermal properties of interplanetary dust[END_REF][START_REF] Mann | Dust near the Sun[END_REF][START_REF] Lasue | Polarimetry of the interplanetary dust cloud[END_REF] 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 [START_REF] Weinberg | Zodiacal light and interplanetary dust[END_REF]. 18 published in [START_REF] Leinert | The 1997 reference of diffuse night sky brightness[END_REF]. 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, λ [START_REF] Levasseur-Regourd | Physical properties of dust in the solar system: relevance of a computational approach and of measurements under microgravity conditions[END_REF][START_REF] Lumme | Light scattering by dense clusters of spheres[END_REF]. 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 [START_REF] Weinberg | Negative polarization in the zodiacal light[END_REF], from sounding rockets [START_REF] Wolstencroft | Observations of the zodiacal light from a sounding rocket[END_REF] as well as from balloons [START_REF] Frey | Photometry of the zodiacal light with the balloon-borne telescope THISBE[END_REF], 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 ( [START_REF] Belskaya | Optical Polarimetry of Small Solar System Bodies: From Asteroids to Debris Disks[END_REF]).

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 [START_REF] Leinert | The 1997 reference of diffuse night sky brightness[END_REF]. This could be due partly to contamination by the thermal dust emission, which starts at about 3 µm [START_REF] Berriman | COBE DIRBE nearinfrared polarimetry of the zodiacal light: Initial results[END_REF]. 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 [START_REF] Belskaya | Refining the asteroid taxonomy by polarimetric observations[END_REF] and cometary dust clouds (see e.g. Levasseur-Regourd et al., 1996;[START_REF] Hadamcik | evolution of comet C/1995 O1 (Hale-Bopp) by imaging polarimetric observations[END_REF] further more precise observations may indicate potential genetic links to the dust sources.

Solar F-corona observations

Closer to the Sun, the zodiacal cloud starts to merge with the solar corona.

The brightness of the solar corona originates mainly from three physical processes: i. Thomson scattering by free electrons (K-corona, major component at < 4 solar radii), ii. line emission from ions (L-corona), and iii. light scattering from dust particles (F-corona). The F-corona emission dominates at the largest distances and can be separated from the other emissions by spectrometric or polarimetric analyses.

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 [START_REF] Mann | Dust near the Sun[END_REF]. 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 [START_REF] Levasseur-Regourd | Optical and thermal properties of interplanetary dust[END_REF] references therein).

Color trends are not evidenced from measurements at wavelengths below 1 µm [START_REF] Kimura | Influence of dust shape and material composition on the solar F-corona[END_REF]. 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 [START_REF] Mann | The influence of circumsolar dust on the whitelight corona-study of the visual F-corona brightness[END_REF]. 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 [START_REF] Blackwell | The distribution of dust in interplanetary space[END_REF]. Polarimetric observations of the Solar corona [START_REF] Skomorovsky | White-light observations and polarimetric analysis of the solar corona during the eclipse of 1 August 2008[END_REF][START_REF] Burkepile | Polarization observations of the total solar eclipse of August 21[END_REF] and data obtained by Clementine during solar occulations by the Moon [START_REF] Hahn | Clementine observations of the zodiacal light and the dust content of the inner solar system[END_REF] have confirmed the variations measured in previous works (see reviews by [START_REF] Levasseur-Regourd | Optical and thermal properties of interplanetary dust[END_REF][START_REF] Mann | Dust near the Sun[END_REF].

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. [START_REF] Mishchenko | Scattering, absorption, and emission of light by small particles[END_REF]. 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 [START_REF] Wolstencroft | Observations of the zodiacal light from a sounding rocket[END_REF][START_REF] Staude | Circular polarisation measurements of the zodiacal light[END_REF][START_REF] Wolstencroft | Circular polarization of the nightsky radiation[END_REF]. 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 • [START_REF] Wolstencroft | Circular polarization of the nightsky radiation[END_REF]. In the case of the zodiacal cloud, one expects the major processes to explain a possible circular polarization to be either: i.particle alignment [START_REF] Lazarian | Grain Alignment: Role of Radiative Torques and Paramagnetic Relaxation[END_REF][START_REF] Lasue | Polarimetry of the interplanetary dust cloud[END_REF], 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 [START_REF] Rosenbush | Circular polarization in comets: Observations of comet C/1999 S4 (LINEAR) and tentative interpretation[END_REF][START_REF] Kolokolova | Clues to cometary circular polarization from studying the magnetic field in the vicinity of the nucleus of comet 67P/Churyumov-Gerasimenko[END_REF]. 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. [START_REF] Dermott | Origin of the solar system dust bands discovered by IRAS[END_REF][START_REF] Sykes | IRAS observations of extended zodiacal structures[END_REF][START_REF] Reach | The structured zodiacal light: IRAS, COBE, and ISO observations[END_REF] 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. [START_REF] Ipatov | Dynamical zodiacal cloud models constrained by high resolution spectroscopy of the zodiacal light[END_REF] 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 [START_REF] Berriman | COBE DIRBE nearinfrared polarimetry of the zodiacal light: Initial results[END_REF], and as summarized in (Leinert et al., 1998, and references therein).

The spectral gradient can be calculated to be S = 8.5 ± 1.0%/100 nm at 460 nm [START_REF] Yang | Origin of interplanetary dust through optical properties of zodiacal light[END_REF]. Combining this value with the albedo of A = 0.06 ± 0.01 measured at the same wavelength in the direction of the gegenschein by [START_REF] Ishiguro | High-resolution imaging of the gegenschein and the geometric albedo of interplanetary dust[END_REF] shows that such particles scattering properties are most consistent with the properties of primitive small bodies such as comets or D-type asteroids [START_REF] Yang | Origin of interplanetary dust through optical properties of zodiacal light[END_REF].

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. [START_REF] Leinert | The 1997 reference of diffuse night sky brightness[END_REF].

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 [START_REF] Fixsen | The zodiacal emission spectrum as determined by COBE and its implications[END_REF], which is consistent with the size distribution variations measured in situ [START_REF] Grün | Collisional balance of the meteoritic complex[END_REF]. Comparison of the model with the data by [START_REF] Fixsen | The zodiacal emission spectrum as determined by COBE and its implications[END_REF] 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 [START_REF] Reach | The midinfrared spectrum of the zodiacal and exozodiacal light[END_REF] 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 [START_REF] Ootsubo | Mid-Infrared spectrum of the zodiacal light observed with AKARI/IRC[END_REF].

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 [START_REF] Ade | Planck 2013 results[END_REF].

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 dust particles and their short lifetimes [START_REF] Campbell-Brown | High resolution radiant distribution and orbits of sporadic radar meteoroids[END_REF][START_REF] Nesvorný | Dynamical model for the zodiacal cloud and sporadic meteors[END_REF][START_REF] Yang | Evolution of cometary dust particles to the orbit of the Earth: Particle size, shape, and mutual collisions[END_REF], it is possible to invert the light integrated over the line-of-sight to deduce local scattering properties of the dust particles. Local rigorous inversion and local inversions through mathematical methods are, together with their main results, summarized in [START_REF] Levasseur-Regourd | Optical and thermal properties of interplanetary dust[END_REF], section 4.

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 [START_REF] Dumont | Phase function and polarization curve of interplanetary scatterers from zodiacal light photopolarimetry[END_REF]. 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 [START_REF] Schuerman | Inverting the zodiacal light brightness integral[END_REF]. The local linear polarization is calculated to be about 0.30 ± 0.03 [START_REF] Levasseur-Regourd | A comparison between polarimetric properties of cometary dust and interplanetary dust particles[END_REF]. [START_REF] Dumont | Zodiacal light gathered along the line of sight: The vicinity of the terrestrial orbit studied with photopolarimetry and with Doppler spectrometry[END_REF]; [START_REF] Dumont | Zodiacal Light gathered along the line of sight: Retrieval of the local scattering coefficient from photometric surveys of the ecliptic plane[END_REF] 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 [START_REF] Renard | Properties of interplanetary dust from infrared and optical observations. II. Brightness, polarization, temperature, albedo and their dependence on the elevation above the ecliptic[END_REF] and at other positions in the solar system [START_REF] Schuerman | Inverting the zodiacal light brightness integral[END_REF]. 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, 1996;[START_REF] Levasseur-Regourd | Optical and thermal properties of interplanetary dust[END_REF]. 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 [START_REF] Hanner | On the definition of albedo and application to irregular particles[END_REF].

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 [START_REF] Levasseur-Regourd | Optical and thermal properties of interplanetary dust[END_REF].

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 [START_REF] Leinert | Observations of zodiacal light from Helios 1 and 2[END_REF]. 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) [START_REF] Levasseur-Regourd | The zodiacal cloud complex[END_REF]. 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 [START_REF] Hanner | On the albedo of the interplanetary dust[END_REF].

The degree of linear polarization, P , variation includes data obtained from a compilation by [START_REF] Fechtig | Interplanetary dust and zodiacal light[END_REF], data inverted by Levasseur-Regourd (1996), data from the Dumont & Sanchez (1975b) survey and data observed close to the Sun inverted by Mann (1998). Numerical simulations from [START_REF] Lasue | Inferring the interplanetary dust properties-from remote observations and simulations[END_REF] and laboratory measurements from [START_REF] Hadamcik | Analogues of interplanetary dust particles to interpret the zodiacal light polarization[END_REF] 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 [START_REF] Mann | Observation and analysis of the F-corona brightness[END_REF].

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 [START_REF] Muinonen | Light scattering by solar system dust: the opposition effect and the reversal of polarization[END_REF].

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 [START_REF] Levasseur-Regourd | Physical properties of dust in the solar system: relevance of a computational approach and of measurements under microgravity conditions[END_REF][START_REF] Lumme | Light scattering by dense clusters of spheres[END_REF]. 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' [START_REF] Umow | Chromatische depolarisation durch lichtzerstreuung[END_REF]. From the inverted data, a geometric albedo value of 0.15 ± 0.08 can be calculated at 1 au [START_REF] Levasseur-Regourd | Zodiacal light, certitudes and questions[END_REF], which is consistent within the error bars with the value given in Table 1 and also consistent with the more recent analysis done by [START_REF] Ishiguro | High-resolution imaging of the gegenschein and the geometric albedo of interplanetary dust[END_REF].

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 [START_REF] Renard | Properties of interplanetary dust from infrared and optical observations. II. Brightness, polarization, temperature, albedo and their dependence on the elevation above the ecliptic[END_REF]. 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. The intensity and polarization phase curves can be compared to numerical and experimental simulations in order to retrieve the likely physical properties of the scatterers.

A classical assumption, based on astronomical observations, is that interplanetary, cometary and interstellar dust particles consist of silicates, carbonaceous compounds and possibly ices (see e.g. [START_REF] Greenberg | From interstellar dust to comets-A unification of observational constraints[END_REF][START_REF] Ehrenfreund | From interstellar material to cometary particles and molecules[END_REF]. The optical indices of such materials would present both the possibility for high absorption (carbonaceous compounds and dirty ices) and for low absorption of light (pure ice and some silicates). Light scattering by irregular solid particles, especially micrometer-sized fluffy aggregates of submicrometer grains made of such materials are typically consistent with the intensity and polarization phase curves observed from interplanetary dust [START_REF] Giese | Large fluffy particles-A possible explanation of the optical of interplanetary dust[END_REF][START_REF] Haudebourg | Theoretical polarimetric responses of fractal aggregates, in relation with experimental studies of dust in the solar system[END_REF][START_REF] Nakamura | Optical properties of fluffy aggregates as analogue of interplanetary dust particles[END_REF]. 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 [START_REF] Lasue | Inferring the interplanetary dust properties-from remote observations and simulations[END_REF][START_REF] Kolokolova | Effects of electromagnetic interaction in the polarization of light scattered by cometary and other types of cosmic dust[END_REF].

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 [START_REF] Lasue | Inferring the interplanetary dust properties-from remote observations and simulations[END_REF], that have been confirmed by laboratory simulations [START_REF] Hadamcik | Analogues of interplanetary dust particles to interpret the zodiacal light polarization[END_REF], 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 [START_REF] Mann | Dust near the Sun[END_REF].

As interplanetary dust particles spiral inwards and vanish in the Sun, with a typical lifetime of a couple of million years at 1 au [START_REF] Burns | Radiation forces on small particles in the solar system[END_REF], 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 [START_REF] Whipple | A comet model. I. The acceleration of Comet Encke[END_REF][START_REF] Dermott | Origin of the solar system dust bands discovered by IRAS[END_REF]. Other minor contributions come from the giant planets environment, the Kuiper Belt, and from interstellar dust particles (see e.g. [START_REF] Grün | Discovery of Jovian dust streams and interstellar grains by the Ulysses spacecraft[END_REF]. 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 [START_REF] Hahn | Clementine observations of the zodiacal light and the dust content of the inner solar system[END_REF]. 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% [START_REF] Nesvorný | Cometary origin of the zodiacal cloud and carbonaceous micrometeorites. Implications for hot debris disks[END_REF][START_REF] Nesvorný | Dynamical model for the zodiacal cloud and sporadic meteors[END_REF]. 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 [START_REF] Yang | Evolution of cometary dust particles to the orbit of the Earth: Particle size, shape, and mutual collisions[END_REF]. Consistent results were also obtained from the study of the Doppler shift from the Fraunhofer lines of the scattered solar light [START_REF] Ipatov | Dynamical zodiacal cloud models constrained by high resolution spectroscopy of the zodiacal light[END_REF]). An improved model of the IRAS and COBE thermal emission observations by [START_REF] Rowan-Robinson | An improved model for the infrared emission from the zodiacal dust cloud: cometary, asteroidal and interstellar dust[END_REF] 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 [START_REF] Carrillo-Sánchez | Sources of cosmic dust in the Earth's atmosphere[END_REF].

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 [START_REF] Jenniskens | Meteor showers and their parent comets[END_REF]. 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) [START_REF] Flynn | Composition of cosmic dust: sources and implications for the early solar system[END_REF][START_REF] Koschny | Interplanetary Dust, Meteoroids, Meteors and Meteorites[END_REF]. 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 [START_REF] Duprat | Extreme deuterium excesses in ultracarbonaceous micrometeorites from central[END_REF].

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 [START_REF] Flynn | Composition of cosmic dust: sources and implications for the early solar system[END_REF]. 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. [START_REF] Levasseur-Regourd | Cometary dust[END_REF]. Most notably, the presence of both fluffy aggregates of grains and compact particles have been evidenced by both Stardust and Rosetta space missions [START_REF] Hörz | Impact on Stardust: Implications for comet 81P/Wild 2 dust[END_REF][START_REF] Güttler | Synthesis of the morphological description of cometary dust at comet 67P/Churyumov-Gerasimenko[END_REF]. 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 [START_REF] Sandford | Organics Captured from Comet 81P/Wild 2 by the Stardust Spacecraft[END_REF][START_REF] Fray | High-molecular-weight organic matter in the particles of comet 67P/Churyumov-Gerasimenko[END_REF]. Finally, the scattered light phase curves present a similar behavior between ejected cometary dust and interplanetary dust particles analogues (Levasseur-Regourd et al., 2019).

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 [START_REF] Levasseur-Regourd | Linking studies of micrometeoroids, zodiacal dust, cometary dust and circumstellar disks[END_REF].

Conclusions and future observations

In this work, we attempted to give an overview of the properties of the zodiacal light and its link with the physical properties of the interplanetary dust particles that evolve in the Solar System and their sources. 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 [START_REF] Leinert | The 1997 reference of diffuse night sky brightness[END_REF][START_REF] Levasseur-Regourd & Lasue | Figure 7: Variation of the local near-ecliptic brightness (top) and linear degree of polarization, P[END_REF] 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 [START_REF] Abell | LSST science book, version 2.0[END_REF], 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 [START_REF] Valls-Gabaud | The MESSIER surveyor: unveiling the ultra-low surface brightness universe[END_REF][START_REF] Muslimov | Fast, wide-field and distortionfree telescope with curved detectors for surveys at ultralow surface brightness[END_REF]. 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 [START_REF] Gater | Comet mission given green light by European Space Agency[END_REF]. 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 [START_REF] Weinberg | Brightness and Polarization of the Zodiacal Light: Results of Fixed-Position Observations from Skylab[END_REF]; rocket Astro 7 (λe -λ = 25 • ) [START_REF] Pitz | Colour and polarization of the zodiacal light from the ultraviolet to the near infrared[END_REF]; Helios (λe -λ = 90 • ) [START_REF] Leinert | Helios zodiacal light measurements-a tabulated summary[END_REF]; COBE (λe -λ = 90 • ) [START_REF] Berriman | COBE DIRBE nearinfrared polarimetry of the zodiacal light: Initial results[END_REF]. 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.

Figure 1 :

 1 Figure 1: zodiacal light (ESO/Y.Beletsky)

Figure 3 :

 3 Figure 3: Brightness and polarization maps Meanwhile, several space mission observations were also obtained with the added advantage of not being contaminated by the atmospheric night-glow, but limited in viewing directions and mission time span. The characteristics of the dust were assessed in different zones of the Solar System thanks to their mobility. The D2A satellite around the Earth measured the zodiacal light in a plane perpendicular to the Sun-Earth direction. Its results confirmed the validity of ground-based measurements, and gave some evidence for the presence of interplanetary dust heterogeneities[START_REF] Levasseur | Satellite observations of intensity variations of the zodiacal light[END_REF]. The Helios 1 and

Fig. 3

 3 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 Table18published in[START_REF] Leinert | The 1997 reference of diffuse night sky brightness[END_REF]. The error in the mea-

Figure 5 :

 5 Figure 5: Spectral gradient for polarization

Figure 6 :

 6 Figure 6: Observation geometry for zodiacal light inversion.

Figure 7 :

 7 Figure 7: Local brightness and polarimetric data

Figure 8 :

 8 Figure 8: Local brightness and polarimetric data

  The zodiacal light observations provide us with an overall picture of the interplanetary dust environment: its mass flux distribution; chemical, material and structural properties; and the importance of different cometary and asteroidal sources in populating the inner Solar System with dust. This information has implications for the determination of the meteoroid flux deposited over the planetary surfaces.

Figure 1 :

 1 Figure 1: Zodiacal light seen from Paranal, Chile. The bright planet at the top center is Jupiter (ESO/Y.Beletsky).

Figure 2 :

 2 Figure 2: Zodiacal light geometry of observation in helio-ecliptic coordinates.

Figure 4 :

 4 Figure4: Most recent map of the zodiacal light brightness (10 -8 W m -2 sr -1 µm -1 ). Brightness levels are comparable to those presented in the upper part of Fig.3, while the resolution is significantly improved. The outer circle represents the ecliptic, the point at 0 • is the Sun, the central point is the ecliptic pole, and the near 180 • region is the Gegenschein. Based on data from[START_REF] Buffington | Measurements and an empirical model of the Zodiacal brightness as observed by the Solar Mass Ejection Imager (SMEI)[END_REF] and adapted from[START_REF] Levasseur-Regourd & Lasue | Figure 7: Variation of the local near-ecliptic brightness (top) and linear degree of polarization, P[END_REF].

Figure 5 :

 5 Figure 5: Wavelength variation of the linear polarization, P , in different directions of the sky. Skylab at the North Celestial Pole[START_REF] Weinberg | Brightness and Polarization of the Zodiacal Light: Results of Fixed-Position Observations from Skylab[END_REF]; rocket Astro 7 (λe -λ = 25 • )[START_REF] Pitz | Colour and polarization of the zodiacal light from the ultraviolet to the near infrared[END_REF]; Helios (λe -λ = 90 • )[START_REF] Leinert | Helios zodiacal light measurements-a tabulated summary[END_REF]; COBE (λe -λ = 90 • )[START_REF] Berriman | COBE DIRBE nearinfrared polarimetry of the zodiacal light: Initial results[END_REF].

Figure 6 :

 6 Figure 6: Geometry of observations for inversion of the zodiacal light, for lines of sight in the symmetry plane (near the ecliptic) of the zodiacal cloud (Figure adapted from Renardet 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.

Figure 8 :

 8 Figure 8: Phase curves of the local near-ecliptic brightness (top) and linear degree of polarization, P , (in percent bottom), at 1.5 au, adapted from Dumont & Sanchez (1975b); Fechtig et al. (1981); Levasseur-Regourd et al. (1990) and references within.