On the use of mean groundwater age, life expectancy and capture probability for defining aquifer vulnerability and time-of-travel zones for source water protection

Protection and sustainability of water supply wells requires the assessment of vulnerability to contamination and the delineation of well capture zones. Capture zones, or more generally, time-of-travel zones corresponding to specific contaminant travel times, are most commonly delineated using advective particle tracking. More recently, the capture probability approach has been used in which a probability of capture of P=1 is assigned to the well and the growth of a probability-of-capture plume is tracked backward in time using an advective-dispersive transport model. This approach accounts for uncertainty due to local-scale heterogeneities through the use of macrodispersion. In this paper, we develop an alternative approach to capture zone delineation by applying the concept of mean life expectancy E (time remaining before being captured by the well), and we show how life expectancy E is related to capture probability P. Either approach can be used to delineate time-of-travel zones corresponding to specific travel times, as well as the ultimate capture zone. The related concept of mean groundwater age A (time since recharge) can also be applied in the context of defining the vulnerability of a pumped aquifer. In the same way as capture probability, mean life expectancy and groundwater age account for local-scale uncertainty or unresolved heterogeneities through macrodispersion, which standard particle tracking neglects. The approach is tested on 2D and 3D idealized systems, as well as on several watershed-scale well fields within the Regional Municipality of Waterloo, Ontario, Canada.     

Probabilistic study of well capture zones distribution at the Lauswiesen field site

The delineation of well capture zones is of utmost environmental and engineering relevance as pumping wells are commonly used both for drinking water supply needs, where protection zones have to be defined, and for investigation and remediation of contaminated aquifers. We analyze the probabilistic nature of well capture zones within the well field located at the "Lauswiesen" experimental site. The test site is part of an alluvial heterogeneous aquifer located in the Neckar river valley, close to the city of Tübingen in South-West Germany. We explore the effect of different conceptual models of the structure of aquifer heterogeneities on the delineation of three-dimensional probabilistic well catchment and time-related capture zones, in the presence of migration of conservative solutes. The aquifer is modeled as a three-dimensional, doubly stochastic composite medium, where distributions of geo-materials and hydraulic properties are uncertain. We study the relative importance of uncertain facies geometry and uncertain hydraulic conductivity and porosity on predictions of catchment and solute time of travel to the pumping well by focusing on cases in which (1) the facies distribution is random, but the hydraulic properties of each material are fixed, and (2) both facies geometry and material properties vary stochastically. The problem is tackled within a conditional numerical Monte Carlo framework. Results are provided in terms of probabilistic demarcations of the three-dimensional well catchment and time-related capture zones. Our findings suggest that the uncertainty associated with the prediction of the location of the outer boundary of well catchment at the "Lauswiesen" site is significantly affected by the conceptual model adopted to incorporate the heterogeneous nature of the aquifer domain in a predictive framework. Taking into account randomness of both lithofacies distribution and materials hydraulic conductivity allows recognizing the existence of preferential flow paths that influence the extent of the well catchment and the solute travel time distribution at the site.     

Transport of reactive colloids and contaminants in groundwater: effect of nonlinear kinetic interactions

Transport of reactive colloids in groundwater may enhance the transport of contaminants in groundwater. Often, the interpretation of results of transport experiments is not a simple task as both reactions of colloids with the solid matrix and reactions of contaminants with the solid matrix and mobile and immobile colloids may be time dependent and nonlinear. Further colloid transport properties may differ from solute transport properties. In this paper, a one-dimensional model for coupled and contaminant in a porous medium (COLTRAP) is presented together with simulation results. Calculated breakthrough curves (BTC's) during contamination and decontamination show systematically the effect of nonlinear and kinetic interactions on contaminant transport in the presence of reactive colloids, and the effect of colloid transport properties that differ from solute transport properties. It is shown that in case of linear kinetic reactions, the rate of exchange of mobile and immobile colloids have a large impact on the shape of BTC's even if the solid matrix is saturated with respect to colloids. BTC's during the contamination and decontamination phase have identical shapes in this case. Moreover, the slow reactions of contaminants and colloids may lead to unretarded breakthrough of contaminants. Independent of reaction rates, nonlinear reactions lead to BTC's that are steeper during contamination than in the linear case. A characteristic aspect of nonlinear sorption is that shapes of BTC's differ during the contamination and decontamination phase. It has been observed that shapes of some of the simulated adsorption and desorption curves are similar as shapes found in experiments reported in literature. This stresses the importance of incorporating both kinetics and nonlinearity in models for coupled colloid and contaminant transport and the capability of COLTRAP to interpret experimental results. Finally, to figure out whether nonlinear processes play a role, it is very important to consider both contamination and decontamination in transport experiments.     

A probability model for evaluating building contamination from an environmental event

Asbestos dust and bioaerosol sampling data from suspected contaminated zones in buildings allowed development of an environmental data evaluation protocol based on the differences in frequency of detection of a target contaminant between zones of comparison. Under the assumption that the two test zones of comparison are similar, application of population proportion probability calculates the significance of observed differences in contaminant levels. This was used to determine whether levels of asbestos dust contamination detected after a fire were likely the result of smoke-borne contamination, or were caused by pre-existing/background conditions. Bioaerosol sampling from several sites was also used to develop the population proportion probability protocol. In this case, significant differences in indoor air contamination relative to the ambient conditions were identified that were consistent with the visual observations of contamination. Implicit in this type of probability analysis is a definition of "contamination" based on significant differences in contaminant levels relative to a control zone. Detection of a suspect contaminant can be assessed as to possible sources(s) as well as the contribution made by pre-existing (i.e., background) conditions, provided the test and control zones are subjected to the same sampling and analytical methods.     

A Probability Model for Evaluating Building Contamination from an Environmental Event

ABSTRACT

Asbestos dust and bioaerosol sampling data from suspected contaminated zones in buildings allowed development of an environmental data evaluation protocol based on the differences in frequency of detection of a target contaminant between zones of comparison. Under the assumption that the two test zones of comparison are similar, application of population proportion probability calculates the significance of observed differences in contaminant levels. This was used to determine whether levels of asbestos dust contamination detected after a fire were likely the result of smoke-borne contamination, or were caused by pre-existing/background conditions.

Bioaerosol sampling from several sites was also used to develop the population proportion probability protocol. In this case, significant differences in indoor air contamination relative to the ambient conditions were identified that were consistent with the visual observations of contamination. Implicit in this type of probability analysis is a definition of “contamination” based on significant differences in contaminant levels relative to a control zone. Detection of a suspect contaminant can be assessed as to possible sources(s) as well as the contribution made by pre-existing (i.e., background) conditions, provided the test and control zones are subjected to the same sampling and analytical methods.     

Assessment of metal pollution based on multivariate statistical modeling of 'hot spot' sediments from the Black Sea

The paper deals with application of different statistical methods like cluster and principal components analysis (PCA), partial least squares (PLSs) modeling. These approaches are an efficient tool in achieving better understanding about the contamination of two gulf regions in Black Sea. As objects of the study, a collection of marine sediment samples from Varna and Bourgas "hot spots" gulf areas are used. In the present case the use of cluster and PCA make it possible to separate three zones of the marine environment with different levels of pollution by interpretation of the sediment analysis (Bourgas gulf, Varna gulf and lake buffer zone). Further, the extraction of four latent factors offers a specific interpretation of the possible pollution sources and separates natural from anthropogenic factors, the latter originating from contamination by chemical, oil refinery and steel-work enterprises. Finally, the PLSs modeling gives a better opportunity in predicting contaminant concentration on tracer (or tracers) element as compared to the one-dimensional approach of the baseline models. The results of the study are important not only in local aspect as they allow quick response in finding solutions and decision making but also in broader sense as a useful environmetrical methodology.     

Travel time and trajectory moments of conservative solutes in two-dimensional convergent flows

We address advective transport of a solute traveling toward a single pumping well in a two-dimensional randomly heterogeneous aquifer. The two random variables of interest are the trajectory followed by an individual particle from the injection point to the well location and the particle travel time under steady-state conditions. Our main objective is to derive the predictors of trajectory and travel time and the associated uncertainty, in terms of their first two statistical moments (mean and variance). We consider a solute that undergoes mass transfer between a mobile and an immobile zone. Based on Lawrence et al. [Lawrence, A.E., Sánchez-Vila, X., Rubin, Y., 2002. Conditional moments of the breakthrough curves of kinetically sorbing solute in heterogeneous porous media using multirate mass transfer models for sorption and desorption. Water Resour. Res. 38 (11), 1248, doi:10.1029/2001WR001006.], travel time moments can be written in terms of those of a conservative solute times a deterministic quantity. Moreover, the moments of solute particles trajectory do not depend on mass transfer processes. The resulting mean and variance of travel time and trajectory for a conservative species can be written as functions of the first, second moments and cross-moments of trajectory and velocity components. The equations are developed from a consistent second order expansion in sigmaY (standard deviation of the natural logarithm of hydraulic conductivity). Our solution can be completely integrated with the moment equations of groundwater flow of Guadagnini and Neuman [Guadagnini, A., Neuman, S.P., 1999a. Nonlocal and localized analyses of conditional mean steady state flow in bounded, randomly non uniform domains 1. Theory and computational approach. Water Resour. Res. 35(10), 2999-3018.,Guadagnini, A., Neuman, S.P., 1999b. Nonlocal and localized analyses of conditional mean steady state flow in bounded, randomly non uniform domains 2. Computational examples. Water Resour. Res. 35(10), 3019-3039.], it is free of distributional assumptions regarding the log conductivity field, and formally includes conditioning. We present analytical expressions for the unconditional case by making use of the results of Riva et al. [Riva, M., Guadagnini, A., Neuman, S.P., Franzetti, S., 2001. Radial flow in a bounded randomly heterogeneous aquifer. Transport in Porous Media 45, 139-193.]. The quality of the solution is supported by numerical Monte Carlo simulations. Potential uses of this work include the determination of aquifer reclamation time by means of a single pumping well, and the demarcation of the region potentially affected by the presence of a contaminant in the proximity of a well, whenever the aquifer is very thin and Dupuit-Forchheimer assumption holds.     

Description of hydrogeologic heterogeneity and evaluation of radionuclide transport at an underground nuclear test

Realistic models of contaminant transport in groundwater demand detailed characterization of the spatial distribution of subsurface hydraulic properties, while at the same time programmatic constraints may limit collection of pertinent hydraulic data. Fortunately, alternate forms of data can be used to improve characterization of spatial variability. We utilize a methodology that augments sparse hydraulic information (hard data) with more widely available hydrogeologic information to generate equiprobable maps of hydrogeologic properties that incorporate patterns of connected permeable zones. Geophysical and lithologic logs are used to identify hydrogeologic categories and to condition stochastic simulations using Sequential Indicator Simulation (SIS). The resulting maps are populated with hydraulic conductivity values using field data and Sequential Gaussian Simulation (SGS). Maps of subsurface hydrogeologic heterogeneity are generated for the purpose of examining groundwater flow and transport processes at the Faultless underground nuclear test, Central Nevada Test Area (CNTA), through large-scale, three-dimensional numerical modeling. The maps provide the basis for simulation of groundwater flow, while transport of radionuclides from the nuclear cavity is modeled using particle tracking methods. Sensitivity analyses focus on model parameters that are most likely to reduce the long travel times observed in the base case. The methods employed in this study have improved our understanding of the spatial distribution of preferential flowpaths at this site and provided the critical foundation on which to build models of groundwater flow and transport. The results emphasize that the impacts of uncertainty in hydraulic and chemical parameters are dependent on the radioactive decay of specific species, with rapid decay magnifying the effects of parameters that change travel time.     

Estimating spatially-variable first-order rate constants in groundwater reactive transport systems

Numerical reactive transport models are often used as tools to assess aquifers contaminated with reactive groundwater solutes as well as investigating mitigation scenarios. The ability to accurately simulate the fate and transport of solutes, however, is often impeded by a lack of information regarding the parameters that define chemical reactions. In this study, we employ a steady-state Ensemble Kalman Filter (EnKF), a data assimilation algorithm, to provide improved estimates of a spatially-variable first-order rate constant λ through assimilation of solute concentration measurement data into reactive transport simulation results. The methodology is applied in a steady-state, synthetic aquifer system in which a contaminant is leached to the saturated zone and undergoes first-order decay. Multiple sources of uncertainty are investigated, including hydraulic conductivity of the aquifer and the statistical parameters that define the spatial structure of the parameter field. For the latter scenario, an iterative method is employed to identify the statistical mean of λ of the reference system. Results from all simulations show that the filter scheme is successful in conditioning the λ ensemble to the reference λ field. Sensitivity analyses demonstrate that the estimation of the λ values is dependent on the number of concentration measurements assimilated, the locations from which the measurement data are collected, the error assigned to the measurement values, and the correlation length of the λ fields.     

Analytical solution to transport in three-dimensional heterogeneous well capture zones

Solute transport is investigated in a heterogeneous aquifer for combined natural-gradient and well flows. The heterogeneity is associated with the spatially varying hydraulic conductivity K(x, y, z), which is modelled as a log-normal stationary-random function. As such, the conductivity distribution is characterized by four parameters: the arithmetic mean K(A), the variance sigma(Y)(2) (Y=lnK), the horizontal integral scale I of the axisymmetric log-conductivity autocorrelation and the anisotropy ratio e=I(v)/I (I(v) is the vertical integral scale). The well fully penetrates an aquifer of constant thickness B and has given constant discharge QB, while the background aquifer flow is driven by an uniform mean head-gradient, - J. Therefore, for a medium of homogeneous conductivity K(A), the steady-state capture zone has a width 2L=Q/(K(A)|J|) far from the well (herein the term capture zone is used to refer both to the zone from which water is captured by a pumping well and the zone that captures fluid from an injecting well). The main aim is to determine the mean concentration as a function of time in fluid recovered by a pumping well or in a control volume of the aquifer that captures fluid from an injecting well. Relatively simple solutions to these complex problems are achieved by adopting a few assumptions: a thick aquifer B>I(v) of large horizontal extent (so that boundary effects may be neglected), weak heterogeneity sigma(Y)(2)<1, a highly anisotropic formation e<0.2 and neglect of pore-scale dispersion. Transport is analyzed to the first-order in sigma(Y)(2) in terms of the travel time of particles moving from or towards the well along the steady streamlines within the capture zone. Travel-time mean and variance to any point are computed by two quadratures for an exponential log-conductivity two-point covariance. Spreading is reflected by the variance value, which increases with sigma(Y)(2) and I/L. For illustration, the procedure is applied to two particular cases. In the first one, a well continuously injects water at constant concentration. The mean concentration as function of time for different values of the controlling parameters sigma(Y)(2) and I/L is determined within control volumes surrounding the well or in piezometers. In the second case, a solute plume, initially occupying a finite volume Omega(0), is drawn towards a pumping well. The expected solute recovery by the well as a function of time is determined in terms of the previous controlling parameters as well as the location and extent of Omega(0). The methodology is tested against a full three-dimensional simulation of a multi-well forced-gradient flow field test ([Lemke, L., W.B. II, Abriola, L., Goovaerts, P., 2004. Matching solute breakthrough with deterministic and stochastic aquifer models. Ground Water 42], SGS simulations). Although the flow and transport conditions are more complex than the ones pertinent to capture zones in uniform background flow, it was found that after proper adaptation the methodology led to results for the breakthrough curve in good agreement with a full three-dimensional simulation of flow and transport.     

Modeling of flow and contaminant transport in coupled stream-aquifer systems

This paper describes an integrated approach for modeling flow and contaminant transport in hydraulically connected stream-aquifer systems. The code, FTSTREAM, extended the capabilities of the ground-water model, FTWORK, to incorporate chemical fate and transport in streams. Flow in the stream network is modeled as an unsteady, spatially varying flow, while transport modeling is based on a one-dimensional advection-dispersion equation. In addition to sorption and decay during transport in ground water, the model incorporates volatilization, settling and decay during transport in surface water. The interaction between surface water and ground water is accommodated by a leakage term and is implemented in the model using an iterative Picard-type procedure to ensure mass conservation across the interface between the two systems. The modeling approach is used to simulate contaminant transport in the Mad River, Ohio, which is hydraulically connected to a buried valley aquifer of sand and gravel outwash. The river is a receiving stream in the upstream part of the modeled area. Downstream, heavy pumping from a municipal well field causes the river to become a loosing stream. Induced infiltration from the river is responsible for a considerable portion of the well yield. The flow and transport model, developed for this study, simulates coupling between flow in the aquifer and the river. Hypothetical sources of contamination are introduced at selected locations in the upstream portion of the aquifer. The model is then used to simulate the expected transport in both the aquifer and the stream. A series of simulations elucidates the role of the river in facilitating the transport of the hypothetical contaminants in ground water and surface water. Effect of sorption, retardation and volatilization on contaminant transport is also examined for the case of the volatile organic compounds.     

New methodology to investigate potential contaminant mass fluxes at the stream-aquifer interface by combining integral pumping tests and streambed temperatures

The spatial pattern and magnitude of mass fluxes at the stream-aquifer interface have important implications for the fate and transport of contaminants in river basins. Integral pumping tests were performed to quantify average concentrations of chlorinated benzenes in an unconfined aquifer partially penetrated by a stream. Four pumping wells were operated simultaneously for a time period of 5 days and sampled for contaminant concentrations. Streambed temperatures were mapped at multiple depths along a 60m long stream reach to identify the spatial patterns of groundwater discharge and to quantify water fluxes at the stream-aquifer interface. The combined interpretation of the results showed average potential contaminant mass fluxes from the aquifer to the stream of 272microgm(-2)d(-1) MCB and 71microgm(-2)d(-1) DCB, respectively. This methodology combines a large-scale assessment of aquifer contamination with a high-resolution survey of groundwater discharge zones to estimate contaminant mass fluxes between aquifer and stream.     

Optimization-Based Atmospheric Plume Source Identification

This work applies optimization and an Eulerian inversion approach presented by Bagtzoglou and Baun in 2005 in order to reconstruct contaminant plume time histories and to identify the likely source of atmospheric contamination using data from a real test site for the first time. Present-day distribution of an atmospheric contaminant plume as well as data points reflecting the plume history allow the reconstruction and provide the plume velocity, distribution, and probable source. The method was tested to a hypothetical case and with data from the Forest Atmosphere Transfer and Storage (FACTS) experiment in the Duke experimental forest site. In the scenarios presented herein, as well as in numerous cases tested for verification purposes, the model conserved mass, successfully located the peak of the plume, and managed to capture the motion of the plume well but underestimated the contaminant peak.     

Redistribution of contaminants by a fluctuating water table in a micro-porous, double-porosity aquifer: field observations and model simulations

Large seasonal fluctuations of the water table are characteristic of aquifers with a low specific yield, including those fractured, double-porosity aquifers that have significant matrix porosity containing virtually immobile porewater, such as the Chalk of northern Europe. Where these aquifers are contaminated, a strong relationship between water table elevation and contaminant concentration in groundwater is commonly observed, of significance to the assessment, monitoring, and remediation of contaminated groundwater. To examine the processes governing contaminant redistribution by a fluctuating water table within the 'seasonally unsaturated zone', or SUZ, profiles of porewater solute concentrations have been established at a contaminated site in southern England. These profiles document the contaminant distribution in porewater of the Chalk matrix over the SUZ at a greater level of detail than recorded previously. A novel double-porosity solute transport code has been developed to simulate the evolution of the SUZ matrix porewater contaminant profiles, given a fluctuating water table, when the groundwater is initially contaminated and the SUZ is initially free of contamination. The model is simply characterised by: the matrix-fracture porosity ratio, the matrix block geometry, and a characteristic diffusion time. De-saturation and re-saturation of fractures is handled by a new approximation method. Contaminant accumulates in the upper levels of the SUZ, where it is less accessible to mobile groundwater, and acts as a persistent secondary source of contamination once the original source of contamination has been removed or has become depleted. The 'SUZ process' first attenuates the progress of contaminants in groundwater, and subsequently controls the slow release of contamination back to the mobile groundwater, thus prolonging the duration of groundwater contamination by many years. The SUZ process should operate in any fractured, micro-porous lithology e.g. fractured clays and mudstones, making this approach widely applicable.     

Influence of indoor transport and mixing time scales on the performance of sensor systems for characterizing contaminant releases

Optimizing real-time sensor systems to detect and identify relevant characteristics of an indoor contaminant event is a challenging task. The interpretation of incoming sensor data is confounded by uncertainties in building operation, in the forces driving contaminant transport, and in the physical parameters governing transport. In addition, simulation tools used by the sensor interpretation algorithm introduce modeling uncertainties. This paper explores how the time scales inherent in contaminant transport influence the information that can be extracted from real-time sensor data. In particular, we identify three time scales (within room mixing, room-to-room transport, and removal from the building) and study how they affect the ability of a Bayesian Monte Carlo (BMC) sensor interpretation algorithm to identify the release location and release mass from a set of experimental data, recorded in a multi-floor building. The research shows that some limitations in the BMC approach do not depend on details of the models or the algorithmic implementation, but rather on the physics of contaminant transport. These inherent constraints have implications for the design of sensor systems.     

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.     

Assessing performance and closure for soil vapor extraction: integrating vapor discharge and impact to groundwater quality

Soil vapor extraction (SVE) is typically effective for removal of volatile contaminants from higher-permeability portions of the vadose zone. However, contamination in lower-permeability zones can persist due to mass transfer processes that limit the removal effectiveness. After SVE has been operated for a period of time and the remaining contamination is primarily located in lower-permeability zones, the remedy performance needs to be evaluated to determine whether the SVE system should be optimized, terminated, or transitioned to another technology to replace or augment SVE. Numerical modeling of vapor-phase contaminant transport was used to investigate the correlation between measured vapor-phase mass discharge, MF(r), from a persistent, vadose-zone contaminant source and the resulting groundwater contaminant concentrations. This relationship was shown to be linear, and was used to directly assess SVE remediation progress over time and to determine the level of remediation in the vadose zone necessary to protect groundwater. Although site properties and source characteristics must be specified to establish a unique relation between MF(r) and the groundwater contaminant concentration, this correlation provides insight into SVE performance and support for decisions to optimize or terminate the SVE operation or to transition to another type of treatment.     

Freezing of a clayey silt contaminated with an organic solvent

In a study of water migration characteristics and organic contaminant transfer mechanisms in a freezing fine-grained saturated soil, a series of one-dimensional freezing tests were conducted on a clayey silt contaminated with a miscible, non-reactive organic compound, 1-propanol, at various concentrations. The experimental results indicate that the frost heave behaviour and solute rejection mechanisms of a soil contaminated with 1-propanol is similar to that of the same soil contaminated with sodium-chloride salt. It was found that 1-propanol is rejected from the pore for rates of cooling smaller than 4 ± 1°C/day. Diffusion appears to control contaminant redistribution in the unfrozen soil. Finally, there has been no contaminant redistribution in the frozen soil for periods up to 245 hours.     

Using the multiple cell balance method to solve the problem of two-dimensional groundwater flow and contaminant transport with nonequilibrium sorption

The problem of large-scale contamination of groundwater by relatively low levels of organic contaminants is most frequently addressed by extracting and treating the impacted groundwater. This pump-and-treat strategy is often unsuccessful because of difficulties encountered in recovering the contaminants from relatively immobile zones within the porous medium. These zones can exist at the particle scale, as intraparticle or intra-aggregate porosity, and at the larger scales, as low-permeability layers or lenses interspersed in substantially more permeable layers. This work focuses on achieving an efficient numerical solution to a system of groundwater flow and contaminant transport equations that sufficiently captures the dynamics of slow desorption in a two-dimensional porous medium. The conceptual model and governing equations are presented. A numerical method for solving the governing equations, the upstream-weighted, multiple cell balance (UMCB) method, is proposed. The UMCB algorithm has been employed previously for the case of solute transport with equilibrium sorption, and is extended here to the nonequilibrium case. The approach employs a finite-element basis function and a finite-difference local mass balance, and is designed to reduce computational and storage requirements, while minimizing the mass balance error. The computational grid is formed by division of the flow domain into triangular elements. An invented node at the center of each element divides the element into three subtriangular regions. By linking the center of each triangular element and the mid-point of each elemental side, a multiangular region, referred to as an exclusive subdomain, is defined. The discretized system of governing equations is derived from the integral form that describes the mass balance in the exclusive subdomain of each node. The paper details the application of the numerical method, and demonstrates that the method is reasonably accurate and computationally efficient for a two-dimensional domain subject to nonequilibrium sorption.     

Root establishment of perennial ryegrass (L. perenne) in diesel contaminated subsurface soil layers

The efficiency of rhizosphere biodegradation of petroleum hydrocarbons heterogeneously distributed in soils is dependent on the ability of plant roots to prospect into contaminated zones. Rhizobox experiments were conducted to study the influence of diesel contaminated layers on the spatial distribution and the development of the roots of perennial ryegrass. Root distribution and root and shoot development were monitored over time. The final root and above ground biomass and the final TPH concentration were determined. The spatial distribution of the contaminant as well as the irrigation method used affected root distribution, plant development and TPH degradation and therefore ryegrass remediation potential. The results show that roots colonise fully uncontaminated soil and grow preferentially between zones of contamination. Conversely, when no immediate uncontaminated soil is available, roots grow through contaminated zones in order to prospect for uncontaminated soil.