We present two-temperature, radiative general relativistic magnetohydrodynamic simulations of Magnetically Arrested Discs (MAD) that launch powerful relativistic jets. The mass accretion rates of our simulations are scaled to match the luminosity of the accretion flow around the supermassive black hole in M87. We consider two sub-grid prescriptions for electron heating: one based on a Landau-damped turbulent cascade, and the other based on heating from trans-relativistic magnetic reconnection. The simulations produce jets with power on the order of the observed value for M87. Both simulations produce spectra that are consistent with observations of M87 in the radio, millimetre, and submillimetre. Furthermore, the predicted image core-shifts in both models at frequencies between 15 and 86 GHz are consistent with observations. At 43 and 86 GHz, both simulations produce wide opening angle jets consistent with very long baseline interferometry images. Both models produce 230 GHz images with distinct black hole shadows that are resolvable by the Event Horizon Telescope (EHT), although at a viewing angle of 17°, the 230 GHz images are too large to match EHT observations from 2009 and 2012. The 230 GHz images from the simulations are dynamic on time-scales of months to years, suggesting that repeated EHT observations may be able to detect the motion of rotating magnetic fields at the event horizon.
