Spatially explicit simulations of gene flow within complex landscapes could help forecast the responses of populations to global and anthropological changes. Simulating how past climate change shaped intraspecific genetic variation can provide a validation of models in anticipation of their use to predict future changes. We review simulation models that provide inferences on population genetic structure. Existing simulation models generally integrate complex demographic and genetic processes but are less focused on the landscape dynamics. In contrast to previous approaches integrating detailed demographic and genetic processes and only secondarily landscape dynamics, we present a model based on parsimonious biological mechanisms combining habitat suitability and cellular processes, applicable to complex landscapes. The simulation model takes as input (a) the species dispersal capacities as the main biological parameter, (b) the species habitat suitability, and (c) the landscape structure, modulating dispersal. Our model emphasizes the role of landscape features and their temporal dynamics in generating genetic differentiation among populations within species. We illustrate our model on caribou/reindeer populations sampled across the entire species distribution range in the Northern Hemisphere. We show that simulations over the past 21 kyr predict a population genetic structure that matches empirical data. This approach looking at the impact of historical landscape dynamics on intraspecific structure can be used to forecast population structure under climate change scenarios and evaluate how species range shifts might induce erosion of genetic variation within species.