Chromium transport, oxidation, Nieznany

background image

Ž

.

Journal of Contaminant Hydrology 49 2001 311–334

www.elsevier.comrlocaterjconhyd

Chromium transport, oxidation, and adsorption in

manganese-coated sand

Hillol Guha

a

, James E. Saiers

b,)

, Scott Brooks

c

, Phil Jardine

c

,

Krishnaswamy Jayachandran

d

a

Hemispheric Center for EnÕironmental Technology, Florida International UniÕersity, 10555 West Flagler

Street CEAS 2100, Miami, FL 33174, USA

b

School of EnÕironmental Studies, Yale UniÕersity, 370 Prospect Street, New HaÕen, CT 06511, USA

c

EnÕironmental Sciences DiÕision, Oak Ridge National Laboratory, Oak Ridge, TN 06831, USA

d

Department of EnÕironmental Studies, Florida International UniÕersity, Miami, FL 33199, USA

Received 26 April 2000; received in revised form 25 October 2000; accepted 15 December 2000

Abstract

We examine how the processes of advection, dispersion, oxidation–reduction, and adsorption

Ž

.

combine to affect the transport of chromium through columns packed with pyrolusite b-MnO -

2

Ž

.

Ž

.

coated sand. We find that b-MnO effectively oxidizes Cr III to Cr VI and that the extent of

2

Ž

.

oxidation is sensitive to changes in pH, pore water velocity, and influent concentrations of Cr III .

Ž

.

Cr III oxidation rates, although initially high, decline well before the supply of b-MnO

is

2

Ž

.

Ž

.

depleted, suggesting that a reaction product inhibits the conversion of Cr III to Cr VI . Rate-limited

Ž

.

reactions govern the weak adsorption of each chromium species, with Cr III adsorption varying

Ž

.

directly with pH and Cr VI adsorption varying inversely with pH. The breakthrough data on

Ž .

chromium transport can be matched closely by calculations of a simple model that accounts for 1

Ž

.

Ž

.

Ž .

advective-dispersive transport of Cr III , Cr VI , and dissolved oxygen, 2 first-order kinetics

Ž .

adsorption of the reduced and oxidized chromium species, and 3 nonlinear rate-limited oxidation

Ž

.

Ž

.

of Cr III to Cr VI . Our work supplements the limited database on the transport of redox-sensitive
metals in porous media and provides a means for quantifying the coupled processes that contribute
to this transport. q 2001 Elsevier Science B.V. All rights reserved.

Keywords: Solute transport; Geochemical models; Oxidation; Kinetics; Chromium

)

Corresponding author. Tel.: q1-203-432-5121; fax: q1-203-432-3929.

Ž

.

E-mail address: james.saiers@yale.edu J.E. Saiers .

0169-7722r01r$ - see front matter q 2001 Elsevier Science B.V. All rights reserved.

Ž

.

PII: S 0 1 6 9 - 7 7 2 2 0 0 0 0 1 9 9 - 6

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

312

1. Introduction

Chromium, released from liquid and solid waste sources, has polluted groundwaters

Ž

.

throughout North America and Europe Fetter, 1993; Westbrook, 1983; Hartford, 1983 .
A combination of geochemical processes, including oxidation–reduction, adsorption,
and precipitation, govern the subsurface transport and distribution of this toxic metal.
Manganese-oxides, which occur commonly as coatings on mineral grains and are the

Ž

.

only known naturally occurring oxidants of chromium, catalyze the oxidation of Cr III

Ž

. Ž

to Cr VI

Barlett and James, 1979; Johnson and Xyla, 1991; Manceau and Charlet,

.

1992 . In most soil–water systems, the oxidized chromium species is relatively soluble
and does not sorb strongly to geologic materials; as a result, it is capable of traveling in

Ž

.

a nearly conservative fashion Davis et al., 1993 . The reduced chromium species, on the
other hand, is characterized by low aqueous solubility, and precipitation and adsorption

Ž

reduce its pore water concentrations and retard its overall rate of migration Leckie et

.

al., 1984; Sass and Rai, 1987 .

Recognition of the critical link between the oxidation state and the mobility of

chromium has motivated research into the environmental controls on the oxidative

Ž

.

transformation of Cr III by manganese oxides. Results of batch experiments demon-

Ž

.

strate that solution pH and ionic composition influence the kinetics and extent of Cr III

Ž

.

oxidation by manganese-containing minerals, such as pyrolusite b-MnO , hausmannite

2

Ž

.

Ž

. Ž

Mn O , and birnessite d-MnO

Eary and Rai, 1987; Fendorf et al., 1993; Chung et

2

4

2

.

Ž

.

al., 1994 . The oxidation reaction is nonlinear: Cr III oxidation rates, although initially
rapid, decline precipitously well before the concentrations of reactants are depleted
Ž

.

Ž

.

Fendorf and Zasoski, 1992 . This inhibition in Cr III oxidation has been attributed to

poisoning of the manganese surface by products of the oxidation–reduction reaction
Ž

.

Amacher and Baker, 1982; Eary and Rai, 1987 , as well as to occlusion of the oxidizing

Ž

.

Ž

.

sites by precipitation of Cr OH , a redox-stable solid phase Fendorf et al., 1992 .

3

The kinetics of the chromium redox reactions, which ultimately control the distribu-

Ž

.

Ž

.

tion of Cr III and Cr VI , complicate chromium transport because each species exhibits
a unique sorptive response to changes in physico-chemical conditions. The oxidized
chromium species exists in anionic form in natural waters, and its adsorption to oxide
surfaces is readily reversible and varies inversely with pH and concentrations of

Ž

.

Ž

.

competing anions Mesuere and Fish, 1992; Anderson et al., 1994 . In contrast, Cr III is
cationic, its adsorption tends to increase with pH, and, owing to coordination in

Ž

inner-sphere complexes, its desorption from oxide surfaces is slow Charlet and Manceau,

.

1992; Fendorf et al., 1994; Fendorf and Sparks, 1994 .

Although chromium oxidation and adsorption has been investigated rather extensively

in batch experiments, little effort has been devoted to examining how these geochemical
reactions combine with hydrologic processes to influence chromium migration. In the
work reported here, we examine the transport of chromium through columns packed
with b-MnO -coated quartz sand. We develop a simple mathematical model for

2

chromium transport and reaction, and we compare calculations of this model to data
collected from the column experiments. The experimental results, together with the
mathematical modeling, illuminate how changes in physical and chemical conditions

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

313

affect the kinetics of chromium oxidation and adsorption and define the role of these
mass transfer processes in the fate of chromium within MnO -containing porous media.

2

2. Laboratory methods

We measured the transport of chromium under water-saturated conditions in labora-

tory columns packed homogeneously with b-MnO -coated sand. We conducted six

2

Ž

.

column experiments with influent solutions containing Cr III and six column experi-

Ž

.

ments with influent solutions containing Cr VI . In these experiments, we tested the
sensitivity of chromium breakthrough to changes in pH, flow rate, and influent chromium

Ž

.

concentration Table 1 . All experiments were conducted in duplicate.

Ž

.

A mineralogically pure quartz sand Unimin with a diameter between 150 and 250

mm was coated with b-MnO and used in the column experiments. We coated the sand

2

Ž

.

according to the procedure described by Jardine and Taylor 1995 . This procedure

Ž

.

involved adding diluted Mn NO

to acid-washed quartz sand, drying the sand in an

3 2

oven at 430 K for 3 days, and washing the coated sand with 5 mM CaCl until the loss

2

of b-MnO

to the solution was negligible. X-ray diffraction analysis of precipitated

2

material prepared in this fashion, but without the sand, confirmed that the Mn-oxide
phase was nearly 100% b-MnO . The concentration of Mn on the sand, determined by

2

digesting the coating with 40% nitric acid and measuring the extracted manganese

Ž

.

concentrations by inductively coupled plasma emission spectroscopy ICPES , was
0.045 mmol Mnrg sand.

Ž

.

Glass chromatography columns Kontes with an internal diameter of 1 cm were

packed with b-MnO -coated sand and used in the chromium displacement experiments.

2

These columns were oriented vertically and sealed at the top and bottom with Teflon
end-fittings. The columns were packed to a height of 4.2 cm by slowly pouring 4.5 g of
b-MnO -coated sand into 4 ml of 5 mM CaCl

standing in the columns. The pore

2

2

volumes of the sand packs measured 1.6 cm

3

, corresponding to a porosity of 0.48.

Table 1

Ž .

Ž

.

Žw

Ž

.x .

Ž

.

Measurements of average pore water velocity Õ , influent Cr III concentration

Cr III

, influent Cr VI

0

Žw

Ž

.x .

Ž

.

concentration

Cr VI

, and duration of the input pulse TP

0

w

Ž

.x

w

Ž

.x

Ž

Experiment

pH

Õ

Cr III

Cr VI

TP pore

0

0

Ž

.

Ž

.

Ž

.

.

Number

cmrh

mM

mM

volumes

1

3

10.6

0.2

0

17.6

2

3

10.6

0.4

0

17.8

3

3

101.2

0.4

0

168.5

4

4

10.4

0.2

0

25.2

5

4

10.4

0.4

0

17.5

6

4

100.8

0.4

0

120.6

7

3

10.6

0

0.2

17.6

8

3

10.6

0

0.4

17.8

9

3

101.2

0

0.4

24.1

10

4

10.4

0

0.2

12.6

11

4

10.4

0

0.4

17.5

12

4

100.8

0

0.4

120.6

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

314

In order to equilibrate the columns, a 5 mM CaCl solution, adjusted to pH 3 or to

2

pH 4 by addition of HCl, was injected into the columns until the pH of the effluent
matched the pH of the influent. The CaCl solution was maintained in equilibrium with

2

the atmosphere and had a dissolved O

concentration of 0.26 mM. A variable-flow

2

Ž

.

syringe pump

Model M365, Orion Research , positioned at the inlet end of the

columns, controlled the upward flow of the CaCl solution through the columns. Flow

2

was temporarily stopped once the columns had equilibrated, whereupon a solution

Ž

.

Ž

.

containing 5 mM CaCl and either Cr III or Cr VI was introduced to the bottom of the

2

Ž

Ž

.

Ž

.

.

columns at a constant rate see Table 1 for influent Cr III and Cr VI concentrations .

Ž

.

Ž

.

The Cr III and Cr VI solutions were prepared from CrCl

salts and K CrO salts,

3

2

4

respectively, had a dissolved O concentration of 0.26 mM, and were adjusted to pH 3

2

or pH 4 with HCl. Injection of the chromium solution was terminated once chromium

Ž

concentrations in the effluent reached steady-state levels see Table 1 for influent pulse

.

duration , and a water of identical pH and CaCl

concentration, but containing no

2

chromium, was applied to the bottom of the columns until effluent concentrations
returned to baseline levels.

Effluent samples were collected from the top of the columns in Pyrex test tubes by

Ž

.

use of an automated fraction collector Eldex Laboratories, CA . Aqueous concentra-

Ž

.

Ž

tions of Cr VI were determined by a modified s-diphenyl carbazide procedure Barlett

.

and James, 1979 , aqueous concentrations of total Cr and Mn were determined by

Ž

.

ICPES, and aqueous concentrations of Cr III

were calculated from the difference

Ž

.

between concentrations of total chromium and Cr VI . Dissolved oxygen levels in the

Ž

column effluent were monitored inline with a microelectrode oxygen probe Microelec-

.

trodes, Bedford, NH . Analysis of column experiments conducted without the sand

Ž

.

Ž

.

demonstrated that dissolved oxygen did not oxidize Cr III to Cr VI and that the

Ž

.

Ž

.

column apparatus did not adsorb detectable quantities of Cr III or Cr VI from solution.

3. Theory

3.1. Model description

We quantify the factors that control the fate of chromium by comparing results of the

laboratory experiments to calculations of a numerical model that accounts for transport,
oxidation-reduction reactions, and adsorption. This model shares a structure similar to a

Ž

.

transport model developed by Saiers et al. 2000 , but has been modified to reflect the
stoichiometry unique to the chromium oxidation–reduction reactions. The model solves

Ž

.

Ž

.

advection–dispersion equations for the movement of Cr III , Cr VI , and dissolved
oxygen. We assume that first-order rate laws govern the adsorption of each chromium
species and that dissolved oxygen transport is unaffected by sorptive mass transfer.

We couple the equations for transport and adsorption to nonlinear equations for

Ž

.

rate-limited oxidation–reduction reactions. We assume that oxidation of Cr III gener-

Ž

.

Ž

.

ates Cr VI and MnOOH, a Mn III surface-bound precipitate:

Cr

3q

q

3MnO s q 4H O ™ HCrO

y

q

3MnOOH s q 4H

q

1

Ž .

Ž .

Ž .

2

2

4

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

315

This overall reaction is supported by results of recent surface spectroscopic studies that

Ž

.

Ž .

show Cr III , as well as Co II EDTA and oxalate, reduce MnO to a surface-associated

2

Ž

.

Ž

.

Mn III species Banerjee and Nesbitt, 1999a,b; Fendorf et al., 1999 . An alternate
reaction can be written with Mn

2q

, instead of MnOOH, as the end product of b-MnO

2

reduction. Nevertheless, we believe that such a reaction does not contribute significantly

Ž

.

to Cr III oxidation in our study because, across the range of conditions tested, the molar

2q

Ž

.

mass of Mn

eluted from the columns did not exceed 5% of the molar mass of Cr III

Ž

2q

oxidized within the columns.

Mn

does not adsorb to b-MnO

in appreciable

2

Ž

.

2q

quantities at low pH Fendorf et al., 1992 , so the eluted mass of Mn

closely

2q

.

approximates the total mass of Mn

generated within the columns.

The low levels of Mn

2q

observed in our experiments also serve as evidence against

Ž

.

Ž

.

the relative importance of MnOOH as a secondary oxidant of Cr III . Although Mn III

Ž

.

Ž

. Ž

.

Ž

.

phases do oxidize Cr III to Cr VI

Johnson and Xyla, 1991 , Cr III oxidation by

MnOOH, which generates Mn

2q

as a reaction product, must proceed at substantially

Ž

.

Ž .

slower rates than Cr III oxidation by b-MnO , as described by reaction 1 . Because

2

Ž

.

redox reactions with MnOOH do not produce Cr VI at concentrations comparable to

Ž

.

redox reactions with b-MnO , we assume b-MnO represents the sole oxidant of Cr III

2

2

in our experimental system.

We hypothesize that b-MnO , depleted in the oxidation reduction reaction with

2

Ž

.

Cr III , can be regenerated in the presence of oxygenated pore water. The model
equations for this process are based on the overall reaction

MnOOH q 0.5O q H

q

™ MnO q H O

2

Ž .

2

2

2

Once regenerated, the solid phase b-MnO is free to participate in oxidation–reduction

2

Ž

.

Ž .

reactions with Cr III . A reaction similar in form to reaction 2 was proposed by Jardine

Ž

.

Ž .

and Taylor

1995

to account for prolonged oxidation of Co II EDTA in column

Ž

.

experiments with b-MnO -coated sand. Saito 1951 demonstrated that oxygen adsorp-

2

Ž

.

tion on Mn III oxides resulted in the transformation to a more catalytically active MnO

2

Ž .

solid. If reaction 2 plays a key role in the redox dynamics of our experiments, we

Ž

.

should observe sustained oxidation of Cr III , as well as a decline in concentrations of
dissolved oxygen.

3.2. Model equations

The one-dimensional form of the advection-dispersion equation quantifies the spatial

and temporal changes in the concentrations of reduced and oxidized species of chromium
within the laboratory columns:

2

E Cr III

r E Cr III

E

Cr III

E Cr III

Ž

.

Ž

.

Ž

.

Ž

.

S

q

s

D

y Õ

q

Q1

3

Ž .

Cr ŽIII.

2

Et

Q

Et

E x

E x

2

E Cr VI

r E Cr VI

E

Cr VI

E Cr VI

Ž

.

Ž

.

Ž

.

Ž

.

S

q

s

D

y Õ

q

Q1

4

Ž .

Cr Ž VI.

2

Et

Q

Et

E x

E x

w

Ž

.x

w

Ž

.x

where Cr III

and Cr VI

are the aqueous concentrations of the reduced and oxidized

w

Ž

. x

w

Ž

. x

species, respectively, Cr III

and Cr VI

are the sorbed concentrations of each

S

S

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

316

species, r is the bulk density, Q is the porosity, t is time, x is the coordinate parallel to
flow, D is the hydrodynamic dispersion coefficient, and Õ is the average pore water
velocity. Q1

and Q1

are sourcersink terms that describe the rates of species

CrŽIII.

CrŽVI.

concentration change resulting from the kinetics of the oxidation–reduction reaction

Ž .

given by reaction 1 . The sourcersink terms, which are defined below, are positive if
the species is generated in the reaction and negative if the species is consumed in the
reaction.

The equations for chromium transport are coupled with kinetics expressions for

sorptive solute mass transfer. Various forms of linear and nonlinear rate laws have been

Ž

adopted to describe chromium adsorption to geologic solids Selim and Amacher, 1988;

.

Jardine et al., 1999 . We quantify the temporal change in sorbed concentrations as the
difference between first-order adsorption and desorption rates:

E Cr III

k Q

Ž

.

S

f 1

s

k

Cr III

y

Cr III

Ž

.

Ž

.

S

b1

ž

/

Et

k

r

b

1

s

k

K

Cr III

y

Cr III

5

Ž

.

Ž

.

Ž .

Ž

.

S

b1

d1

E Cr VI

k Q

Ž

.

S

f 2

s

k

Cr VI

y

Cr VI

Ž

.

Ž

.

S

b 2

ž

/

Et

k

r

b 2

s

k

K

Cr VI

y

Cr VI

6

Ž

.

Ž

.

Ž .

Ž

.

S

b 2

d 2

Ž

.

Ž

.

where k

and k

are adsorption rate coefficients for Cr III and Cr VI , respectively,

f1

f2

Ž

.

Ž

.

and k

and k

are desorption rate coefficients for Cr III and Cr VI , respectively. The

b1

b2

Ž

.

Ž

.

equilibrium distribution coefficient, given by K

for Cr III and K

for Cr VI ,

d1

d2

incorporates the ratio of the adsorption rate coefficient to the desorption rate coefficient.

Ž .

Ž .

An equation similar in form to Eqs. 3 and 4 describes the movement of dissolved

oxygen through the sand columns:

w

x

2

w

x

w

x

E O

E

O

E O

2

2

2

s

D

y Õ

q

Q2

7

Ž .

O

2

2

Et

E x

E x

w

x

where O

is the aqueous-phase oxygen concentration and Q2

quantifies temporal

2

O

2

changes in oxygen concentration resulting from the kinetics of the oxidation–reduction

Ž .

reaction expressed by reaction 2 .

The solid-phase concentrations of b-MnO and MnOOH are not affected directly by

2

advective-dispersive transport processes, but vary only in response to redox reactions.
The total change in b-MnO concentrations is equal to the sum of Q1

, the change

2

MnO

2

Ž

.

in b-MnO concentrations that arises from Cr III oxidation, and Q2

, the change in

2

MnO

2

b-MnO concentrations that arises from reoxidation of MnOOH by dissolved oxygen:

2

w

x

r E MnO

2

s

Q1

q

Q2

8

Ž .

Mn O

MnO

2

2

Q

Et

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

317

w

x

where MnO

is the solid-phase concentration of b-MnO . Variations in solid-phase

2

2

MnOOH occur concomitantly and are inversely proportional to changes in b-MnO

2

concentrations. As such, the mass balance equation for MnOOH is written as

w

x

r E MnOOH

s

Q1

q

Q2

9

Ž .

Mn OOH

MnOOH

Q

Et

w

x

where MnOOH is the solid-phase concentration of MnOOH, Q1

is the rate that

Mn OOH

Ž

.

MnOOH is generated from the reduction of b-MnO by Cr III and Q2

is the rate

2

MnOOH

that MnOOH is consumed through its oxidation by O .

2

Ž . Ž .

Eqs. 3 – 9 represent the most general form of the mass balance equations for

chromium fate and transport. Solution of these equations requires specification of the
rate laws that define the source-sink terms. We derive the mathematical relationship

Ž .

between Q1

, Q1

, Q1

, and Q1

from the stoichiometry of reaction 1

CrŽIII.

CrŽVI.

MnO

MnOOH

2

and the mathematical relationship between Q2

, Q2

, and Q2

from the

Mn OOH

O

MnO

2

2

Ž .

stoichiometry of reaction 2 , such that

1

1

y

Q1

s y

Q1

s

Q1

s

Q1

10

Ž

.

Cr ŽIII.

MnO

Cr Ž VI.

MnOOH

2

3

3

and

y

Q2

s y

2Q2

s

Q2

11

Ž

.

Mn OOH

O

MnO

2

2

The rate of a chemical reaction is usually proportional to the concentrations of

Ž

.

reactants raised to some power Chang, 1984 ; therefore, the sourcersink terms of Eqs.
Ž

.

Ž

.

10 and 11 can be written in terms of the reactant species concentrations of reactions

Ž .

Ž .

1 and 2 :

1

1

a

g

w

x

y

Q1

s y

Q1

s

Q1

s

Q1

s

k

Cr III

MnO

b

Ž

.

Cr ŽIII.

MnO

Cr Ž VI.

MnOOH

R 1

2

2

3

3

12

Ž

.

and

k

l

w

x w

x

y

Q2

s y

2Q2

s

Q2

s

k

MnOOH

O

13

Ž

.

Mn OOH

O

MnO

R 2

2

2

2

where k

and k

are kinetic rate constants, a , g , k , and l are exponential constants

R1

R2

Ž

.

Ž

.

that define the reaction order, and b is a function defined below that describes Cr III
oxidation inhibition. The exponential constants can be integers, non-integers, or even

Ž

.

equal to zero Chang, 1984 . We cannot determine the values of the exponents a priori
because these parameters are generally unrelated to the stoichiometry of the overall

Ž

.

reactions see Sposito, 1994 . To reduce the number of adjustable parameters in our
model, we assume that the exponents equal one, and thus the oxidation–reduction
reactions are first order with respect to each of the reactants.

We incorporate a function for oxidation inhibition because several studies have

Ž

.

demonstrated that rates of Cr III oxidation decline, and even cease, well before the

Ž

.

supply of MnO becomes limiting Amacher and Baker, 1982; Eary and Rai, 1987 . To

2

Ž

.

our knowledge, three different explanations for Cr III oxidation inhibition have been

Ž

.

2q

Ž

.

proposed. Amacher and Baker 1982 speculated that Mn

competed with Cr III for

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

318

Ž

.

sorption sites on d-MnO and reduced the amount of Cr III available to participate in

2

Ž

.

Ž

.

the oxidation reaction. Eary and Rai 1987 , on the other hand, theorized that Cr VI

Ž

.

adsorption to b-MnO drove the sharp decline in Cr III oxidation rates by limiting the

2

Ž

.

amount of Cr III able to contact oxidizing sites on the b-MnO . In later work, Fendorf

2

Ž

.

2q

Ž

.

Ž

.

et al. 1992, 1993 proved that neither Mn

nor Cr VI impeded Cr III oxidation, but

Ž

.

Ž

.

that precipitation of Cr OH

occluded the MnO surface and arrested Cr III oxidation.

3

2

This mechanism, although important under a range of conditions, should not contribute

Ž

.

to Cr III oxidation inhibition in our experiments because chromium precipitated only

Ž

.

when the pH exceeded 4 Fendorf et al., 1992 .

We hypothesize that oxidation rates slow because MnOOH, a reaction product with

low redox reactivity compared to b-MnO , forms a physical barrier that blocks access of

2

Ž

.

Cr III

to the b-MnO

oxidizing sites. A similar inhibition mechanism has been

2

Ž

.

Ž .

proposed by Jardine and Taylor 1995 to explain rapidly declining rates of Co II EDTA
oxidation in column experiments with b-MnO -coated sand. We quantify the relation-

2

Ž

.

ship between Cr III redox kinetics and MnOOH accumulation with a Langmuir-type

Ž

.

function that depicts a linear decline in Cr III oxidation rates with increasing MnOOH
concentrations:

w

x

G MnOOH

b s 1 y

14

Ž

.

w

x

MnO

2 init

w

x

where MnO

is the initial solid-phase concentration of b-MnO . The slope of this

2 init

2

linear response is dictated by the constant, G . This dimensionless parameter is equiva-
lent to the ratio of the molar mass of b-MnO blocked from reaction to the molar mass

2

Ž

.

of MnOOH precipitated. The Langmuirian blocking function of the form of Eq. 14 was
originally developed to quantify surface exclusion of crystalline lattice sites by molecules
and also has been used to describe blocking dynamics associated with colloid deposition

Ž

.

in granular porous media Langmuir, 1918; Johnson and Elimelech, 1995 .

Ž

.

Ž

. Ž

Substitution of the relationships given by Eqs.

12

and

13

for the case of

.

Ž . Ž . Ž . Ž .

a s g s k s l s 1 into the mass balance equations given by Eqs. 3 , 4 , 7 , 8 , and
Ž .

9 yields

2

E Cr III

r E Cr III

E

Cr III

E Cr III

Ž

.

Ž

.

Ž

.

Ž

.

S

q

s

D

y Õ

2

Et

Q

Et

E x

E x

w

x

y

k

Cr III

MnO

b

15

Ž

.

Ž

.

R 1

2

2

E Cr VI

r E Cr VI

E

Cr VI

E Cr VI

Ž

.

Ž

.

Ž

.

Ž

.

S

q

s

D

y Õ

2

Et

Q

Et

E x

E x

w

x

q

k

Cr III

MnO

b

16

Ž

.

Ž

.

R 1

2

w

x

2

w

x

w

x

E O

E

O

E O

k

2

2

2

R 2

w

x w

x

s

D

y Õ

y

MnOOH O

17

Ž

.

2

2

Et

E x

2

E x

w

x

r E MnO

2

w

x

w

x w

x

s y

3k

Cr III

MnO

b q k

MnOOH O

18

Ž

.

Ž

.

R 1

2

R 2

2

Q

Et

w

x

r E MnOOH

w

x

w

x w

x

s

3k

Cr III

MnO

b y k

MnOOH O

19

Ž

.

Ž

.

R 1

2

R 2

2

Q

Et

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

319

Ž

. Ž

.

Ž .

Ž .

Ž

.

Eqs. 15 – 19 , together with Eqs. 5 , 6 , and 14 represent the complete set of
equations for the transport and fate of reduced and oxidized forms of chromium within
the b-MnO -coated sand. These equations are solved numerically using a finite-dif-

2

ference method with a predictor–corrector time-stepping scheme for a semi-infinite
column with first-type influent boundary conditions for the aqueous species. Zero initial

Ž

.

Ž

.

concentrations are assumed for Cr III , Cr VI , and MnOOH, while the initial concentra-
tions of 0.26 mM and 0.045 mmolrg are assumed for O and b-MnO , respectively.

2

2

3.3. Model sensitiÕity analysis

We conducted simulations to examine the response of model-calculated chromium

breakthrough to changes in k

, k

, and G . For this analysis, we set the adsorption and

R1

R2

desorption rate coefficients to zero in order to isolate the effects of redox reactions on

Ž

.

Ž

.

the temporal variability in Cr III and Cr VI concentrations. The influent concentrations

Ž

.

Ž

.

Fig. 1. Sensitivity of calculated concentrations of Cr III and Cr VI to variation in k

. Values used for k

:

R1

R1

Ž

.

y

1

y

1

Ž .

3

y

1

y

1

A 10 g mmol

h

and B 1=10 g mmol

h

. Both k

and G were assigned values of zero for the

R2

two simulations.

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

320

Ž

.

Ž

.

Ž

.

Fig. 2. Sensitivity of calculated concentrations of Cr III and Cr VI to variation in G . Values used for G : A

Ž .

Ž .

3

y

1

y

1

20, B 100, and C 350. The modeled results were obtained for k

s

1=10 g mmol

h

and k

s

0.

R1

R2

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

321

Ž

.

Ž

.

Fig. 3. Sensitivity of calculated concentrations of Cr III , Cr VI , and dissolved oxygen to variation in k

.

R2

Ž

.

4

y

1

y

1

Ž .

4

y

1

y

1

Ž .

5

y

1

Values used for k

: A 1=10

g mmol

h

, B 8=10

g mmol

h

, and C 5=10 g mmol

R2

h

y

1

. The modeled results were obtained for k

s

1=10

3

g mmol

y

1

h

y

1

and G s 350.

R1

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

322

Ž

.

Ž

.

of dissolved oxygen, Cr III , and Cr VI were set at 0.26, 0.201, and 0 mM, respec-
tively, and the average pore water velocity was fixed at 10.6 cmrh.

In the first set of model simulations, we varied k

by two orders of magnitude while

R1

Ž .

holding k

and G at zero. Reaction 1 exerts sole control on redox dynamics when

R2

both k

and G equal zero, and, as a result, increases in k

promote generation and

R2

R1

subsequent breakthrough of progressively higher concentrations of the oxidized species
Ž

.

Ž

.

Ž

.

Fig. 1A and B . The percentage of influent Cr III converted to Cr VI increases from

18% for the simulation in which k

s

10 g mmol

y

1

h

y

1

to 100% for the simulation in

R1

3

y

1

y

1

Ž

.

which k

s

1 = 10 g mmol

h

Fig. 1A and B .

R1

The next sequence of simulations was devoted to examining the relationship between

redox kinetics and solid-phase concentrations of MnOOH, assuming that accumulation

Ž

.

of MnOOH interferes with Cr III oxidation. We adjusted G from 20 to 350, while
holding k

at 1 = 10

3

g mmol

y

1

h

y

1

and k

at zero. For the simulation with G s 20,

R1

R2

Ž

.

Ž

Ž

.

Cr III

oxidation rates are initially high

i.e., large concentrations of Cr VI

are

.

Ž

.

generated , but begin to decrease slowly after nine pore volumes Fig. 2A . The onset of

Ž

.

the decline in the rates of Cr III oxidation occurs sooner and the rapidity of the decline

Ž

.

increases as G increases in value from 20 to 350 Fig. 2A, B, and C . These simulations
demonstrate that the sensitivity of oxidation kinetics to changes in MnOOH concentra-
tions increases with the magnitude of G ; or, in other words, for a given MnOOH

Ž

.

concentration, the level of Cr III oxidation inhibition increases with G .

We examined the sensitivity of the model solution to k

, the rate coefficient for

R2

MnOOH reoxidation, by increasing the magnitude of this parameter from 1 = 10

4

g

mmol

y

1

h

y

1

to 5 = 10

5

g mmol

y

1

h

y

1

. For each of these simulations, the values of

k

and G were set at 1 = 10

3

g mmol

y

1

h

y

1

and 350, respectively. The rate of

R1

b-MnO

restoration increases with the magnitude of k

, so the long-term oxidative

2

R2

Ž

potential of the column is greater for higher values of this parameter Fig. 3A, B, and

.

C . Consider the simulation for the lowest value of k

tested. In this case, b-MnO

R2

2

Ž

.

regeneration is slow and low rates of Cr III oxidation are maintained under steady-state

Ž

.

conditions Fig. 3A . For the highest value of k

tested, on the other hand, MnOOH is

R2

Ž

.

reoxidized to b-MnO nearly as rapidly as it forms; thus, complete conversion of Cr III

2

Ž

.

Ž

.

to Cr VI occurs until the supply of oxygen becomes limiting Fig. 3C . Dissolved

Ž

Ž ..

oxygen catalyzes the oxidation of MnOOH to b-MnO

see reaction 2 . As a result,

2

progressively greater concentrations of dissolved oxygen are consumed as the magnitude

Ž

of k

and, correspondingly, the rate of MnOOH oxidation increases Fig. 3A, B, and

R2

.

C .

4. Results

4.1. Response of chromium transport to changes in pH, flow rate, and influent
concentration

Ž

.

Results of the column experiments with Cr III -containing influent solutions demon-

Ž

.

Ž

.

strate that b-MnO effectively catalyzes the oxidation of Cr III to Cr VI . We find that

2

changes in pore water pH, average pore water velocity, and influent concentration

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

323

Ž

.

control the characteristics of Cr III oxidation. Oxidation rates are higher and thus

Ž

.

Ž

.

effluent concentrations of Cr VI are greater at pH 4 relative to pH 3 Fig. 4 . Fendorf

Ž

.

Ž

.

and Zasoski 1992 reported a qualitatively similar relationship between pH and Cr III

Ž

.

oxidation rates in batch experiments and attributed the behavior to increases in Cr III
adsorption by b-MnO with increasing pH.

2

Ž

.

An increase in pore water velocity decreases the mass of Cr III oxidized within the

Ž

.

Ž

.

columns Fig. 4 . In experiments at pH 4, for example, b-MnO oxidizes 41% of Cr III

2

Ž

.

applied to the columns at a pore water velocity of 10.4 cmrh, but only 7% of Cr III
applied to the columns at a pore water velocity of 100.8 cmrh. The dependence of

Ž

.

Cr III

oxidation on flow rate demonstrates that the redox reaction is kinetically

Ž

.

controlled; or, in other words, the time-scale for Cr III oxidation is on the order of the
time-scale for advective transport.

Chromium transport and reaction also responds to changes in the influent concentra-

Ž

.

tion of Cr III . A two-fold increase in influent concentration at a given pH produces
lower relative peak and steady-state concentrations of the oxidized chromium species
Ž

.

Fig. 4 . We infer from the dependence of the dimensionless breakthrough curves on

Ž

.

influent concentration that a nonlinear rate-law governs the kinetics of Cr III oxidation.

Ž

.

For each experimental treatment, application of the Cr III solutions produces a spike

Ž

.

Ž

.

Ž

.

in effluent Cr VI concentrations Fig. 4 . The rapid increase in Cr VI concentrations,
especially apparent between 1 and 2.5 pore volumes in the pH 4 experiments, reflects

Ž

.

Ž

.

fast rates of MnO -catalyzed Cr III oxidation. With continued Cr III injection, effluent

2

Ž

.

concentrations of Cr VI decline from their peak levels, signaling a proportionate decline

Ž

.

Ž

.

in Cr III oxidation rates. The magnitude of the decline in Cr VI concentrations varies

Ž

.

as a function of pH, flow rate, and influent Cr III concentration, but is generally greater
for the pH 4 relative to pH 3 experiments. Mass balance calculations indicate that,
depending on the experimental treatment, 94–98% of the b-MnO within the columns

2

remains unreacted when oxidation rates begin to decline. This observation is at least

Ž

.

qualitatively consistent with our hypothesis that an Mn III precipitate forms over the

Ž

.

redox-reactive surface and impedes Cr III oxidation.

Ž

.

Ž

.

Although conversion of Cr III to Cr VI is inhibited, it does not cease during the

Ž

.

course of the experiments, as effluent Cr VI concentrations eventually stabilize at levels
well above zero. We suspect that regeneration of MnO by oxidation of MnOOH is

2

Ž

.

responsible for sustaining the low levels of Cr III oxidation throughout the course of
the experiment. Oxidation of MnOOH to MnO consumes dissolved oxygen. Our results

2

show that dissolved oxygen concentrations decrease concomitantly with the break-

Ž

. Ž

.

Ž

.

through of Cr VI

Fig. 4 . Conditions that promote greater Cr III oxidation and greater

MnOOH production—low flow rates and high pore water pH—also promote a more
precipitous decline in oxygen levels.

Oxidation–reduction reactions clearly play an important role in the fate of chromium;

however, sorptive mass transfer reactions exert a relatively small influence on the

Ž

.

Ž

.

overall movement of chromium. Experiments in which Cr VI is used in place of Cr III

Ž

.

in the influent solution show that Cr VI retardation is small across the range of

Ž

.

conditions tested and that Cr VI transport approximates that of bromide, a conservative

Ž

.

Ž

tracer Fig. 5 . In a similar fashion, the breakthrough curves for total chromium i.e.,

Ž

.

Ž

..

Ž

.

Cr III q Cr VI , measured in experiments with Cr III -containing influents, nearly

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

324

Ž

.

Ž

.

Fig. 4. Measured breakthrough curves symbols and model-calculated breakthrough curves lines for the

Ž

.

displacement of Cr III influent solutions. The three experiments are distinguished on the basis of pH, influent

Ž

.

Žw

Ž

.x .

Ž .

Cr III concentration

Cr III

, and average pore water velocity Õ . The asterisks represent concentrations

0

Ž

.

Ž

.

of Cr III , the open circles represent concentrations of Cr VI , and the squares represent concentrations of
dissolved oxygen.

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

325

Ž

.

Ž

.

Fig. 5. Measured breakthrough curves asterisks and model-calculated breakthrough curves solid lines for

Ž

.

the displacement of Cr VI influent solutions. The treatments are distinguished on the basis of pH, influent

Ž

.

Žw

Ž

.x .

Ž .

Ž

.

Cr VI concentration Cr VI

, and average pore water velocity Õ . Comparison of Cr VI concentrations to

0

Ž

.

Ž

.

measured bromide concentrations open circles demonstrates that Cr VI travels in a nearly conservative
fashion.

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

326

Ž

.

match the breakthrough curves for bromide, indicating that Cr III also exhibits a weak

Ž

.

affinity for the b-MnO -coated sand Fig. 6 .

2

4.2. Comparison of model calculations and experimental results

Analysis of the breakthrough data, in coordination with the results of the sensitivity

analysis, suggests that advective–dispersive transport, rate-limited redox reactions, and,
to a lesser extent, sorptive mass transfer control the fate of chromium within the sand
columns. Here we attempt to define the response of these coupled processes to changes
in experimental conditions by comparing the breakthrough data to model calculations.

Seven parameters for oxidation–reduction and adsorption must be specified in order

Ž

.

to simulate transport in the experiments with the Cr III

influent solutions. These

parameters can be divided into three groups: k

, k

, and G quantify oxidation–reduc-

R1

R2

Ž

.

Ž

.

tion kinetics, K

and k

describe Cr III adsorption, and K

and k

describe Cr VI

d1

b1

d2

b2

adsorption.

We estimated values of k

, k

, and G , as well as values of K

and k , directly

R1

R2

d1

b1

Ž

.

from the column experiments performed with the Cr III -containing influent solutions.
We employ a Levenberg–Marquadt least-squares algorithm to find the values of these
parameters that minimize an objective function, defined as the sum of the squared
residuals between modeled and observed concentrations:

d

J s

Cr III

y

Cr III

q

Cr VI

y

Cr VI

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

Ý

exp

exp

i

i

i

i

is1

2

w

x

w

x

q

O

y

O

20

Ž

.

Ž

.

exp

2

2

i

i

w

Ž

.x

w

Ž

.x

w

x

where d equals the number of observation times, Cr III

, Cr VI

, and O

exp

exp

2 exp

i

i

i

are experimental concentrations of each species, as measured at the ith observation time,

w

Ž

.x

w

Ž

.x

w

x

and Cr III

, Cr VI

, O

are the modeled concentrations of each species, as

i

i

2 i

calculated for the ith observation time. By defining the objective function in terms of all
three aqueous species, we obtain unique solutions to the inverse problem. Estimates of
K

and k

required for these inverse simulations were determined from analysis of

d2

b2

Ž .

Ž .

separate experiments. That is, we fit solutions of Eqs. 4 and 6 to the breakthrough

Ž

.

Ž

data measured in experiments in which Cr VI was used in the influent solutions. See

.

Fig. 5 for model fits.

Comparison of calculated concentrations to those measured in the experiments

Ž

.

conducted with the Cr III influent solutions suggests that model accounts for the key

Ž

.

mechanisms influencing chromium transport Fig. 4 . The model captures most of the
temporal variability in effluent concentrations during the initial phases of the experi-

Ž

.

ments, when Cr VI concentrations increase rapidly, peak, and then decline toward
steady-state levels. Results of the sensitivity analysis illustrate that MnO -mediated

2

Ž

.

Ž

.

Cr III oxidation causes the initial increase in effluent Cr VI concentrations, while

Ž

.

Ž

.

accumulation of MnOOH inhibits Cr III oxidation and promotes the decline in Cr VI

Ž

.

concentrations from peak levels. The good description of the spike in Cr VI concentra-

Ž

.

tions suggests that the kinetics of MnO -mediated Cr III oxidation can be approximated

2

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

327

Ž

.

Fig. 6. Breakthrough curves from experiments conducted with Cr III -containing influent solutions. The

Ž

w

Ž

.x

w

Ž

.x.

asterisks denote total chromium concentrations i.e., Cr III q Cr VI

and the open circles denote bromide

concentrations.

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

328

Ž

.

Ž

.

Ž

.

by Eqs. 12 and 14 , which relate rates of Cr VI generation to concentrations of
surface-associated MnOOH.

The model also reproduces the magnitude of chromium concentrations under steady-

state conditions and, at the same time, reproduces with reasonable success the temporal

Ž

.

variability in dissolved oxygen concentrations Fig. 4 . In development of the model, we
have assumed that oxidation of surface-bound MnOOH maintains the long-term oxida-
tive potential of the porous medium and consumes dissolved oxygen. The ability of the

Ž

.

model to simulate simultaneously oxygen concentrations and steady-state Cr VI con-

Ž .

centrations suggests that stoichiometry of reaction 2 and the rate law derived from it
Ž

Ž

..

i.e., Eq. 13

are appropriate for quantifying the kinetics of MnOOH reoxidation.

Best-fit estimates of k

exhibit small variability between experiments conducted at

R1

pH 3, ranging from 25.5 g mmol

y

1

h

y

1

to 42.0 g mmol

y

1

h

y

1

. At pH 4, optimal values

of k

show greater sensitivity to changes in influent concentration and pore water

R1

Ž

.

velocity, but still vary by less than a factor of 3.5 Table 2 . During the stage of rapid

Ž

.

Ž

.

Cr III oxidation i.e., t - 4 pore volumes , the solid-phase concentration of MnO

2

remains essentially constant and the concentrations of MnOOH remain negligibly small.

w

x

Ž

.

For these conditions, the product of k

and MnO

appearing in Eq. 12 defines a

R1

2 init

first-order rate coefficient, and the reciprocal of this first-order rate coefficient quantifies

Ž

.

the time scale for Cr III oxidation, such that

1

1

TS

s

s

21

Ž

.

X

R 1

w

x

k

MnO

k

R 1

2

R 1

init

X

Ž

.

where k

is a first-order rate coefficient for Cr III oxidation. Comparison of TS

R1

R1

Ž

.

values demonstrates that Cr III oxidation rates increase with pH, or, in other words, the

Ž

.

Ž

.

time scales for Cr III oxidation to Cr VI are shorter at pH 4 than at pH 3. Calculations
of TS

for the pH 3 experiments range from 0.53 to 0.87 h and average 0.68 h, while

R1

calculations of TS

for the pH 4 experiments range from 0.10 to 0.37 h and average

R1

0.20 h.

The magnitude of k

controls the rate of MnOOH conversion to redox-reactive

R2

b-MnO and governs the long-term oxidative capacity of the column. Optimal values of

2

k

are slightly greater at pH 4 compared to pH 3, but are nearly constant for

R2

Ž

.

experiments conducted at the same pH Table 2 . These estimates of k

can be used in

R2

Ž

.

an equation similar in form to Eq. 21 to calculate the time scale for MnOOH oxidation
to b-MnO :

2

Q

1

TS

s

22

Ž

.

R 2

w

x

r k

O

R 2

2

Both the experimental and model-calculated results suggest that oxidation of MnOOH to
b-MnO consumes dissolved oxygen as it is transported along the flow path; therefore,

2

Ž

.

according to Eq. 22 , the magnitude of TS

depends on time and, for a particular time,

R2

will increase with distance along the flow path. Consider, for example, the pH 4

Ž

.

experiment conducted at the low pore water velocity and with the 0.2 mM Cr III
influent solution. At a time corresponding to elution of 10 pore volumes, values of TS

R2

range from 0.5 h at the influent end of the column, where oxygen concentrations equal

background image

()

H.
Guha

et
al.
r
Journal

of
Contaminant

Hydrology

49

2001

311

334

329

Table 2

Ž

.

Fitted parameters with standard errors

4

y

4

y

5a

w

Ž

.x

Number

pH

Õ

Cr III

k

k

=10

G

K

=10

k

K

=10

k

0

R1

R2

d1

b1

d2

b2

y

1

y

1

y

1

y

1

y

1

y

1 a

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

cmrh

mM

g h

mmol

g h

mmol

lrg

h

lrg

h

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

1

3

10.6

0.2

42.0 5.9

1.62 0.26

299.0 5.1

1.29 0.25

14.5 3.9

9.33 1.62

11.0 1.9

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

2

3

10.6

0.4

25.5 3.2

1.42 0.21

250.5 36.1

0.90 0.16

38.3 8.2

6.93 1.33

25.1 3.7

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

3

3

101.2

0.4

32.1 9.1

1.63 0.22

294.0 46.0

0.61 0.12

7.1 1.1

6.64 1.47

3.1 0.1

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

4

4

10.4

0.2

158.2 29.2

2.10 0.27

226.0 29.0

3.12 0.44

36.4 8.9

5.42 1.31

6.2 0.4

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

5

4

10.4

0.4

60.6 9.9

2.26 0.35

109.5 15.5

3.07 0.61

39.8 4.6

5.28 0.98

6.1 0.5

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

Ž

.

6

4

100.8

0.4

217.2 37.2

2.00 0.37

83.8 5.2

1.94 0.25

11.9 2.0

0.76 0.09

6.9 0.2

a

Ž .

Ž .

Ž

.

K

and k

were determined by fitting solutions of Eqs. 4 and 6 to breakthrough data measured in experiments with the Cr VI -containing influent solutions.

d2

b2

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

330

0.26 mM to1.1 h at the effluent end of the column, where oxygen concentrations equal
0.13 mM. Based on the minimum and maximum oxygen concentrations calculated for
each experiment and on the corresponding best-fit estimates of k

, calculations of TS

R2

R2

span between 0.7 and 1.1 h at pH 3 and between 0.5 and 1.7 h at pH 4.

The parameter, G , is equal to the ratio of the molar mass of b-MnO blocked from

2

reaction to the molar mass of MnOOH precipitated on the surface. Best-fit estimates of
this parameter are smaller, but show greater variation with changes in influent concentra-

Ž

.

tion and pore water velocity, at pH 4 relative to pH 3 Table 2 . In all cases, optimal
values of this parameter exceed 80, signifying that small accumulations of MnOOH are
capable of occluding a large mass of b-MnO . As discussed below, the high blocking

2

efficiency of MnOOH reflects the manner in which the b-MnO is distributed on the

2

surfaces of the sand grains.

Both the reduced and oxidized chromium species adsorb weakly to the b-MnO -coated

2

sand. Damkohler numbers, calculated as the ratio of the time-scale for advective

Ž

.

transport i.e., LrÕ, where L equals column length to the time-scale for sorptive mass

Ž

.

transfer

i.e., 1rk

and 1rk

, are less than 100, indicating that adsorption is

b1

b2

kinetically controlled. Analysis of fitted values of the partition coefficients demonstrate

Ž

.

Ž

.

that Cr III adsorption varies directly with pH, while Cr VI adsorption varies inversely

Ž

.

Ž

.

Ž

.

with pH Table 2 . Calculations of the Cr III retardation factor

R s 1 q r K rQ

1

d1

range from 1.2 to 1.4 at pH 3 and from 1.5 to 1.8 at pH 4, whereas calculations of the

Ž

.

Ž

.

Cr VI retardation factor R s 1 q r K rQ

vary between 1.2 and 1.3 at pH 3 and

2

d2

between 1.0 and 1.1 at pH 4. The low values of the R and R suggest that, even if the

1

2

residence times were great enough to achieve equilibrium, the effect of adsorption on
chromium transport would be small.

5. Discussion

With a relatively simple model, we are capable of describing how multiple oxida-

Ž

.

Ž

.

tion–reduction reactions combine to influence the breakthrough of Cr III , Cr VI , and
dissolved oxygen. The oxidation–reduction reactions are described with three adjustable
parameters—k

, k

, and G . Values of these parameters estimated from the pH 3

R1

R2

experiments are nearly constant, varying by less than a one-standard error range. At pH
4, best-fit estimates of k

are essentially constant, but optimal values of k

and G

R2

R1

depend on both flow rate and influent concentration. This variation in parameter values
with changes in physical conditions has been reported without exception in studies
where kinetics models have been tested against data from column experiments and has
been commonly attributed to an overly simplistic representation of the mass transfer

Ž

process van Genucthen and Wierenga, 1977; Rao et al., 1980; Nkedi-Kizza et al., 1983;

.

Skopp, 1986; Chen and Wagenet, 1997; Saiers and Tao, 2000 . The equations for
oxidation–reduction kinetics that govern our model, like the equations of other models
for redox-sensitive metals transport, are derived from overall chemical reactions.
Therefore, these rate laws do not explicitly account for the elementary steps of the
oxidation–reduction process, such as diffusional transport across the boundary layer,
adsorption, intermediate electron transfer steps, and desorption. Given the complexity of
the oxidation–reduction process, development and parameterization of mechanistically

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

331

precise kinetics models with complete predictive capability has proven to be exceedingly
difficult. Although our model does not capture the sequence of individual steps that

Ž

.

contribute to Cr III transformation, its application does enable us to test hypotheses
regarding the identity of the overall oxidation–reduction reactions, quantify the rates by
which these reactions proceed, and define the effects that interactions between redox
reactions have on the distribution and transport of reduced and oxidized chromium
species.

Ž

.

In general, the temporal variability in the kinetics of Cr III

oxidation can be

Ž

.

characterized in terms of three stages: a stage of relatively rapid conversion of Cr III to

Ž

.

Cr VI , a short interval of declining oxidation rates, followed by stabilization of
transformation rates at low, but nonzero levels. During the first stage, the fresh b-MnO

2

Ž

.

surface catalyzes Cr III oxidation, and the time scales for this reaction are on the order
of minutes. In natural environments, characterized by low flow rates and long flow

Ž

.

paths, the Cr III oxidation reaction will proceed essentially instantaneously compared to

Ž

.

Ž

.

the solute residence time; therefore, the extent of Cr III transformation to Cr VI will be
limited by the availability of unreacted MnO , but not by the rate of reaction.

2

We assume, in accordance with recent surface spectroscopic studies, that solid-phase

Ž

.

Ž

.

MnOOH forms concomitantly with the conversion of Cr III to Cr VI . Accumulation of

Ž

.

this precipitate promotes the onset of the second stage of Cr III oxidation kinetics,

Ž

.

Ž

.

signaled by decreasing rates of Cr III

transformation to Cr VI . We successfully

Ž

.

modeled the relationship between Cr III oxidation rates and MnOOH concentrations
with a Langmuir-type blocking function. Inspection of the best-fit values of G indicates
that 1 mg of MnOOH is capable of covering over 80 mg of b-MnO . The effectiveness

2

of MnOOH in excluding b-MnO from reaction can be understood by considering the

2

physical characteristics of the b-MnO

coating. Analysis of the sand by scanning

2

electron microscopy reveals that the b-MnO exists as aggregates. During the experi-

2

ments, reductive transformation of the outermost layers of the aggregates generates a
rind of MnOOH that effectively shields the bulk of the b-MnO from reaction with the

2

Ž

.

Cr III . Because MnOOH is an insoluble precipitate, it is not swept from the surface by
the advecting pore fluid. Consequently, the b-MnO that is associated with the aggre-

2

gate interiors is never exposed.

Regeneration of the b-MnO surface through MnOOH oxidation arrests the decline

2

Ž

.

Ž

.

in Cr III oxidation rates and prevents effluent Cr VI concentrations from dropping to

Ž

.

zero. The stabilization of Cr VI concentrations at non-zero levels marks the beginning
of the third stage in oxidation kinetics. During this stage, the rate of b-MnO depletion

2

Ž

.

due to Cr III oxidation equals the rate of b-MnO

regeneration due to oxidation of

2

MnOOH by dissolved oxygen. The kinetics of MnOOH oxidation to b-MnO

are

2

described well with a second-order equation that expresses the rate of reaction as a
function of the product of the concentrations of solid-phase MnOOH and dissolved
oxygen. Quantification of this oxidation reaction is critical to defining the fate of
chromium in subsurface environments, as our results suggest that sustained production

Ž

.

of Cr VI , the more mobile of the two chromium species, ultimately depends on the rate
of conversion of MnOOH to b-MnO . The oxidative capacity of aquifer materials, then,

2

is not controlled solely by the mass of MnO , but is a complex function of both MnO

2

2

availability and dissolved oxygen levels.

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

332

In a general sense, the oxidation–reduction kinetics measured in these experiments

are consistent with those observed in studies performed with other redox-sensitive
contaminants and solid-phase oxidants. Experiments on CoEDTA transport through

Ž

.

MnO -coated sand Jardine and Taylor, 1995 , for example, show the same three stages

2

of oxidation kinetics revealed in our experiments: fast initial oxidation, a brief period of
declining oxidation rates, and then a prolonged period of low-level CoEDTA oxidation.
Breakthrough data on the transport of CoEDTA through ferrihydrite-coated sand express

Ž

.

these oxidation–reduction dynamics as well Brooks et al., 1996 . In addition, inhibition

Ž .

of Fe II oxidation has been observed in batch studies and has been linked to reaction

Ž

.

products that poison the redox-reactive surface e.g., Roden and Urrutia, 1999 . Based
on these studies, we conclude that the overall structure of our mathematical model has
applicability outside of the chromium-MnO system. It is clear that continued advances

2

in our understanding of the transport of redox-sensitive metals relies on additional
experimental work aimed at elucidating the key mechanisms that govern mass transfer
kinetics and on tests of mathematical models that appropriately combine quantitative
treatments of mass transfer kinetics with expressions for transport.

Acknowledgements

The Environmental Management Science Program of the US Department of Energy

sponsored this research. We thank two anonymous reviewers for their thorough review
of the manuscript.

References

Amacher, M.C., Baker, D.E., 1982. Redox reactions involving chromium, plutonium, and manganese in soils,

Final Rept. DE-AS08-77DPO4515, Institute for research on land and water research, Pennsylvania State
University, University Park, PA.

Ž

.

Anderson, L.D., Kent, D.B., Davis, J.A., 1994. Batch experiments characterizing the reduction of Cr VI using

suboxic material from a mildly reducing sand and gravel aquifer. Environ. Sci. Technol. 28, 178–185.

Ž

.

Banerjee, D., Nesbitt, H.W., 1999a. Oxidation of aqueous Cr III at birnessite surfaces: constraints on reaction

mechanism. Geochim. Cosmochim. Acta 63, 1671–1687.

Banerjee, D., Nesbitt, H.W., 1999b. XPS Study of reductive dissolution of birnessite by oxalate: rates and

mechanistic aspects of dissolution and redox processes. Geochim. Cosmochim. Acta 63, 3025–3038.

Barlett, R.J., James, B.J., 1979. Behavior of chromium in soils: III. Oxidation. J. Environ. Qual. 8, 31–35.
Brooks, S.C., Taylor, D.L., Jardine, P.M., 1996. Reactive transport of EDTA-complexed cobalt in the presence

of ferrihydrite. Geochim. Cosmochim. Acta 60, 1899–1908.

Chang, R., 1984. Chemistry. Random House, New York, 824 pp.

Ž

.

Charlet, L., Manceau, A., 1992. X-ray absorption spectroscopic study of the sorption of Cr III at the

oxide-water interface: II. Adsorption, coprecipitation, and surface precipitation on hydrous ferric oxide. J.
Colloid Interface Sci. 148, 443–458.

Chen, W., Wagenet, R.J., 1997. Description of atrazine transport in soil with heterogeneous nonequilibrium

sorption. Soil Sci. Soc. Am. J. 61, 360–371.

Ž

.

Chung, J., Zasoski, R.J., Lim, S., 1994. Kinetics of chromium III oxidation by various manganese oxides.

Agric. Chem. Biotechnol. 37, 414–420.

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

333

Davis, J.A., Kent, D.B., Rea, B.A., Maest, A.S., Garabedian, S.P., 1993. Influence of redox environment and

aqueous speciation on metal transport in groundwater: preliminary result of tracer injection studies. In:

Ž

.

Allen, H.E., Perdue, E.M., Brown, D.S. Eds. , Metals in Groundwater. Lewis, Boca Raton, FL, pp.
223–273.

Ž

.

Ž

.

Eary, L.E., Rai, D., 1987. Kinetics of chromium III oxidation to chromium VI by reaction with manganese

dioxide. Environ. Sci. Technol. 21, 1187–1193.

Ž

.

Fendorf, S.E., Sparks, D.L., 1994. Mechanisms of chromium III sorption on silica: 2. Effect of reaction

conditions. Environ. Sci. Technol. 28, 290–297.

Ž

.

Fendorf, S.E., Zasoski, R.J., 1992. Chromium III oxidation by d-MnO : 1. Characterization. Environ. Sci.

2

Technol. 26, 79–85.

Ž

.

Fendorf, S.E., Fendorf, M., Sparks, D.L., Gronsky, R., 1992. Inhibitory mechanisms of Cr III oxidation by

d-MnO . J. Colloid Interface Sci. 153, 37–54.

2

Ž

.

Fendorf, S.E., Zasoski, R.J., Burau, R.G., 1993. Competing metal ion influences on chromium III oxidation

by birnessite. Soil Sci. Soc. Am. J. 57, 1508–1515.

Ž

.

Fendorf, S.E., Lamble, G.M., Stapleton, M.G., Kelley, M.J., Sparks, D.L., 1994. Mechanisms of chromium III

Ž

.

sorption on silica: 1. Cr III

surface structure derived by extended X-ray absorption fine structure

spectroscopy. Environ. Sci. Technol. 28, 284–289.

Fendorf, S., Jardine, P.M., Patterson, R.R., Taylor, D.L., Brooks, S.C., 1999. Pyrolusite surface transforma-

Ž .

2y

tions measured in real-time during the reactive transport of Co II EDTA

. Geochim. Cosmochim. Acta

63, 3049–3057.

Fetter, C.W., 1993. Contaminant Hydrogeology. Prentice Hall, Canada, 500 pp.

Ž

.

Hartford, W., 1983. Chromium chemicals. In: Grayson, M. Ed. , 3rd edn. Kirk–Othmer Encyclopedia of

Chemical Technology vol. 6. Wiley Interscience, New York, pp. 83–120.

Ž .

Jardine, P.M., Taylor, D.L., 1995. Kinetics and mechanisms of Co II EDTA oxidation by pyrolusite. Geochim.

Cosmochim. Acta 59, 4193–4203.

Jardine, P.M., Fendorf, S.F., Mayes, M.A., Larsen, L., Brooks, S.C., Bailey, B., 1999. Fate and transport of

chromium in undisturbed heterogeneous soil. Environ. Sci. Technol. 33, 2939–2944.

Johnson, P.R., Elimelech, M., 1995. Dynamics of colloid deposition in porous media: blocking based on

random sequential adsorption. Langmuir 11, 801–812.

Ž

.

Ž

.

Johnson, C.A., Xyla, A.G., 1991. The oxidation of chromium III to chromium VI on the surface of

Ž

.

magnetite g-MnOOH . Geochim. Cosmochim. 55, 2861–2866.

Langmuir, I., 1918. J. Am. Chem. Soc. 40, 1361.
Leckie, J.O., Appleton, A.R., Ball, N.B., Hayes, K.F., Honeyman, B.D., 1984. Adsorption removal of trace

elements from fly-ash pond effluents onto iron oxyhydroxide. Final Rep. EPRI-RP-910-1. Elec. Power
Inst., Palo Alto, CA.

Ž

.

Manceau, A., Charlet, L., 1992. X-ray absorption spectroscopic study of the sorption of Cr III at the

Ž

.

oxiderwater interface: I. Molecular mechanism of Cr III oxidation on manganese oxides. J. Colloid
Interface Sci. 148, 425–442.

Mesuere, K., Fish, W., 1992. Chromate and oxalate adsorption on goethite. 2. Surface complexation modeling

of competitive adsorption. Environ. Sci. Technol. 26, 2365–2370.

Nkedi-Kizza, P., Biggar, J.W., van Genuchten, M.Th., Wierenga, P.J., Selim, H.M., Davidson, J.M., Nielsen,

D.R., 1983. Modeling tritium and chloride 36 transport through an aggregated oxisol. Water Resour. Res.
19, 691–700.

Rao, P.S.C., Rolston, D.E., Jessup, R.E., Davidson, J.M., 1980. Solute transport in aggregated in porous

media: theoretical and experimental evaluation. Soil Sci. Soc. Am. J. 44, 1139–1146.

Ž

.

Roden, R.R., Urrutia, M.M., 1999. Ferrous iron removal promotes microbial reduction of crystalline iron III

oxides. Environ. Sci. Technol. 33, 1847–1853.

Saiers, J.E., Tao, G., 2000. Evaluation of continuous distribution models for rate-limited solute adsorption to

geological materials. Water Resour. Res. 36, 1627–1639.

Saiers, J.E., Guha, H., Jardine, P.M., Brooks, S.C., 2000. Development and evaluation of a mathematical

model for the transport and oxidation–reduction of CoEDTA. Water Resour. Res. 36, 3151–3165.

Saito, H., 1951. The manganese oxide catalysts. II. Catalytic oxidation of carbon monoxide by several kinds of

oxides of manganese. Nippon Kagaku Zasshi 59, 333–336.

background image

(

)

H. Guha et al.r Journal of Contaminant Hydrology 49 2001 311–334

334

Ž

.

Ž

.

Sass, B.M., Rai, D., 1987. Solubility of amorphous chromium III –iron III hydroxide solid solution. Inorg.

Chem. 26, 2228–2232.

Selim, H.M., Amacher, M.C., 1988. A second-order kinetic approach for modeling solute retention and

transport in soils. Water Resour. Res. 24, 2061–2075.

Skopp, J., 1986. Analysis of time-dependent chemical processes in soils. J. Environ. Qual. 15, 205–213.
Sposito, G., 1994. Chemical Equilibria and Kinetics in Soils. Oxford Univ. Press, New York, 268 pp.
van Genucthen, M.Th., Wierenga, P.J., 1977. Mass transfer studies in sorbing porous media: II. Experimental

evaluation with tritium. Soil Sci. Soc. Am J. 41, 273–278.

Ž

.

Westbrook, J., 1983. Chromium and chromium alloys. In: Grayson, M.

Ed. , 3rd edn. Kirk–Othmer

Encyclopedia of Chemical Technology, vol. 6. Wiley Interscience, New York, pp. 54–82.


Wyszukiwarka

Podobne podstrony:
11 transportid 12256 Nieznany (2)
problematyka transportu samocho Nieznany
1112 3 zl di 6 transport wewnet Nieznany (2)
Bezpieczenstwo transportu samoc Nieznany (2)
96 Rogalski Transport kolejowy Nieznany (2)
O publicznym transporcie zbioro Nieznany
Kopia Mrpd Bud sieci transp id Nieznany
116 USTAWA o transporcie kolej Nieznany
11 transportid 12256 Nieznany (2)
12 Uzytkowanie urzadzen transpo Nieznany (2)
barometr logistyki i transportu Nieznany
kolm srodki transportu 4 id 239 Nieznany
kolm srodki transportu 1 id 239 Nieznany
instrukcja transport ryzyko id Nieznany

więcej podobnych podstron