We proved the wellposedness of variable-order mobile-immobile time-fractional diffusion equations and the regularity of their solutions. Optimal-order finite element approximation was presented and analyzed. Numerical experiments were carried out for demonstration.