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