We present three epochs of VLBA observations of Galactic HI absorption toward the quasar 3C~138 with resolutions of 20 mas (~ 10 AU). This analysis includes VLBA data from observations in 1999 and 2002 along with a reexamination of 1995 VLBA data. Improved data reduction and imaging techniques have led to an order of magnitude improvement in sensitivity compared to previous work. With these new data we confirm the previously detected milliarcsecond scale spatial variations in the HI opacity at the level of Delta(tau_{max}) =0.50 \pm 0.05. The typical size scale of the optical depth variations is ~ 50 mas or 25 AU. In addition, for the first time we see clear evidence for temporal variations in the HI opacity over the seven year time span of our three epochs of data. We also attempted to detect the magnetic field strength in the HI gas using the Zeeman effect. From this analysis we have been able to place a 3 sigma upper limit on the magnetic field strength per pixel of ~45 muG. We have also been able to calculate for the first time the plane of sky covering fraction of the small scale HI gas of ~10%. We also find that the line widths of the milliarcsecond sizescale HI features are comparable to those determined from previous single dish measurements toward 3C~138, suggesting that the opacity variations cannot be due to changes in the HI spin temperature. From these results we favor a density enhancement interpretation for the small scale HI structures, although these enhancements appear to be of short duration and are unlikely to be in equilibrium.