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

Re: an interesting water

>I ran this file through PHREEQC to get a CO2 number. I wanted to plot the
water on a Phase diagram which uses pCO2 , to find out if We are in the
calcite or the dolomite field. at 25 deg the file was not converging, then
I realized that the field temperature was 15 degs C. It runs at 15 and at
every temperature up to and including 20. At 21 degrees it no longer
converges. I thought this was curious. I am wondering if increasing the
allowable iterations will allow convergence. Or do I recalculate the phase
diagram for 15 deg C?

Hi Jorie,

The problem you are running into is related to alkalinity. At pH 12.3,
there is a lot of OH- alkalinity, nearly equal to the alkalinity that you
input (903 mg/L). As you raise the temperature, the OH- plus smaller
contributions from Al, Fe, NH3, and OH- complexes is greater than the input
alkalinity. The program is trying to find the right amount of carbon for
the given alkalinity and by the time you get to 21 degrees, there can be no
carbon contribution to the alkalinity with the given pH and total
alkalinity. The program just fails.

I'm not sure what the best way to go about this is. If you enter TIC [C(4)]
instead of alkalinity, you will not have a problem, but you may not know
TIC. You could set it by calcite equilibrium or fix the pCO2, but that may
not be appropriate either. Anyway, the problem is that the alkalinity and
the pH become incompatible at higher T. You can write or give me a call if
you want to talk about it.

Finally, you can look directly at saturation indices for calcite and
dolomite to see whether they should form, rather than going through some of
the assumptions needed for the activity-activity diagram.


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

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/msg00176.html
Last modified: $Date: 2005-09-13 21:04:21 -0600 (Tue, 13 Sep 2005) $
Visitor number 3469 since Jan 22, 1998.