SciELO - Scientific Electronic Library Online

vol.40 issue2Water quality of Flag Boshielo Dam, Olifants River, South Africa: Historical trends and the impact of droughtShear rheological properties of fresh human faeces with different moisture content author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Water SA

On-line version ISSN 0378-4738

Water SA vol.40 n.2 Pretoria Apr. 2014




Finite difference simulation of biological chromium (VI) reduction in aquifer media columns



Phalazane J Mtimunye; Evans MN Chirwa

Water Utilisation Division, Department of Chemical Engineering, University of Pretoria, Pretoria 0002, South Africa





A mechanistic mathematical model was developed that successfully traced the Cr(VI) concentration profiles inside porous aquifer media columns. The model was thereafter used to calculate Cr(VI) removal rate for a range of Cr(VI) loadings. Internal concentration profiles were modelled against data collected from intermediate sample ports along the length of the test columns. For the first time, the performance of a simulated barrier was evaluated internally in porous media using a finite difference approach. Parameters in the model were optimised at transient-state and under near steady-state conditions with respect to biomass and effluent Cr(VI) concentration respectively. The best fitting model from this study followed non-competitive inhibition kinetics for Cr(VI) removal with the best fitting steady-state parameters: Cr(VI) reduction rate coefficient, k = 5.2x108 l-mg"1-h"1; Cr(VI) threshold inhibition concentration, C = 50 mg-l-1; and a semi-empirical reaction order, n = 2. The model results showed that post-barrier infusion of biomass into the clean aquifer downstream of the barrier could be limited by depletion of the substrates within the barrier. The model when fully developed will be used in desktop evaluation of proposed in situ biological barrier systems before implementation in actual aquifer systems.

Keywords: in situ bioremediation, porous media reactors, Cr(VI) reduction, finite difference model




C Cr(VI) concentration (ML-3)

Cr Cr(VI) toxicity threshold coefficient (ML-3)

Fc molar flow rate of Cr(VI) (MT-1)

k reaction rate coefficient (L3M-1T-1)

Q inflow rate (L3T-1)

Δc rate of change of Cr(VI) (ML-3T-1)

W mass of the aquifer soil particles (M)

ρ density of aquifer soil particles (ML-3)



Industrial activities such as leather tanning, electroplating, wood preservation/processing, and alloy preparation have led to widespread chromium contamination of soils and ground-water resources around the world (Allen et al., 1998; Shi et al., 2009). The improper disposal of Cr(VI)-containing waste from these industries and its subsequent mobility in ground-water aquifers is a subject of paramount concern (Chirwa and Molokwane, 2011; Blowes, 2002). Aquifer systems below and around industrial and mining sites are especially vulnerable to contamination. Contamination of aquifers has resulted in elevated Cr(VI) levels, sometimes rendering the groundwater unsuitable as a potable water supply source (Egwuonwu et al., 2011; DWAF, 2005; Francisco et al., 2002).

Cr(VI) has been classified by the USEPA since 1978 as a Class A human carcinogen based on its acute toxicity at high concentration and chronic and subchronic toxicity under low exposure conditions (Federal Register, 2004). Deleterious effects resulting from exposure to Cr(VI) include skin dermatitis, irritation of the mucous membrane, and DNA damage or mutations (Cheung and Gu, 2007; Caglieri et al., 2006). Cr(VI) has also been known to hinder seed germination and cause morphological changes in plants (James et al., 1984; Towill, 1978). So far, contamination of groundwater resources around mining areas has been shown to exceed the internationally acceptable exposure limit of 0.05 mg-£-1 in a significant number of groundwater resources (Molokwane et al., 2008; Kobya, 2004; Baral and Engelken, 2002).

Reduction of Cr(VI) to C(III) is favoured since Cr(III) is 1 000 times less toxic and more stable than Cr(VI) in aquatic environments under natural conditions (pH 7, 15-25°C temperature) (Boni and Sbaffoni, 2009). Currently, Cr(VI) contamination in groundwater systems is conventionally treated using the pump-and-treat method which involves the extraction of contaminated water from the aquifer, treatment above ground, and injection of the treated water back into the aquifer (Kurniawan et al., 2006; Watts, 1998). Experience at sites where pump-and-treat remediation of Cr-contaminated groundwater is underway suggests that, although it is feasible to remove high levels of Cr from the subsurface, it becomes more difficult to remove the remaining Cr as concentration decreases (Wittbrodt and Palmer, 1992). Additionally, the method is environmental intrusive in that chemical agents used in the reduction process and hydroxide precipitates of reduced chromium produce toxic sludge that may be difficult to remediate (Gonzalez et al., 2003; Blowes, 2002).

Recent studies have shown the possibility of using locally isolated organisms to treat Cr(VI) contaminated groundwater and/or soil. The treatment can be implemented ex situ utilising a pump-and-treat approach or in situ as a permeable reactive barrier (Molokwane and Chirwa, 2009; Tseng and Bielefeldt, 2002; Turick et al., 1998). In the study by Molokwane and co-workers (2008, 2009), microorganisms isolated from dried digested sludge (from sand drying beds) successfully reduced and immobilised Cr(VI) in extracted aquifer columns in the absence of any added external carbon sources. This demonstrated that a remediation strategy could be implemented using microorganisms without the risk of contaminating the aquifers with external carbon sources and nutrients. The study also demonstrated that locally isolated microorganisms could be adapted to operate at higher Cr(VI) loading conditions thereby eliminating the need for using foreign isolates of bacteria. However, in spite of these new findings, the activity of the organisms inside the columns was not evaluated.

In this study, the model of Cr(VI) removal along the reactor columns is evaluated in detail using intermediate sampling ports installed along the reactor columns. This enabled us to plot the Cr(VI) concentration profiles along the length of the reactor at any given time. The simulation of internal concentration profiles enabled us to evaluate the fundamental physical processes within the soil columns. Without the internal analysis, only a black-box approach may be employed allowing prediction only in a narrow range of operating conditions (Zabiri and Mazuki, 2010; Gomm et al., 1996).



Culture and media

Cr(VI) reducing organisms were collected from sand drying beds of sludge digester effluent at Brits Wastewater Treatment Works (Brits, North West Province, South Africa). The wastewater treatment plant has received periodic loadings of Cr(VI)-contaminated effluent since 1996, thus some of the microorganisms in the sludge were expected to be Cr(VI) tolerant and candidates for Cr(VI) reduction. Isolates capable of reducing Cr(VI) were purified and sub-cultured under micro-aerobic conditions. All the culturable Cr(VI)-reducing species from the original sludge samples were represented as shown in Table 1.



Bacteria was sub-cultured directly from single species culture slants into sterile Luria-Bettani or nutrient broth amended with 75 mg-£-1 Cr(VI) and then incubated at 30±2°C for approximately 24 h. The overnight-grown mixed culture of bacteria was then harvested by centrifugation and washed twice with 0.85% sterile NaCl. Finally, the cells were suspended in sterilised modified basal mineral medium (BMM) prepared according to Roslev et al. (1998) and the suspension was used to inoculate the columns.

Reactor system

Aquifer columns extracted from the actual contaminated site (Brits, North West Province, South Africa) were installed in the laboratory as continuous-flow columns (Fig. 1). The extracted columns were encased in tight fitting Plexiglas cylinders (60 cm long, 5 cm internal diameter). Five equally spaced longitudinal sampling ports were drilled along each column. The columns were capped on both ends with PVC caps with a provision for influent and effluent sampling ports. Before closing the columns on both ends, one of the columns (R3) was sterilised by autoclaving for 20 min at 121°C to serve as an abiotic control. The operational conditions of the experimental columns are shown in Fig. 1.

The aquifer soil columns were then installed vertically on a mounting board. Reservoirs feeding the reactors were interconnected with an overflow to maintain a constant hydraulic head of 76 cm in all columns. The applied hydraulic head corresponded to the average summer water table above the sample collection depth at the target test site. The reactors (R2) and (R4) were each seeded with 100 m£ suspension containing 8 χ 1010 cells of mixed-culture at a very slow rate of 0.28 m£-min-1 for 48 h, enough time to allow uniform distribution and adhesion of cells to soil particles. Cr(VI) loading in the reactor columns was simulated by gravity feeding as in the case of open aquifers at the actual contaminated site. Microcosm reactors were operated as packed-bed reactors at different influent Cr(VI) concentrations (20-50 mg-£-1) corresponding to the seasonal change in concentration at the site. Samples withdrawn from the equally-spaced sampling ports and from the effluent port at the bottom were centrifuged at 6 000 r-min-1 (2 820 χ g) for 10 min to remove cells and soil particles before analysis.

Since the system is being developed for application in a nutrient-deficient groundwater environment, some reactors were operated without added organic carbon source to test performance under autotrophic conditions and some were operated with carbon sources from decaying roots and leachate from decaying vegetation on the surface. Organic sources were simulated by a layer of sawdust (2 cm) to continuously leach organics into the soil column. The leachate was characterised by GC/MS and quantified by TOC Analyser to check for biodegradable organic content from the feed stream.

Analytical methods

Sample pre-concentration and loading

Volatile organics were mobilised from soil and water samples using a Turbo Matrix 40 Headspace (Perkin-Elmer, Connecticut, USA) equipped with a trap. The organics were extracted by boiling at the headspace oven temperature of 80°C for 40 min at a vial pressure of 206.84 kPa (30 psi) and desorbed onto the GC column at 103.42 kPa (15 psi).

Characterisation of organics

Organics extracted in the headspace unit above were desorbed from the trap into the PerkinElmer GC/MS system comprising of a Clarus 600 GC, and a Clarus 600T Mass Spectrometer (MS) (Perkin Elmer, Connecticut, USA) using helium as carrier gas. The GC was equipped with a Perkin Elmer Elite - 5MS capillary column (30 m χ 0.25 mm ID χ 0.5 µηι fixed phase). Injection port temperature was set at 250°C. The column oven program was initiated at 40°C and ramped at 5°C per minute to 150°C and at 15°C per minute to 250°C with a carrier gas split flow of 3 m£-min-1 and a flow rate of 1 m£-min-1. The GC detectors were bypassed in this method. The MS method was operated using electron ionisation (EI+) mode, centroid data setting, with scan duration of 0.3 s and inter-scan delay of 0.2 s. Degradation intermediates were identified by searching for closest matches in the Wiley Mass Spectral Libraries.

Quantification of organics

The total carbon content in the organic leachate from natural organic carbon sources was analysed using the Total Organic Carbon (TOC) Analyser (Model TOC-VWP, Shimadzu Corporation, Kyoto, Japan). The TOC analyser was calibrated by dissolving different proportions of a 1 000 mg-£-1 potassium hydrogen phthalate stock solution in concentrations ranging from 0-5 mg-£-1 in a 100 m£ volumetric flask prior analysing for total carbon.

Cr(VI) and total Cr

Cr(VI) was measured using a UV/Vis spectrophotometer (WPA, Light Wave II, Labotech, South Africa) at a wavelength of 540 nm (10 mm light path) after acidification and reaction with 1,5-diphenyl carbazide to produce a purple-pink colour (APHA, 2005). Total Cr was measured at a wavelength of 359.9 nm using a Varian AA-1275 Series Atomic Adsorption Spectrophotometer (AAS) (Varian, Palo Alto, CA (USA)) equipped with a 3 mA chromium hollow cathode lamp. Cr(III) was determined as the difference between total Cr and Cr(VI) concentration.

Total and viable biomass

Total biomass was measured as a function of volatile solids (VS) in the reactor. Volatile solids were determined as the loss in dry weight of samples of the known volume after igniting the thoroughly dried sample at 600°C in furnace. A conversion factor f = 0.8 was used to convert VS to total biomass (Reynolds and Richards, 1996). Viable cells were determined using the heterotrophic (pour) plate method and colony counts as described in the Standard Methods for the Examination of Water and Wastewater (APHA, 2005), with the colonies grown on Luria-Bettani (LB) and Plate Count (PC) agar. Samples of known mass collected from the sample ports were suspended in 5 m£ sterile BMM and subjected to serial dilution. CFU counts from the top three plates were averaged and the CFU count was converted to mass concentration using a calibration conversion factor of 1.833 χ 10-10 mg-cell-1 (with R2 = 0.997).

Model development

Diffusion process

The diffusivity of Cr(VI) in water was estimated using the Nernst-Haskell equation (Longsworth, 1972):


Dc = diffusion coefficient of Cr(VI) in water (cm2-s-1)

T = absolute temperature (K)

R = universal gas constant (J-mol-1-K)

F = Faraday's constant (C-g-1-mol)

γi = ionic valence (charge) of Cr(VI) (as CrO42-)

λ = electrolytic conductance (cm2-ohm-1-C-1)

The conductance was computed from an equivalent conductance Λ = 1 000 K7C., where K = specific conductance (1 ohm-1-cm-1) estimated from tables of physical properties (Perry, 1997), and Cc = electrolyte charge concentration (equivalent charge-cm-3). A value of Dc = 1.05 χ 10-5 cm2-s-1 was found using a K value of 4.5 χ 10-6 ohm-cm-1 determined from standard 1 molar Cr(VI) solutions in water at 25oC (Garrels and Christ, 1965).

Mass transport of all dissolved species u in water follows Fick's second law of molecular diffusion:


ui = concentration of dissolved species (ML3)

Diw = diffusivity of dissolved species in water (L2Tl)

t = time (T)

z is the spatial coordinate (L) for a 1-dimensional diffusion model in the direction of z

Hydraulic characteristics

The void volume available for free flow of water in the soil column decreased as cells grew around the soil particles. This resulted in a decrease in the effective cross-sectional area [A(t)] and constrained flow with increasing time. The initial value of A(t) before biomass attachment is computed as:


= π(ID)2. h/4, the empty-column volume (L3)

Vs = volume occupied by impermeable solids (L3)

ID = internal diameter of the reactor column (L)

h = height of the packed-bed (L)

The available cross-sectional area A(t) affects the interstitial fluid velocity u (LT1) and the stagnant liquid layer around the soil particles Lw. The thickness of the stagnant liquid layer is estimated using the Frossling's correlation (Frossling, 1938):


Diw = diffusion coefficient of dissolved species in water (L2T-1)

u = Qc/A(t), the interfacial discharge rate (LT -1)

Qc = Q in+ QR, the combined flow across the bulk liquid zone

(L3T-1), where Qin = feed flow rate (L3T-1) and QR = recycle flow rate (L3T-1)'"

dp = average particle size (L)

ν = kinematic viscosity (LT-1)

A(t) = effective cross-sectional area of the reactor (L2)

Reactor mass balance and kinetics

Removal of dissolved species over time across the bulk liquid of the soil column is represented by a set of ordinary differential equations consisting of the following aspects in its mass balance equation:



C = effluent Cr(VI) concentration (ML-3)

V = volume of the reactor (L3)

Cin = influent Cr(VI) concentration (ML-3)

Q (Au) = influent flow rate (L3T-1)

A = cross sectional area of reactor (L2)

u = velocity of the flow (LT-1)

t = time (T)

Reduction by microbial kinetics (reaction process)

The associated rate kinetics, which were based on non-competitively inhibited rate kinetics, were proposed by Molokwane (2010). In Molokwane's work, a finite Cr(VI) reduction capacity was assumed in immobilised cells. The acclimated cells demonstrated the capability to function normally until a toxicity threshold concentration of approximately 100 mg.-1 Cr(VI).


K = dimensionless Cr(VI) inhibition factor (MM-1)

Kc = half-velocity coefficient (ML-3)

Xa = concentration of Cr(VI) reducing cells (ML -3)

Cr= inhibition threshold concentration (ML -3)

Cin = influent Cr(VI) concentration (ML-3)

Rc = Cr(VI) reduction capacity of cells [g Cr(VI) reduced/g immobilised cells] (MM-1)

Y = cell yield coefficient (MM-1)

kms = maximum specific substrate utilisation rate coefficient (T--1)

kmc = maximum specific Cr(VI) reduction rate (T-1)

kd = cell death rate coefficient (T-1)

X = viable cell concentration (ML-3)

Due to space limitations in the reactor, cells may only grow to a certain maximum concentration [Xmax ]. The time at which the cells reach the maximum allowable concentration is dependent on initial cells, Cr(VI) toxicity and hydraulic loading rate. These conditions allow us to simplify the cell growth term (Eq. 7) to the following logistic function:


X = viable cell concentration (ML-3) at any time t (T)

Xmax = maximum attainable viable cell concentration (ML-3) in the microcosm

Xo = initial viable cell concentration (ML 3), t0 = logistic interval (T)

b = pitch (dimensionless).

Diffusion process

Mass transfer within the attached cell layer is described by Fick's law for dispersion as presented in (Eq.(2)). The contaminant flux across the stagnant layer to the biofilm is a function of the contaminant dispersion coefficient and concentration and is represented as:


Diw = dispersion coefficient of Cr(VI) in water (L2T-l)

dC/dz = Cr(VI) concentration gradient (ML-3L-1)

LW = thickness of stagnant layer (L)

CB= bulk liquid Cr(VI) concentration (ML-3)

jC = flux of dissolved species into the biofilm (ML-2T-1)

CS = Cr(VI) surface concentration (ML-3)

NB: In most mass transfer-limited reactions CB >> Cs, therefore Cs is negligible.

Adsorption process


kad = adsorption rate coefficient (T--1)

Ceq = equilibrium concentration at surface area (ML-3)

C = Cr(VI) concentration at any time (ML-3)

qc = rate of Cr(VI) removal by adsorption (T-1)

The adsorption and diffusion processes are facilitated by the terms (q) and (j), respectively, and the reaction process is represented by the term (r) which represents reaction by biomass (X). The term (q) can approach equilibrium easily whereas the terms (r) and (j) depend on the active biomass. The overall mass balance for dissolved species (Eqs (5) to (10)) and biomass is given by the equation below:


jc = flux of dissolved species into the biofilm (ML-2T-1)

rc = rate of removal of dissolved species by suspended cells in the bulk liquid (ML-3T-1)

λ = cell detachment rate coefficient (T-1) - a function of the interfacial velocity u (LT-1)

kd = cell death rate coefficient (T-1)

Af = available surface area (L2)

V = the bulk liquid volume (L3)

XB = biomass density in the bulk liquid (ML-3)

Xf = biomass density in the biofilm zone (ML-3)

Lf = thickness of biofilm layer (L).

Quasi-steady state kinetics

Steady-state conditions were evaluated with the soil column as a plug flow system with incremental loading per mass of reactor:


Fc = molar loading rate of Cr(VI) (MT-1)

W = mass of aquifer soil particles (M)

yc = reaction rate (ML-3T-1)

Under quasi-steady-state conditions biomass, adsorption and diffusion process in Eqs( 9) and (10) reached equilibrium. Therefore, the local reaction rate at a particular point can be simplified in terms of the Cr(VI) reduction capacity of individual cells. The reaction was originally first order but was slowed down by the toxicity of the reactant; a mixed kinetic may be represented as follows:

(Levenspiel, 1999)


C = Cr(VI) concentration along the reactor column (ML-3)

Cr = Cr(VI) toxicity threshold coefficient (ML-3)

k = reaction rate coefficient (L3M-lT-l)

n = semi-empirical dimensionless variable (MM-1)

The mass (W) and loading rate (Fc) are defined as ρAƒL and QinC, respectively, so that the differential term can be expressed in terms of dC/dL and the whole reaction term can be simplified by rearranging terms:


C = Cr(VI) concentration (ML-3) at any point (L) along the reactor column

k = reaction rate coefficient (L3M-lTl)

Qin = hydraulic loading rate (L3T1)

ρ = density of aquifer soil particles (ML-3)

Af = available surface area (L2)

L = length of the reactor (L)

The parameters ρ, Qin, and Cr are constant whereas the surface area Af varies with time as described in Eq. (3).

Assumptions governing the model in Eq.(15) are as follows:

There is no mixing in the axial direction, this implies that molecular and/or turbulent mass dispersion is negligible in flow direction

Uniform properties in the direction perpendicular to the flow (flow is 1-dimensional)

The net growth of bacteria is zero near steady-state, i.e., Cr(VI) reduction is by resting cells

The system is a homogenous catalytic system



Cr(VI) reduction kinetics in columns

Spatial and temporal variation of Cr(VI) in the columns inoculated with mixed-culture of bacteria from the dried sludge samples collected within the vicinity of the contaminated site was evaluated under both carbon source and non-carbon source conditions. Complete Cr(VI) reduction was observed in the inoculated carbon source reactor after 102 h of operation while a Cr(VI) removal efficiency of 69.5% was observed in the inoculated non-carbon source reactor at the same time of operation (Fig. 2A-B). These results suggest that the presence of the carbon source greatly enhances or stimulates the performance of Cr(VI) reduction in a biological system.



Additionally, data collected from equally-spaced longitudinal sampling ports of columns inoculated with a mixed-culture of bacteria demonstrated significant Cr(VI) reduction with maximum observable rate of approximately 2 mg--1.h-1 between 50 and 75 h, with a removal efficiency of about 70% (Fig. 2A-B). These results demonstrate that the effluent Cr(VI) concentration in the columns inoculated with mixed-culture of bacteria did not stabilise until a quasi-steady state was achieved. The insignificant Cr(VI) reduction (< 5%) along the column at quasi-steady state, i.e. between 75 and 102 h may be associated with Cr(VI) inhibitory effects on Cr(VI)-reducing bacteria, loss of Cr(VI) reduction capacity due to blockage of pores with Cr(III) precipitate, and/or the saturation of physical chemical process-adsorption and biosorption over time.

The performance of Cr(VI) reduction between the native soil culture column amended with sawdust and the sterile-control column was also evaluated. Insignificant Cr(VI) reduction was observed in both columns with a removal efficiency of about 10% for the native soil culture column and only 2% for the sterile-control column after 102 h of operation (Fig. 3A-B). Pronounced Cr(VI) reduction observed in columns inoculated with a mixed-culture of bacteria, as compared to the native soil culture column, indicated the prominence of in situ bioremediation using target Cr(VI)-reducing species. The insignificant Cr(VI) reduction observed in the native soil culture column amended with sawdust as carbon source indicates that the native species present in the soil samples may not be Cr(VI) reducers but rather Cr(VI) resistors, as they remained persistent in the column even after long periods of exposure to Cr(VI) loadings. Additionally, insignificant Cr(VI) reduction in a sterile-control (heat-killed cells) reactor over time indicated that the abiotic processes are negligible.



Evaluation of internal processes

Parameter optimisation and model evaluation

The effluent Cr(VI) concentration across the reactor at transient-state was evaluated using the model defined in Eq. (11). The kinetic parameters were obtained by performing a nonlinear regression analysis using AQUASIM 2.0 (a computer program for identification and simulation of aquatic systems). The reliability of the model was demonstrated by near-constant optimisable parameters under different loading conditions, as shown in Table 2A. Conversely, spatial modelling of Cr(VI) concentration profile across the reactor was evaluated using the model developed under near steady-state conditions (Eq. 15).



The model was initially tested with data from the experimental run of 50 mg-£-1 Cr(VI) in the inoculated columns. The experimental run at this concentration was treated as a control due to insignificant Cr(VI) reduction (<10%) observed under these loading conditions after operation. The kinetic parameters for the model developed under near steady-state conditions were obtained using a computer program for solving numerical problems (Octave 3.0). Trial values of the optimisable kinetic parameters were initially adapted from the transient-state kinetic parameters and then adjusted as a result of low levels of biomass under near steady-state conditions. The optimum kinetic parameters (k and n) obtained from the initially tested Cr(VI) concentration under quasi-steady-state conditions were then used to simulate the experimental data under different loading conditions in both carbon source and non-carbon source columns. The model successfully captured the trends of Cr(VI) response profiles under quasi-steady-state conditions for different loading conditions. The reliability of the model was demonstrated by near-constant optimisable parameters under different loading conditions, as shown in Table 2B.



Performance evaluation under near steady-state conditions

The quasi-steady-state model was evaluated in the column inoculated with mixed-culture of bacteria under both carbon and non-carbon source conditions (Fig. 4A-C). Cr(VI) toxicity threshold; Cr in the model was assumed to be 50 mg-£-1 based on the experimental data obtained at this concentration. Fig. 4A demonstrates that the inhibitory effects of Cr(VI) removal in a non-carbon source column occurs at a faster rate than in the carbon source reactor. This indicates that the presence of organic compounds in the carbon source reactor enhances Cr(VI) reduction activity by stimulating microbial cell activity. It is also shown in Table 2A and Table 2B that minor adjustments of the kinetic parameter, k, were achieved under non-carbon source conditions as a result of slower reaction rate in the absence of carbon source. The toxicity term in the model, suggests that:



The above term implies that the greater the value of (n) in the biological system, the slower the rate at which Cr(VI) toxicity in the cells is approached. As a result of slower Cr(VI) toxic-ity rates achieved in cells at higher (n) values, the subsequent inhibitory effects of Cr(VI) reduction by Cr(VI)-reducing culture in the system are also approached at a slower rate.

Evaluation of organics in the sawdust

In order to increase the rate of Cr(VI) reduction in the inoculated column, sawdust which simulated the carbon source from decaying vegetation within the contaminated site was introduced in the column. The results from TOC analyser showed that at the initial time of incubation (tg), the total amount of organic carbon in the leachate sample was 268.7 mg-£-1and then decreased to 198.3 mg-£-1 after 12 days of incubation, suggesting that sawdust contain certain organic compounds which may be biodegradable. Additionally, GC/MS also confirmed that certain organic compounds in the sawdust leachate, such as 1-ethyl-3-vinyl, mercaptoacetic acid, Frenchol, and Bicyclo 2.2.1 heptane of molecular weight 161 g-mol-1, 73 g-mol-1, 81 g-mol-1, and 95 g-mol-1 may be partially degradable, as they were observed to be slowly disappearing in a sawdust leachate over time (Fig. 5A-B).

Microbial culture dynamics

Characterisation of initially inoculated microcosm bacteria

The mixed-culture of bacteria used for column inoculation was isolated under micro-aerobic conditions from the dried sludge. The predominant species under these conditions were the Gram-positive Bacilli. A partial sequence of 16S rRNA matched the Bacillus groups - Bacillus cereus ATCC 10987, Bacillus cereus 213 16S, Bacillus thuringiensis, Bacillus mycoides, and 2 Microbacterium groups - Microbacterium foliorum and Microbacterium sp. S15-M4.

Characterisation of microcosm bacteria after of operation

After 7 weeks of microcosm operation under oxygen-stressed conditions, a microbial community shift was expected in the non-sterile soil reactors inoculated with mixed-culture of bacteria from the dried sludge (R2 and R4), as the soil samples themselves contain a wide range of soil-dwelling species of bacteria (native soil culture). The direct measurement of viable cell concentration was challenging in the reactor since the majority of microbial species were not culturable using the conventional methods. The microbial dynamics monitored by the 16S rRNA fingerprinting method showed a decrease in culturable species. After exposure to Cr(VI) for 7 weeks the only serotypes which persisted, either through resilience against toxicity or through adaptation to changing conditions in the soil columns inoculated with mixed-culture of bacteria, were Enterobacter sp., B. cereus, and B. thuringiensis (Fig. 6). Both these Bacilli species are well known Cr(VI)-reducing species (Faisal and Hasnain, 2010; Francisco et al., 2002).



Microcosm experiments show that it is possible to introduce microbial cultures, both pure and mixed-cultures of bacteria, into the environment in a controlled manner to successfully achieve Cr(VI) reduction in flowing water. Compared to the indigenous soil culture which insignificantly removed Cr(VI) in the system, even in the presence of the added carbon source, the mixed-culture of bacteria from the dried sludge samples collected within the vicinity of a contaminated site was able to sustain high Cr(VI) removal at 40 mg.-1, which is the highest groundwater Cr(VI) concentration at the study site. Although the microbial culture showed good Cr(VI) reduction capacity, the application of a biological barrier system for in situ biore-mediation of Cr(VI) remained a challenge due to foreseeable blockage by Cr(OH)4 species in the barrier zone. The model developed in this study under quasi-steady-state conditions, which is incorporated with flow characteristics and inactivation term (C), successfully predicted the effluent Cr(VI) concentration profiles in the column inoculated with a mixed-culture of bacteria for different loading conditions, with a R2 value of 95%. Kinetic parameters obtained in the finite difference model showed highest removal rate in the column inoculated with a mixed-culture of bacteria in the presence of non-toxic carbon source, sawdust, than in the inoculated column without carbon source. These results indicate that the rate of Cr(VI) reduction in the column inoculated with target Cr(VI)-reducing species from the dried sludge samples may be limited by the substrate concentration in the system, which stimulates or encourages microbial activity. Modification of the model for application in engineering biological systems will facilitate final scale-up and operation of microbial barriers in the field.



The research was funded partially through the National Research Foundation (NRF) Incentive Funding for Rated Researchers, Grant No. IFR2010042900080 awarded to Evans MN Chirwa of the University of Pretoria.



ALLEN H, GARRISON A and LUTHER IG (1998) Metals in Soil Surface. Sleeping Beer Press, New York, NY.17.         [ Links ]

APHA (2005) Standard Methods for the Examination of Water and Wastewater (21st edn.). In: Eaton AD, Clesceri LS, Rice EW, Greenberg AE and Franson MAH (eds.). American Public Health Association, American Water Works Association, Water Environment Federation, USA.         [ Links ]

BARAL A and ENGELKEN RD (2002) Chromium-based regulations and greening in the metal finishing industries in the USA. Environ. Sci. Polic. 5 (2) 121-133.         [ Links ]

BONI MR and SBAFFONI S (2009) The potential of compost-based biobarriers for Cr(VI) removal from contaminated groundwater: Column test. J. Hazardous Mater. 166 (2-3) 1087-1095.         [ Links ]

BLOWES D (2002) Environmental chemistry: tracking hexavalent chromium in groundwater. Science 295 (5562) 2024-2025.         [ Links ]

CAGLIERI A, GOLDONI M, ACAMPA O, ANDREOLI R, VETTORI MV, CARRADI M, APOSTOLI P and MUTTI A (2006) The effect of inhaled chromium on different exhaled breath condensate bio-markers among chrome-plating workers. Environ. Health Perspect. 114 (4) 542-546.         [ Links ]

CHEUNG KH and GU JD (2007) Mechanism of hexavalent chromium detoxification by microorganisms and bioremediation application potential: a review. Int. Biodeter. Biodegr. 59 (1) 8-15.         [ Links ]

CHIRWA EMN and MOLOKWANE PE (2011) Biological Cr(VI) reduction: Microbial diversity, kinetics and biotechnological solutions to pollution. In: Sofo A (ed.) Biodiversity. InTech Online Publishers, United Kingdom.         [ Links ]

DWAF (DEPARTMENT OF WATER AFFAIRS AND FORESTRY, SOUTH AFRICA) (2005) Pollution of chrome-6 in the Crocodile River between Brits and the Roodekopjies Dam. Press Release by the Department of Water Affairs and Forestry, 23 December 2005. DWAF, Pretoria.         [ Links ]

EGWUONWU GN, OLABODE VO, BUKAR PH, OKOLO VN and ODUNZE AC (2011) Characterization of topsoil and groundwater at leather industrial area, Challawa, Kano, Northern Nigeria. The Pacific J. Sci. Technol. 12 (1) 628-641.         [ Links ]

FAISAL M and HASNAIN S (2010) Detoxification of Cr(VI) by Bacillus cereus S-6. Res. J. Microbiol. 7 (5) 651-656.         [ Links ]

FEDERAL REGISTER (2004) Occupational Safety and Health Administration. Occupational Exposure to Hexavalent Chromium. 69 Federal Register 59404, 4 October 2004.         [ Links ]

FRANCISCO R, ALPOIM MC and MORAIS PV (2002) Diversity of chromium-resistant and chromium-reducing bacteria in a chromium-contaminated activated sludge. J. Appl. Microbiol. 92 (5) 837-843.         [ Links ]

FROSSLING N (1938) Of [Over] the volatilization of falling drops. Gerlands Beitreage Geophys. 52 170-216 (in German).         [ Links ]

GARRELS RM and CHRIST CL (1965) Solutions, Minerals and Equilibria. Harper and Row Publishers, New York. 403-435.         [ Links ]

GONZALEZ CF, ACKERLEY DF, PARK CH and MATIN A (2003) A soluble flavoprotein contributes to chromate reduction and tolerance by Pseudomonas putida. Acta Biotechnol. 23 (2-3) 233-239.         [ Links ]

GOMM JBD and WILLIAMS JTE (1996) Development of a Neural-Network Predictive Controller. Liverpool John Moores University, Liverpool.         [ Links ]

JAMES BR and BARTLETT RJ (1984) Plant-soil interaction of chromium. J. Environ. Qual. 13 (1) 67-70.         [ Links ]

KOBYA M (2004) Removal of Cr(VI) from the aqueous solutions by adsorption onto hazelnut shell activated carbon: Kinetic and equilibrium studies. Bioresour. Technol. 91 (3) 317-321.         [ Links ]

KURNIAWA TA, CHAN GYS, LO WH and BABEL S (2006) Physico-chemical treatment technique for wastewater laden with heavy metals. Chem. Eng. J. 118 (1-2) 83-98.         [ Links ]

LEVENSPIEL O (1999) Chemical Reaction Engineering (3rd edn.). John Wiley & Sons, New York. 668 pp.         [ Links ]

LONGSWORTH LG (1972) Diffusion in liquids. In: Gray DW (ed.) American Institute of Physics Handbook. McGraw-Hill Book Company, New York.         [ Links ]

MOLOKWANE PE (2010) Simulation of In Situ Bioremediation of Cr(VI) in Groundwater. PhD thesis, University of Pretoria, Pretoria, South Africa. URL:         [ Links ]

MOLOKWANE PE and CHIRWA EMN (2009) Microbial culture dynamics and chromium (VI) removal in packed-column microcosm reactors. Water Sci. Technol. 60 (2) 381-388.         [ Links ]

MOLOKWANE PE, NKALAMBAYAUSI-CHIRWA EM and MELI KC (2008) Chromium (VI) reduction in activated sludge bacteria exposed to high chromium loading: Brits culture (South Africa). Water Res. 42 (17) 4538-4548.         [ Links ]

REYNOLDS TD and RICHARDS PA (1996) Unit Operations and Processes in Environmental Engineering (2nd edn.). PWS Publishing Company, New York.         [ Links ]

ROSLEV P, MADSEN PL, THYME JB and HENRIKSEN K (1998) Degradation of phthalate and di-(2-ethylhexyl) phthalate by indigenous and inoculated microorganisms in sludge-amended soil. Appl. Environ. Microbiol. 64 (12) 4711-4719.         [ Links ]

SHI WY, SHAO HB, LI H, SHAO M and DU S (2009) Progress in the remediation of hazardous heavy metal-polluted soils by natural zeolite. J. Hazardous Mater. 170 (1) 1-6.         [ Links ]

TEKERLEKOPOULOU AG, TSIAMIS G, DERMOU E, SIOZIOUS S, BOURTZIS K and VAYENAS DV (2010) The effect of carbon source on microbial community structure and Cr(VI) reduction rate. Biotechnol. Bioeng. 107 (3) 478-487.         [ Links ]

TSENG JK and BIELEFELDT AR (2002) Low temperature chromium(VI) biotransformation in soil with varying electron acceptors. J. Environ. Qual. 31 (6) 1831-1841.         [ Links ]

TOWILL LE, SHRINER CR and DRURY JS (1978) Reviews of the Environmental Effects of Pollutants: Chromium. National Academy Press, Cincinnati.         [ Links ]

TURICK CE, GRAVES C and APEL WA (1998) Bioremediation potential of Cr(VI)-contaminated soil using indigenous microorganisms. Bioremed. J. 2 (1) 1-6.         [ Links ]

WATTS RJ (1998) Hazardous Wastes: Sources, Pathways, Receptors. John Wiley & Sons, Inc., New York. 123 pp.         [ Links ]

WITTBRODT PR and PALMER CD (1992) Limitations to pump-and-treat remediation of a chromium contaminated site. Aquifer restoration: pump-and-treat and the alternatives. National Groundwater Association Convention, 30 September - 2 October 1992, Las Vegas.         [ Links ]

ZABIRI H and MAZUKI N (2010) A black-box approach in modeling valve stiction. Int. J. Eng. Nat. Sci. 4 (1) 9-16.         [ Links ]

ZAKARIA ZA, ZAKARIA Z, SURIF S and AHMAD WA (2007) Biological detoxification of Cr(VI) using wood-husk immobilized Acinetobacter haemolyticus. J. Hazardous Mater. 148 (1-2) 164-171.         [ Links ]



Phalazane J Mtimunye
+27 (12) 420-5894; fax: +27 (12) 362-5089

Received 8 May 2013
Accepted in revised form 3 March 2014