Primaries and multiples are imaged jointly using least-squares reverse time migration. Joint least-squares reverse time migration (LSRTM) uses the multiple energy as signal while removing the cross talk in the image. A class of multiply scattered events can be imaged by using the two-way propagator in modeling and migration. Such events can scatter off sharp interfaces in the migration velocity as well as the sea surface. The concept and methodology can be demonstrated with a 2D synthetic Sigsbee2B model and the Deimos ocean-bottom field data set. By including the multiples in imaging, the image quality is improved in areas close to complex salt structures.