“…Various hypotheses have been invoked to explain possible factors that affect the migration and distribution of solutes under unsaturated, heterogeneous conditions, including turbulent mixing due to high rainfall (Havis et al, 1992); solute transfer between mobile and immobile water (De Smedt and Wierenga, 1984); mobile-immobile exchange and hysteresis (Butters et al, 1989;Russo et al, 1989aRusso et al, , b, 2014; lateral mixing due to velocity fluctuations (Russo et al, 1998); isotope effects (Barnes and Allison, 1988;LaBolle et al, 2008;Zhang et al, 2009);variable, state-dependent anisotropy (McCord et al, 1991); non-Gaussian early-time mean tracer plume behavior (Naff, 1990); non-Fickian solute migration at low water contents (Padilla et al, 1999) and for macroscopically homogeneous sand (Bromly and Hinz, 2004); and saturation-dependent dispersivity (Raoof and Hassanizadeh, 2013). In addition, Konikow et al (1997) and Parker and van Genuchten (1984) discuss the importance of boundary condition treatment (e.g., water-solute injection, solute exchange between soil and atmosphere).…”