Climate change has created new and evolving environmental conditions, impacting all species, including hosts and parasitoids. I therefore present integrodifference equation (IDE) models of host--parasitoid systems to model population dynamics in the context of climate-driven shifts in habitats. I describe and analyze two IDE models of host--parasitoid systems to determine criteria for coexistence of the host and parasitoid. Specifically, I determine the critical habitat speed, beyond which the parasitoid cannot survive. By comparing the results from two IDE models, I investigate the impacts of assumptions that reduce the system to a single-species model. I also compare critical speeds predicted by a spatially-implicit difference-equation model with critical speeds determined from numerical simulations of the IDE system. The spatially-implicit model uses approximations for the dominant eigenvalue of an integral operator. The classic methods to approximate the dominant eigenvalue for IDE systems do not perform well for asymmetric kernels, including those that are present in shifting-habitat IDE models. Therefore, I compare several methods for approximating dominant eigenvalues and ultimately conclude that geometric symmetrization and iterated geometric symmetrization give the best estimates of the parasitoid critical speed.