Reverse time migration (RTM) exploits time-reversal symmetry and adjoint duality to focus wavefields and reconstruct subsurface reflectivity, but large surveys remain limited by the cost of forward and backward propagation. We present a Graphics Processing Unit (GPU)-accelerated adjoint RTM workflow for depth imaging of the Chicxulub impact structure using the marine A0/A1 composite profile (1996). The processed stacked section contains 14,172 traces with 6.25 m Common Depth Point (CDP) spacing, 1 ms sampling, and 18 s record length. Forward and adjoint wavefields are computed with a staggered-grid finite-difference scheme (fourth order in space, second in time) and Convolutional Perfectly Matched Layers (CPMLs), which provide stable finite-domain simulations while introducing controlled symmetry breaking through absorption. The solver is verified with the Lamb half-space analytical benchmark and applied through five interpretation-guided velocity/density updates. The final depth image improves reflector continuity and interpretability of crater-scale elements, including post-impact sedimentary fill, melt and breccia units, terrace fault blocks, and deep uplift-related structure. Compute Unified Device Architecture (CUDA) acceleration reduces runtime from ∼32.36 h on a CPU baseline to ∼34.10 min on an RTX 3070 (≈56.9×), enabling practical, reproducible iterative RTM on accessible hardware.
Herrera-Pérez et al. (Wed,) studied this question.