afni.nimh.nih.gov

Try 3dNwarpFuncs. I'm no expert on VBM, but the -bulk option gives you the Jacobian information that I think you want. Warps that are composed similarly by other packages should work the same as long as they are dx,dy,dz volumes.