“…Inverse modeling of hydraulic tomography (HT) estimates spatially distributed fields of hydrogeological parameters, such as hydraulic conductivity, transmissivity and specific storage, using steady or time‐dependent hydraulic head measurements in sequential pumping tests (Butler Jr. et al., 1999; Cardiff et al., 2009; Gottlieb & Dietrich, 1995; Huang et al., 2022; Illman et al., 2008, 2010; Liu et al., 2013; Liu et al., 2014; Tosaka et al., 1993; Yeh & Lee, 2007; Yeh & Liu, 2000; Zha et al., 2018; Zhu & Yeh, 2005). The most widely used approach for solving HT inverse problems is geostatistical approach (GA), including quasilinear GA (Fienen et al., 2008; Kitanidis, 1995) and successive linear estimator (SLE) (Yeh et al., 1995) The major bottleneck of GA is that it requires iterative forward model simulations to evaluate the Jacobian matrix, which are computationally expensive for large‐scale, high‐dimensional models.…”