EM propagation over and radiation from radially varying ground systems are numerically modeled using the compensation theorem integral formulation. The source is taken to be a vertical electric dipole (VED) situated either on or above an azimuthally symmetric plane ground system having an assumed surface impedance profile which can vary arbitrarily in the radial direction.