Nowadays, the space missions employing the global navigation satellite system reflectometry (GNSS-R) are UK TDS-1, NASA CYGNSS, and the Chinese BuFeng-1A/B twin satellites, part of the first Chinese global navigation satellite system reflectometry (GNSS-R) satellite mission. They provide delay-Doppler map (DDM) measurements reflected from the land as well as the ocean. Using land reflected DDMs, several studies have been conducted to retrieve land geophysical parameters, such as soil moisture and biomass. Despite the clear dependence on these parameters, many other parameters impact the DDMs as well, such as topography, surface roughness, etc. The impact of these perturbing factors must be analyzed, and modeled in various conditions. This article presents a comprehensive end-to-end simulator that can generate synthetic DDMs reflected over land. It is an extension of a previously developed simulator validated for ocean applications. This simulator is very generic, and it includes numerous configurable parameters such as arbitrary scattering geometry (transmitter and receiver positions and speeds), arbitrary GPS and Galileo transmitted signals and frequencies, GNSS-R instrument antenna and receiver errors, as well as surface topography, roughness, soil moisture, vegetation cover, etc. The sensitivities of GNSS-R observables with respect to soil moisture and vegetation are obtained and compared to previous experimental results, and synthetic DDMs are compared and validated against TDS-1 ones.