Vlasiator is a kinetic space plasma model that simulates the behavior of plasma, solar wind and magnetic fields in near-Earth space. So far, these simulations have been run without any interaction with the ionosphere, the uppermost layer of Earth's atmosphere. In this paper, we present the new methods that add an ionospheric electrodynamics model to Vlasiator, coupling it with the existing methods and presenting new simulation results of how space plasma and Earth's ionosphere interact.