The mean power of the reflected Global Navigation Satellite System (GNSS) signals acquired by a GNSS-Reflectometry (GNSS-R) receiver can be modeled through the integral bistatic radar equation by weighting the contributions of all scatterers on the surface by the system impulse response. The geophysical properties of the scattering surface affect the magnitude of the reflected navigation signals through the bistatic scattering coefficient which, in case the observed surface is land, is a function of the soil dielectric properties, surface roughness, and vegetation cover. In this paper, the GNSS-R signal simulator developed in the framework of the Land MOnitoring with Navigation signal (LEiMON) Project, supported by European Space Agency, is presented. The simulator is able to predict the power reflected by land, taking as input the system and observation parameters, as well as the land surface parameters. The latter are used to simulate both the coherent and the incoherent scattering, taking advantage of widely used theoretical models of bistatic scattering from bare soils and vegetated surfaces. First, the geometrical formulation is discussed, and then, the problem of polarization mismatch due to real antennas at circular polarization is faced following the polarization synthesis approach. Finally, a comparison with some experimental data collected during the LEiMON campaign is presented. The simulations display the same trend of the experimental data, thus showing that the simulator can be used as an efficient tool for the interpretation of GNSS-R measurements. © 1980-2012 IEEE.