Dual-domain solute transfer and transport processes: evaluation in batch and transport experiments

Soil macropore networks establish a dual-domain transport scenario in which water and solutes are preferentially channeled through soil macropores while slowly diffusing into and out of the bulk soil matrix. The influence of macropore networks on intra-ped solute diffusion and preferential transport in a soil typical of subsurface-drained croplands in the Midwestern United States was studied in batch- and column-scale experiments. In the batch diffusion studies with soil aggregates, the estimated diffusion radius (length) of the soil aggregates corresponded to the half-spacing of the aggregate fissures, suggesting that the intra-ped fissures reduced the diffusion impedance and preferentially allowed solutes to diffuse into the soil matrix. In the column-scale solute transport experiments, the average diffusion radius (estimated from HYDRUS-2D simulations and a first-order diffusive transfer term) was nearly double that of the batch-scale study. This increase may be attributed to a loss of pore continuity and a compounding of the small diffusion impedance through macropores at the larger scale. The column-scale solute transport experiments also suggest that two preferential networks exist in the soil. At and near soil saturation, a primary network of large macropores (possibly root channels and earthworm burrows) dominate advective transport, causing a high degree of physical and sorption nonequilibrium and simultaneous breakthrough of a nonreactive (bromide) and a reactive (alachlor) solute. As the saturation level decreases, the primary network drains, while transport through smaller macropores (possibly intra-ped features) continues, resulting in a reduced degree of nonequilibrium and separation in the breakthrough curves of bromide and alachlor.     

Effects of initial solute distribution on contaminant availability, desorption modeling, and subsurface remediation

Low permeability regions in which solute movement is governed by diffusion reduce the availability of pollutants for remediation and can function as long-term sources of groundwater contamination. The inherent difficulty in understanding mass transfer from these regions of sequestered contamination is further complicated by unknown solute distributions within the low-permeability regions (sequestering regions). When models are calibrated to reproduce temporal histories of solute release from a sequestering region (desorption), the fitted parameter values are used to infer the physical or chemical characteristics of the media; however, the calibrated parameters also reflect the case-specific initial conditions (i.e., the solute distribution within the sequestering region domain at the onset of desorption). This phenomenon is demonstrated using model simulations of solute diffusion from hypothetical solids with characteristics similar to those of the well studied Borden, Ontario aquifer system. Solute release from the solids is simulated using a batch diffusion model under different initial solute distributions within the solids. The results of these model simulations are used to calibrate parameters of a multiple first-order rate desorption model (MRM) to illustrate how the fitted MRM parameters increase or decrease depending on the initial "aging" of the solids. Further numerical simulations are conducted for a one-dimensional flow system under steady-state and variable-rate hydraulic flushing. These simulations show that although aging reduces desorptive mass flux during early stages of flushing, aged sites have greater desorptive mass flux (greater solute availability) than "freshly" contaminated media during the later stages of remediation. Overall, the results demonstrate why the physicochemical meaning of observed desorption rates cannot be accurately deduced without first understanding the initial solute distribution within the media.     

Hexavalent chromium reduction kinetics in rodent stomach contents

Reduction of hexavalent chromium (Cr(VI)) to trivalent chromium (Cr(III)) in the stomach prior to absorption is a well-recognized detoxification process thought to limit the toxicity of ingested Cr(VI). However, administration of high concentrations of Cr(VI) in drinking water cause mouse small intestinal tumors, and quantitative measures of Cr(VI) reduction rate and capacity for rodent stomach contents are needed for interspecies extrapolation using physiologically-based toxicokinetic (PBTK) models. Ex vivo studies using stomach contents of rats and mice were conducted to quantify Cr(VI) reduction rate and capacity for loading rates (1-400 mg Cr(VI) L⁻¹ stomach contents) in the range of recent bioassays. Cr(VI) reduction was measured with speciated isotope dilution mass spectrometry to quantify dynamic Cr(VI) and Cr(III) concentrations in stomach contents at select time points over 1 h. Cr(VI) reduction followed mixed second-order kinetics, dependent upon concentrations of both Cr(VI) and the native reducing agents. Approximately 16 mg Cr(VI)-equivalents of reducing capacity per L of fed stomach contents (containing gastric secretions, saliva, water and food) was found for both species. The second-order rate constants were 0.2 and 0.3 L mg⁻¹ h⁻¹ for mice and rats, respectively. These findings support that, at the doses that caused cancer in the mouse small intestine (?20 mg Cr(VI) L⁻¹ in drinking water), the reducing capacity of stomach contents was likely exceeded. Thus, for extrapolation of target tissue dose in risk assessment, PBTK models are necessary to account for competing kinetic rates including second order capacity-limited reduction of Cr(VI) to Cr(III).     

Modeling and interpreting bioavailability of organic contaminant mixtures in subsurface environments

Bioavailability often controls the fate of organic contaminants in surface and subsurface aquatic environments. Bioavailability can be limited by sorption, mass transfer, and intrinsic biodegradation potential and can be further altered by the presence of other compounds. This paper reviews current perspectives on the processes influencing subsurface contaminant bioavailability, how these processes are modeled, and how the relative role of the various processes can be assessed through bioavailability indices. Although these processes are increasingly well understood, the use of sophisticated models and indices often are precluded by an inability to estimate the many parameters that are associated with complex models. Nonetheless, the proper representation of sorption, mass transfer, biodegradation, and co-solute effects can be critical in predicting bio-attenuation. The influence of these processes on contaminant fate is illustrated with numerical simulations for the simultaneous degradation of toluene (growth substrate) and trichloroethylene (nongrowth cometabolite) in hypothetical, aerobic, solid-water systems. The results show how the relative impacts on contaminant fate of the model's various component processes depends upon system conditions, including co-solute concentrations. Slow biodegradation rates increase the inhibition effects of a cometabolite and suppress the rate enhancement effects of a growth substrate. Irrespective of co-solute effects, contaminant fate is less sensitive to biodegradation processes in systems with strong sorption and slow desorption rates. Bioavailability indices can be used to relate these findings and to help identify appropriate modeling simplifications. In general, however, there remains a need to redefine such indices in order that bioavailability concepts can be better incorporated into site characterization, remediation design, and regulatory oversight.