skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Chemical speciation models based upon the Pitzer activity coefficient equations, including the propagation of uncertainties. III. Seawater from the freezing point to 45 ◦C, including acid-base equilibria
A quantitative understanding of pH, acid-base equilibria, and chemical speciation in natural waters including seawater is needed in applications ranging from global change to environmental and water quality management. In a previous study (Humphreys et al., 2022) we implemented a model of solutions containing the ions of artificial seawater, based upon the use of the Pitzer equations for the calculation of activity coefficients and including, for the first time, the propagation of uncertainties. This was extended (Clegg et al., 2022) to include the Tris buffer solutions that are used to calibrate the seawater total pH scale. Here we apply the same methods to develop a model of solutions containing the ions of standard reference seawater, based upon studies by Millero and co-workers. We compare the predictions of the model to literature data for: the dissociation of dissolved CO2 and bicarbonate ion; boric acid dissociation; saturation with respect to calcite, the ion product of water, and osmotic coefficients of seawater. Estimates of the uncertainty contributions of all thermodynamic equilibrium constants and Pitzer parameters to the variance of the calculated quantity are used to determine which elements of the model need improvement, with the aim of agreeing with properties noted above to within their experimental uncertainty. Further studies are recommended. Comparisons made with several datasets for carbonate system dissociation in seawater suggest which are the most reliable, and identify low salinity waters (S <10) as a region for which dissociation constants of bicarbonate are not yet accurately known. At present, the model is likely to be most useful for the direct calculation of equilibria in natural waters of arbitrary composition, or for adjusting dissociation constants known for seawater media to values for natural waters in which the relative compositions of the major ions are different.  more » « less
Award ID(s):
1744653 1840868
PAR ID:
10423768
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Marine chemistry
Volume:
250
ISSN:
0304-4203
Page Range / eLocation ID:
1 - 21
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We have predicted acid dissociation constants (pKa), octanol–water partition coefficients (KOW), and DMPC lipid membrane–water partition coefficients (Klipid-w) of 150 different eight-carbon-containing poly-/perfluoroalkyl carboxylic acids (C8-PFCAs) utilizing the COnductor-like Screening MOdel for Realistic Solvents (COSMO-RS) theory. Different trends associated with functionalization, degree of fluorination, degree of saturation, degree of chlorination, and branching are discussed on the basis of the predicted values for the partition coefficients. In general, functionalization closest to the carboxylic headgroup had the greatest impact on the value of the predicted physicochemical properties. 
    more » « less
  2. Here, we explicitly define a half-cell reaction approach for pH calculation using the electrode couple comprised of the solid-state chloride ion-selective electrode (Cl-ISE) as the reference electrode and the hydrogen ionselective ion-sensitive field effect transistor (ISFET) of the Honeywell Durafet as the hydrogen ion (H+)-sensitive measuring or working electrode. This new approach splits and isolates the independent responses of the Cl-ISE to the chloride ion (Cl−) (and salinity) and the ISFET to H+ (and pH), and calculates pH directly on the total scale (pHEXT total) in molinity (mol (kg-soln)−1) concentration units. We further apply and compare pHEXT total calculated using the half-cell and the existing complete cell reaction (defined by Martz et al. (2010)) approaches using measurements from two SeapHOx sensors deployed in a test tank. Salinity (on the Practical Salinity Scale) and pH oscillated between 1 and 31 and 6.9 and 8.1, respectively, over a six-day period. In contrast to established Sensor Best Practices, we employ a new calibration method where the calibration of raw pH sensor timeseries are split out as needed according to salinity. When doing this, pHEXT total had root-mean squared errors ranging between ±0.0026 and ±0.0168 pH calculated using both reaction approaches relative to pHtotal of the discrete bottle samples (pHdisc total). Our results further demonstrate the rapid response of the Cl-ISE reference to variable salinity with changes up to ±12 (30 min)−1. Final calculated pHEXT total were ≤±0.012 pH when compared to pHdisc total following salinity dilution or concentration. These results are notably in contrast to those of the few in situ field deployments over similar environmental conditions that demonstrated pHEXT total calculated using the Cl-ISE as the reference electrode had larger uncertainty in nearshore waters. Therefore, additional work beyond the correction of variable temperature and salinity conditions in pH calculation using the Cl-ISE is needed to examine the effects of other external stimuli on in situ electrode response. Furthermore, whereas past work has focused on in situ reference electrode response, greater scrutiny of the ISFET as the H+-sensitive measuring electrode for pH measurement in natural waters is also needed. 
    more » « less
  3. We have predicted acid dissociation constants (pKa), octanol-water partition coefficients (KOW), and DPMC lipid membrane-water partition coefficients (Klipid-w) of 150 different 8-carbon containing poly-/per-fluoroalkyl carboxylic acids (C8-PFCAs) utilizing COMSO-RS theory. Different trends associated with functionalization, degree of fluorination, degree of saturation, degree of chlorination, and branching are discussed based upon the predicted values for the partition coefficients. In general, functionalization closest to the carboxylic head group had the greatest impact on the value of the predicted physicochemical properties. 
    more » « less
  4. Abstract The capacity of aquatic systems to buffer acidification depends on the sum contributions of various chemical species to total alkalinity (TA). Major TA contributors are inorganic, with carbonate and bicarbonate considered the most important. However, growing evidence shows that many rivers, estuaries, and coastal waters contain dissolved organic molecules with charge sites that create organic alkalinity (OrgAlk). This study describes the first comparison of (1) OrgAlk distributions and (2) acid–base properties in contrasting estuary‐plume systems: the Pleasant (Maine, USA) and the St. John (New Brunswick, CA). The substantial concentrations of OrgAlk in each estuary were sometimes not conservative with salinity and typically associated with very low pH. Two approaches to OrgAlk measurement showed consistent differences, indicating acid–base characteristics inconsistent with the TA definition. The OrgAlk fraction of TA ranged from 78% at low salinity to less than 0.4% in the coastal ocean endmember. Modeling of titration data identified three groups of organic charge sites, with mean acid–base dissociation constants (pKa) of 4.2 (± 0.5), 5.9 (± 0.7) and 8.5 (± 0.2). These represented 21% (± 9%), 8% (± 5%), and 71% (± 11%) of titrated organic charge groups. Including OrgAlk, pKa, and titrated organic charge groups in carbonate system calculations improved estimates of pH. However, low and medium salinity, organic‐rich samples demonstrated persistent offsets in calculated pH, even using dissolved inorganic carbon and CO2partial pressure as inputs. These offsets show the ongoing challenge of carbonate system intercomparisons in organic rich systems whereby new techniques and further investigations are needed to fully account for OrgAlk in TA titrations. 
    more » « less
  5. Lead (Pb) solubility is commonly limited by dissolution–precipitation reactions of secondary mineral phases in contaminated soils and water. In the research described here, Pb solubility and free Pb2+ ion activities were measured following the precipitation of Pb minerals from aqueous solutions containing sulfate or carbonate in a 1:5 mole ratio in the absence and presence of phosphate over the pH range 4.0–9.0. Using X-ray diffraction and Fourier-transform infrared spectroscopic analysis, we identified anglesite formed in sulfate-containing solutions at low pH. At higher pH, Pb carbonate and carbonate-sulfate minerals, hydrocerussite and leadhillite, were formed in preference to anglesite. Precipitates formed in the Pb-carbonate systems over the pH range of 6 to 9 were composed of cerussite and hydrocerussite, with the latter favored only at the highest pH investigated. The addition of phosphate into the Pb-sulfate and Pb-carbonate systems resulted in the precipitation of Pb3(PO4)2 and structurally related pyromorphite minerals and prevented Pb sulfate and carbonate mineral formation. Phosphate increased the efficiency of Pb removal from solution and decreased free Pb2+ ion activity, causing over 99.9% of Pb to be precipitated. Free Pb2+ ion activities measured using the ion-selective electrode revealed lower values than predicted from thermodynamic constants, indicating that the precipitated minerals may have lower KSP values than generally reported in thermodynamic databases. Conversely, dissolved Pb was frequently greater than predicted based on a speciation model using accepted thermodynamic constants for Pb ion-pair formation in solution. The tendency of the thermodynamic models to underestimate Pb solubility while overestimating free Pb2+ activity in these systems, at least in the higher pH range, indicates that soluble Pb ion-pair formation constants and KSP values need correction in the models. 
    more » « less