The work presented in this paper has been carried out with the aim of interpreting the data of a GNSS Reflectometer (GNSS-R) over land. The problem involves the analysis of bistatic scattering of the incoming signal collected around the specular direction. This requires to model the coherent component associated to the mean surface but at the same time the diffuse incoherent component due to roughness at wavelength scale. In presence of vegetation, both components will be affected, the former mainly because of the canopy attenuation and the latter for the combined effect of attenuation and volume scattering. The paper reviews the problem and presents the approach followed to develop a simulator of GNSS-R data over land, aiming to support potential applications of GNSS-R for soil moisture and biomass retrieval.