  We use the Very Large Array (VLA) in the A configuration with the Pie Town
(PT) Very Long Baseline Array (VLBA) antenna to spatially resolve the extended
atmosphere of Betelgeuse over multiple epochs at 0.7, 1.3, 2.0, 3.5, and 6.1
cm. The extended atmosphere deviates from circular symmetry at all wavelengths
while at some epochs we find possible evidence for small pockets of gas
significantly cooler than the mean global temperature. We find no evidence for
the recently reported e-MERLIN radio hotspots in any of our multi-epoch VLA/PT
data, despite having sufficient spatial resolution and sensitivity at short
wavelengths, and conclude that these radio hotspots are most likely
interferometric artefacts. The mean gas temperature of the extended atmosphere
has a typical value of 3000 K at 2 $R_{\star}$ and decreases to 1800 K at 6
$R_{\star}$, in broad agreement with the findings of the single epoch study
from Lim et al. (1998). The overall temperature profile of the extended
atmosphere between $2 R_{\star} \lesssim r \lesssim 6 R_{\star}$ can be
described by a power law of the form $T_{\mathrm{gas}}(r) \propto r^{-0.6}$,
with temporal variability of a few 100 K evident at some epochs. Finally, we
present over 12 years of V band photometry, part of which overlaps our
multi-epoch radio data. We find a correlation between the fractional flux
density variability at V band with most radio wavelengths. This correlation is
likely due to shock waves induced by stellar pulsations, which heat the inner
atmosphere and ionize the more extended atmosphere through radiative means.
Stellar pulsations may play an important role in exciting Betelgeuse's extended
atmosphere.
