Electromagnetic wave propagation in the Earth-ionosphere cavity presents an interesting challenge for simulations. Three-dimensional latitude-longitude finite-difference time-domain (FDTD) models accounting for the bathymetry, topography and ionosphere have been developed and applied towards a number of applications previously. However, to date most of these models treat the ionosphere as a simple, isotropic exponential conductivity profile. Only recently has a latitude-longitude FDTD model been developed that treats the ionosphere as a magnetized cold plasma. This opens the door to modeling electromagnetic phenomena at higher frequencies and higher altitudes by accommodating more physics. Further, a geodesic (hexagonal-pentagonal) FDTD model that is more efficient, is easier to implement, and executes faster than latitude-longitude models has been recently developed. In this thesis, the magnetized cold plasma global latitude-longitude algorithm is adapted and implemented for the first time in a geodesic FDTD model of the Earth-ionosphere cavity