Effect of Temperature on Seawater

In this section are presented in numerical form the results for various temperature simulations on seawater.

Seawater composition use in all of the simulations (not only this one) is linked to salinity as the following:

'**************************************************************
' DOE 1994 Standard mean composition of seawater with borates *
'**************************************************************
[Cl-]  = 0.54586 * S/35  mol/kg(solution)
[Na+]  = 0.46906 * S/35  mol/kg(solution)
[Mg2+] = 0.05282 * S/35  mol/kg(solution)
[Ca2+] = 0.01028 * S/35  mol/kg(solution)
[SO4--]= 0.02824 * S/35  mol/kg(solution)
[K+]   = 0.01021 * S/35  mol/kg(solution)
[Br-]  = 0.00084 * S/35  mol/kg(solution)
[Sr2+] = 0.00009 * S/35  mol/kg(solution)
[F-]   = 0.00007 * S/35  mol/kg(solution)
[Btot] = 0.00042 * S/35  B(OH)3 + B(OH)4- mol/kg(solution)

where S is the salinity in g/(kg-soln).  The standard seawater composition as listed may be supersaturated in respect to calcite or aragonite formation. Calcite is the less soluble form of calcium carbonate, so theoretically should precipitate as first. Coral reef is however constituted by aragonite, a fact that shoul be considered.

The terms 'sat. C' and 'sat. A' in the two colums at the right, refer to the oversaturation of the solution, defined as the ratio between the concentrations product [Ca++] [CO3--] and the solubility product. If greater then 1 , the corresponding form of CaCO3 should start to precipitate. In some of the simulations the precipitation of CaCO3 is simulated to a certain extent (e.g. 10%) or to completion (100%).

To examine the equilibrium reactions considered here, just refer to the proper section.

Taking as average value of Dissolved Inorganic carbon (DIC) the value of 2.060, with the aid of SeaWaterCalc is possible to calculate the ppmCO2 which equilibrate at every temperature from 0 to 30C.

The results are presented in graphical form, to point clearly the fact that to-day the 410 ppm value (red dot) is still over the equilibrium value for 17C (average ocean temperature)


Tc = -2 to 30 step 1            ' temperature in Celsius (C)
ppmCO2 = 410      ' parts per million (in volume) of CO2 in the atmosphere
P = 0             ' pressure in atm. P=0 means ambient pressure, 1 atm.
S = 35            ' salinity, in grams of salts in 1 kg of solution
pptF = 0          ' fraction of CaCO3 which actually precipitates 
Cl  = 0.54586     ' Cl-   Mol/kg(solution)
Na  = 0.46906     ' Na+   Mol/kg(solution)
Mg  = 0.05282     ' Mg++  Mol/kg(solution)
Ca  = 0.01028     ' Ca++  Mol/kg(solution)
SO4 = 0.02824     ' SO4-- Mol/kg(solution)
K   = 0.01021     ' K+    Mol/kg(solution)
Br  = 0.00084     ' Br-   Mol/kg(solution)
Sr  = 0.00009     ' Sr++  Mol/kg(solution)
F   = 0.00007     ' F-    Mol/kg(solution)
B   = 0.00042     ' B(OH)3 + B(OH)4- Mol/kg(solution)
save .csv file = 0' (0=no  1 = yes)
---------------------------------------------------------------------------------------------------------------------------------
 *** CO2 SOLUBILITY EQUILIBRIA IN SEAWATER AND CaCO3 FORMATION ***
 Concentrations [..] are in mmol/kg              DIC = [H2CO3] + [HCO3-] + [CO3--]
 Cout>0 means CO2 is outgassed from SeaWater to Atmosphere, the opposite if Cout<0
 +----+------+-----+----+----+------+------+-------+-------+-------+-----+------+-------+------+-------+------+-----+-----+-----+
 |(C)|ppmCO2| pptF|  P |  S |   pH |  pOH |[H2CO3]|[HCO3-]|[CO3--]| DIC |[Ca++]|[CaOH+]|[Mg++]|[MgOH+]| Cout |CaCO3|sat.C|sat.A|
 +----+------+-----+----+----+------+------+-------+-------+-------+-----+------+-------+------+-------+------+-----+-----+-----+
 | -2 |  410 |0.000|   0|35.0| 8.051| 6.354|  0.028|  2.167|0.08615|2.281|10.280| 0.0001|52.805| 0.0146| 0.000|0.000|  2.1|  1.3|
 | -1 |  410 |0.000|   0|35.0| 8.056| 6.299|  0.027|  2.155|0.09052|2.273|10.280| 0.0001|52.804| 0.0161| 0.009|0.000|  2.2|  1.4|
 |  0 |  410 |0.000|   0|35.0| 8.061| 6.246|  0.026|  2.143|0.09500|2.264|10.280| 0.0002|52.802| 0.0177| 0.018|0.000|  2.3|  1.4|
 |  1 |  410 |0.000|   0|35.0| 8.065| 6.192|  0.025|  2.130|0.09959|2.255|10.280| 0.0002|52.800| 0.0195| 0.027|0.000|  2.4|  1.5|
 |  2 |  410 |0.000|   0|35.0| 8.070| 6.140|  0.024|  2.117|0.10428|2.245|10.280| 0.0002|52.799| 0.0214| 0.036|0.000|  2.5|  1.6|
 |  3 |  410 |0.000|   0|35.0| 8.074| 6.088|  0.023|  2.104|0.10906|2.236|10.280| 0.0002|52.796| 0.0235| 0.046|0.000|  2.6|  1.6|
 |  4 |  410 |0.000|   0|35.0| 8.078| 6.037|  0.022|  2.090|0.11393|2.226|10.280| 0.0002|52.794| 0.0257| 0.056|0.000|  2.7|  1.7|
 |  5 |  410 |0.000|   0|35.0| 8.082| 5.986|  0.021|  2.076|0.11888|2.216|10.280| 0.0003|52.792| 0.0282| 0.066|0.000|  2.8|  1.8|
 |  6 |  410 |0.000|   0|35.0| 8.085| 5.937|  0.021|  2.061|0.12390|2.206|10.280| 0.0003|52.789| 0.0307| 0.076|0.000|  3.0|  1.9|
 |  7 |  410 |0.000|   0|35.0| 8.088| 5.887|  0.020|  2.046|0.12898|2.195|10.280| 0.0003|52.786| 0.0335| 0.086|0.000|  3.1|  1.9|
 |  8 |  410 |0.000|   0|35.0| 8.092| 5.839|  0.019|  2.031|0.13413|2.184|10.280| 0.0003|52.783| 0.0365| 0.097|0.000|  3.2|  2.0|
 |  9 |  410 |0.000|   0|35.0| 8.094| 5.791|  0.019|  2.016|0.13932|2.173|10.280| 0.0004|52.780| 0.0397| 0.108|0.000|  3.3|  2.1|
 | 10 |  410 |0.000|   0|35.0| 8.097| 5.744|  0.018|  2.000|0.14455|2.162|10.280| 0.0004|52.777| 0.0431| 0.119|0.000|  3.4|  2.2|
 | 11 |  410 |0.000|   0|35.0| 8.100| 5.697|  0.017|  1.984|0.14981|2.151|10.280| 0.0005|52.773| 0.0468| 0.130|0.000|  3.6|  2.3|
 | 12 |  410 |0.000|   0|35.0| 8.102| 5.651|  0.017|  1.967|0.15510|2.139|10.280| 0.0005|52.769| 0.0506| 0.142|0.000|  3.7|  2.4|
 | 13 |  410 |0.000|   0|35.0| 8.104| 5.605|  0.016|  1.951|0.16039|2.127|10.279| 0.0005|52.765| 0.0548| 0.154|0.000|  3.8|  2.4|
 | 14 |  410 |0.000|   0|35.0| 8.106| 5.560|  0.016|  1.934|0.16569|2.115|10.279| 0.0006|52.761| 0.0591| 0.166|0.000|  3.9|  2.5|
 | 15 |  410 |0.000|   0|35.0| 8.107| 5.516|  0.015|  1.917|0.17097|2.103|10.279| 0.0006|52.756| 0.0638| 0.178|0.000|  4.1|  2.6|
 | 16 |  410 |0.000|   0|35.0| 8.109| 5.472|  0.015|  1.899|0.17624|2.090|10.279| 0.0007|52.751| 0.0687| 0.191|0.000|  4.2|  2.7|
 | 17 |  410 |0.000|   0|35.0| 8.110| 5.429|  0.014|  1.882|0.18148|2.077|10.279| 0.0008|52.746| 0.0739| 0.204|0.000|  4.3|  2.8|
 | 18 |  410 |0.000|   0|35.0| 8.111| 5.386|  0.014|  1.864|0.18668|2.064|10.279| 0.0008|52.741| 0.0794| 0.217|0.000|  4.5|  2.9|
 | 19 |  410 |0.000|   0|35.0| 8.112| 5.344|  0.014|  1.846|0.19183|2.051|10.279| 0.0009|52.735| 0.0852| 0.230|0.000|  4.6|  3.0|
 | 20 |  410 |0.000|   0|35.0| 8.113| 5.303|  0.013|  1.827|0.19693|2.038|10.279| 0.0010|52.729| 0.0913| 0.244|0.000|  4.7|  3.1|
 | 21 |  410 |0.000|   0|35.0| 8.113| 5.262|  0.013|  1.809|0.20195|2.024|10.279| 0.0011|52.722| 0.0977| 0.258|0.000|  4.8|  3.1|
 | 22 |  410 |0.000|   0|35.0| 8.114| 5.222|  0.013|  1.790|0.20689|2.010|10.279| 0.0011|52.716| 0.1045| 0.272|0.000|  5.0|  3.2|
 | 23 |  410 |0.000|   0|35.0| 8.114| 5.182|  0.012|  1.772|0.21174|1.995|10.279| 0.0012|52.708| 0.1116| 0.286|0.000|  5.1|  3.3|
 | 24 |  410 |0.000|   0|35.0| 8.114| 5.142|  0.012|  1.753|0.21649|1.981|10.279| 0.0013|52.701| 0.1190| 0.300|0.000|  5.2|  3.4|
 | 25 |  410 |0.000|   0|35.0| 8.114| 5.104|  0.012|  1.733|0.22114|1.966|10.279| 0.0014|52.693| 0.1267| 0.315|0.000|  5.3|  3.5|
 | 26 |  410 |0.000|   0|35.0| 8.113| 5.065|  0.011|  1.714|0.22566|1.951|10.278| 0.0015|52.685| 0.1348| 0.330|0.000|  5.4|  3.6|
 | 27 |  410 |0.000|   0|35.0| 8.113| 5.028|  0.011|  1.695|0.23006|1.936|10.278| 0.0016|52.677| 0.1433| 0.345|0.000|  5.6|  3.7|
 | 28 |  410 |0.000|   0|35.0| 8.112| 4.990|  0.011|  1.676|0.23433|1.921|10.278| 0.0018|52.668| 0.1521| 0.361|0.000|  5.7|  3.8|
 | 29 |  410 |0.000|   0|35.0| 8.111| 4.954|  0.011|  1.656|0.23845|1.905|10.278| 0.0019|52.659| 0.1612| 0.376|0.000|  5.8|  3.9|
 | 30 |  410 |0.000|   0|35.0| 8.110| 4.917|  0.010|  1.637|0.24241|1.889|10.278| 0.0020|52.649| 0.1708| 0.392|0.000|  5.9|  3.9|
 +----+------+-----+----+----+------+------+-------+-------+-------+-----+------+-------+------+-------+------+-----+-----+-----+
---> END OF PROGRAM !          e-mail and info --> mazzad50@gmail.com