Evaporative mass transfer behavior of a complex immiscible liquid

A series of laboratory experiments was conducted with a multiple-component immiscible liquid, collected from the Picillo Farm Superfund Site in Rhode Island, to examine liquid-vapor mass-transfer behavior. The immiscible liquid, which comprises solvents, oils, pesticides, PCBs, paint sludges, explosives, and other compounds, was characterized using gas chromatography and gas chromatography/mass spectrometry to determine mole fractions of selected constituents. Batch experiments were conducted to evaluate equilibrium phase-partitioning behavior. Two sets of air-stripping column studies were conducted to examine the mass-transfer dynamics of five selected target compounds present in the immiscible-liquid mixture. One set of column experiments was designed to represent a system with free-phase immiscible liquid present; the other was designed to represent a system with a residual phase of immiscible liquid. Initial elution behavior of all target components generally appeared to be ideal for both systems, as the initial vapor-phase concentrations were similar to vapor-phase concentrations measured for the batch experiment and those estimated using Raoult's law (incorporating the immiscible-liquid composition data). Later-stage removal of 1,2-dichlorobenzene appeared to be rate limited for the columns containing free-phase immiscible liquid and no porous medium. Conversely, evaporative mass transfer appeared to be ideal throughout the experiment conducted with immiscible liquid distributed relatively uniformly as a residual phase within a sandy porous medium.     

The effect of local-scale physical heterogeneity and nonlinear,rate-limited sorption/desorption on contaminant transport in porous media

Nonideal transport of contaminants in porous media has often been observed in laboratory characterization studies. It has long been recognized that multiple processes associated with both physical and chemical factors can contribute to this nonideal transport behavior. To fully understand system behavior, it is important to determine the relative contributions of these multiple factors when conducting contaminant transport and fate studies. In this study, the relative contribution of physical-heterogeneity-related processes versus those of nonlinear, rate-limited sorption/desorption to the observed nonideal transport of trichloroethene in an undisturbed aquifer core was determined through a series of miscible-displacement experiments. The results of experiments conducted using the undisturbed core, collected from a Superfund site in Tucson, AZ, were compared to those obtained from experiments conducted using the same aquifer material packed homogeneously. The results indicate that both physical and chemical factors, specifically preferential flow and associated rate-limited diffusive mass-transfer and rate-limited sorption/desorption, respectively, contributed to the nonideal behavior observed for trichloroethene transport in the undisturbed core. A successful prediction of trichloroethene transport in the undisturbed core was made employing a mathematical model incorporating multiple sources of nonideal transport, using independently determined model parameters to account for the multiple factors contributing to the nonideal transport behavior. The simulation results indicate that local-scale physical heterogeneity controlled the nonideal transport behavior of trichloroethene in the undisturbed core, and that nonlinear, rate-limited sorption/desorption were of secondary importance.     

Nonideal transport of contaminants in heterogeneous porous media: 9 - impact of contact time on desorption and elution tailing

A series of miscible-displacement experiments was conducted to examine the impact of sorption contact time on desorption and elution of trichloroethene from a well-characterized soil. A large number of contact times were examined, spanning 1 h to 4 years (∼2 × 10⁶ h). Effluent trichloroethene concentrations were monitored over a range of greater than six orders of magnitude, allowing characterization of potential asymptotic tailing. The results of the column experiments showed that trichloroethene exhibited extensive elution tailing for all experiments. Each increase in contact time resulted in a successive increase in the extent of tailing. In total, the number of pore volumes of water flushing required to reach the analytical detection limit increased from approximately 1000 for the 1-h contact time to almost 9000 for the 4-year contact time. These results indicate that a contact time of less than 1 h produced a sorbed phase that is relatively resistant to desorption, and that a progressive increase in resistance to desorption occurred with increased contact time. A mathematical model incorporating nonlinear, rate-limited sorption/desorption described by a continuous-distribution reaction function was used to successfully simulate the measured data. The nonlinear sorption, the apparent rapid development of desorption resistance, and the progressive increase in resistance with increasing contact time are consistent with behavior anticipated for sorbate interactions with hard-carbon components of the soil.     

Source-zone characterization of a chlorinated-solvent contaminated Superfund site in Tucson, AZ

An extensive site-characterization project was conducted at a large chlorinated-solvent contaminated Superfund site in Tucson, AZ. The project consisted of several components, including traditional site-characterization activities, tracer tests, laboratory experiments conducted with core material collected from the site, and mathematical modeling. The primary focus of the work presented herein is the analysis of induced-gradient contaminant elution tests conducted in a source zone at the site, investigation of the potential occurrence of immiscible liquid in the saturated zone, characterization of the relationship between mass flux reduction and mass removal, and evaluation of the impact of source-zone management on site remediation. The results of the present study, along with those of prior work, indicate that immiscible liquid is likely present in the saturated zone at the site source zones. Extensive tailing and rebound was observed for the contaminant-elution tests, indicating nonideal transport and mass-transfer behavior. The elution data were analyzed with a source-zone-scale mathematical model, and the results indicated that nonideal immiscible-liquid dissolution was the primary cause of the observed behavior. The time-continuous relationship between mass flux reduction and mass removal associated with the plume-scale pump-and-treat operation exhibited an initial large drop in mass flux with minimal mass removed, followed by a period of minimal mass flux reduction and a second period of large reduction. This behavior reflects the impact of both source-zone and aqueous-plume mass removal dynamics. Ultimately, a greater than 90% reduction in mass flux was achieved for a mass removal of approximately 50%. The influence of source-zone management on site remediation was evaluated by conducting two predictive simulations, one for which the source zones were controlled and one for which they were not. A plume-scale model was used to simulate the composite contaminant concentrations associated with groundwater extracted with the pump-and-treat system, which were compared to measured data. The information generated from this study was used to enhance the site conceptual model, help optimize operation of the pump-and-treat system, and evaluate the utility of source-zone remediation.     

Prediction of mass-transfer coefficient for solute transport in porous media

Several previously reported laboratory studies related to transport of solutes through packed columns were utilized to develop predictive relationships for mass-transfer rate coefficient. The data were classified into two groups: those obtained under rate-limited mass transfer between mobile and immobile water regions (physical nonequilibrium conditions), and those derived from rate-limited mass transfer between instantaneous and slow sorption sites (sorption nonequilibrium conditions). The mass-transfer coefficient in all these studies was obtained by fitting breakthrough data to a transport model employing a first-order rate limitations with a "constant" mass-transfer coefficient, independent of flow conditions. This study demonstrated that the mass-transfer coefficient in these models is dependent on system parameters including pore-water velocity, length-scale, retardation coefficient, and particle or aggregate size. Predictive relationships were developed, through regression analysis, relating mass-transfer coefficient to residence time. The developed relationships adequately estimated previously reported field mass-transfer values. Successful simulations of field desorption data reported by Bahr [J. Contam. Hydrol. 4 (1989) 205] further demonstrate the potential applicability of the developed relationships.     

Impact of nonaqueous phase liquid (NAPL) source zone architecture on mass removal mechanisms in strongly layered heterogeneous porous media during soil vapor extraction

An existing multiphase flow simulator was modified in order to determine the effects of four mechanisms on NAPL mass removal in a strongly layered heterogeneous vadose zone during soil vapor extraction (SVE): a) NAPL flow, b) diffusion and dispersion from low permeability zones, c) slow desorption from sediment grains, and d) rate-limited dissolution of trapped NAPL. The impacts of water and NAPL saturation distribution, NAPL-type (i.e., free, residual, or trapped) distribution, and spatial heterogeneity of the permeability field on these mechanisms were evaluated. Two different initial source zone architectures (one with and one without trapped NAPL) were considered and these architectures were used to evaluate seven different SVE scenarios. For all runs, slow diffusion from low permeability zones that gas flow bypassed was a dominant factor for diminished SVE effectiveness at later times. This effect was more significant at high water saturation due to the decrease of gas-phase relative permeability. Transverse dispersion contributed to fast NAPL mass removal from the low permeability layer in both source zone architectures, but longitudinal dispersion did not affect overall mass removal time. Both slow desorption from sediment grains and rate-limited mass transfer from trapped NAPL only marginally affected removal times. However, mass transfer from trapped NAPL did affect mass removal at later time, as well as the NAPL distribution. NAPL flow from low to high permeability zones contributed to faster mass removal from the low permeability layer, and this effect increased when water infiltration was eliminated. These simulations indicate that if trapped NAPL exists in heterogeneous porous media, mass transfer can be improved by delivering gas directly to zones with trapped NAPL and by lowering the water content, which increases the gas relative permeability and changes trapped NAPL to free NAPL.     

Surfactant enhanced recovery of tetrachloroethylene from a porous medium containing low permeability lenses. 2. Numerical simulation

A numerical model of surfactant enhanced solubilization was developed and applied to the simulation of nonaqueous phase liquid recovery in two-dimensional heterogeneous laboratory sand tank systems. Model parameters were derived from independent, small-scale, batch and column experiments. These parameters included viscosity, density, solubilization capacity, surfactant sorption, interfacial tension, permeability, capillary retention functions, and interphase mass transfer correlations. Model predictive capability was assessed for the evaluation of the micellar solubilization of tetrachloroethylene (PCE) in the two-dimensional systems. Predicted effluent concentrations and mass recovery agreed reasonably well with measured values. Accurate prediction of enhanced solubilization behavior in the sand tanks was found to require the incorporation of pore-scale, system-dependent, interphase mass transfer limitations, including an explicit representation of specific interfacial contact area. Predicted effluent concentrations and mass recovery were also found to depend strongly upon the initial NAPL entrapment configuration. Numerical results collectively indicate that enhanced solubilization processes in heterogeneous, laboratory sand tank systems can be successfully simulated using independently measured soil parameters and column-measured mass transfer coefficients, provided that permeability and NAPL distributions are accurately known. This implies that the accuracy of model predictions at the field scale will be constrained by our ability to quantify soil heterogeneity and NAPL distribution.     

Modeling the effects of concentration history on the slow desorption of trichloroethene from a soil at 100% relative humidity

To determine the effects of concentration history on slow sorption processes, desorption kinetic profiles for trichloroethene (TCE) were measured for a soil at 100% relative humidity subject to different exposure concentrations and exposure times. Exposure concentrations ranged from 1% to 80% of the saturation vapor pressure (Ps) for TCE, and exposure times (i.e., time allowed for sorption before desorption begins) ranged from 1 to 96 days. A spherical diffusion model based on a gamma distribution of sorption rates and a gamma distribution of desorption rates was developed and applied to the data. At 80% P/Ps, the entire gamma distributions of sorption and desorption rates were available for TCE. In accordance with a micropore filling mechanism, the fraction of these distributions available for TCE sorption decreased with decreasing P/Ps. Experimental results are consistent with a micropore-filling mechanism, where the amount of slow desorbing mass decreased with decreasing exposure time, and the fraction of slow desorbing sites filled decreased with decreasing exposure concentration. Simulation results suggest that diffusion limits the rates that micropores fill, and that rates of sorption and desorption for soil contaminated at smaller values of P/Ps are, on average, less than those at larger values of P/Ps (i.e., slow desorption rates are a function of exposure concentration). Simulation results also suggest that the model adequately describes the effects of exposure concentration and exposure time on the rates of sorption and desorption, but not on the capacity of the slow sites for TCE. This work is important because contaminant concentrations in the subsurface vary in space and time, and the proposed model represents a new and mechanistically based approach to capture the effects of this heterogeneity on slow desorption.     

Analysis of the long-term behavior of solute transport with nonlinear equilibrium sorption using breakthrough curves and temporal moments

We analyzed the long-term behavior of breakthrough curves (BTCs) and temporal moments of a solute subjected to Freundlich equilibrium sorption (s = kc(n)). For one-dimensional transport in a homogeneous porous medium, we derived a power-law relation between travel time, tau, and solute displacement, chi, with the exponent being equal to the Freundlich n exponent. The mean solute velocity, derived from the first time moment, was found to change as tau(n-1). For n values larger than 0.66, the second time moment could be related to c chi(2/n), where c is a constant. An approach based on the use of a critical concentration was developed to estimate the presence of the asymptotic regime in the tail of the BTC. This approach was tested successfully using numerical case studies. One-dimensional numerical simulations with varying values of k, n and initial mass were run to verify the closed form analytical expressions for the large time behavior of temporal moments and the tailing part of breakthrough curves. Good agreement between the slope of the tailing part of log-log transformed BTCs and the predicted slope using asymptotic theory was found. Asymptotic theory in general underestimated the magnitude of the concentration in the tail. The quality of the estimated concentrations in the tail improved for small values of the dispersivity. Experimental BTCs of uranin and benazolin were analyzed in combination with sorption/desorption batch experiments using asymptotic theory. A good agreement between the value of n parameter derived from desorption experiment with benazolin and the value of the n parameter derived from the tail of the BTC was found.     

Model coupling intraparticle diffusion/sorption, nonlinear sorption, and biodegradation processes

Diffusion, sorption and biodegradation are key processes impacting the efficiency of natural attenuation. While each process has been studied individually, limited information exists on the kinetic coupling of these processes. In this paper, a model is presented that couples nonlinear and nonequilibrium sorption (intraparticle diffusion) with biodegradation kinetics. Initially, these processes are studied independently (i.e., intraparticle diffusion, nonlinear sorption and biodegradation), with appropriate parameters determined from these independent studies. Then, the coupled processes are studied, with an initial data set used to determine biodegradation constants that were subsequently used to successfully predict the behavior of a second data set. The validated model is then used to conduct a sensitivity analysis, which reveals conditions where biodegradation becomes desorption rate-limited. If the chemical is not pre-equilibrated with the soil prior to the onset of biodegradation, then fast sorption will reduce aqueous concentrations and thus biodegradation rates. Another sensitivity analysis demonstrates the importance of including nonlinear sorption in a coupled diffusion/sorption and biodegradation model. While predictions based on linear sorption isotherms agree well with solution concentrations, for the conditions evaluated this approach overestimates the percentage of contaminant biodegraded by as much as 50%. This research demonstrates that nonlinear sorption should be coupled with diffusion/sorption and biodegradation models in order to accurately predict bioremediation and natural attenuation processes. To our knowledge this study is unique in studying nonlinear sorption coupled with intraparticle diffusion and biodegradation kinetics with natural media.     

Sorption and transport of trichloroethylene in caliche soil

Sorption of TCE to the caliche soil exhibited linear isotherm at the high TCE concentrations (Co = 122-1300 mg L⁻¹) but Freundlich isotherm at the low concentration range (1-122 mg L⁻¹). Sorption strength of the carbonate fraction of the soil was about 100-fold lower than the sorption strength of soil organic matter (SOM) in the caliche soil, indicating weak affinity of TCE for the carbonate fraction of the soil. Desorption of TCE from the caliche soil was initially rapid (7.6 × 10⁻⁴ s⁻¹), then continued at a 100-fold slower rate (7.7 × 10⁻⁶ s⁻¹). Predominant calcium carbonate fraction of the soil (96%) was responsible for the fast desorption of TCE while the SOM fraction (0.97%) controlled the rate-limited desorption of TCE. Transport of TCE in the caliche soil was moderately retarded with respect to the water (R = 1.75-2.95). Flow interruption tests in the column experiments indicated that the rate-limited desorption of TCE controlled the non-ideal transport of TCE in the soil. Modeling studies showed that both linear and non-linear nonequilibrium transport models provided reasonably good match to the TCE breakthrough curves (r² = 0.95-0.98). Non-linear sorption had a negligible impact on both the breakthrough curve shape and the values of sorption kinetics parameters at the high TCE concentration (Co = 1300 mg L⁻¹). However, rate-limited sorption/desorption processes dominated at this concentration. For the low TCE concentration case (110 mg L⁻¹), in addition to the rate-limited sorption/desorption, contribution of the non-linear sorption to the values of sorption kinetics became fairly noticeable.     

A method for estimating distributions of mass transfer rate coefficients with application to purging and batch experiments

Mass transfer between aquifer material and groundwater is often modeled as first-order rate-limited sorption or diffusive exchange between mobile zones and immobile zones with idealized geometries. Recent improvements in experimental techniques and advances in our understanding of pore-scale heterogeneity demonstrate that two (or even a few) rate coefficients are insufficient in many cases. Here, we investigate a piece-wise linear model for a continuous distribution of rate coefficients, that has several advantages over previously used ‘statistical' distribution models (with functional form from gamma or lognormal PDF's): (1) distributions of arbitrary, even bimodal, shapes can be represented; (2) linear estimation methods can be applied to determine the distribution from experimental data; (3) the uncertainty in the distribution can be determined for each of its sections; and (4) the relationship between the time scales of available data and those of estimatable mass transfer processes can be investigated. A statistical model refinement algorithm is presented that reduces the number of parameters (sections of the piece-wise linear model) to the admissible minimum. We show that purging experiments allow estimation of a wider zone of the rate distribution than do batch experiments, and hence will provide predictions that are accurate over a wider range of time scales. Finally, in an application to TCE gas-purging desorption data, the piece-wise linear rate-distribution model has a higher probability of being adequate than those using a gamma or lognormal distribution or a single rate coefficient.     

Temperature effects on chlorinated-benzene sorption to hydrophobic surfaces

Sorption and desorption of chlorinated benzenes were investigated in a series of column experiments using porous silica that had phenyl groups bonded to the surface; the mass-fraction organic carbon was 0.016. Both sorption and desorption curves were asymmetrical, but they were mirror images of each other for most experiments, indicating good sorption reversibility. The resulting breakthrough curves were fit to an advection-dispersion mathematical model, with sorption as a first-order, reversible reaction. Significantly greater tailing in the chlorinated-benzene breakthrough curves versus the salt-tracer ones was evidence of slow sorption and desorption. ΔH° values for di-, tri- and tetra-chlorobenzene were 13–21 kJ mol₋₁, indicative of strong van der Waals binding. Despite these small values, slow desorption was attributed to slow binding and release rather than diffusion through the bonded organic phase. Desorption rates decreased in going from di- to tri- to tetra- to penta-chlorobenzene. This decrease was significantly more than the decrease in molecular-diffusion coefficients in the same series, suggesting a chemical rather than a physical rate control. There was less difference in sorption rates through the series, suggesting an inverse relation between partition coefficient and desorption rate. ΔG° values were −17 to −23 kJ mol⁻¹, giving TΔS° values of about 4 kJ mol⁻¹. Thus enthalpic contributions to sorption appear to be of greater importance than entropic contributions.     

Mass-removal and mass-flux-reduction behavior for idealized source zones with hydraulically poorly-accessible immiscible liquid

A series of flow-cell experiments was conducted to investigate aqueous dissolution and mass-removal behavior for systems wherein immiscible liquid was non-uniformly distributed in physically heterogeneous source zones. The study focused specifically on characterizing the relationship between mass flux reduction and mass removal for systems for which immiscible liquid is poorly accessible to flowing water. Two idealized scenarios were examined, one wherein immiscible liquid at residual saturation exists within a lower-permeability unit residing in a higher-permeability matrix, and one wherein immiscible liquid at higher saturation (a pool) exists within a higher-permeability unit adjacent to a lower-permeability unit. The results showed that significant reductions in mass flux occurred at relatively moderate mass-removal fractions for all systems. Conversely, minimal mass flux reduction occurred until a relatively large fraction of mass (>80%) was removed for the control experiment, which was designed to exhibit ideal mass removal. In general, mass flux reduction was observed to follow an approximately one-to-one relationship with mass removal. Two methods for estimating mass-flux-reduction/mass-removal behavior, one based on system-indicator parameters (ganglia-to-pool ratio) and the other a simple mass-removal function, were used to evaluate the measured data. The results of this study illustrate the impact of poorly accessible immiscible liquid on mass-removal and mass-flux processes, and the difficulties posed for estimating mass-flux-reduction/mass-removal behavior.     

Sorption of heterocyclic compounds from a complex mixture of coal-tar compounds on natural clayey till

The sorption and desorption of heterocyclic organic compounds in a complex multisolute system to a natural clayey till was investigated. The composition of the solutes reflect a simplified composition of an aqueous phase in contact with coal tar. Sorption was studied for two ratios (s:l) of clayey till (solid) to aqueous phase (liquid). The effect of the complex mixture of solutes on sorption of the four heterocyclic compounds: benzofuran, dibenzofuran, benzothiophene, and dibenzothiophene is evaluated by comparison with their sorption measured in single-solute systems. Sorption of the four compounds is affected by the complex mixture, with sorption decreases for all four compounds at high s:l ratio indicating competitive sorption. The effect on sorption of the individual compounds is not related to solubility or hydrophobicity of the compounds. Freundlich-type isotherms are observed for all compounds in the high s:l-ratio experiments, but for the most hydrophobic compounds isotherms are close to linear. The sorption of N-compounds and benzofuran is apparently influenced by cation exchange and dipole–dipole attraction to clay minerals. At high concentrations a dramatic increase in the sorption of the most strongly sorbing compounds is observed in the low s:l-ratio experiment. The dramatic increase in sorption appears to be a result of multimolecular layer sorption or condensation on surfaces in the clayey till at high surface density of organic compounds, and the data are fitted by BET (Brunauer, Emmet, and Teller) type 2 isotherms. The increase may or may not be induced by the presence of N-heterocyclic compounds sorbed by cation exchange and dipole–dipole attraction. The desorption of the compounds was studied for the low s:l ratio where multimolecular layer formation apparently had occurred. Partially irreversible sorption, hysteric Langmuir type desorption with isotherms of very high K_l coefficient, or behaviour reflecting dissolution of a condensed phase is observed.     

Influence of soil properties on heavy metal sequestration by biochar amendment: 2. Copper desorption isotherms

Contaminant desorption constrains the long-term effectiveness of remediation technologies, and is strongly influenced by dynamic non-equilibrium states of environmental and biological media. Information is currently lacking in the influence of biochar and activated carbon amendments on desorption of heavy metal contaminants from soil components. In this study, copper sorption-desorption isotherms were obtained for clay-rich, alkaline San Joaquin soil with significant heavy metal sorption capacity, and eroded, acidic Norfolk sandy loam soil having low capacity to retain copper. Acidic pecan shell-derived activated carbon and basic broiler litter biochar were employed in desorption experiments designed to address both leaching by rainfall and toxicity characteristics. For desorption in synthetic rain water, broiler litter biochar amendment diminished sorption-desorption hysteresis. In acetate buffer (pH 4.9), significant copper leaching was observed, unless acidic activated carbon (pH_pzc = 3.07) was present. Trends observed in soluble phosphorus and zinc concentrations for sorption and desorption equilibria suggested acid dissolution of particulate phases that can result in a concurrent release of copper and other sorbed elements. In contrast, sulfur and potassium became depleted as a result of supernatant replacements only when amended carbon (broiler litter biochar) or soil (San Joaquin) contained appreciable amounts. A positive correlation was observed between the equilibrium aluminum concentration and initial copper concentration in soils amended with acidic activated carbon but not basic biochar, suggesting the importance of cation exchange mechanism, while dissolution of aluminum oxides cannot be ruled out.     

Migration of Cs through quartz sand: experimental results and modeling approaches

We present results from experiments on the migration of ¹³⁷Cs through columns containing quartz sand. Times for ¹³⁷Cs movement through these columns and the quantity of ¹³⁷Cs adsorbed by the sand decreased as the ionic strength of the pore water increased from 0.002 to 0.1 m. The breakthrough curves were characterized by a slow approach towards steady-state concentrations as well as by long tails, indicating that ¹³⁷Cs adsorption to the sand grains was, at least in part, controlled by rate-limited reactions. Various formulations for solute mass transfer were tested for their ability to fit the experimental breakthrough curves. Based on a statistical analysis, a nonlinear, two-site model was identified as the most appropriate for describing the suite of experimental data. Variation in the model parameter that describes the rate of ¹³⁷Cs adsorption to the sand showed no consistent pattern with changes in ionic strength. In contrast, model parameters describing the sorption capacity of the sand grains and the fraction of kinetic sorption sites on the sand decreased with increasing ionic strength. The parameter describing the rate of ¹³⁷Cs desorption varied directly with changes in ionic strength.     

Transport of polar and non-polar volatile compounds in polystyrene foam and oriented strand board

Transport of hexanal and styrene in polystyrene foam (PSF) and oriented strand board (OSB) was characterized. A microbalance was used to measure sorption/desorption kinetics and equilibrium data. While styrene transport in PSF can be described by Fickian diffusion with a symmetrical and reversible sorption/desorption process, hexanal transport in both PSF and OSB exhibited significant hysteresis, with desorption being much slower than sorption. A porous media diffusion model that assumes instantaneous local equilibrium governed by a nonlinear Freundlich isotherm was found to explain the hysteresis in hexanal transport. A new nonlinear sorption and porous diffusion emissions model was, therefore, developed and partially validated using independent chamber data. The results were also compared to the more conventional linear Fickian-diffusion emissions model. While the linear emissions model predicts styrene emissions from PSF with reasonable accuracy, it substantially underestimates the rate of hexanal emissions from OSB. Although further research and more rigorous validation is needed, the new nonlinear emissions model holds promise for predicting emissions of polar VOCs such as hexanal from porous building materials.     

Temporal evolution of DNAPL source and contaminant flux distribution: impacts of source mass depletion

We investigated, using model simulations, the changes occurring in the distribution of dense non-aqueous phase liquid (DNAPL) mass (Sn) within the source zone during depletion through dissolution, and the resulting changes in the contaminant flux distribution (J) at the source control plane (CP). Two numerical codes (ISCO3D and T2VOC) were used to simulate selected scenarios of DNAPL dissolution and transport in three-dimensional, heterogeneous, spatially correlated, random permeability fields with emplaced sources. Data from the model simulations were interpreted based on population statistics (mean, standard deviation, coefficient of variation) and spatial statistics (centroid, second moments, variograms). The mean and standard deviation of the Sn and J distributions decreased with source mass depletion by dissolution. The decrease in mean and standard deviation was proportional for the J distribution resulting in a constant coefficient of variation (CV), while for the Sn distribution, the mean decreased faster than the standard deviation. The spatial distributions exhibited similar behavior as the population distribution, i.e., the CP flux distribution was more stable (defined by temporally constant second moments and range of variograms) than the Sn distribution. These observations appeared to be independent of the heterogeneity of the permeability (k) field (variance of the log permeability field=1 and 2.45), correlation structure (positive vs. negative correlation between the k and Sn domains) and the DNAPL dissolution model (equilibrium vs. rate-limited), for the cases studied. Analysis of data from a flux monitoring field study (Hill Air Force Base, Utah) at a DNAPL source CP before and after source remediation also revealed temporal invariance of the contaminant flux distribution. These modeling and field observations suggest that the temporal evolution of the contaminant flux distribution can be estimated if the initial distribution is known. However, the findings are preliminary and broader implications to sampling strategies for remediation performance assessment need to be evaluated in additional modeling and experimental studies.