[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

RE: Phreeqc modeling of kinetic Fe(OH)3 reductive dissolution

David: thanks very much for your quick reply. I was able to successfully
implement the input file you provided!!! However, I would in fact like to
try to depict the Fe(OH)3 reductive dissolution reaction kinetically, i.e.
in a manner analogous to how Postma and Appelo (2000) modeled MnO2 reduction
by Fe2+, according to a standard surface-area controlled reaction mechanism.
Could you show me how to set this one up? And then field an important
question: within this kinetically-controlled framework is it still possible
to model surface complexation by the diminishing Fe(OH)3? I really
appreciate your assistance, and will by all means acknowledge it in the
paper I'm working on.


-----Original Message-----
From: David L Parkhurst [mailto:dlpark@xxxxxxxx]
Sent: Wednesday, February 13, 2002 12:46 PM
To: Roden, Eric E
In-Reply-To: <OF4165DB05.EF86737E-ON87256B5F.0071110D@xxxxxxxx>
Subject: Re: Phreeqc modeling of kinetic Fe(OH)3 reductive dissolution

There are other ways to set this up, but it is simplest if you consider the
consumption of organic carbon as the kinetic reaction. Then all you must do
is add the correct amount of CH2O and let equilibrium do the rest. Signs
may be confusing, but positive "moles" times positive coefficient of CH2O
is positive, which means CH2O is added to the solution. Fe(OH)3(a)
dissolves to equilibrium in response to the addition of CH2O. Carbon from
CH2O must end up as C(IV)(carbonate species) and/or C(-IV) (methane),
because those are the only aqueous species for carbon that are defined to
the program. With Fe(OH)3(a) present, the carbon ends up in carbonate
(oxidized) and Fe(OH)3 is reduced because of thermodynamics.

Another reason to set it up this way is that the number of surface
complexation sites can be related (proportional) to the number of moles of
an EQUILIBRIUM_PHASE, so as you develop this model, you can allow for
surface sites that diminish as Fe(OH)3 dissolves.

     pH   7
     Na   10
     C    10 charge          # note total C ends up at 1.2 mmol for charge
balance at pH 7
     Fe(OH)3(a)     0.0  .01
     -formula  CH2O 1
     -steps 86400 172800 345600 691200 1382400 # 1, 2, 4, 8, 16 days
10 m_feoh3 = EQUI("Fe(OH)3(a)")
20 k = .05/(24*3600)    # constant /sec
30 rate = k*m_feoh3
40 moles = rate*TIME    # TIME in sec
50 SAVE moles

David Parkhurst (dlpark@xxxxxxxx)
U.S. Geological Survey
Box 25046, MS 413
Denver Federal Center
Denver, CO 80225


                    "Roden, Eric

                    E"                   To:     "'dlpark@xxxxxxxx'"
                    <Eric.Roden@x        cc:

In-Reply-To: <OF4165DB05.EF86737E-ON87256B5F.0071110D@xxxxxxxx>
                    nl.gov>              Subject:     Phreeqc modeling of
kinetic Fe(OH)3 reductive        


                    11:20 AM



Dear Dr. Parkhurst:

I wonder if I might trouble you with a question regarding the use of
(version 1.5.08 for Windows) for a paper I'm working on to be submitted to
Applied Geochemistry. This is actually the same problem, more-or-less,
I had written to you about sometime ago, as I was just starting to try my
hand at mixed equilibrium-kinetic geochemical modeling. I wish to simulate
changes in pH and mineral precipitation accompanying kinetic reductive
dissolution of hydrous ferric oxide [HFO, assumed to be represented by
Fe(OH)3] by dissmilatory Fe(III)-reducing bacteria in HCO3-buffered medium.
An approximate set of starting conditions for such a simulation is as

pH 7
10 mM NaHCO3
10 mmol/L Fe(OH)3

We have good evidence that bacterial HFO reduction follows first-order
kinetics, with a rate constant of ca. 0.05/d in our normal culture systems.
Hence, the initial idea would be to predict changes in pH and dissolved
inorganic carbon speciation assuming that HFO reductive dissolution took
place according the following reaction:

Fe(OH)3(s) + 0.25 CH2O + 1.75H+  =>  Fe2+(aq) + 0.25HCO3- + 2.5H2O

in which the labile organic carbon (CH2O) is present in excess and
does not need to be tracked during the simulation, and in which Fe(OH)3
consumption is depicted as a first-order reaction

RFe(OH)3(s) = -k[Fe(OH)3(s)]

Alternatively, the reaction could be modeled as a surface-area controlled
process according to a more standard formulation such as

RFe(OH)3 = -k(m/m0)^n

Ultimately I would like to include other heterogeneous reactions and their
impact on pH changes during HFO reduction, namely: (i) H+ complexation by
Fe(OH)3 surfaces, whose abundance would be declining as reductive
dissolution takes place, e.g. as a first approximation in direct proportion
to molar/mass concentration and an assumed surface area of 600 m2/g; (ii)
precipitation of FeCO3(s) (siderite) and/or Fe(OH)2.

I have made some initial attempts to set-up a Phreeqc simulation of the
simple case where only reductive dissolution of Fe(OH)3 takes place, but I
am uncertain about how to depict the parallel consumption of the mineral
phase and H+, together with production of Fe2+(aq) and HCO3-.

I would be very grateful if you could provide assistance in setting-up this
problem. If I can get this to work, I can think of a great many uses for
Phreeqc in understanding and interpreting experimental studies of bacterial
Fe(III) oxide reduction and associated aqueous/solid-phase geochemical

Best regards, Eric

Eric E. Roden
Current address (01/07/02 to 05/10/02):
Analytical Microbiology
Pacific Northwest National Laboratory
900 Battelle Blvd, Mail Stop P7-50
Richland, WA   99352
(509) 373-1043 (office)
(509) 376-5154 (lab)
(509) 376-1321 (fax)
(509) 627-0118 (home)

Permanent address:
The University of Alabama
Department of Biological Sciences
A122 Bevill Bldg 7th Ave
Tuscaloosa, AL 35487-0206
(205) 348-0556 (office)
(205) 348-1813 (lab)
(205) 348-1403 (fax)
(205) 349-1134 (home)

Project Home Page
Complete Water Resources Division Software
USGS Home Page
Water Resources Division Home Page
NRP Home Page
Help Page
USGS Privacy Statement       

Please note that some U.S. Geological Survey (USGS) information accessed through this page may be preliminary in nature and presented prior to final review and approval by the Director of the USGS. This information is provided with the understanding that it is not guaranteed to be correct or complete and conclusions drawn from such information are the sole responsibility of the user.

Any use of trade, product, or firm names in this publication is for descriptive purposes only and does not imply endorsement by the U.S. Government. 

The URL of this page is: https://wwwbrr.cr.usgs.gov/projects/GWC_coupled/phreeqc/mail/msg00389.html
Last modified: $Date: 2005-09-13 21:04:21 -0600 (Tue, 13 Sep 2005) $
Visitor number 4109 since Jan 22, 1998.