In this chapter, I propose an efficient multisource least-squares reverse time migration algorithm for a VTI medium. With blended sources processing, many conventionally acquired shot gathers are phase-encoded and blended together to form supergathers to reduce the computational cost and I/O burden of migration (Romero et al., 2000; Dai and Schuster, 2009,2010b; Krebs et al., 2009). However, blended sources processing introduces crosstalk noise, which needs to be removed from the final migration images. Our synthetic results demonstrate that LSRTM can mitigate both migration artifacts and crosstalk noise introduced by phase encoding, balance the amplitudes of reflectors, and improve the spatial resolution of the image. Moreover, the efficiency of multisource LSRTM can be significantly higher than conventional RTM depending on the number of shots encoded in one supergather, the number of migration operations at every iteration, and the number of iterations needed for an image of acceptable quality.