Minnaert and two-term phase function Hapke models are used to describe the photometric properties of the martian surface using HRSC (High Resolution Stereo Camera) multi-angular observations acquired along the ongoing Mars-Express mission. These observations can be pieced together to derive integrated phase functions over a wide range of phase angles. The photometric diversity at 675 nm, as seen from orbit, of the martian surface properties across Gusev is depicted with seven units. Three photometric units widespread across the flanks of Apollinaris Patera flank and the floor of Gusev Crater are identified as having high single scattering albedo with rather forward scattering properties, low to intermediate macroscopic roughness and porous or not compacted powdered surface state as indicated by the opposition parameters. Another unit has the highest single scattering albedo, the smoothest surface in terms of macroscopic roughness, associated with an extremely forward scattering behavior. The opposition parameters are consistent with the presence of transparent particles in the surface powder layer. The distribution of this unit appears quite intermittent across the crater and does not seem to indicate any relationship with a given morphological structure. It may correspond to sparse areas where the structure of the surface dust layer is the most preserved. The most pronounced photometric changes are observed in three units associated with the low-albedo features corresponding to dark wind streaks. These units have a low single scattering albedo, are the most backscattering surfaces across Gusev, have a high surface roughness and present variable surface states as shown by the opposition parameters estimates, consistent with the occurrence of large grains organized in more or less packed layers. Clear differences are seen among these units in terms of opposition effect. While one exhibits typical characteristics for the opposition effect, another appears more unusual in terms of lobe width and the third suggests the occurrence of a packed/compressed/narrow size distribution powder surface. The opposition effect thus appears to play a significant role suggesting that the surface state optical properties across Gusev are strongly influenced by the porosity and packing characteristics or grain size distribution of the upper layer of the martian soil. The mapping aspect of the investigation is quite useful to get a better sense of the meaning of the observed photometric variations. Indeed, the Hapke modeling suggests that surface organization (surface roughness, packing state) is more important than the simple physical characterization of the intrinsic optical properties of the constitutive particles. Given the overall spatial patterns derived from the photometric analysis, the variations, at least for the western and central part of Gusev Crater, are likely partly driven by the prevailing wind regimes, considered to be oriented north-northwest/south-southeast and disturbing the very upper surface layer. The present photometric results agree with independent investigations based on thermal inertia, reflectance spectroscopy, in situ photometric and microscopic imaging and support the idea of a thin layer of fine-grained dust, being stripped off in the low albedo units to reveal a dark basaltic substrate comprising coarse-grained materials.