“…To interpret the heterogeneity pattern of aquifer‐aquitard systems, geostatistical and stochastic inverse methods were developed (Gutjahr & Wilson, 1989; Harvey & Gorelick, 1995; Hoeksema & Kitanidis, 1984, 1989; Kitanidis, 1995; Kitanidis & Vomvoris, 1983; Rubin & Dagan, 1987; Yeh et al., 1995, 1996), in which geostatistical parameters of aquifer systems were introduced into the inversion, and the most likely spatial distribution of hydraulic properties and their uncertainty are estimated for accurate prediction of drawdown responses. To date, HT has been developed as a mature technique to more accurately map spatial heterogeneity of hydraulic properties through synthetic studies (Cardiff & Barrash, 2011; Gottlieb & Dietrich, 1995; Luo et al., 2020; Tso et al., 2015; Yeh & Liu, 2000; Zhu & Yeh, 2005, 2006), laboratory experiments (Berg & Illman, 2011a; Illman et al., 2012, 2008, 2007, 2010; S. Liu et al., 2002; X. Liu et al., 2007; Luo et al., 2017; Zhao & Illman, 2021), field studies (Berg & Illman, 2011b, 2013, 2015; Bohling et al., 2007; Brauchler et al., 2013; Cardiff et al., 2013; Huang et al., 2011; Illman et al., 2009; Klepikova et al., 2020; Tiedeman & Barrash, 2020; Zha et al., 2016; Zhao & Illman, 2017, 2018, 2022) and basin‐scale aquifer characterization (Daranond et al., 2020; F. Liu, Yeh, Wang, Hao, et al., 2020; F. Liu, Yeh, Wang, Song, et al., 2020; Luo et al., 2022; Xia et al., 2023; Yeh et al., 2009).…”