“…However, the current understanding of this approach mainly comes from the analyses and applications of the first-order functional derivative (the gradient), for example, Brossier et al (2015); Li and Alkhalifah (2020); Ma and Hale (2013); Xu et al (2012Xu et al ( ), (2019. It has been observed that numerous iterations in the gradient-based RWI workflow are required to converge toward an appropriate background model (Wang et al, 2018;Xu et al, 2019), which makes it computationally demanding for real scale applications. Although the class of second-order optimization methods that exploits the importance of the inverse Hessian operator is well known in the numerical optimization community and has been widely applied to the FWI problems (e.g., Akcelik et al, 2003;Chen et al, 2007;Epanomeritakis et al, 2008;Fichtner & Trampert, 2011;Métivier et al, 2017;Operto et al, 2013;Pan et al, 2017;Pratt et al, 1998;Santosa & Symes, 1988), the applications of Hessian-based method in the context of RWI has still not been fully investigated.…”