Bioremediation With Oxygen-Releasing Solids: A Mathematical Model

The use of oxygen releasing solids in passive wells has become of considerable interest in the bioremediation of groundwater contaminated with hydrocarbons and other biodegradable contaminants. Labor, operating, and maintenance costs of this technology are quite competitive with more conventional techniques in many instances. Because of the small rates of mass transport by transverse dispersion, however, a pin-stripe pattern of alternating contaminant and oxygen distribution is expected downgradient of a line of wells which is unacceptable if it extends beyond the point of compliance downgradient. A mathematical model is presented which permits determination of the maximum well spacing in an array of wells which will yield effective remediation within a specified distance of a plume of contaminated groundwater passing through the array. Transverse dispersion of both oxygen and the contaminant is assumed, and biodegradation is modeled by means of Monod kinetics. Longitudinal numerical dispersion is minimized by the use of a asymmetrical upwind algorithm for modeling advection. Modeling results are presented showing the dependence of model behavior on input parameters.     

Temporal Moment Analysis of Multi-Species Radionuclide Transport in a Coupled Fracture-Skin-Matrix System with a Variable Fracture Aperture

This paper presents an analysis using temporal moments to study multi-species radionuclide transport along a single fracture with variable fracture aperture in a fracture-skin-matrix system. In the present study, a decay chain having three elements is considered and transport of each member in the decay chain is modeled by solving a coupled system of partial differential equations for fracture, fracture-skin, and rock matrix using explicit finite difference method. Having obtained the concentration distribution, lower order temporal moments of radionuclide distribution are computed to analyze the effective velocity and macro-dispersion coefficient of radionuclides in the fracture. In the present study, effect of varying fracture aperture on transport characteristics of radionuclides is also analyzed. It is found that the fracture aperture variation profile has significant impact on radionuclide distribution along the fracture. Sensitivity analysis is carried out in a fracture-skin-matrix system with sinusoidal fracture aperture to study the effect of various parameters like fracture aperture thickness, fracture-skin porosity, fracture-skin diffusion coefficient, flow rate, radioactive decay constant, and Freundlich sorption isotherm exponent on transport characteristics of radionuclide. The results suggest that the above mentioned parameters significantly influence concentration, distribution, mobility, and effective macro-dispersion coefficient of radionuclides along the fracture.     

Respirometric Evaluation by Graphical Analysis for Microbial Systems

Analysis of respirometric graphs for system diagnosis has recently drawn much attention due to the nature of the test being a fast method with reliable oxygen uptake (Ou) data. In this study, a new graphical analysis procedure is proposed, using the oxygen uptake rate (OUR) vs. Ou respirogram because of its distinct characteristics that provides easy visual inspection. The OUR plot is typically divided into two respiration phases, and each one of these phases can be simulated by a previously derived analytical equation. Based on the 2-phase model, five characteristic parameters (CPs) were derived for the OUR respirogram. Mathematically the 5 CPs can be related to the system parameters (SPs), including the four Monod kinetic parameters and two operating conditions (initial substrate and biomass concentration). Sensitivity analyses were then conducted to assess the 5 CPs for the level of effect caused by system parameters (SPs). This study also conducted respirometric tests under batch operation to evaluate the applicability of the graphical analysis method. Based on the 5 CPs estimated from the OUR respirogram of the tests, it can be concluded that respirometric tests normally experience the problems associated with initial lag, non-uniform seeding, and inadequate mixing.     

Use of microbial activity measurements for monitoring RBC biofilms

Fixed biomass technologies, such as rotating biological contactors or biodiscs, have been applied for wastewater depuration both in large and medium–small-sized agglomerations. Biofilm’s composition and microorganism activity are essential parameters for the successful operation and control of these systems. Biofilm’s thickness and total dry weight have been widely used for biofilm’s characterization but, actually, are not sufficient to describe biofilm activity. In fact, biofilm’s activity is not proportional to the quantity of fixed biomass, but increases with the thickness of biofilm up to a determined level called the “active thickness”. Above this level, the diffusion of nutrients through the film becomes a limiting factor. A stable, thin, and active biofilm thus offers numerous advantages in water and wastewater treatment. Different parameters have been used to evaluate biofilm’s activity. The specific oxygen uptake rate, INT-dehydrogenase activity, and the ratio active/total cells have been applied for that purpose. These methods are not only simple and rapid but also sensitive, precise, and representative. The results obtained confirm the potential of the microbial activity measurements studied for an accurate biofilm’s characterization and biomass activity estimation in fundamental research and for the practical operation and control of fixed biomass depuration systems.     

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 Conceptual Nitrate Transport Model and Its Application at Different Scales

Groundwater and surface water pollution from nonpoint sources are problems of increasing concern. Nitrate transport is a crucial part in the whole process of nitrate loss. Previous methods for simulating the nitrate transport process usually involved solving complex advection–dispersion equations and required large amount of parameters and input data, which makes a model application expensive or invalid in terms of computation and data availability. The aim of this research is to develop a simple conceptual method to simulate the nitrate transport. In this paper, a simple conceptual method utilizing the distribution of nitrate loss against cumulative drainage to simulate the nitrate transport process is introduced. The concept and construction of this method are presented, and its applicability at different spatial scales is demonstrated. Results of model application at laboratory, pilot, and catchment scales show that this method has a good capability to simulate the transport of nitrate. Impact of spatial scales on parameters’ values is also discussed. Compared to conventional methods, this method requires fewer parameters, less computing time, and less input data.     

Kinetic parameters for the removal of lead and chromium from wastewater using activated carbon developed from fertilizer waste material

The waste slurry generated in fertilizer plants in India has been converted into a cheap carbonaceous adsorbent material. The prepared adsorbent has been characterised and used for the removal of lead and chromium metals. The kinetics of adsorption and the extent of adsorption at equilibrium are dependent on the physical and chemical characteristics of the adsorbate, adsorbent and experimental system. Results of laboratory scale studies conducted to delineate the effect of such parameters on the kinetics of adsorption of metal ions are reported. Parameters evaluated include: hydronium ion concentration, temperature, initial adsorbate concentration, size of adsorbent, and amount of adsorbent. On the basis of these studies the various physical parameters such as effective diffusion coefficient, activation energies and entropy of activation are evaluated, as these provide some information regarding the mechanistic aspects. Mass transfer coefficient values suggest a rapid transport of the adsorbate from bulk to solid phase.     

Estimating overall persistence and long-range transport potential of persistent organic pollutants: a comparison of seven multimedia mass balance models and atmospheric transport models

Two different approaches to modeling the environmental fate of organic chemicals have been developed in recent years. The first approach is applied in multimedia box models, calculating average concentrations in homogeneous boxes which represent the different environmental media, based on intermedia partitioning, transport, and degradation processes. In the second approach, used in atmospheric transport models, the spatially and temporally variable atmospheric dynamics form the basis for calculating the environmental distribution of chemicals, from which also exchange processes to other environmental media are modeled. The main goal of the present study was to investigate if the multimedia mass balance models CliMoChem, SimpleBox, EVn-BETR, G-CIEMS, OECD Tool and the atmospheric transport models MSCE-POP and ADEPT predict the same rankings of the overall persistence (P(ov)) and long-range transport potential (LRTP) of POPs, and to explain differences and similarities between the rankings by the mass distributions and inter-compartment mass flows. The study was performed for a group of 14 reference chemicals. For P(ov), the models yield consistent results, owing to the large influence of phase partitioning parameters and degradation rate constants, which are used similarly by all models. Concerning LRTP, there are larger differences between the models than for P(ov), due to different LRTP calculation methods and spatial model resolutions. Between atmospheric transport models and multimedia fate models, no large differences in mass distributions and inter-compartment flows can be recognized. Deviations in mass flows are mainly caused by the geometrical design of the models.     

A Stage-Structured Pest-Natural Enemy Dynamics with Holling Type-II Interaction and Maturation Delay for Pest Species

In this paper, a pest-natural enemy dynamics with maturation delay for pest species is proposed. Here, pest is categorized into two stages, namely immature and mature and the natural enemy only harvest mature pest population with Holling type-II interaction. The boundedness and positivity of the solution of the system are presented. The asymptotic behavior of the system is studied and analyzed for all feasible equilibrium points. Sensitivity analysis of the system at interior equilibrium point for the system parameters is performed, and respective sensitive indices of the variables are identified. Finally, numerical simulations are presented to support our analytic results.

     

Mathematical Modeling of Column and Field Dense Nonaqueous Phase Liquid Tracer Tests

Mathematical models for the simulation of dense nonaqueous phase liquid tracer tests (DTTs) in laboratory columns and in the field are developed and examined. The DTT technique is a means of estimating the quantity of dense nonaqueous phase liquid (DNAPL) in a domain of interest in an aquifer. The two-dimensional field DTT model uses the Method of Principal Directions and an asymmetrical upwind algorithm for describing advective transport. Both models include diffusion transport of tracer into and from low-permeability porous structures such as clay lenses, as well as the mass transport kinetics of partitioning tracer to and from the DNAPL droplets. The dependence of the effluent tracer concentration curves on the parameters of the models is explored, and conclusions are drawn regarding the applicability of, and several possible problems with, the DTT technique. Model results indicate that the DTT performs well at locating distributed droplets of DNAPL, but is unlikely to be useful in the assessment of pooled DNAPL.     

Calibration and verification of a dissolved oxygen management model for a highly polluted river with extreme flow variations in Pakistan

A dissolved oxygen (DO) model is calibrated and verified for a highly polluted River Ravi with large flow variations. The model calibration is done under medium flow conditions (431.5 m³/s), whereas the model verification is done using the data collected during low flow conditions (52.6 m³/s). Biokinetic rate coefficients for carbonaceous biochemical oxygen demand (CBOD) and nitrogenous biochemical oxygen demand (NBOD) (i.e, K _cr and K _n ) are determined through the measured CBOD and ammonia river profiles. The calculated values of K _cr and K _n are 0.36 day^?1 and 0.34 day^?1, respectively. The close agreement between the DO model results and the field values shows that the verified model can be used to develop DO management strategies for the River Ravi. The biokinetic coefficients are known to vary with degree of treatment (DOT) and therefore need to be adjusted for a rational water quality management model. The effect of this variation on level of treatment has been evaluated by using the verified model to attain a DO standard of 4 mg/L in the river using the biokinetic rate coefficients as determined during the model calibration and verification process. The required DOT in this case is found to be 96 %, whereas the DOT is 86 % if adjusted biokinetic rate coefficients are used to reflect the effect of wastewater treatment. The cost of wastewater treatment is known to increase exponentially as the removal efficiency increases; therefore, the use of appropriate biokinetic coefficients to manage the water quality in rivers is important.     

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).     

Almond shell activated carbon: adsorbent and catalytic support in the phenol degradation

In this work, two technologies are studied for the removal of phenol from aqueous solution: dynamic adsorption onto activated carbon and photocatalysis. Almond shell activated carbon (ASAC) was used as adsorbent and catalytic support in the phenol degradation process. The prepared catalyst by deposition of anatase TiO₂ on the surface of activated carbon was characterized by scanning electron microscopy, sorption of nitrogen, X-ray diffraction, Fourier transform infrared (FT-IR) spectroscopy, and pH_ZPC point of zero charge. In the continuous adsorption experiments, the effects of flow rate, bed height, and solution temperature on the breakthrough curves have been studied. The breakthrough curves were favorably described by the Yoon–Nelson model. The photocatalytic degradation of phenol has been investigated at room temperature using TiO₂-coated activated carbon as photocatalyst (TiO₂/ASAC). The degradation reaction was optimized with respect to the phenol concentration and catalyst amount. The kinetics of disappearance of the organic pollutant followed an apparent first-order rate. The findings demonstrated the applicability of ASAC for the adsorptive and catalytic treatment of phenol.     

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.     

Sensitivity of Simulated Conservation Practice Effectiveness to Representation of Field and In-Stream Processes in the Little River Watershed

Evaluating the effectiveness of conservation practices (CPs) is an important step to achieving efficient and successful water quality management. Watershed-scale simulation models can provide useful and convenient tools for this evaluation, but simulated conservation practice effectiveness should be responsive to parameter values used to represent the practices in the modeling. The objectives of this study were to (1) assess the impacts of a set of conservation practices on hydrology and water quality of a watershed and (2) evaluate the sensitivity of Soil and Water Assessment Tool (SWAT) modeling outputs and simulated conservation practice effectiveness to parameters. The modeling study was conducted in an agricultural watershed, the subwatershed K (16.9 km²) of the Little River Experimental watershed located in the South Atlantic Coastal Plain of the USA. Sensitivity analysis showed that hydrologic response unit (HRU) and watershed-scale simulations for water quality were most sensitive to CN and FILTERW parameters. Load reduction rates as a function of increased aerial coverage of the conservation practices were greatest for total phosphorus (TP), followed by sediment and total nitrogen (TN). The results indicated that conservation practices would have a limited impact on stream flow volume but could have a significant impact on sediment and TP loads within this region. Watershed-scale TN and TP loads were also sensitive to an in-stream nutrient transformation process represented using the QUAL2E algorithm in SWAT. The study clearly demonstrated the most sensitive model parameters and the optimal conservation practices for this watershed.     

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.     

Sorption and desorption of benzo(a)pyrene in aquatic systems

The sorptive behaviour of polycyclic aromatic hydrocarbons (PAHs) is critical to controlling their transport, fates and effects in the environment. Experiments are described which detail the behaviour of a model compound (benzo(a)pyrene) under simulated aquatic conditions. The kinetics of sorption were comparable to those of other PAH compounds. The sorption equilibrium was extensively studied and found to be affected by several key parameters, notably sediment concentration. The sorption coefficient decreased substantially with the sediment concentration, from 9580 ml g(-1) at a sediment concentration of 0.067 g l(-1) to 1,110 ml g(-1) at a sediment concentration of 9.8 g l(-1). The results are consistent with previous reports and often explained by the presence of colloids. In this paper the dry weight concentration of colloids was determined and used for deriving the true sorption coefficient, which is up to an order of magnitude higher than the observed partition coefficient. The sorption of benzo(a)pyrene was also dependent on some of the particle properties, and the sorption coefficient was found to increase with the organic carbon content and specific surface area of sediment particles. The desorption of benzo(a)pyrene from sediment was shown to be relatively rapid, with implications for the potential remobilisation of benzo(a)pyrene and similar compounds.     

Effect of sublethal concentrations of phenol on some enzyme activities and blood sugar level of carp (Cyprinus carpio L.)

To measure the phenol induced tissue necrosis and its adverse effects on nervous system, GOT, GPT, LDH, AChE enzyme activities and blood sugar level detections were used under aquarium conditions and in field experiments.Under aquarium conditions all experimental animals survived 5 ppm phenol concentration even at 96 hr exposure time. However the increase of plasma GOT, GPT, and LDH activities reflected serious tissue damage.In field experiments the plasma GOT. GPT, and LDH activities three times elevated parallel with the increasing phenol contamination of the water.The measurement of our biochemical parameters appeared to be a sensitive biomonitoring indicators to assess exposure of fish to antropogenic agents in the natural environment.     

垃圾渗出液中酚类物质的气相色谱分析

建立了一种用气相色谱法准确、快速测定垃圾渗出液中酚类物质的新方法。外标法定量分析 4 -甲基苯酚及苯酚。采用填装了 80～ 1 0 0目 Chromosorb WAW DCMS的玻璃填充柱 ,担体上涂渍了 2 % OV-1 7,2 .5% QF-1的固定液 ,分离酚类物质效果较好 ,线性范围为 1 0 0～ 1 50 0 mg/ L,相关系数为 0 .9996,方法回收率为 94 .7%。     

Transient modeling of trickling filters for biological ammonia removal

A dynamic model for nitrifying trickling filters is developed, based on material balances in the biofilm and the bulk liquid. The model predicts the profile of ammonia as a function of the operating parameters (volumetric flow rate and feed ammonia concentration) and the biofilm thickness as a function of filter depth and time. This revised version was published online in July 2006 with corrections to the Cover Date.