Modeling the flow and leachate transport in the vadose and saturated zones of a municipal landfill

The purpose of the present study was to evaluate the groundwater contamination due to the construction and operation of the municipal landfill of Amari, Rethymno, Crete. The groundwater flow and leachate transport in the vadose and saturated zones were studied and simulated, using three different models: the one-dimensional groundwater flow and contaminant transport model for the vadose zone Pesticide Root Zone Model (PRZM-3), the Geographical Information System (GIS) Argus ONE and the three-dimensional groundwater flow and contaminant transport model Princeton Transport Code (PTC). The simulation time was 30 and 20 years, and the results obtained, according to the models and the existing hydrogeological conditions, were very encouraging and reassuring about the groundwater quality of the broad region.     

Associations between distance lags,groundwater velocities,and detection efficiencies in groundwater monitoring networks

Effects of distance lags between landfills and monitoring wells on contaminant detection capability were quantified in several groundwater velocity settings. Detection efficiency calculations were made with and without imposing a time limit on contaminant travel. In general, longer distance lags yieldedhigher detection efficiencies. However, detection efficienciesdecreased as monitoring wells approached a buffer zone boundaryimposing a maximum permissible contaminant transport distance.Imposing a time limit on contaminant travel substantially reduced detection efficiency in low velocity settings, especiallyat longer distance lags. Time limits were less significant in high velocity settings where contaminants more quickly reachedmonitoring wells. Detection efficiencies also decreased as velocity increased, but decreases were minor once the velocityreached a threshold value.     

Numerical Modeling on the Effect of Dissolved Oxygen on Nitrogen Transformation and Transport in Unsaturated Porous System

Nitrogen pollution in groundwater resulting from wastewater application to land is a common problem, and it causes a major threat to groundwater-based drinking water supplies. In this study, a numerical model is developed to study the nitrogen species transport and transformation in unsaturated porous media. Further, a new mass transfer module for dissolved oxygen (DO) is incorporated in the one-dimensional numerical model for nitrogen species transport to describe the fate and transport of nitrogen species, dissolved oxygen, dissolved organic carbon (DOC), and biomass. The spatial and temporal variation of dissolved oxygen is incorporated in the model through the mass transfer from gaseous phase to water phase in an unsaturated porous system. The numerical results of the water flow model and single species and multispecies transport model in an unsaturated zone developed for this purpose have been validated with the available analytical/numerical solution. The developed model is applied in clay loam, silt, and sand soils to analyze the transport behavior of nitrogen species under unsaturated condition. The numerical results suggest that the high rate of oxygen mass transfer from the air phase to the water phase positively increases the dissolved oxygen in the applied wastewater and enhances the nitrification process. Because of this high oxygen mass transfer, the nitrate nitrogen concentration significantly increases in the unsaturated zone and the same is transported to a larger depth at higher simulation period. On the other hand, the low rate of oxygen mass transfer implicitly enhances the denitrification process and finally reduces the nitrate nitrogen concentration in the unsaturated zone. The numerical results also show that the nitrate nitrogen transport is rapid in sandy soil when compared with clay loam and silty soils under high oxygen mass transfer rate. In essence, the high oxygen mass transfer rate significantly increases the nitrate nitrogen in the unsaturated zone, especially at a greater depth at larger time levels and eventually affects the groundwater quality.     

Colloid-Associated Groundwater Contaminant Transport in Homogeneous Saturated Porous Media: Mathematical and Numerical Modeling

During the past two decades, significant efforts have been made to study contaminant transport in the presence of colloids. Several researchers reported that colloidal particles could enhance the migration of contaminants in groundwater by reducing retardation factor. When the colloidal particles are present in the aquifer, the subsurface system can be considered as a three-phase system with two solid phases and an aqueous phase. The interaction between contaminants, colloids, and solid matrix should be considered in assessing the fate and transport of the contaminant in the groundwater flow system. In this study, a one-dimensional numerical model is developed by employing a fully implicit finite difference method. This model is based on mass balance equations and mass partition mechanisms between the carriers and solid matrix, as well as between the carriers and contaminants in a saturated homogeneous porous medium. This phenomenon is presented by two approaches: equilibrium approach and fully kinetic first-order approach. The formulation of the model can be simplified by employing equilibrium partitioning of particles. However, contaminant transport can be predicted more accurately in realistic situations by kinetic modeling. To test the sensitivity of the model, the effect of the various chemical and physical coefficients on the migration of contaminant was investigated. The results of numerical modeling matched favorably with experimental data reported in the literature.     

Linear alkylbenzene sulfonate (LAS) in water of Lake Dianchi—spatial and seasonal variation, and kinetics of biodegradation

Agroforesty systems, which are recommended as a management option to lower the shallow groundwater level and to reuse saline subsurface drainage waters from the tile-drained croplands in the drainage-impacted areas of Jan Joaquin Valley of California, have resulted in excessive boron buildup in the soil root zone. To assess the efficacy of the long-term impacts of soil boron buildup in agroforesty systems, a mathematical model was developed to simulate non-conservative boron transport. The developed dynamic two-dimensional finite element model simulates water flow and boron transport in saturated–unsaturated soil system, including boron sorption and boron uptake by root-water extraction processes. The simulation of two different observed field data sets by the developed model is satisfactory, with mean absolute error of 1.5 mg/L and relative error of 6.5%. Application of the model to three different soils shows that boron adsorption is higher in silt loam soil than that in sandy loam and clay loam soils. This result agrees with the laboratory experimental observations. The results of the sensitivity analysis indicate that boron uptake by root-water extraction process influences the boron concentration distribution along the root zone. Also, absorption coefficient and maximum adsorptive capacity of a soil for boron are found to be sensitive parameters.     

Two-dimensional finite elements model for selenium transport in saturated and unsaturated zones

A two-dimensional finite element model was developed to simulate species of selenium transport in two dimensions in both saturated and unsaturated soil zones. The model considers water, selenate, selenite, and selenomethionine uptake by plants. It also considers adsorption and desorption, oxidation and reduction, volatilization, and chemical and biological transformations of selenate, selenite, and selenomethionine. In addition to simulating water flow, selenate, selenite, and selenomethionine transport, the model also simulates organic and gaseous selenium transport. The developed model was applied to simulate two different observed field data. The simulation of the observed data was satisfactory, with mean absolute error of 48.5 μg/l and mean relative error of 8.9%.     

Temporal moments analysis of preferential solute transport in soils

Temporal moments analysis of solute breakthrough curves is used to investigate the preferential leaching of chloride, nitrate and phosphate through an Australian soil. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals, and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and non-equilibrium between soil constituents, water and solutes. Therefore simple models are required to accurately characterise solute transport in natural and agricultural soils under non-equilibrium conditions. A multiple sample percolation system, consisting of 25 individual collection wells was constructed to study the effects of localised soil heterogeneities on the transport of nutrients (NO₃ ^–, Cl^–, PO₄ ^{3 –}) in the vadose zone of an agricultural soil predominantly dominated by clay. Using data collected from the multiple sample percolation experiments, this paper compares and contrasts the performance of temporal moments analysis with two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR) and a two-region preferential flow model (TRM) suitable for modelling preferential transport. The values for solute transport parameters predicted by temporal moments analysis were in excellent agreement with experimental data and results from ADR and TRM. It is concluded that temporal moments analysis when applied with other physical models such as the ADR and TRM, provide an excellent means of obtaining values for important solute transport parameters and gaining insight of preferential flow. These results have significant ramifications for modelling solute transport and predicting nutrient loadings.     

Augmenting Groundwater Monitoring Networks near Landfills with Slurry Cutoff Walls

This study investigated the use of slurry cutoff walls in conjunction with monitoring wells to detect contaminant releases from a solid waste landfill. The 50 m wide by 75 m long landfill was oriented oblique to regional groundwater flow in a shallow sand aquifer. Computer models calculated flow fields and the detection capability of six monitoring networks, four including a 1 m wide by 50 m long cutoff wall at various positions along the landfill's downgradient boundaries and upgradient of the landfill. Wells were positioned to take advantage of convergent flow induced downgradient of the cutoff walls. A five-well network with no cutoff wall detected 81% of contaminant plumes originating within the landfill's footprint before they reached a buffer zone boundary located 50 m from the landfill's downgradient corner. By comparison, detection efficiencies of networks augmented with cutoff walls ranged from 81 to 100%. The most efficient network detected 100% of contaminant releases with four wells, with a centrally located, downgradient cutoff wall. In general, cutoff walls increased detection efficiency by delaying transport of contaminant plumes to the buffer zone boundary, thereby allowing them to increase in size, and by inducing convergent flow at downgradient areas, thereby funneling contaminant plumes toward monitoring wells. However, increases in detection efficiency were too small to offset construction costs for cutoff walls. A 100% detection efficiency was also attained by an eight-well network with no cutoff wall, at approximately one-third the cost of the most efficient wall-augmented network.     

Simulating Nitrate Leaching Profiles in a Highly Permeable Vadose Zone

An approach is developed to simulate leaching of a dissolved chemical constituent in the vadose zone of an aquifer. Specifically, nitrate loading at the water table for different water table depths, for a range of aquifer permeability values, and for different cases of heterogeneity of the aquifer, are considered. Models from the literature are first used to derive soil–water characteristic curves (water retention and hydraulic conductivity) from a grain size distribution curve for unsaturated conditions. Given infiltration from the surface, the initial conditions for the chemical concentration, and the water content profile, leaching of the chemical in the vadose zone is simulated as a function of both time and depth. The methodology is illustrated for a permeable aquifer. Simulations are undertaken using a finite element code for saturated and unsaturated flow. Different scenarios are simulated depending on the heterogeneity of the aquifer and the depth of the water table. Modeling results show that in the example case studied, nitrate concentration loading at the water table does not depend strongly on the position of the water table, but rather on the material properties of the aquifer. The contribution of this endeavor resides in the methodology which allows a prediction of nitrate leaching using only the grain size property of the aquifer. It allows practitioners to obtain a first assessment of leaching with limited data.     

Saturated and Unsaturated Water Flow in Inclined Porous Media

This paper considers the two-dimensional saturated and unsaturated flow of water through inclined porous media, namely a waste dump or hill slope. Since the partial differential equation governing this water flow transforms from being parabolic to elliptic as the water flow varies from unsaturated to saturated, an iterative, finite differencing scheme is used to develop a numerical solution. The model can be used to investigate the effects that hill slope angle, depth of soil cover and hilltop width have on water accumulation in the dump and the time required for saturation to occur at different areas in the dump domain. The accuracy and reliability of the computer based solution is tested for two different boundary conditions – (1) no flow on all boundaries (i.e., the internal redistribution of soil moisture to steady state) and (2) a constant rainfall flux on the dump surface. Numerical studies then show the effects of changing the hill slope angle, depth of layer, and dump geometry on the flow characteristics in the dump.     

Continuous ‘Passive’ flow-proportional monitoring of drainage using a new modified Sutro weir (MSW) unit

In view of their crucial role in water and solute transport, enhanced monitoring of agricultural subsurface drain tile systems is important for adequate water quality management. However, existing monitoring techniques for flow and contaminant loads from tile drains are expensive and labour intensive. The aim of this study was to develop a cost-effective and simple method for monitoring loads from tile drains. The Flowcap is a modified Sutro weir (MSW) unit that can be attached to the outlet of tile drains. It is capable of registering total flow, contaminant loads and flow-averaged concentrations. The MSW builds on a modern passive sampling technique that responds to hydraulic pressure and measures average concentrations over time (days to months) for various substances. Mounting the samplers in the MSW allowed a flow-proportional part of the drainage to be sampled. Laboratory testing yielded high linear correlation between the accumulated sampler flow, q _total, and accumulated drainage flow, Q _total (r ²?>?0.96). The slope of these correlations was used to calculate the total drainage discharge from the sampled volume, and therefore contaminant load. A calibration of the MSW under controlled laboratory condition was needed before interpretation of the monitoring results was possible. The MSW does not require a shed, electricity, or maintenance. This enables large-scale monitoring of contaminant loads via tile drains, which can improve contaminant transport models and yield valuable information for the selection and evaluation of mitigation options to improve water quality. Results from this type of monitoring can provide data for the evaluation and optimisation of best management practices in agriculture in order to produce the highest yield without water quality and recipient surface waters being compromised.     

Treatment of Urban Areas within a Regional Transport Model of Sulphur and Nitrogen Oxides

The air pollution transport model UGEM (The University of Greenwich Evaluation Model) has been developed to evaluate medium-range transport and deposition of sulphur and oxidised nitrogen from all types of sources of emissions in the UK and to estimate their average annual deposition and concentrations across the UK. The model has been tested for its predictions against the available measurements.This study was focused on a possibility of applying the UGEM model to the assessment of air quality on a local scale. One parameter in the model is crucial, the local deposition fraction. The effect of this parameter on quality of the model predictions has been studied for different scales of UGEM output, such as the whole territory of the UK, a rural region and an urban area.The results of the study show that the magnitude of the local deposition fraction should be different for each grid square to reach the best agreement of predictions of concentrations with measurements. Applying a local value of the parameter to each grid square will improve the model predictions of the concentrations in urban areas in particular and will not affect the quality of model predictions of the wet deposition.     

Highlighting the implications of selenium dispersion from disposal of Kahota Industrial Triangle area, Islamabad, Pakistan using three-dimension solute transport model

This study highlights the implications of selenium (Se) dispersion in groundwater flow regimes of Kahota Industrial Triangle area located adjacent to the Soan River, Islamabad. Initially, a regional groundwater 3-D flow model has been developed, calibrated to the known observed heads of 24 water wells, verified, and confirmed that convergence has actually arrived to satisfy the steady state condition. Later, the transient simulation was carried out adding in the known recharge, storage factor, porosity, and observed drawdown matched with the simulated drawdown that appears to fall in close agreement with a difference of 0.25 m. As such the steady state groundwater model has facilitated to understand the mechanism of groundwater flow regimes in reference to the implications of selenium dispersion from disposal of Kahota Industrial Triangle area. Thirty-five water samples were collected mainly from the industrial water wells for the evaluation of heavy metals. Selenium being the major contributor of pollution has been short listed to monitor its dispersion using a solute transport model modular three-dimensional transport model (MT3D). Chemical parameters related to selenium characteristics including horizontal and vertical transverse dispersivity/longitudinal dispersivity, effective molecular diffusion coefficient and bulk density of the porous medium of aquifers have been used in MT3D contaminant transport model. MT3D is run for 30 years in steady state condition. As usual first run did not produce the exact field conditions. Therefore, the contaminant transport model is calibrated against the 32 values of observed selenium concentrations in boreholes by minor adjustments in the chemical parameter values. The final calibration has been achieved with residual value of 3.88 × 10^???5 Kg/m³. Seven hypothetical observation wells are used to monitor the selenium concentrations over a long-term period of time.     

Managing Uncertainty in Risk-Based Corrective Action Design: Global Sensitivity Analysis of Contaminant Fate and Exposure Models Used in the Dose Assessment

A variance-based global sensitivity analysis (GSA) was applied to the dose assessment model used in the risk-based corrective action methodology of environmental risk analysis to identify key sources of variability and uncertainty and quantify the relative contribution of these sources to the variance of estimated dose. GSA was performed applying extended Fourier amplitude sensitivity test technique. The soil-to-air contaminant transport pathway within an inhalation exposure scenario was addressed. Three persistent semi-volatile carcinogenic chemicals, including polychlorinated biphenyls, benzo(a)pyrene, and 2,3,7,8-tetrachlorodibenzo-p-dioxin, were chosen as contaminants of concern.     

Modelling of spatial contaminant probabilities of occurrence of chlorinated hydrocarbons in an urban aquifer

In this study, a 3D urban groundwater model is presented which serves for calculation of multispecies contaminant transport in the subsurface on the regional scale. The total model consists of two submodels, the groundwater flow and reactive transport model, and is validated against field data. The model equations are solved applying finite element method. A sensitivity analysis is carried out to perform parameter identification of flow, transport and reaction processes. Coming from the latter, stochastic variation of flow, transport, and reaction input parameters and Monte Carlo simulation are used in calculating probabilities of pollutant occurrence in the domain. These probabilities could be part of determining future spots of contamination and their measure of damages. Application and validation is exemplarily shown for a contaminated site in Braunschweig (Germany), where a vast plume of chlorinated ethenes pollutes the groundwater. With respect to field application, the methods used for modelling reveal feasible and helpful tools to assess natural attenuation (MNA) and the risk that might be reduced by remediation actions.     

A method for monitoring ground water quality near waste storage facilities

A practical optimization approach developed in this paper derives effective monitoring configurations for detecting contaminants in ground water. The approach integrates numerical simulation of contaminant transport and mathematical programming. Well sites identified by the methodology can be monitored to establish the occurrence of a contaminant release before a plume migrates to a regulatory compliance boundary. Monitoring sites are established along several horizons located between the downgradient margin of a contaminant source and a compliance boundary. A horizon can form an effective line of defense against contaminant migration to the compliance boundary if it is spanned (covered) by a sufficient number of sites to yield a well spacing that is equal to or less than a maximum value established by numerical modeling. The objective function of the integer programming model formulation expresses the goals of: (1) covering a maximum number of siting horizons, and (2) allocating wells to the single most effective horizon. The latter is determined from well spacing requirements and the width of the zone of potential contaminant migration traversed by the horizon. The methodology employs a highly tractable linear programming model formulation, and the user is not required to predefine a set of potential well sites. These attributes can facilitate its implementation in practice.     

An Experimental Methodology for Monitoring Contaminant Transport Through Geotechnical Centrifuge Models

In this paper, an attempt has been made to highlight an experimental methodology for monitoring contaminant transport through locally available silty soil and commercially available clay in geotechnical centrifuge models, for different compaction states. Use of multiple depth sensors to determine depth distribution of sodium chloride in the soil column has been detailed. The obtained results have been compared with argentometric method. To validate the centrifuge modelling, modelling of models has been used. The test setup developed can simulate contaminant transport mechanisms through the soil mass, which is approximately 10 m deep, over a period of 600 days. R _e and P _e are found to be N times higher in the centrifuge models. These numbers are found to be several orders less than unity. This indicates that laminar flow prevails and the dominating Cl⁻ transport mechanism in centrifuge is diffusion. The study also highlights the fact that the geotechnical centrifuge modelling can be used as a viable alternative to field scale experimentation.     

Assessment of phenol infiltration resilience in soil media by HYDRUS-1D transport model for a waste discharge site

The movement of contaminants through soil imparts a variety of geo-environmental problem inclusive of lithospheric pollution. Near-surface aquifers are often vulnerable to contamination from surface source if overlying soil possesses poor resilience or contaminant attenuation capacity. The prediction of contaminant transport through soil is urged to protect groundwater from sources of pollutants. Using field simulation through column experiments and mathematical modeling like HYDRUS-1D, assessment of soil resilience and movement of contaminants through the subsurface to reach aquifers can be predicted. An outfall site of effluents of a coke oven plant comprising of alarming concentration of phenol (4–12.2 mg/L) have been considered for studying groundwater condition and quality, in situ soil characterization, and effluent characterization. Hydrogeological feature suggests the presence of near-surface aquifers at the effluent discharge site. Analysis of groundwater of nearby locality reveals the phenol concentration (0.11–0.75 mg/L) exceeded the prescribed limit of WHO specification (0.002 mg/L). The in situ soil, used in column experiment, possess higher saturated hydraulic conductivity (K _S ?=?5.25?×?10^?4 cm/s). The soil containing 47 % silt, 11 % clay, and 1.54 % organic carbon content was found to be a poor absorber of phenol (24 mg/kg). The linear phenol adsorption isotherm model showed the best fit (R ²?=?0.977, RMSE?=?1.057) to the test results. Column experiments revealed that the phenol removal percent and the length of the mass transfer zone increased with increasing bed heights. The overall phenol adsorption efficiency was found to be 42–49 %. Breakthrough curves (BTCs) predicted by HYDRUS-1D model appears to be close fitting with the BTCs derived from the column experiments. The phenol BTC predicted by the HYDRUS-1D model for 1.2 m depth subsurface soil, i.e., up to the depth of groundwater in the study area, showed that the exhaustion point was reached within 12 days of elapsed time. This clearly demonstrated poor attenuation capacity of the soil to retard migration of phenol to the groundwater from the surface outfall site. Suitable liner, based on these data, may be designed to inhibit subsurface transport of phenol and thereby to protect precious groundwater from contamination.     

Geophysical and hydrological data assimilation to monitor water content dynamics in the rocky unsaturated zone

In recent years, geophysics is increasingly used to study the flow and transport processes in the vadose zone. Particularly, when the vadose zone is made up of rocks, it is difficult to install sensors in the subsurface to measure hydrological state variables directly. In these cases, the electrical resistivity tomography (ERT) represents a useful tool to monitor the hydrodynamics of the infiltration and to estimate hydraulic parameters and state variables, such as hydraulic conductivity and water content. We propose an integrated approach aimed at predicting water content dynamics in calcarenite, a sedimentary carbonatic porous rock. The uncoupled hydrogeophysical approach proposed consists in 4D ERT monitoring conducted during an infiltrometer test under falling head conditions. Capacitance probes were installed to measure water content at different depths to validate the estimations derived from ERT. A numerical procedure, based on a data assimilation technique, was accomplished by combining the model (i.e., Richards’ equation) with the observations in order to provide reliable water content estimations. We have used a new data assimilation method that is easy to implement, based on the ensemble Kalman filter coupled with Brownian bridges. This approach is particularly suitable for strongly non-linear models, such as Richards’ equation, in order to take into account both the model uncertainty and the observation errors. The proposed data assimilation approach was tested for the first time on field data. A reasonable agreement was found between observations and predictions confirming the ability of the integrated approach to predict water content dynamics in the rocky subsoil.