A system, comprising: a seismic imaging engine configured to create a seismic image of multiples and a seismic image of primaries based on seismic data, wherein the seismic data includes primaries and multiples, wherein the seismic image of multiples includes causal crosstalk and anti-causal crosstalk, and the seismic image of primaries includes causal crosstalk, and wherein: causal crosstalk comprises a false indication of an up-going wavefield being in phase with a down-going wavefield that corresponds to a location deeper than a location where two different wavefields are in phase; and anti-causal crosstalk comprises a false indication of the up-going wavefield being in phase with a down-going wavefield that corresponds to a location that is shallower than the location where the two different wavefields are in phase; and a prediction engine configured to: predict causal crosstalk (448, 649, 962) based on the seismic data by: migration of a down-going wavefield to a subsurface image point (236, 336, 536) as a source wavefield (237, 337, 437) of the seismic data; and migration of up-going multiples (346, 446, 546, 646) to the subsurface image point (236, 336, 536) as a receiver wavefield of the seismic data; and predict anti-causal crosstalk (552, 753, 1064) based on the seismic data by: migration of the down-going multiples to the subsurface image point (236, 336, 536) as a source wavefield (237, 337, 437) of the seismic data; and migration of up-going primaries to the subsurface image point (236, 336, 536) as a receiver wavefield of the seismic data; and an attenuation engine configured to attenuate the predicted causal crosstalk (448, 649, 962) and the predicted anti-causal crosstalk (552, 753, 1064) from the seismic image of multiples, and to attenuate the predicted causal crosstalk from the seismic image of primaries, wherein the seismic imaging engine is configured to combine the causal crosstalk and anti-causal crosstalk attenuated seismic image of the multiples with a causal crosstalk attenuated seismic image of primaries to obtain a crosstalk attenuated seismic image of the primaries and the multiples. a seismic imaging engine configured to create a seismic image of multiples and a seismic image of primaries based on seismic data, wherein the seismic data includes primaries and multiples, wherein the seismic image of multiples includes causal crosstalk and anti-causal crosstalk, and the seismic image of primaries includes causal crosstalk, and wherein: causal crosstalk comprises a false indication of an up-going wavefield being in phase with a down-going wavefield that corresponds to a location deeper than a location where two different wavefields are in phase; and anti-causal crosstalk comprises a false indication of the up-going wavefield being in phase with a down-going wavefield that corresponds to a location that is shallower than the location where the two different wavefields are in phase; and a prediction engine configured to: predict causal crosstalk (448, 649, 962) based on the seismic data by: migration of a down-going wavefield to a subsurface image point (236, 336, 536) as a source wavefield (237, 337, 437) of the seismic data; and migration of up-going multiples (346, 446, 546, 646) to the subsurface image point (236, 336, 536) as a receiver wavefield of the seismic data; and predict anti-causal crosstalk (552, 753, 1064) based on the seismic data by: migration of the down-going multiples to the subsurface image point (236, 336, 536) as a source wavefield (237, 337, 437) of the seismic data; and migration of up-going primaries to the subsurface image point (236, 336, 536) as a receiver wavefield of the seismic data; and predict causal crosstalk (448, 649, 962) based on the seismic data by: migration of a down-going wavefield to a subsurface image point (236, 336, 536) as a source wavefield (237, 337, 437) of the seismic data; and migration of up-going multiples (346, 446, 546, 646) to the subsurface image point (236, 336, 536) as a receiver wavefield of the seismic data; and migration of a down-going wavefield to a subsurface image point (236, 336, 536) as a source wavefield (237, 337, 437) of the seismic data; and migration of up-going multiples (346, 446, 546, 646) to the subsurface image point (236, 336, 536) as a receiver wavefield of the seismic data; and predict anti-causal crosstalk (552, 753, 1064) based on the seismic data by: migration of the down-going multiples to the subsurface image point (236, 336, 536) as a source wavefield (237, 337, 437) of the seismic data; and migration of up-going primaries to the subsurface image point (236, 336, 536) as a receiver wavefield of the seismic data; and migration of the down-going multiples to the subsurface image point (236, 336, 536) as a source wavefield (237, 337, 437) of the seismic data; and migration of up-going primaries to the subsurface image point (236, 336, 536) as a receiver wavefield of the seismic data; and an attenuation engine configured to attenuate the predicted causal crosstalk (448, 649, 962) and the predicted anti-causal crosstalk (552, 753, 1064) from the seismic image of multiples, and to attenuate the predicted causal crosstalk from the seismic image of primaries, wherein the seismic imaging engine is configured to combine the causal crosstalk and anti-causal crosstalk attenuated seismic image of the multiples with a causal crosstalk attenuated seismic image of primaries to obtain a crosstalk attenuated seismic image of the primaries and the multiples.