“…In this context, the variability of the hydraulic properties typically leads to preferential flow channels and lowpermeability areas where contaminants can be temporarily trapped by rate-limited mass transfer [e.g., Gomez-Hernandez and Wen, 1998;Zinn and Harvey, 2003;Bianchi et al, 2011]. The formation of these fast flow channels is typically associated with the presence of well-connected, highly permeable geological bodies or structures that can concentrate flow and solute transport [e.g., Knudby and Carrera, 2005; Incorporating hydrogeological uncertainty in human health predictions has been a topic of intense research in the past [e.g., Andričević and Cvetković, 1996;de Barros and Rubin, 2008;Cvetković and Molin, 2012;Rodak and Silliman, 2011;Andričević et al, 2012;Siirila and Maxwell, 2012;Atchley et al, 2013;de Barros and Fiori, 2014]. Probabilistic risk models allow one to determine the likelihood of risk exceeding a given regulatory threshold value [Tartakovsky, 2007], to delineate the spatial distribution of a plume for monitoring adaptation or intensification [James and Gorelick, 1994;Smalley et al, 2000;Maxwell et al, 2007;Fernandez-Garcia et al, 2012] and to better allocate characterization efforts to reduce the overall uncertainty of a given environmental performance metric [e.g., de Barros et al, 2009].…”