In this paper, we study the joint inversion method for acoustic and electromagnetic data based on the contrast source inversion method and the cross-gradient constraint. During the iteration process, gradients of the cost function and the cross-gradient regularization term are calculated separately and then coupled together. Compared with separate inversion using acoustic or electromagnetic data, joint inversion can improve the accuracy of reconstruction through the correlation between different measured data. Numerical results show good performance of the reconstruction.