Joint total variation functional, its anisotropic generalization, and its application to joint inversion of seismic and EM data is discussed. Algorithm of 2-D joint inversion of seismic first-Arrival traveltimes and MT impedances is presented. In numerical experiment with synthetic data it is demonstrated, that the discussed algorithm improves both seismic and EM models in comparison with single-domain inversions.