Effects of large-scale heterogeneity and temporally varying hydrologic processes on estimating immobile pore space: A mesoscale-laboratory experimental and numerical modeling investigation

Facilitated transport of toluene and naphthalene with humic acid in high- and low-permeability systems: Role of ionic strength and cationic type

The occurrence of colloids on pollutants transport in groundwater has attracted more attention. However, the research on the regulation mechanism of colloids on combined pollutants transport in heterogeneous aquifers is limited. In this study, a series of tank experiments were conducted to systematically investigate the effects of ionic strength, and cation type on humic acid (HA) facilitated transport of toluene (TOL), and naphthalene (NAP) in high- and low-permeability systems. The results showed that HA facilitated pollutants transport in low Na+ solution. In Ca2+ solution, the presence of HA hindered pollutants transport, and the inhibition increased with the increase of ionic strength. Both in Na+ solution and low Ca2+ solution, the influence of heterogeneous structure on pollutant transport played a dominant role, and TOL and NAP had a greater transport potential in the high permeability zone (HPZ) due to the preferential flow. Whereas, deposition of HA aggregates, and electrostatic attractive interaction had negative effects on transport than groundwater flow in high Ca2+ solution. Pollutants were prone to accumulate at the bottom of the HPZ, and the top of the low permeability zone (LPZ). These new findings provide insights into the mechanism of colloids influence on the pollutants transport in heterogenous aquifer.

Evaluation of Dual Domain Mass Transfer in Porous Media at the Pore Scale

Dual-porosity models are often used to describe solute transport in heterogeneous media, but the parameters within these models (e.g., immobile porosity and mobile/immobile exchange rate coefficients) are difficult to identify experimentally or relate to measurable quantities. Here, we performed synthetic, pore-scale millifluidics simulations that coupled fluid flow, solute transport, and electrical resistivity (ER). A conductive-tracer test and the associated geoelectrical signatures were simulated for four flow rates in two distinct pore-scale model scenarios: one with intergranular porosity, and a second with an intragranular porosity also defined. With these models, we explore how the effective characteristic-length scale estimated from a best-fit dual-domain mass transfer (DDMT) model compares to geometric aspects of the flow field. In both model scenarios we find that: (1) mobile domains and immobile domains develop even in a system that is explicitly defined with one domain; (2) the ratio of immobile to mobile porosity is larger at faster flow rates as is the mass-transfer rate; and (3) a comparison of length scales associated with the mass-transfer rate (Lα ) and those associated with calculation of the Peclet number (LPe ) show LPe is commonly larger than Lα . These results suggest that estimated immobile porosities from a DDMT model are not only a function of physically mobile or immobile pore space, but also are a function of the average linear pore-water velocity and physical obstructions to flow, which can drive the development of immobile porosity even in single-porosity domains.