Hapke proposed a convenient and widely used analytical model to describe the spectro-photometry of granular materials. Using a compilation of the published data, Hapke (2012, Icarus, 221, 1079-1083) recently studied the relationship of b and c for natural examples and proposed the hockey stick relation (excluding b>0.5 and c>0.5). For the moment, there is no theoretical explanation for this relationship. One goal of this article is to study a possible bias due to the retrieval method. We expand here an innovative Bayesian inversion method in order to study into detail the uncertainties of retrieved parameters. On Emission Phase Function (EPF) data, we demonstrate that the uncertainties of the retrieved parameters follow the same hockey stick relation, suggesting that this relation is due to the fact that b and c are coupled parameters in the Hapke model instead of a natural phenomena. Nevertheless, the data used in the Hapke (2012) compilation generally are full Bidirectional Reflectance Diffusion Function (BRDF) that are shown not to be subject to this artifact. Moreover, the Bayesian method is a good tool to test if the sampling geometry is sufficient to constrain the parameters (single scattering albedo, surface roughness, b, c, opposition effect). We performed sensitivity tests by mimicking various surface scattering properties and various single image-like/disk resolved image, EPF-like and BRDF-like geometric sampling conditions. The second goal of this article is to estimate the favorable geometric conditions for an accurate estimation of photometric parameters in order to provide new constraints for future observation campaigns and instrumentations.

The analysis of the surface texture from the particle (grain size, shape and internal structure) to its organization (surface roughness) provides information on the geological processes. CRISM multi-angular observations (varied emission angles) allow to characterize the surface scattering behavior which depends on the composition but also the material physical properties (e.g., grain size, shape, internal structure, the surface roughness). After an atmospheric correction by the Multi-angle Approach for Retrieval of the Surface Reflectance from CRISM Observations, the surface reflectances at different geometries are analyzed by inverting the Hapke photometric model depending on the single scattering albedo, the 2-term phase function, the macroscopic roughness and the 2-term opposition effects. Surface photometric maps are created to observe the spatial variations of surface scattering properties as a function of geological units at the CRISM spatial resolution (200m/pixel). An application at the Mars Exploration Rover (MER) landing sites located at Gusev Crater and Meridiani Planum where orbital and in situ observations are available, is presented. Complementary orbital observations (e.g. CRISM spectra, THermal EMission Imaging System, High Resolution Imaging Science Experiment images) are used for interpreting the estimated Hapke photometric parameters in terms of physical properties. The in situ observations are used as ground truth to validate the interpretations. Varied scattering properties are observed inside a CRISM observation (5x10km) suggesting that the surfaces are controlled by local geological processes (e.g. volcanic resurfacing, aeolian and impact processes) rather than regional or global. Consistent results with the in situ observations are observed thus validating the approach and the use of photometry for the characterization of Martian surface physical properties.

The present article proposes an approach to analyze the photometric properties of the surface materials from multi-angle observations acquired by the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) on-board the Mars Reconnaissance Orbiter. We estimate photometric parameters using Hapke model in a Bayesian inversion framework. This work also represents a validation of the atmospheric correction provided by the Multi-angle Approach for Retrieval of Surface Reflectance from CRISM Observations (MARS-ReCO) proposed in the companion article.The latter algorithm retrieves photometric curves of surface materials in reflectance units after removing the aerosol contribution. This validation is done by comparing the estimated photometric parameters to those obtained from in situ measurements by Panoramic Camera instrument at the Mars Exploration Rover (MER)-Spirit and MER-Opportunity landing sites. Consistent photometric parameters with those from in situ measurements are found, demonstrating that MARS-ReCO gives access to accurate surface reflectance. Moreover the assumption of a non-Lambertian surface as included in MARS-ReCO is shown to be significantly more precise to estimate surface photometric properties from space in comparison to methods based on a Lambertian surface assumption. In the future, the presented method will allow us to map from orbit the surface bidirectional reflectance and the related photometric parameters in order to characterize the Martian surface.