“…For instance, the over‐simplification of the 3D elastic domain as a homogeneous or layered half‐space, a conventional practice in coseismic slip inversion exercises, could result in over‐/under‐prediction of slip‐induced displacement when the approximated crust substantially differs from the actual environment. Thus, the resolved fault‐slip distribution and fit to the geodetic observations could be improved by incorporating more realistic lithological configurations into forward models linking fault dislocations to observed surface deformation (Hearn & Bürgmann, 2005; Kyriakopoulos et al., 2013; Tung, Fielding, et al., 2019; Tung, Katzenstein, et al., 2019; Tung & Masterlark, 2016, 2018b; 2018c; Tung, Masterlark, & Lo, 2018; Williams & Wallace, 2015). In addition to modeling the fault kinematics, the near‐fault elastic variabilities should also be considered for accurately calculating the geodetic slip moment (and moment magnitude) and Coulomb stress changes (Das et al., 2019; King et al., 1994; Langenbruch & Shapiro, 2014; Tung, Fielding, et al., 2019; Tung, Katzenstein, et al., 2019).…”