Experiments on vertical transverse mixing in a large-scale heterogeneous model aquifer

Vertical transverse mixing is known to be a controlling factor in natural attenuation of extended biodegradable plumes originating from continuously emitting sources. We perform conservative and reactive tracer tests in a quasi two-dimensional 14 m long sand box in order to quantify vertical mixing in heterogeneous media. The filling mimics natural sediments including a distribution of different hydro-facies, made of different sand mixtures, and micro-structures within the sand lenses. We quantify the concentration distribution of the conservative tracer by the analysis of digital images taken at steady state during the tracer-dye experiment. Heterogeneity causes plume meandering, leading to distorted concentration profiles. Without knowledge about the velocity distribution, it is not possible to determine meaningful vertical dispersion coefficients from the concentration profiles. Using the stream-line pattern resulting from an inverse model of previous experiments in the sand box, we can correct for the plume meandering. The resulting vertical dispersion coefficient is approximately approximately 4 x 10(-)(9) m(2)/s. We observe no distinct increase in the vertical dispersion coefficient with increasing travel distance, indicating that heterogeneity has hardly any impact on vertical transverse mixing. In the reactive tracer test, we continuously inject an alkaline solution over a certain height into the domain that is occupied otherwise by an acidic solution. The outline of the alkaline plume is visualized by adding a pH indicator into both solutions. From the height and length of the reactive plume, we estimate a transverse dispersion coefficient of approximately 3 x 10(-)(9) m(2)/s. Overall, the vertical transverse dispersion coefficients are less than an order of magnitude larger than pore diffusion coefficients and hardly increase due to heterogeneity. Thus, we conclude for the assessment of natural attenuation that reactive plumes might become very large if they are controlled by vertical dispersive mixing.     

Choice of dispersion coefficients in reactive transport calculations on smoothed fields

Local dispersion dominates the mixing of compounds that are introduced separately into the subsurface and do not partition into any other than the aqueous phase. Thus, reactions between these compounds are controlled by dispersive mixing if they are not limited by kinetics. I quantify longitudinal dispersive mixing by the longitudinal effective dispersion coefficient of a conservative tracer introduced by a point-like injection [Water Resour. Res. 36 (12) (2000) 3591-3604]. In the upscaling of mixing-controlled reactive transport, I apply the mean velocity and the effective dispersion coefficient to the macroscopic transport calculations, whereas the reactive parameters on the macro-scale are identical to those on the local scale. The applicability of the approach is demonstrated for the transport of compounds undergoing a second-order irreversible bimolecular reaction. Ten realizations of a two-dimensional heterogeneous log-conductivity field are considered. Using the effective dispersion parameters, the overall mass balance is met in the ensemble average, whereas solute spreading is underestimated. I assess the lack of spreading by the difference between the expected macrodispersion and effective dispersion coefficients. I extend the approach to simulations on log-conductivity fields obtained by kriging of regularly spaced conductivity measurements. These fields contain the large-scale features of the true fields but do not resolve the small-scale variability. For the calculations on the kriged fields, the corresponding conditional covariance is substituted into the analytical expressions of effective dispersion, yielding a correction effective dispersion coefficient. The comparison between simulations on the fully resolved fields and on the kriged fields indicates that the approach is valid for wide plumes meeting the ergodicity condition. The high variability of mixing on small scales unresolved by kriging, however, leads to severe uncertainty when mixing-controlled reactions are predicted for narrow plumes.     

Model-based evaluation of controlled-release systems in the remediation of dissolved plumes in groundwater

Controlled-release, semi-passive reactive barrier systems have been recently developed as a long-term treatment option for controlling the spread of contaminant plumes in groundwater. This paper describes a new computer code, and applies it to study coupled processes of solute release, reaction, and mass transport in an in situ remediation scheme using the controlled release of potassium permanganate. Confidence with the modeling approach was developed by model verifications and simulating results of a pilot-scale test-cell experiment. Sensitivity analyses indicated the possibilities of treatment inefficiencies due to inability of transverse dispersion to mix the permanganate (MnO(4)(-)) within the zone of reaction, fluctuations in source strength due to variations in flow velocity, and the small length of treatment zone due to strong soil utilization of MnO(4)(-). Although problems associated with the fluctuating source strength and strong soil utilization can be addressed by optimizing the release rate, the inefficiency of transverse dispersion to create mixing could pose a serious limitation. Through a series of model simulations, a system of injection/withdrawal wells in a doublet arrangement was developed to facilitate lateral spreading and mixing of MnO(4)(-). A well-mixed, stable MnO(4)(-) zone with predetermined size (DxL=8m x 2m) and concentration ranges (1.5-20 mg l(-1)) was created by four 1-day injection/withdrawal pumping periods over 24 d. This type of mixing zone may persist for many years with periodic well mixing and replacements of exhausted controlled-release forms. Coupled use of the generalized code with field hydrologic data will help to optimize the design and operation of controlled-release systems in practice.     

Numerical experiments and field results on the size of steady state plumes

Contaminated groundwater poses a serious risk for drinking water supplies. Under certain conditions, however, groundwater contamination remains restricted to a tolerable extent because of natural attenuation processes. We present an innovative approach to evaluate the size of these so-called steady-state plumes by 2-D and 1-D modelling in homogeneous aquifers. If longitudinal mixing is negligible, scenarios can be modelled in a simplified way using a 1-D domain vertical to the direction of flow. We analysed the sensitivity of the plume length with respect to biodegradation kinetics, flow velocity, transverse vertical dispersivity alphat, the source and aquifer geometry and reaction stoichiometry. Our findings indicate that for many readily biodegradable compounds transverse-dispersive mixing rather than reaction kinetics is the limiting factor for natural attenuation. Therefore, if alphat, aquifer and source geometry and concentrations of electron acceptors and donors are known, the length of the steady state contaminant plume can be predicted. The approach is validated under field conditions for an ammonium plume at a former landfill site in SW Germany.     

A controlled field experiment on groundwater contamination by a multicomponent DNAPL: creation of the emplaced-source and overview of dissolved plume development

A unique field experiment has been undertaken at the CFB Borden research site to investigate the development of dissolved chlorinated solvent plumes from a residual dense non-aqueous phase liquid (DNAPL) source. The "emplaced-source" tracer test methodology involved a controlled emplacement of a block-shaped source of sand containing chlorinated solvents below the water table. The gradual dissolution of this residual DNAPL solvent source under natural aquifer conditions caused dissolved solvent plumes of trichloromethane (TCM), trichloroethene (TCE) and perchloroethene (PCE) to continuously develop down gradient. Source dissolution and 3-D plume development were successfully monitored via 173 multilevel samplers over a 475-day tracer test period prior to site remediation research being initiated. Detailed groundwater level and hydraulic conductivity data were collected. Development of plumes with concentrations spanning 1-700,000 micrograms/1 is described and key processes controlling their migration identified. Plumes were observed to be narrow due to the weakness of transverse dispersion processes and long due to advection and significant longitudinal dispersion, very limited sorptive retardation and negligible, if any, attenuation due to biodegradation or abiotic reaction. TCM was shown to be essentially conservative, TCE very nearly conservative and PCE, consistent with its greater hydrophobicity, more retarded yet having a greater mobility than observed in previous Borden field tests. The absence of biodegradation was ascribed to the prevailing aerobic conditions and lack of any additional biodegradable carbon substrates. The transient groundwater flow regime caused significant transverse lateral plume movement, plume asymmetry and was likely responsible for most of the, albeit limited, transverse horizontal plume spreading. In agreement with the widespread incidence of extensive TCE and PCE plumes throughout the industrialized world, the experiment indicates such solvent plumes are likely to be highly mobile and persistent, at least in aquifers that are aerobic and have low sorption potential (low foc content).     

Transverse dispersion of non-reactive tracers in porous media: a new nonlinear relationship to predict dispersion coefficients

Assessing the potential of natural attenuation in groundwater relies on the ability to predict and quantify the processes that occur in contaminant plumes. Transverse dispersion is a significant mass transfer mechanism for mixing of electron acceptors and donors and thus may control the lengths of steady state plumes. Laboratory experiments were carried out using a 2-dimensional acrylic glass tank filled with glass beads, quartz sand and field site material as porous media. Flow velocities and grain sizes were varied in order to cover a large range of Peclet numbers including typical field scenarios. The laboratory study was extended by a comprehensive literature search to compare the new results with earlier work. As a result we propose a new empirical relationship for prediction of transverse dispersion coefficients (Dt) which is based on the Peclet number (Pe). This new relationship indicates a nonlinear dependency on the flow velocity (nu a) and grain size (d), namely a relative decrease of the dispersion coefficient with increasing flow velocity in relatively fast flowing water: Dt/Daq=Dp/Daq+0.28(Pe)0.72 (with Pe=nu a d/Daq; Daq and Dp denote the aqueous and pore diffusion coefficients, resp.).     

Transverse vertical dispersion in groundwater and the capillary fringe

Transverse dispersion is the most relevant process in mass transfer of contaminants across the capillary fringe (both directions), dilution of contaminants, and mixing of electron acceptors and electron donors in biodegrading groundwater plumes. This paper gives an overview on literature values of transverse vertical dispersivities alpha(tv) measured at different flow velocities and compares them to results from well-controlled laboratory-tank experiments on mass transfer of trichloroethene (TCE) across the capillary fringe. The measured values of transverse vertical dispersion in the capillary fringe region were larger than in fully saturated media, which is credited to enhanced tortuosity of the flow paths due to entrapped air within the capillary fringe. In all cases, the values observed for alpha(tv) were < 1 mm. The new measurements and the literature values indicate that alpha(tv) apparently declines with increasing flow velocity. The latter is attributed to incomplete diffusive mixing at the pore scale (pore throats). A simple conceptual model, based on the mean square displacement and the pore size accounting for only partial diffusive mixing at increasing flow velocities, shows very good agreement with measured and published data.     

Estimation of local scale dispersion from local breakthrough curves during a tracer test in a heterogeneous aquifer: the Lagrangian approach

The local scale dispersion tensor, Dd, is a controlling parameter for the dilution of concentrations in a solute plume that is displaced by groundwater flow in a heterogeneous aquifer. In this paper, we estimate the local scale dispersion from time series or breakthrough curves, BTCs, of Br concentrations that were measured at several points in a fluvial aquifer during a natural gradient tracer test at Krauthausen. Locally measured BTCs were characterized by equivalent convection dispersion parameters: equivalent velocity, v(eq)(x) and expected equivalent dispersivity, [lambda(eq)(x)]. A Lagrangian framework was used to approximately predict these equivalent parameters in terms of the spatial covariance of log(e) transformed conductivity and the local scale dispersion coefficient. The approximate Lagrangian theory illustrates that [lambda(eq)(x)] increases with increasing travel distance and is much larger than the local scale dispersivity, lambda(d). A sensitivity analysis indicates that [lambda(eq)(x)] is predominantly determined by the transverse component of the local scale dispersion and by the correlation scale of the hydraulic conductivity in the transverse to flow direction whereas it is relatively insensitive to the longitudinal component of the local scale dispersion. By comparing predicted [lambda(eq)(x)] for a range of Dd values with [lambda(eq)(x)] obtained from locally measured BTCs, the transverse component of Dd, DdT, was estimated. The estimated transverse local scale dispersivity, lambda(dT) = DdT/U1 (U1 = mean advection velocity) is in the order of 10(1)-10(2) mm, which is relatively large but realistic for the fluvial gravel sediments at Krauthausen.     

Concentration field and turbulent fluxes during the mixing of two buoyant plumes

In this work an experimental study of mixing of two identical plumes, carried out in a turbulent neutral boundary layer generated in a wind tunnel, is presented. Measurements have been performed with fast flame ionisation detectors (FFIDs) and a two-component Laser-Doppler Anemometer system. Results allow the study of both the average and the fluctuating concentration field, including the turbulent vertical and longitudinal mass fluxes, in single plumes and during the interaction of two identical plumes. This information gives insight into the details of the mixing phase of the two plumes. Results of trajectories and additional rise (due to plume interactions) have been compared with previous measurements carried out in laminar cross-flows, showing similar behaviour. Concentration distributions in plume cross-sections in turbulent cross-flows differ from those measured in laminar cross-flows. Average vertical and longitudinal velocity measurements into the plume core show the strength of the shielding effect of the upwind plume and some details of interaction between the counter-rotating vortex pairs (CVPs). For large values of the alignment angle φ, between the line joining the stacks and the cross-flow, an average negative vertical velocity is measured in the middle of the plume even if its centre of mass is rising. This downward velocity is induced by the slow interaction of the CVPs and generates a vertical stretching of the plume and negative rise enhancement. Vertical turbulent fluxes change sign on the plume centreline and are of opposite sign with respect to the longitudinal turbulent fluxes. Results indicate a good linearity between vertical turbulent fluxes and concentration gradients, with different proportionality for the top and bottom parts of the plume (especially in the near field) indicating that dispersion could be described by a gradient-transfer model.     

Theoretical analysis of deviations from local equilibrium during sorbing solute transport through idealized stratified aquifers

This paper studies the spreading characteristics of reactive solute plumes in idealized stratified aquifers. The aquifer consists of two layers having different permeabilities with flow parallel to the stratification. The solute is assumed to adsorb onto the aquifer solids according to a first-order reversible kinetic rate law; the adsorption parameters are spatially uniform. We use the Aris moment method to examine analytically the time evolution of the lower-order spatial moments of the depth-averaged contaminant plume for an instantaneous input of mass. The results demonstrate that sorption kinetics cause the total dissolved mass and average velocity of the contaminant plume to decrease with increasing travel time. The plume variance is shown to depend upon three factors: intra-layer longitudinal dispersion, intra-layer kinetics, and vertical averaging. The results indicate that the relative importance of sorption kinetics diminishes as the permeability contrast between the layers increases. We present a simple criterion that can be used to assess the applicability of the local equilibrium assumption in idealized stratified systems.     

Sorption of naphthoic acids and quinoline compounds to estuarine sediment

The sorption of 16 ionizable organic compounds (IOCs) to an estuarine sediment was measured in synthetic estuarine water as a function of IOC concentration (1-100 microM) at fixed ionic strength (0.4 M), pH (7.6), and sediment concentration (0.018 g sediment kg(-1) suspension). Of the 16 IOCs, 11 were naphthoic acids and five were quinoline compounds. The linear sorption distribution coefficient (Kd) was used to correlate sorption to IOC physicochemical and molecular characteristics. With respect to naphthoic acid, sorption increased with the addition of ortho-substituent groups and with increasing chain length of the 1-acid group, and the greatest increase occurred with ortho-hydroxyl, carbonyl, and carboxyl groups. With respect to quinoline, sorption decreased with substituent group addition (except for nitro group) and with additional heterocyclic N atoms. For the naphthoic acids, log Kd exhibited a positive correlation with water solubility (log Sw) indicative of sorption primarily to mineral surfaces under the solution chemistry. For the quinoline compounds, log Kd exhibited a negative correlation with log Sw and a positive correlation with n-octanol/water partition coefficient (log K(OW)) indicative of sorption primarily to organic matter. For both compounds, poor or no correlations were established between log Kd and acid dissociation constant (pKa1), and between log Kd and a variety of molecular connectivity indexes. The results from this study demonstrate that the sorption of IOCs differ depending on their backbone structure and may differ between parent compound and ionizable degradation product.     

Characterization of metal adsorption kinetic properties in batch and fixed-bed reactors

Copper adsorption kinetic properties in batch and fixed-bed reactors were studied in this paper. The isothermal adsorption experiments showed that the copper adsorption capacity of a granular activated carbon (Filtrasorb 200) increased when ionic strength was higher. The presence of EDTA diminished the adsorption. An intraparticle diffusion model and a fixed-bed model were successfully used to describe the batch kinetic and fixed-bed operation behaviors. The kinetics became faster when the solution pH was not controlled, implying that the surface precipitation caused some metal uptake. The external mass transfer coefficient, the diffusivity and the dispersion coefficient were obtained from the modeling. It was found that both external mass transfer and dispersion coefficients increased when the flow rate was higher. Finally effects of kinetic parameters on simulation of fixed-bed operation were conducted.     

A controlled field experiment on groundwater contamination by a multicomponent DNAPL: dissolved-plume retardation

A natural gradient emplaced-source (ES) controlled field experiment was conducted at the Borden aquifer research site, Ontario, to study the transport of dissolved plumes emanating from residual dense nonaqueous-phase liquid (DNAPL) source zones. The specific objective of the work presented here is to determine the effects of solute and co-solute concentrations on sorption and retardation of dissolved chlorinated solvent-contaminant plumes. The ES field experiment comprised a controlled emplacement of a residual multicomponent DNAPL below the groundwater table and intensive monitoring of dissolved-phase plumes of trichloromethane (TCM), trichloroethylene (TCE), and perchloroethylene (PCE) plumes continuously generated in the aquifer down gradient from gradual source dissolution. Estimates of plume retardation (and dispersion) were obtained from 3-D numerical simulations that incorporated transient source input and flow regimes monitored during the test. PCE, the most retarded solute, surprisingly exhibited a retardation factor approximately 3 times lower than observed in a previous Borden tracer test by Mackay et al. [Water Resour. Res. 22 (1986) 2017] conducted approximately 150 m away. Also, an absence of temporal trend in PCE retardation contrasted with the previous Borden test. Supporting laboratory studies on ES site core indicated that sorption was nonlinear and competitive, i.e. reduced sorption of PCE was observed in the presence of TCE. Consideration of the effects of relatively high co-solute (TCE) concentration (competitive sorption) in addition to PCE concentration effects (nonlinear sorption) was necessary to yield laboratory-based PCE retardation estimates consistent with the field plume values. Concentration- and co-solute-based sorption and retardation analysis was also applied to the previous low-concentration pulse injection test of Mackay et al. [Water Resour. Res. 22 (1986) 2017] and was able to successfully predict the temporal field retardation trends observed in that test. While it is acknowledged that other "nonideal transport" effects may contribute, our analysis predicts differences in the PCE retardation magnitude and trend between the two experiments that are consistent with field observations based on the marked solute concentration differences that resulted from contrasting source conditions. Solute and co-solute concentration effects have heretofore received little attention, but may have wide significance in aquifers contaminated by point-source pollutants because many plumes contain mixed solutes over wide concentration ranges in strata that are likely subject to nonlinear sorption.     

Coagulation of highly turbid suspensions using magnesium hydroxide: effects of slow mixing conditions

Laboratory experiments were carried out to study the effects of slow mixing conditions on magnesium hydroxide floc size and strength and to determine the turbidity and total suspended solid (TSS) removal efficiencies during coagulation of highly turbid suspensions. A highly turbid kaolin clay suspension (1,213?±?36 nephelometric turbidity units (NTU)) was alkalized to pH 10.5 using a 5 M NaOH solution; liquid bittern (LB) equivalent to 536 mg/L of Mg²⁺ was added as a coagulant, and the suspension was then subjected to previously optimized fast mixing conditions of 100 rpm and 60 s. Slow mixing speed (20, 30, 40, and 50 rpm) and time (10, 20, and 30 min) were then varied, while the temperature was maintained at 20.7?±?1 °C. The standard practice for coagulation-flocculation jar test ASTM D2035-13 (2013) was followed in all experiments. Relative floc size was monitored using an optical measuring device, photometric dispersion analyzer (PDA 2000). Larger and more shear resistant flocs were obtained at 20 rpm for both 20- and 30-min slow mixing times; however, given the shorter duration for the former, the 20-min slow mixing time was considered to be more energy efficient. For slow mixing camp number (Gt) values in the range of 8,400–90,000, it was found that the mixing speed affected floc size and strength more than the time. Higher-turbidity removal efficiencies were achieved at 20 and 30 rpm, while TSS removal efficiency was higher for the 50-rpm slow mixing speed. Extended slow mixing time of 30 min yielded better turbidity and TSS removal efficiencies at the slower speeds.     

Applied dispersion modelling based on meteorological scaling parameters

A method for calculating the dispersion of plumes in the atmospheric boundary layer is presented. The method is easy to use on a routine basis. The inputs to the method are fundamental meteorological parameters, which act as distinct scaling parameters for the turbulence. The atmospheric boundary layer is divided into a number of regimes. For each scaling regime we suggest models for the dispersion in the vertical direction. The models directly give the crosswind-integrated concentrations at the ground, x_y, for nonbuoyant releases from a continuous point source. Generally the vertical concentration profile is proposed to be other than Gaussian. The lateral concentration profile is always assumed to be Gaussian, and models for determining the lateral spread σ_y are proposed. The method is limited to horizontally homogeneous conditions and travel distances less than 10km. The method is evaluated against independent tracer experiments over land. The overall agreement between measurements and predictions is very good and better than that found with the traditional Gaussian plume model.     

A new Lagrangian particle model for the simulation of dense gas dispersion

A Lagrangian stochastic model (MicroSpray), able to simulate the airborne dispersion in complex terrain and in presence of obstacles, was modified to simulate the dispersion of dense gas clouds. This is accomplished by taking into account the following processes: negative buoyancy, gravity spreading and the particle's reflection at the bottom computational boundary. Elevated and ground level sources, continuous and instantaneous emissions, time varying sources, plumes with initial momentum (horizontal, vertical or oblique in any direction), plumes without initial momentum are considered. MicroSpray is part of the model system MSS, which also includes the diagnostic MicroSwift model for the reconstruction of the 3-D wind field in presence of obstacles and orography. To evaluate the MSS ability to simulate the dispersion of heavy gases, its simulation performances are compared in detail to two field experiments (Thorney Island and Kit Fox) and to a chlorine railway accident (Macdona). Then, a comprehensive analysis considering several experiments of the Modelers Data Archive is presented. The statistical analysis on the overall available data reveals that the performance of the new MicroSpray version for dense gas releases is generally reliable. For instance, the agreement between concentration predictions and observations is within a factor of two in the 72% up to 99% of the occurrences for the case studies considered. The values of other performance measures, such as correlation coefficient, geometric mean bias and geometric variance, mostly set in the ranges indicated as good-model performances in the specialized literature.     

Reaction-enhanced mass transfer and transport from non-aqueous phase liquid source zones

The objective of the following research is to theoretically quantify the enhancement of interphase mass transfer of dissolved non-aqueous phase liquid (NAPL) compounds from the non-aqueous phase to the aqueous phase and the enhancement of dispersive mass transport from a NAPL zone due to destruction of dissolved NAPL compounds. For relatively slow reaction rates, such as for permanganate and perchloroethene (PCE), local-scale mass transfer enhancement is expected to be small. Dispersive mass transport with reaction from a horizontal NAPL zone can be quantified using equations derived for a mathematically equivalent falling film reactor system. In contrast to local-scale interphase mass transfer, dispersive mass transport from NAPL zones may be significantly increased by reaction. Enhancement factors due to destruction of the NAPL compound(s) are mainly dependent on NAPL solubility and oxidant concentration and to a lesser extent on reaction rate, stoichiometry, and transverse dispersion coefficients. Higher NAPL solubility and/or lower oxidant concentration reduces the maximum expected enhancement factor. Reaction enhancement factors for mass transport from NAPL zones are expected to be in the range of 5-50 for permanganate and chlorinated solvents. Theoretical results suggest that assuming instantaneous reaction rates may be appropriate for dispersive mass transport from NAPL zones.     

Variability of aquifer sorption properties in a field experiment on groundwater transport of organic solutes: Methods and preliminary results

Methods to characterize the organic solute sorption distribution coefficient, organic carbon content, and specific surface area of aquifer solids from the site of a field experiment on solute transport in groundwater were refined for application to small subsamples of 10-cm depth increments taken from 5-cm diameter cores. Initial results indicate that the average sorption characteristics of the Borden aquifer do not vary appreciably along the trajectory of the solute plumes. However, the sorption distribution coefficient of tetrachloroethylene varied over nearly an order of magnitude among 10-cm depth increments in one core sample. Preliminary evidence suggests that the sorption distribution coefficients for four halogenated organic solutes vary proportionally among core strata. However, the distribution coefficients for sorption of tetrachloroethylene on various depth increments are not well correlated with either organic carbon content or specific surface area, suggesting that as yet unidentified mineral phases may play a significant role in sorption of such solutes by the sandy aquifer solids.     

Modeling depth-variant and domain-specific sorption and biodegradation in dual-permeability media

A dual-permeability model (S_1D_DUAL) was developed to simulate the transport of land-applied pesticides in macroporous media. In this model, one flow domain was represented by the bulk matrix and the other by the preferential flow domain (PFD) where water and chemicals move at faster rates. The model assumed the validity of Darcian flow and the advective-dispersive solute transport in each of the two domains with inter-domain transfer of water and solutes due to pressure and concentration gradients. It was conceptualized that sorption and biodegradation rates vary with soil depth as well as in each of the two flow domains. In addition to equilibrium sorption, kinetic sorption was simulated in the PFD. Simulations were conducted to evaluate the combined effects of preferential flow, depth- and domain-variant sorption, and degradation on leaching of two pesticides: one with strong sorption potential (trifluralin) and the other with weak sorption potential (atrazine). Simulation results for a test case showed that water flux in the PFD was three times more than in the matrix for selected storm events. When equilibrium sorption was considered, the simulated profile of trifluralin in each domain was similar; however, the atrazine profile was deeper in the PFD than in the bulk matrix under episodic storm events. With an assumption of negligible sorption in the PFD, both the atrazine and the trifluralin profiles moved twice deeper into the PFD. The simulated concentrations of the chemicals were several orders higher in the PFD than in the matrix, even at deeper depths. The volume fraction of the macropores and the sorption and biodegradation properties of the chemicals could also affect the amount of pesticides leaving the root zone. For an intense storm event, slow sorption reaction rates in the PFD produced higher breakthrough concentrations of atrazine at the bottom of the simulated soil profile, thus posing the risk for breakthrough of chemicals from the root zone.     

Plume Dispersion in a Mountainous River Valley during Spring

The dispersion of hot plumes emitted from a smelter complex located In the Columbia River Valley, British Columbia, was evaluated under stable and neutral conditions during two mornings In spring. Spatial measurements of SO₂ and temperature within the plume were obtained by immersion probing using fast response helicopter and automobile mounted Instrumentation. In addition, meteorological measurements of vertical wind and temperature profiles at, and downwind from, the smelter were obtained from minisonde balloon releases. With weak down-valley winds, it was found that the plume axis elevations were generally lower during both stable and neutral conditions than would be predicted by Briggs plume-rise formulae. In contrast, plume dispersion, although confined in the horizontal by the steep valley walls during both stability regimes, was significantly enhanced by exceptionally good lateral mixing, particularly close to the source.