Solute transport in unsaturated porous media is of interest in many engineering and environmen- tal applications. The interplay between small-scale, local forces and the porous microstructure exerts a strong control on the transport of fluids and solutes at the larger, macroscopic scales. Heterogeneity in pore geometry is intrinsic to natural materials across a large range of scales. This multiscale nature, and the intricate links between two-phase flow and solute transport, remain far from well understood, by and large. Here, we use high-resolution direct simulation to quantify solute mixing and dispersion behavior within correlated porous media during drainage under an unfavorable viscosity ratio. Through analysis of flow and transport at multiple realizations, we find that increasing spatial correlations in pore sizes increase the size of the required Representative Elementary Volume (REV). We show that increasing the correlation length enhances solute dis- persivity through its impact on the spatial distribution of low-velocity (diffusion-dominated) and high-velocity (advection-dominated) regions. Fluid saturation is shown to directly affect diffusive mass flux among high- and low-velocity zones. Another indirect effect of correlated heterogene- ity on solute transport is through its control of the drainage patterns via repeated alteration in the connectivity of flowing pathways. Our findings improve quantitative understanding of solute mixing and dispersion under two-phase conditions, highly relevant to some of our most urgent environmental problems. 

Solute transport in partially-saturated porous media plays a key role in multiple applications across scales, from the migration of nutrients and contaminants in soils to geological energy storage and recovery. Our understanding of transport in unsaturated porous media remains limited compared to the well-studied saturated case. The focus of this review is the non-reactive transport driven by the displacement of immiscible fluids, where the fluid-fluid interface acts as a barrier that limits the solute to a single fluid phase. State-of-the-art pore-scale models are described, with a critical analysis of the gaps and challenges. A numerical example is provided to demonstrate the acute sensitivity of solute transport prediction to minute, inevitable uncertainties in the spatial distribution of the fluids' velocities and interface configuration associated with the multiphase flow modeling.