Modeling EDTA-facilitated cadmium migration in high- and low-permeability systems using MODFLOW and RT3D.

Cadmium (Cd) has impacted groundwater resources and can pose a serious threat to human health and the environment. Its fate in groundwater is complex and challenging to predict, as it is affected by adsorption to sediments, complexation with aqueous phase ligands, and variations in hydraulic conductivity. In this study, a 2D reactive transport model based on MODFLOW and RT3D is used to simulate published experimental results of cadmium migration without and with EDTA present in a flow cell containing high- and low-permeability zones (i.e., HPZs and LPZs). The model is then extended to conceptual flow cells with more complex LPZ configurations. Simulation results generally match the experimental data well, and analysis of experimental and simulated Cd effluent concentration profiles shows that EDTA enhances Cd removal from LPZs relative to water alone. Simulation results indicate that faster Cd removal is due to EDTA complexation with adsorbed Cd in LPZs, which enhances its solubilization and subsequent back diffusion. Lastly, simulation results show that with increasing LPZ heterogeneity more Cd is retained in flow cells, and EDTA is more effective in enhancing Cd removal relative to water alone; these results are attributed to more LPZ-HPZ interfaces that enhance Cd mass transfer into LPZs during contamination, and enhance EDTA mass transfer into LPZs to promote cleanup. Overall, the results highlight the promise of using EDTA to remove Cd from heterogeneous sites, but caution is advised due to model simplicity and lack of consideration of changes in solution pH, redox potential, or competing cations.

Using MODFLOW and RT3D to simulate diffusion and reaction without discretizing low permeability zones.

Low permeability zones (LPZs) are major sources of groundwater contamination after active remediation to remove pollutants in adjacent high permeability zones (HPZs). Slow back diffusion from LPZs to HPZs can extend management of polluted sites by decades. Numerical models are often used to simulate back diffusion, estimate cleanup times, and develop site management strategies. Sharp concentration gradients of pollutants are present at the interface between HPZs and LPZs, and hence accurate simulation requires fine grid sizes resulting in high computational burden. Since the MODFLOW family of codes is widely used in practice, we develop a new approach for modeling pollutant back diffusion using MODFLOW/RT3D that eliminates the need for fine discretization of the LPZ. Instead, the LPZ is treated as an impermeable region in MODFLOW, while in RT3D the LPZ is conceptualized as a series of immobile zones coupled with a mobile zone at the HPZ/LPZ interface. Finite volume discretization of diffusion and reaction within the LPZ is then modeled as mass transfer and reaction among several immobile species. This results in a simulation domain with significantly fewer grid cells compared to that required if all LPZs are discretized, providing potential for improved computational efficiency. Cases, including a layer of HPZ over an LPZ, a thin/thick lens of LPZ embedded in HPZ, and multiple lens of LPZs embedded in HPZ are tested by the new approach for tracer and reactive scenarios.

Transport and Retention of Concentrated Oil-in-Water Emulsions in Porous Media.

Oil-in-water emulsions are routinely used in subsurface remediation. In these applications, high oil loadings present a challenge to remedial design as mechanistic insights into transport and retention of concentrated emulsions is limited. Column experiments were designed to examine emulsion transport and retention over a range of input concentrations (1.3-23% wt). Droplet breakthrough and retention data from low concentration experiments were successfully described by existing particle transport models. These models, however, failed to capture droplet transport in more concentrated systems. At high oil fraction, breakthrough curves exhibited an early fall at the end of the emulsion pulse and extending tailing. Irrespective of input concentration, all retention profiles displayed hyper-exponential behavior. Here, we extended existing model formulations to include the additional mixing processes occurring when at high oil concentrations-with focus on the influence of deposited mass and viscous instabilities. The resulting model was parametrized with low concentration data and can successfully predict concentrated emulsion transport and retention. The role of retained mass and viscous instabilities on mixing conditions can also be applied more broadly to systems with temporal or spatially variant water saturation or when viscosity contrasts exist between fluids.