“…The data-domain LSRTM method seeks to iteratively update the subsurface reflectivity by minimizing the residual between the simulated data and observed data (Dai and Schuster, 2013;Zhang et al, 2015;Wang et al, 2017). It has been extended to elastic (Feng and Schuster, 2017;Ren et al, 2017;Gu et al, 2018), viscoacoustic (Dutta and Schuster, 2014;Sun et al, 2016;Chen et al, 2017; and anisotropic cases (Qu et al, 2017;Mu et al, 2020) due to its superiority in balancing amplitude, suppressing artifacts, and improving image resolution. However, limited by large computing costs (Dai and Schuster, 2013), the sensitivity to migration velocity (Tan and Huang, 2014;Li et al, 2017), and the mismatch of amplitudes (Zhang et al, 2015), conventional LSRTM is not extensively used in large-scale field data processing.…”