Macro-scale effective constitutive relationships for two-phase flow processes in heterogeneous porous media with emphasis on the relative permeability-saturation relationship

Macro-scale simulations often play an important role in the assessment and remediation of contamination by dense non-aqueous phase liquids (DNAPLs) in the subsurface. Effective parameters for the macro scale are required for these simulations in order to avoid a detailed discretisation of the geological structures. Starting from the observed influence of heterogeneities on multiphase flow processes at the macro scale, we present an upscaling procedure from the local to the macro scale for the derivation of constitutive relationships for multiphase flow processes. The approach is based on the assumption of an equilibrium of (capillary) forces, which allows the application of a percolation model. This results in saturation distributions for different capillary pressures. Averaging these distributions gives rise to a macroscopic capillary pressure-saturation relationship. For the saturation distribution, relative permeabilities and effective conductivities are computed depending on the structure and the flow direction. These are averaged with the help of the renormalisation method. The evolving relative permeability-saturation relationship for the macro scale shows a saturation-dependent anisotropy and pronounced residual saturations of the nonwetting phase (which were not assumed for the local scale). The anisotropy reflects the underlying structure of the considered system that needs not to be known in detail.     

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.     

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.     

A multi-flowpath model for the interpretation of immiscible displacement experiments in heterogeneous soil columns

This work focuses on the phenomenon of the immiscible two-phase flow of water and oil in saturated heterogeneous soil columns. The goal is to develop a fast and reliable method for quantifying soil heterogeneities for incorporation into the relevant capillary pressure and relative permeability functions. Such data are commonly used as input data in simulators of contaminant transport in the subsurface. Rate-controlled drainage experiments are performed on undisturbed soil columns and the transient response of the axial distribution of water saturation is determined from electrical measurements. The transient responses of the axial distribution of water saturation and total pressure drop are fitted with the multi-flowpath model (MFPM) where the pore space is regarded as a system of parallel paths of different permeability. The MFPM enables us to quantify soil heterogeneity at two scales: the micro-scale parameters describe on average the effects of pore network heterogeneities on the two-phase flow pattern; the macro-scale parameters indicate the variability of permeability at the scale of interconnected pore networks. The capillary pressure curve is consistent with that measured with mercury intrusion porosimetry over the low pressure range. The oil relative permeability increases sharply at a very low oil saturation (< 10^− 3) and tends to a high end value. The water relative permeability decreases abruptly at a low oil saturation (~ 0.1), whereas the irreducible wetting phase saturation is quite high. The foregoing characteristics of the two-phase flow properties are associated with critical (preferential) flowpaths that comprise a very small percentage of the total pore volume, control the overall hydraulic conductivity, and are consistent with the very broad range of pore-length scales usually probed in soil porous matrix.     

Three-dimensional visualization and quantification of non-aqueous phase liquid volumes in natural porous media using a medical X-ray Computed Tomography scanner

This study demonstrates the capabilities of a typical medical X-ray Computed Tomography (CT) scanner to non-destructively quantify non-aqueous phase liquid (NAPL) volumes, saturation levels, and three-dimensional spatial distributions in packed soil columns. Columns packed with homogeneous sand, heterogeneous sand, or natural soil, were saturated with water and injected with known quantities of gasoline or tetrachloroethene and scanned. A methodology based on image subtraction was implemented for computing soil porosity and NAPL volumes in each 0.35 mm x 0.35 mm x 1 mm voxel of the columns. Elimination of sample positioning errors and instrument drift artifacts was essential for obtaining reliable estimates of above parameters. The CT data-derived total NAPL volume was in agreement with the measured NAPL volumes injected into the columns. CT data-derived NAPL volume is subject to a 2.6% error for PCE and a 15.5% error for gasoline, at average NAPL saturations as low as 5%, and is mainly due to instrument noise. Non-uniform distributions of NAPL due to preferential flow, and accumulation of NAPL above finer-grained layers could be observed from the data on 3-D distributions of NAPL volume fractions.     

Measurement and analysis of non-Fickian dispersion in heterogeneous porous media

Contaminant breakthrough behavior in a variety of heterogeneous porous media was measured in laboratory experiments, and evaluated in terms of both the classical advection-dispersion equation (ADE) and the continuous time random walk (CTRW) framework. Heterogeneity can give rise to non-Fickian transport patterns, which are distinguished by "anomalous" early arrival and late time tails in breakthrough curves. Experiments were conducted in two mid-scale laboratory flow cells packed with clean, sieved sand of specified grain sizes. Three sets of experiments were performed, using a "homogeneous" packing, a randomly heterogeneous packing using sand of two grain sizes, and an exponentially correlated structure using sand of three grain sizes. Concentrations of sodium chloride tracer were monitored at the inflow reservoir and measured at the outflow reservoir. Breakthrough curves were then analyzed by comparison to fitted solutions from the ADE and CTRW formulations. In all three systems, including the "homogeneous" one, subtle yet measurable differences between Fickian and non-Fickian transport were observed. Quantitative analysis demonstrated that the CTRW theory characterized the full shape of the breakthrough curves far more effectively than the ADE.     

Humic acid enhanced remediation of an emplaced diesel source in groundwater. 2. Numerical model development and application

A pilot scale experiment for humic acid-enhanced remediation of diesel fuel, described in Part 1 of this series, is numerically simulated in three dimensions. Groundwater flow, enhanced solubilization of the diesel source, and reactive transport of the dissolved contaminants and humic acid carrier are solved with a finite element Galerkin approach. The model (BIONAPL) is calibrated by comparing observed and simulated concentrations of seven diesel fuel components (BTEX and methyl-, dimethyl- and trimethylnaphthalene) over a 1500-day monitoring period. Data from supporting bench scale tests were used to estimate contaminant-carrier binding coefficients and to simulate two-site sorption of the carrier to the aquifer sand. The model accurately reproduced the humic acid-induced 10-fold increase in apparent solubility of trimethylnaphthalene. Solubility increases on the order of 2-5 were simulated for methylnaphthalene and dimethylnaphthalene, respectively. Under the experimental and simulated conditions, the residual 500-ml diesel source was almost completely dissolved and degraded within 5 years. Without humic acid flushing, the simulations show complete source dissolution would take about six times longer.     

Microbial physiology-based model of ethanol metabolism in subsurface sediments

Dense non-aqueous phase liquids (DNAPLs) present in the subsurface may contain surface active compounds that impact DNAPL migration and distribution. While a number of studies have revealed the role surface active compounds play in altering the wettability of quartz sand, few have considered the implications for other minerals common to contaminated sites. This study extends understanding of DNAPL/surfactant wettability to iron oxide surfaces. Specifically, quartz and iron oxide-coated sands in a tetrachloroethene (PCE)/water system containing the organic base (an organic molecule that acts as a base) dodecylamine (DDA) were compared at a variety of scales. Wettability of the minerals' surfaces, and the impact of wettability on capillary resistance to DNAPL entry, were assessed as a function of pH through: (i) advancing and receding contact angles, (ii) primary drainage capillary pressure-saturation experiments, and (iii) small, two-dimensional, flow cell experiments. The work revealed that, at neutral pH and under identical boundary capillary pressures, DNAPL invaded quartz sand but not iron oxide-coated sand; however, at low pH, DNAPL invaded both sands equally. These differences were demonstrated to be due to wettability alterations associated with the strength of attractive forces between DDA and the mineral surface, dictated by the isolectric point of the minerals and system pH. Observed differences in DNAPL invasion behavior were consistent with measured intrinsic contact angles and P(c)-S relationships, the latter requiring scaling by the operative contact angle inside the porous medium for a meaningful comparison. This study suggests that the distribution of minerals (and, more specifically, their isoelectric points), as well as the aqueous phase pH at a given site, may have a significant impact on the DNAPL source zone architecture.     

Infiltration and redistribution of LNAPL into unsaturated layered porous media

Enhanced understanding of light non-aqueous phase liquid (LNAPL) infiltration into heterogeneous porous media is important for the effective design of remediation strategies. We used a 2-D experimental facility that allows for visual observation of LNAPL contours in order to study LNAPL redistribution in a layered porous medium. The layers are situated in the unsaturated zone near the watertable and they are inclined to be able to observe the effect of discontinuities in capillary forces and relative permeabilities. Two experiments were performed. The first experiment consisted of LNAPL infiltration into a fine sand matrix with a coarse sand layer, and the second experiment consisted of a coarse sand matrix and a fine sand layer. The numerical multi-phase flow model STOMP was validated with regard to the experimental results. This model is able to adequately reproduce the experimental LNAPL contours. Numerical sensitivity analysis was also performed. The capillarity contrast between sands was found to be the main controlling factor determining the final LNAPL distribution.     

Evaluation of dynamic subgrid-scale models in large-eddy simulations of neutral turbulent flow over a two-dimensional sinusoidal hill

Large-eddy simulation (LES) is used to simulate neutral turbulent boundary-layer flow over a rough two-dimensional sinusoidal hill. Three different subgrid-scale (SGS) models are tested: (a) the standard Smagorinsky model with a wall-matching function, (b) the Lagrangian dynamic model, and (c) the recently developed scale-dependent Lagrangian dynamic model [Stoll, R., Porté-Agel, F., 2006. Dynamic subgrid-scale models for momentum and scalar fluxes in large-eddy simulation of neutrally stratified atmospheric boundary layers over heterogeneous terrain. Water Resources Research 42, W01409. doi:10.1029/2005WR003989]. The simulation results obtained with the different models are compared with turbulence statistics obtained from experiments conducted in the meteorological wind tunnel of the AES (Atmospheric Environment Service, Canada) [Gong, W., Taylor, P.A., Dörnbrack, A., 1996. Turbulent boundary-layer flow over fixed aerodynamically rough two-dimensional sinusoidal waves. Journal of Fluid Mechanics 312, 1–37]. We find that the scale-dependent dynamic model is able to account, without any tuning, for the local changes in the eddy-viscosity model coefficient. It can also capture the scale dependence of the coefficient associated with regions of the flow with strong mean shear and flow anisotropy. As a result, the scale-dependent dynamic model yields results that are more realistic than the ones obtained with the scale-invariant Lagrangian dynamic model.     

Characterization of pore scale NAPL morphology in homogeneous sands as a function of grain size and NAPL dissolution

In this study, we investigate pore scale morphology of nonaqueous phase liquids (NAPLs) trapped in different pore sizes using tracer techniques. Specific interfacial area and saturation of NAPL trapped in homogeneous sands were measured using the interfacial and partitioning tracer techniques. The observed NAPL-water interfacial areas increased in a log-linear fashion with decreasing sand grain size, but showed no clear trend with residual NAPL saturation formed in the various grain sizes. The measured values were used to calculate the NAPL morphology index, which characterizes the spatial NAPL distribution within the pore space. The NAPL morphology indices, increased exponentially with decreasing grain size, indicating that the NAPL becomes smaller, but more blobs. For a fixed grain size, the specific interfacial area and saturation of the NAPL were measured following changes caused by dissolution using alcohol. The observed interfacial areas showed a decrease linearly as a function of the NAPL saturation.     

The influence of field-scale heterogeneity on the infiltration and entrapment of dense nonaqueous phase liquids in saturated formations

Numerical simulations are used for the systematic exploration of the migration and entrapment of dense nonaqueous phase liquids (DNAPLs) in heterogeneous formations. Ensembles of realizations of random, spatially correlated permeability fields are generated and employed in model simulations of a spill event. Statistical techniques are then used to quantify the sensitivity of model predictions to input parameters, thereby identifying the parameters or processes that may be of primary importance in the determination of organic liquid distributions in heterogeneous systems. Results of the study indicate that the most critical factors in modeling organic entrapment include the spill release rate, reliable estimates of the mean, variance, and vertical correlation scale of the formation permeability, and an accurate representation of the correlation between the capillary pressure–saturation function and the permeability. In contrast, the hydraulic gradient and cross-correlation of residual saturations with permeabilities are found to have only minor influence on organic liquid distributions in such heterogeneous formations.     

Relationships between direction of wind flow and fine particle sulfur concentrations within and upwind of St Louis,MO

The relationships are considered between monthly and quarterly means of the fine particle sulfur (S) concentrations and wind flow direction, period of day and season of the year. The measurements used are those obtained at selected urban and at rural monitoring stations in the St Louis area during the Regional Air Pollution Study in 1975, 1976 and 1977. Higher mean fine particle S concentrations are observed with wind flows from the E compared to the W and from the NE and SE quadrants compared to the NW quadrant. Substantially higher fine particle S concentrations are obtained with wind flows from the E compared to the W even when conditions are selected so that the values of temperature, solar radiation intensity and wind speed are within the same restricted ranges. A consistent increase in the fine particle S concentrations occurs through the late morning and afternoon with decrease in the evening and especially in the early morning during spring and summer months with wind flows from the E.The contributions are estimated for local scale and regional scale processes to the observed fine particle S concentrations. Local scale processes include those involving atmospheric formation and primary emissions each contribute 0.6−1.0 and 0.6 μg m⁻³ of the fine particle S. Regional scale processes account for the greater part of the observed concentrations especially when the wind flows are from the SE or SW. Regional scale episodes involving passage of warm high pressure systems to the E of St Louis with accumulation of precursors made especially significant contributions to formation of fine particle S.The atmospheric gas phase and liquid phase chemical reactions contributing to the formation of fine particle S are discussed. Emphasis is placed on the effects of chemistry on the seasonal variations in concentrations of fine particle S.     

2,4-D Mineralization in soil profiles of a cultivated hummocky landscape in Manitoba, Canada

The objective of this study was to quantify 2,4-D (2,4-dichlorophenoxyacetic acid) mineralization in soil profiles characteristic of hummocky, calcareous-soil landscapes in western Canada. Twenty-five soil cores (8 cm inner diameter, 50 to 125 cm length) were collected along a 360 m transect running west to east in an agricultural field and then segmented by soil-landscape position (upper slopes, mid slopes, lower slopes and depressions) and soil horizon (A, B, and C horizons). In the A horizon, 2,4-D mineralization commenced instantaneously and the mineralization rate followed first-order kinetics. In both the B and C horizons, 2,4-D mineralization only commenced after a lag period of typically 5 to 7 days and the mineralization rate was biphasic. In the A horizon, 2,4-D mineralization parameters including the first-order mineralization rate constant (k(1)), the growth-linked mineralization rate constant (k(2)) and total 2,4-D mineralization at the end of the experiment at 56 days, were most strongly correlated to parameters describing 2,4-D sorption by soil, but were also adequately correlated to soil organic carbon content, soil pH, and carbonate content. In both B and C horizons, there was no significant correlation between 2,4-D mineralization and 2,4-D sorption parameters, and the correlation between soil properties and 2,4-D mineralization parameters was very poor. The k(1) significantly decreased in sequence of A horizon (0.113% day(-1)) > B horizon (0.024% day(-1)) = C horizon (0.026% day(-1)) and in each soil horizon was greater than k(2). Total 2,4-D mineralization at 56 days also significantly decreased in sequence of A horizon (42%) > B horizon (31%) = C horizon (27%). In the A horizon, slope position had little influence on k(1) or k(2), except that k(1) was significantly greater in upper slopes (0.170% day(-1)) than in lower slopes (0.080% day(-1)). Neither k(1) nor k(2) was significantly influenced by slope position in the B or C horizons. Total 2,4-D mineralization at 56 days was not influenced by slope positions in any horizon. Our results suggest that, when predicting 2,4-D transport at the field scale, pesticide fate models should consider the strong differences in 2,4-D mineralization between surface and subsurface horizons. This suggests that 2,4-D mineralization is best predicted using a model that has the ability to describe a range of non-linear mineralization curves. We also conclude that the horizontal variations in 2,4-D mineralization at the field scale will be difficult to consider in predictions of 2,4-D transport at the field scale because, within each horizon, 2,4-D mineralization was highly variable across the twenty-five soil cores, and this variability was poorly correlated to soil properties or soil-landscape position.     

Effective parameters for two-phase flow in a porous medium with periodic heterogeneities

Computational simulations of two-phase flow in porous media are used to investigate the feasibility of replacing a porous medium containing heterogeneities with an equivalent homogeneous medium. Simulations are performed for the case of infiltration of a dense nonaqueous phase liquid (DNAPL) in a water-saturated, heterogeneous porous medium. For two specific porous media, with periodic and rather simple heterogeneity patterns, the existence of a representative elementary volume (REV) is studied. Upscaled intrinsic permeabilities and upscaled nonlinear constitutive relationships for two-phase flow systems are numerically calculated and the effects of heterogeneities are evaluated. Upscaled capillary pressure-saturation curves for drainage are found to be distinctly different from the lower-scale curves for individual regions of heterogeneity. Irreducible water saturation for the homogenized medium is found to be much larger than the corresponding lower-scale values. Numerical simulations for both heterogeneous and homogeneous representations of the considered porous media are carried out. Although the homogenized model simulates the spreading behavior of DNAPL reasonably well, it still fails to match completely the results form the heterogeneous simulations. This seems to be due, in part, to the nonlinearities inherent to multiphase flow systems. Although we have focussed on a periodic heterogeneous medium in this study, our methodology is applicable to other forms of heterogeneous media. In particular, the procedure for identification of a REV, and associated upscaled constitutive relations, can be used for randomly heterogeneous or layered media as well.     

Infiltration of PCE in a system containing spatial wettability variations

A two-dimensional infiltration experiment was conducted to investigate and quantify the effect of spatial wettability variations on DNAPL migration and entrapment in saturated sands. Experimental observations of tetrachloroethylene (PCE) infiltration showed that organic-wet sand lenses acted as very effective capillary barriers, retaining PCE and inhibiting its downward migration. A multiphase numerical simulator was used to model this sand box experiment. The simulator incorporates wettability-modified van Genuchten and Brooks-Corey capillary pressure/saturation relationships as well as Burdine and Mualem relative permeability relationships. PCE mass distributions, estimated by image analysis of digital photographs taken during the infiltration event, were compared to simulation results. Although both relative permeability models were qualitatively able to predict the PCE retention in the organic-wet layers, simulations with the Mualem model failed to capture the observed rate of PCE migration. A traditional multiphase simulator, incorporating water-wet capillary retention relations, failed to predict both PCE pathways and retention behavior. This study illustrates the potential influence of subsurface wettability variations on DNAPL migration and entrapment and supports the use of modified capillary relations in conjunction with the Burdine model in multiphase flow simulators.     

自制填料与常用湿地填料除磷能力的比较

用石英砂、石灰、水泥和铝粉等材料制作人工湿地填料,并对自制填料、粉煤灰块和钢渣等6种填料的等温吸附特性进行研究,发现自制填料在较低与较高浓度P溶液中吸附率都最高,在P浓度为(10～50 mg/L)时,平均吸附率达96.8%,其次是粉煤灰和钢渣,分别为87.6%和85.4%;填料对P的理论饱和吸附量分别为:自制填料(2 413.6 mg/kg)>粉煤灰块(1 605.8 mg/kg)>钢渣(1 277.5 mg/kg)>碎砖(451.6 mg/kg)>碎石(182.5 mg/kg)>砾石(18.4 mg/kg);填料中金属矿物成分含量高、比表面积大,是除P效果好的原因。     

Analysis of a mesoscale infiltration and water seepage test in unsaturated fractured rock: Spatial variabilities and discrete fracture patterns

A mesoscale (21 m in flow distance) infiltration and seepage test was recently conducted in a deep, unsaturated fractured rock system at the crossover point of two underground tunnels. Water was released from a 3 mx4 m infiltration plot on the floor of an alcove in the upper tunnel, and seepage was collected from the ceiling of a niche in the lower tunnel. Significant temporal and (particularly) spatial variabilities were observed in both measured infiltration and seepage rates. To analyze the test results, a three-dimensional unsaturated flow model was used. A column-based scheme was developed to capture heterogeneous hydraulic properties reflected by these spatial variabilities observed. Fracture permeability and van Genuchten alpha parameter [van Genuchten, M.T., 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 44, 892-898] were calibrated for each rock column in the upper and lower hydrogeologic units in the test bed. The calibrated fracture properties for the infiltration and seepage zone enabled a good match between simulated and measured (spatially varying) seepage rates. The numerical model was also able to capture the general trend of the highly transient seepage processes through a discrete fracture network. The calibrated properties and measured infiltration/seepage rates were further compared with mapped discrete fracture patterns at the top and bottom boundaries. The measured infiltration rates and calibrated fracture permeability of the upper unit were found to be partially controlled by the fracture patterns on the infiltration plot (as indicated by their positive correlations with fracture density). However, no correlation could be established between measured seepage rates and density of fractures mapped on the niche ceiling. This lack of correlation indicates the complexity of (preferential) unsaturated flow within the discrete fracture network. This also indicates that continuum-based modeling of unsaturated flow in fractured rock at mesoscale or a larger scale is not necessarily conditional explicitly on discrete fracture patterns.     

The impact of methanogenesis on flow and transport in coarse sand

The effects of biofilm growth and methane gas generation on water flow in porous media were investigated in an anaerobic two-dimensional sand-filled cell. Inoculation of the lower portion of the cell with a methanogenic culture and addition of methanol to the bottom of the cell led to biomass growth and formation of a gas phase. Biomass distributions in the water and on the sand in the cell were measured by protein analysis. The biofilm distribution on sand was observed by confocal laser scanning microscopy. The formation, migration, distribution and saturation of gases in the cell were visualized by the charge-coupled device (CCD) camera. The effects of biofilm and gas generation on water flow were separated by performing one tracer test in the presence of both biofilm and a gas phase and a second tracer test after removal of the gas phase through water flushing. The results of tracer tests demonstrated that flow and transport in the two-dimensional cell were significantly affected by both gas generation and biofilm growth. Gas generated at the bottom of the cell in the biologically active zone moved upwards in discrete fingers, so that gas phase saturations (gas-filled fraction of void space) in the biologically active zone at the bottom of the cell did not exceed 40-50%, while gas accumulation at the top of the cell produced gas phase saturations as high as 80%. The greatest reductions in water phase permeability, based on measurements of reductions in water phase saturations, occurred near the top of the box as a result of the gas accumulation. In contrast the greatest reductions in permeability due to biofilm growth, based on measurements of biofilm thickness, occurred in the most biologically active zone at the bottom of the cell, where gas phase saturations were approximately 40-50%, but permeability reductions due to biofilm growth were estimated to be 80-95%.     

Effects of a mixture of tetracyclines to Lemna gibba and Myriophyllum sibiricum evaluated in aquatic microcosms

The impact of a mixture of oxytetracycline, chlortetracycline, tetracycline and doxycycline on Myriophyllum sibiricum and Lemna gibba was investigated using fifteen 12,000-L microcosms (k=5, n=3). Significant concentration-response relationships were only found for M. sibiricum, where dry mass was 69, 47, 30, and 7% of controls at respective treatment concentrations of 0.080, 0.218, 0.668, and 2.289 micromol/L. Somatic endpoints were strongly and negatively correlated with percent light transmission, except plant length, which was positively correlated. Treated microcosms experienced a reduction in the percent of surface irradiance penetrating the water column as high as 99.8% at a depth of 70 cm, relative to controls. Position relative to the water column was likely responsible for the differential effects observed between floating (L. gibba) and submerged (M. sibiricum) species of macrophytes. A hazard quotient assessment of the lowest EC10 value indicated significant risk, exceeding the critical HQ value, but not the lowest EC25 value.