In this paper we show that Sun-viewing images obtained by the Mars Science Laboratory (MSL) Navigation Cameras (Navcam) can be used for retrieving the dust optical depth and constrain the aerosol physical properties at Gale Crater by evaluating the sky brightness as a function of the scattering angle. We have used 65 Sun-pointing images covering a period of almost three Martian years, from MSL mission sol 21 to sol 1646 (MY 31 to 33). Radiometric calibration and geometric reduction were performed on MSL Navcam raw image data records to provide the observed sky radiance as a function of the scattering angle for the near-Sun region (scattering angle from 4{\deg} to 30{\deg}). These curves were fitted with a multiple scattering radiative transfer model for a plane-parallel Martian atmosphere model using the discrete ordinates method. Modelled sky brightness curves were generated as a function of two parameters: the aerosol particle size distribution effective radius and the dust column optical depth at the surface. A retrieval scheme was implemented for deriving the parameters that generated the best fitting curve under a least-square error criterion. The obtained results present a good agreement with previous work, showing the seasonal dependence of both dust column optical depth and the effectiveparticle radius.

Dust is the main driver of Mars' atmospheric variability. The determination of Martian dust aerosol properties is of high relevance for radiative modelling and calculating its weather forcing. In particular, the light scattering behaviour at intermediate and large scattering angles can provide valuable information regarding the airborne dust particle shape. The angular distribution of sky brightness observed by the Mars Science Laboratory engineering cameras (Navcam and Hazcam) is used here to characterise the atmospheric dust single scattering phase function and to constrain the shape of the particles. An iterative radiative transfer based retrieval method was implemented in order to determine the aerosol modelling parameters which best reproduce the observed sky radiance as a function of the scattering angle in the solar almucantar plane. The aerosol models considered in this study for retrieving dust radiative properties were an analytical three term Double Henyey-Greenstein phase function, T-matrix calculations for cylindrical particles with different diameter-to-length aspect ratios and experimental phase functions from laboratory measurements of several Martian dust analogue samples. Results of this study returned mean DHG phase function parameter values in line with Wolff et al. (2009). Although differences were observed during the low opacity aphelion season (lower forward scattering values, presence of a peak in the backward region) compared to the rest of the year, no clear evidences of seasonal behaviour or interannual variability were derived. The obtained average D/L aspect ratios for T-matrix calculated cylindrical particles were 0.70{\pm}0.20 and 1.90{\pm}0.20, and the best fitting Martian dust analogue corresponded to the basalt sample.