Brenden,
If you use 3dMVM as shown in the Frontier paper, the effect estimates from the basis functions are treated as the levels of a repeated-measures factor, and by default they are assumed to have a spherical variance-covariance structure. However, if you are really concerned about the possibility of sphericity violation, you can turn on the option of sphericity correction: -SC.
Gang