Faults are characterized as damaged material that accommodate localized deformation of rocks (Ben-Zion, 2008). The deformation of fault zone rocks is associated with earthquake generation and rupture process (Perrin et al., 2016;Thakur et al., 2020). The fault material with reduced seismic velocity and altered rheological properties can also amplify ground shaking and influence the migration of hydrocarbons and fluids (Caine et al., 1996;Spudich & Olsen, 2001). Thus, mapping the location and properties of faults is critical for understanding earthquake process and assessing seismic hazard. One common method of mapping faults is the observation of exhumed faults in the field (e.g., Collettini et al., 2009;Faulkner et al., 2003;Mitchell & Faulkner, 2009), which utilizes slices through the fault outcrops. Fault zone drilling projects can extend the examination of fault structure to greater depths and be used to monitor long-term changes in physical properties (e.g., Hickman et al., 2004;Hung et al., 2009). These methods provide precise measurements at single points of observation but require considerable labor and resources. Seismological methods can help develop a more complete picture of subsurface fault characteristics. Earthquake locations and focal mechanisms shed light on fault locations and structural complexities (Ross et al., 2017;Wang & Zhan, 2020). Seismic tomography can produce images of seismic velocity and attenuation near a fault zone (e.g., Allam et al., 2014;Liu et al., 2021;Wang et al., 2019). Fault zone trapped waves recorded by the sensors within the fault zones can be used to model fault zone geometries and properties in detail (e.g.