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.     

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.     

复杂岩溶条件下锰矿尾矿库地下水溶质运移特征数值模拟研究

受构造、地层等复杂地质情况综合影响的岩溶地区地下水,其溶质运移及流场模拟向来是水文地质、环境地质研究中的难点.以贵州省松桃县某锰矿尾矿库为例,在充分分析研究区周边水文地质、环境地质条件的基础上,建立地质模型,概化边界条件,利用有限元软件FEFLOW进行流场拟合.在此基础上,模拟特征污染物锰及氨氮的弥散情况.模拟结果表明...     

Modeling nonlinearly transient zero-order bioreduction and transport of Cr(VI) in groundwater

Cr(VI) contamination of soil and groundwater is considered a major environmental concern. Bioreduction of Cr(VI) to Cr(III) can be considered a potentially effective technology in remediating Cr(VI) contaminated sites. Shewanella oneidensis MR-1 (MR-1) is one of the bacteria capable of reducing Cr(VI) to Cr(III) under anaerobic conditions. The kinetics of Cr(VI) reduction by MR-1 is defined by the dual-enzyme kinetic model which is nonlinear, transient, and zero-order. Existing transport models are not designed to simulate such reaction kinetics. The objective of this paper is to present a Petrov–Galerkin finite element model (PGFEM) to simulate transport and bioreduction of Cr(VI), by MR-1, in groundwater. The model developed is unconditionally stable and provides oscillation free accurate results for a wide range of Peclet number (Pn) and Courant number (Cn).     

Composite use of numerical groundwater flow modeling and geoinformatics techniques for monitoring Indus Basin aquifer, Pakistan

The integration of the Geographic Information System (GIS) with groundwater modeling and satellite remote sensing capabilities has provided an efficient way of analyzing and monitoring groundwater behavior and its associated land conditions. A 3-dimensional finite element model (Feflow) has been used for regional groundwater flow modeling of Upper Chaj Doab in Indus Basin, Pakistan. The approach of using GIS techniques that partially fulfill the data requirements and define the parameters of existing hydrologic models was adopted. The numerical groundwater flow model is developed to configure the groundwater equipotential surface, hydraulic head gradient, and estimation of the groundwater budget of the aquifer. GIS is used for spatial database development, integration with a remote sensing, and numerical groundwater flow modeling capabilities. The thematic layers of soils, land use, hydrology, infrastructure, and climate were developed using GIS. The Arcview GIS software is used as additive tool to develop supportive data for numerical groundwater flow modeling and integration and presentation of image processing and modeling results. The groundwater flow model was calibrated to simulate future changes in piezometric heads from the period 2006 to 2020. Different scenarios were developed to study the impact of extreme climatic conditions (drought/flood) and variable groundwater abstraction on the regional groundwater system. The model results indicated a significant response in watertable due to external influential factors. The developed model provides an effective tool for evaluating better management options for monitoring future groundwater development in the study area.     

Modelling of lindane transport in groundwater of metropolitan city Vadodara,Gujarat, India

Migration pattern of organochloro pesticide lindane has been studied in groundwater of metropolitan city Vadodara. Groundwater flow was simulated using the groundwater flow model constructed up to a depth of 60 m considering a three-layer structure with grid size of 40?×?40?×?40 m³. The general groundwater flow direction is from northeast to south and southwest. The river Vishwamitri and river Jambua form natural hydrologic boundary. The constant head in the north and south end of the study area is taken as another boundary condition in the model. The hydraulic head distribution in the multilayer aquifer has been computed from the visual MODFLOW groundwater flow model. TDS has been computed though MT3D mass transport model starting with a background concentration of 500 mg/l and using a porosity value of 0.3. Simulated TDS values from the model matches well with the observed data. Model MT3D was run for lindane pesticide with a background concentration of 0.5?μg/l. The predictions of the mass transport model for next 50 years indicate that advancement of containment of plume size in the aquifer system both spatially and depth wise as a result of increasing level of pesticide in river Vishwamitri. The restoration of the aquifer system may take a very long time as seen from slow improvement in the groundwater quality from the predicted scenarios, thereby, indicating alarming situation of groundwater quality deterioration in different layers. It is recommended that all the industries operating in the region should install efficient effluent treatment plants to abate the pollution problem.     

Modeling Solute Transport in Volcanic Ash Soils with Cation Exchange and Anion Retardation

Groundwater under agricultural lands is often contaminated by nitrate. In southern Japan, aquifers are covered by volcanic ash soils that can leach nitrate to the groundwater. In this study, column experiments using two volcanic ash soils (Kuroboku and Akahoya) were carried out. A mixed solution of KNO₃ and K₂SO₄ was used in the leaching experiments. Based on the experiments, a reactive transport model was developed using the Constrained Interpolation Profile method for ion transport and used to calculate chemical equilibrium for the cation exchange reactions. Anion adsorption was modeled by retardation in the numerical model. The developed simulation model results were compared to results obtained by the reactive transport model PHAST. The developed model was shown to quite well reproduce general characteristics of the experimental results. Also, the developed model results agreed well with results from PHAST. Slight discrepancy between observed and calculated breakthrough curves was probably caused by ignoring the kinetic reaction in the model calculations.     

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.     

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.     

A groundwater flow model for water resources management in the Ismarida plain,North Greece

This paper presents the development of a regional flow simulation model of the stream–aquifer system of Ismarida plain, northeastern Greece. It quantifies the water budget for this aquifer system and describes the components of groundwater and the characteristics of this system on the basis of results of a 3-year field study. The semiconfined aquifer system of Ismarida Lake plain consists of unconsolidated deltaic clastic sediments, is hydraulically connected with Vosvozis River, and covers an area of 46.75 km². The annual precipitation ranges in the study area from 270 to 876 mm. Eighty-seven irrigation wells are densely located and have been widely used for agricultural development. Groundwater flow in this aquifer was simulated with MODFLOW. Model calibration was done with observed water levels, and match was excellent. To evaluate the impacts of the current pumping schedule and propose solutions, four management scenarios were formulated and tested with the model. Based on model results, the simulated groundwater budget indicates that there must be approximately 33% decrease of withdrawals to stop the dramatic decline of groundwater levels. The application of these scenarios shows that aquifer discharge to the nearby river would be very low after a 20-year period.     

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.     

Performance of semipermeable membrane devices for sampling of organic contaminants in groundwater

Lipid-filled semipermeable membrane devices (SPMDs) are receiving increasing attention as passive, in situ samplers for the assessment of environmental pollutant exposure. Although SPMDs have been successfully used in a variety of field studies in surface waters, only a few studies have addressed their characteristics as groundwater samplers. In this study, the performance of the SPMDs for monitoring organic contaminants in groundwater was evaluated in a pilot field application in an area severely contaminated by chemical waste, especially by chlorinated hydrocarbons. The spatial distribution of hydrophobic groundwater contaminants was assessed using a combination of passive sampling with SPMDs and non-target semiquantitative GC-MS analysis. More than 100 contaminants were identified and semiquantitatively determined in SPMD samples. Along the 6 field sites under investigation, a large concentration gradient was observed, which confirms a very limited mobility of hydrophobic substances in dissolved form in the aquifer. The in situ extraction potential of the SPMD is limited by groundwater flow, when the exchange volume of well water during an exposure is lower than the SPMD clearance volume for the analytes. This study demonstrates that SPMDs present a useful tool for sampling and analyzing of groundwater polluted with complex mixtures of hydrophobic chemicals and provides guidance for further development of passive sampling technology for groundwater.     

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.     

Future groundwater extraction scenarios for an aquifer in a semiarid environment: case study of Guadalupe Valley Aquifer,Baja California,Northwest Mexico

Semiarid northwestern Mexico presents a growing water demand produced by agricultural and domestic requirements during the last two decades. The community of Guadalupe Valley and the city of Ensenada rely on groundwater pumping from the local aquifer as its sole source of water supply. This dependency has resulted in an imbalance between groundwater pumpage and natural recharge. A two-dimensional groundwater flow model was applied to the Guadalupe Valley Aquifer, which was calibrated and validated for the period 1984–2005. The model analysis verified that groundwater levels in the region are subject to steep declines due to decades of intensive groundwater exploitation for agricultural and domestic purposes. The calibrated model was used to assess the effects of different water management scenarios for the period 2007–2025. If the base case (status quo) scenario continues, groundwater levels are in a continuous drawdown trend. Some wells would run dry by August 2017, and water demand may not be met without incurring in an overdraft. The optimistic scenario implies the achievement of the mean groundwater recharge and discharge. Groundwater level depletion could be stopped and restored. The sustainable scenario implies the reduction of current extraction (up to about 50 %), when groundwater level depletion could be stopped. A reduction in current extraction mitigates water stress in the aquifer but cannot solely reverse declining water tables across the region. The combination of reduced current extraction and an implemented alternative solution (such as groundwater artificial recharge), provides the most effective measure to stabilize and reverse declining groundwater levels while meeting water demands in the region.     

Assessment and Prediction of Contaminant Transport and Migration at a Florida Superfund Site

An important issue in the application of mathematical models to describe or predict the fate of solutes in soil and groundwater systems is providing the necessary data to support the spatial, temporal, and hydrogeologic model parameter requirements. The larger, more complex the model and scope of application, the more likely that significant amounts of data of several types are required. A management tool is required that allows the model and the modeler to interact with various needed databases, cope with numerous database management issues, efficiently achieve model parameterization and provide visualization of model output. The three-dimensional Groundwater Modeling System (GMS) that contains a state-of-the-art flow and solute transport model and a graphical user interface for data manipulation and analysis was applied to a Superfund site in Florida to demonstrate its capabilities for predicting solute migration.     

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.     

Application of drastic model and GIS: for assessing vulnerability in hard rock granitic aquifer

Geographic information system (GIS) has become one of the leading tools in the field of hydrogeological science that helps in assessing, monitoring, and conserving groundwater resources. Groundwater is a finite resource, which is being overexploited due to increase in demand over the years leading to decrease in its potentiality. In the present study, DRASTIC model has been used to prepare groundwater vulnerable zone in hard rock aquifer of granitic terrain. The main objective is to determine susceptible zone for groundwater pollution by integrating hydrogeological layers in GIS environment. The layers such as depth of aquifer, recharge, aquifer yield, soil type, topography, vadose zone, and transmissivity are incorporated in the DRASTIC model. The final output of the map shows that around 60% of the area falls under low to no risk of pollution zone. The high risk of pollution zones are mostly present towards the margin of southeastern periphery. The lower part of the basin as well as small area on northern side falls under moderate risk of pollution zone. For the assessment of groundwater pollution zone, 24 groundwater samples have been collected from different vulnerable zones. The chemical analysis of sample shows that the southeastern margin of basin has relatively high concentration of nitrate as compared to other parts of the basin. It is present in high pollution zone as well as moderate pollution zone. The present model can be used for assessment and management of groundwater.     

Chloride migration in groundwater for a tannery belt in Southern India

Groundwater in a tannery belt in Southern India is being polluted by the discharge of untreated effluents from 80 operating tanneries. Total dissolved solids and chloride (Cl⁻) measurements in open wells in the tannery cluster vary from 27,686 to 39,100 and 12,000 to 13,652 mg/l, respectively. A mass transport model was constructed using Visual MODFLOW Premium 4.4 software to investigate the chloride migration in an area of 75.56 km². Input to the chloride migration model was a groundwater flow model that considered steady and transient conditions. This model was calibrated with field observations; and sensitivity analysis was carried out whereby model parameters, viz., conductivity, dispersivity, and source concentration were altered slightly, and the effect on calibration statistics was evaluated. Results indicated that hydraulic conductivity played a more sensitive role than did dispersivity. The Cl⁻ migration was mainly through advection rather than dispersion. It was found that even if the pollutant load reduced to 50% of the present level, the Cl⁻ concentration in groundwater, even after 6 years, would not be reduced to the permissible limit of drinking water in the tannery belt.     

Modeling photocatalytic destruction of the organic contaminants in continuous stir flow reactor

Thin-film technique is becoming an industry standard for the preparation of TiO₂-based photocatalyst for organic destruction. The catalyst provides several advantages over the conventional powder TiO₂ in the treatment of wastewater and groundwater. In this study, a continuous stir flow reactor model is developed capable of describing the photocatalytic process. The model incorporates the following fundamental mechanisms: adsorption, diffusion, liquid-film transfer, UV attenuation, and photocatalytic reaction. All of the simulation results indicate that there exists a highly nonlinear relation between each of these parameters and the destruction rate. Various incident light intensities also are incorporated to simulate the energy efficiency. The simulations illustrate that the photocatalytic model can be used to elucidate the effect of process variables. It is also possible to “custom-design” a catalyst for the treatment of a particular waste stream.     

Numerical modelling of the environmental impact of landfill leachate leakage on groundwater quality – a field application

Waste disposal facilities are mainly responsible for the gradual quality degradation of subsurface freshwater reservoirs. The main objective of this work is to identify the groundwater contamination risk due to potential leachate leakage and seepage beneath the municipal landfill of the City of Patras in Greece. A groundwater and leachate mass transport model of the underlying aquifer was developed for this purpose. The derived simulation results indicate that, depending on the permeability of the soil at the location of leakage, the contamination risk for the groundwater can be high. In order to quantify the magnitude and the extent of the leakage, a risk assessment model for the leachate contaminant plume was developed. The risk assessment analysis shows that the municipal drinking wells are under high risk of contamination.