Abstract. The detection of Global Navigation Satellite System (GNSS) signals that are reflected off the surface, together with the reception of direct GNSS signals offers a unique opportunity to monitor water level variations over land and ocean. The time delay between the reception of the direct and the reflected signal gives access to the altitude of the receiver over the reflecting surface. The field of view of the receiver is highly dependent on both the orbits of the GNSS satellites and the configuration of the study site geometries. A simulator has been developed to determine the accurate location of the reflection points on the surface by modelling the trajectories of GNSS electromagnetic waves that are reflected on the surface of the Earth. Only the geometric problem have been considered using a specular reflection assumption. The orbit of the GNSS constellations satellite (mainly GPS, GLONASS and Galileo), and the position of a fixed receiver are used as input. Three different simulation modes are proposed depending on the choice of the Earth surface (local sphere or ellipsoid) and the consideration of topography likely to cause masking effects. Atmospheric delay effects derived from adaptive mapping functions are also taken into account. This simulator was developed to determine where the GNSS-R receivers should be located to monitor efficiently a given study area. In this study, two test sites were considered. The first one at the top of the Cordouan lighthouse (45°35'11'' N; 1°10'24'' W; 65 m) and the second one in the shore of the Geneva lake (46°24'30'' N; 6°43'6'' E, with a 50 m receiver height). This site is hidden by mountains in the South (altitude up to 2000 m), and overlooking the lake in the North (altitude of 370 m). For this second test site configuration, reflections occur until 560 m from the receiver. The geometric differences between the positions of the specular reflection points obtained considering the Earth as a sphere or as an ellipsoid were found to be on average 44 cm for satellites elevation angle greater than 10° and 1 m for satellite elevation angle between 5° and 10°. The simulations highlight the importance of the DEM integration: differences with and without integrating the DEM were found to be about 3.80 m with the minimum elevation angle equal to 5° and 1.4 m with the minimum elevation angle set to 10°. The correction of the tropospheric effects on the signal leads to geometric differences about 24 m maximum for a 50 m receiver height whereas the maximum is 43 cm for a 5 m receiver height. These errors deeply increase with the receiver height. By setting it to 300 m, the geometric errors reach 103 m for satellite elevation angle lower than 10°. The tests performed with the simulator presented in this paper highlight the importance of the choice of the Earth representation and also the non-negligible effect of the troposphere on the specular reflection points positions. Various outputs (time-varying reflection point coordinates, satellites positions and ground paths, wave trajectories, Fresnel first surfaces, etc.) are provided either as text or KML files for a convenient use.