Hi Michael,
there re lots of possibilities there. One is this (as I did not have your functions I just made something up):
ImageCompose[GeoGraphics[GeoDisk[{0, 0}, Quantity[1000, "Miles"]]], ParametricPlot[{Cos[t], Sin[t]}, {t, 0, 2 Pi},
PlotStyle -> Directive[Black, Thick], Axes -> False, PlotRange -> {{-2, 2}, {-2, 2}}]]
Cheers,
Marco