A method for determining displacement between and aligning a first seismic image and a second seismic image, comprising: receiving, via a processor, the first seismic image and the second seismic image, wherein the first seismic image and the second seismic image are representative of the same one or more geological formations within a subsurface region of the Earth at different times; characterized by generating, via the processor, a first scaled image comprising a smaller-size and lower-resolution compared to the first seismic image and a second scaled image comprising a smaller-size and lower-resolution compared to the second seismic image; obtaining, via the processor, a scaled displacement volume between the first scaled seismic image and the second scaled image, wherein obtaining the scaled displacement volume comprises obtaining a scaled displacement volume by comparing the first scaled image and the second scaled image and refining the scaled displacement volume by continuously comparing the first scaled image and a warped second scaled image produced by applying the scaled displacement volume to the second scaled image and updating the scaled displacement volume; determining, via the processor, a displacement volume based on the scaled displacement volume and a scaling function used to generate the first and second scaled images based on the first and second seismic images, wherein the displacement volume is an upscaled displacement volume, and wherein the displacement volume indicates changes in the subsurface region of the Earth between when the first seismic image was acquired and when the second seismic image was acquired; and generating, via the processor, a third seismic image by applying the displacement volume to the second seismic image, wherein one or more seismic events depicted in the third seismic image are aligned with one or more seismic events depicted in the first seismic image. receiving, via a processor, the first seismic image and the second seismic image, wherein the first seismic image and the second seismic image are representative of the same one or more geological formations within a subsurface region of the Earth at different times; characterized by generating, via the processor, a first scaled image comprising a smaller-size and lower-resolution compared to the first seismic image and a second scaled image comprising a smaller-size and lower-resolution compared to the second seismic image; obtaining, via the processor, a scaled displacement volume between the first scaled seismic image and the second scaled image, wherein obtaining the scaled displacement volume comprises obtaining a scaled displacement volume by comparing the first scaled image and the second scaled image and refining the scaled displacement volume by continuously comparing the first scaled image and a warped second scaled image produced by applying the scaled displacement volume to the second scaled image and updating the scaled displacement volume; determining, via the processor, a displacement volume based on the scaled displacement volume and a scaling function used to generate the first and second scaled images based on the first and second seismic images, wherein the displacement volume is an upscaled displacement volume, and wherein the displacement volume indicates changes in the subsurface region of the Earth between when the first seismic image was acquired and when the second seismic image was acquired; and generating, via the processor, a third seismic image by applying the displacement volume to the second seismic image, wherein one or more seismic events depicted in the third seismic image are aligned with one or more seismic events depicted in the first seismic image.