  The accreting millisecond pulsar XTE J1814-338 exhibits oscillations at the
known spin frequency during Type I X-ray bursts. The properties of the burst
oscillations reflect the nature of the thermal asymmetry on the stellar
surface. We present an analysis of the variability of the burst oscillations of
this source, focusing on three characteristics: fractional amplitude, harmonic
content and frequency. Fractional amplitude and harmonic content constrain the
size, shape and position of the emitting region, whilst variations in frequency
indicate motion of the emitting region on the neutron star surface. We examine
both long-term variability over the course of the outburst, and short-term
variability during the bursts. For most of the bursts, fractional amplitude is
consistent with that of the accretion pulsations, implying a low degree of fuel
spread. There is however a population of bursts whose fractional amplitudes are
substantially lower, implying a higher degree of fuel spread, possibly forced
by the explosive burning front of a precursor burst. For the first harmonic,
substantial differences between the burst and accretion pulsations suggest that
hotspot geometry is not the only mechanism giving rise to harmonic content in
the latter. Fractional amplitude variability during the bursts is low; we
cannot rule out the hypothesis that the fractional amplitude remains constant
for bursts that do not exhibit photospheric radius expansion (PRE). There are
no significant variations in frequency in any of the bursts except for the one
burst that exhibits PRE. This burst exhibits a highly significant but small
($\approx 0.1$Hz) drop in frequency in the burst rise. The timescale of the
frequency shift is slower than simple burning layer expansion models predict,
suggesting that other mechanisms may be at work.
