ENVS335: Global Carbon Cycle
Hello, dear friend, you can consult us at any time if you have any questions, add WeChat: daixieit
ENVS335: Global Carbon Cycle
Week 4 workshop: Introduction to CO2Sys
Feburary 2023
What is CO2Sys: Performing carbon chemistry calculations by hand is time consuming and prone to human error. In 2007, two scientists developed a set of spreadsheets to calculate carbon chemistry in seawater. The spreadsheet was developed into a software called CO2Sys, which is available for a suite of platforms including Microsoft Excel, Matlab, Python etc. In this workshop, you will access and download the Microsoft Excel version and learn how to input data and perform calculations. You will need to use CO2Sys to complete assignment 1.
Step 1. Access and download CO2Sys for excel
a. Navigate to the following website:
https://cdiac.ess-dive.lbl.gov/ftp/oceans/co2sys/CO2SYS_calc_XLS_v2.3/
b. Download the file named CO2Sys_v2.3.xls and save to your local hard drive or elsewhere.
CO2 sys v2.3 is used to calculate carbon chemistry in seawater using a series of macros. If you want more information or get stuck, please visit the document on Canvas titled CO2SYS_Macro Instructions.pdf.
Step 2. Open CO2Sys v2.3
(a) Double click to open the excel file. You will be asked to disable or enable macros. You must ENABLE MACROS for the programme to work. If you do not enable macros, the calculations will not work. If you make a mistake and disable macros, simply close the file and reopen it again.
Step 3. Navigating CO2Sys v2.3
The excel file contains three worksheets:
1. INFO: provides details of the software development, version and values used in the software, including references.
2. INPUT: provides information on the source of constants used in the calculations
3. DATA: this is the most useful sheet as it allows the user to INPUT information and after a series of calculations, provides a comprehensive suite of results.
(a) Please check that you have the following settings selected:
1. INFO: you do not need to change anything on this worksheet
2. INPUT: ensure you have the following selected
Set of constants |
In column A |
select K1, K2 from Millero 2010 |
KHSO4 |
In column B |
select Dickson 1990 |
KHF |
In column C |
select Dickson and Riley 1979 |
pH Scale |
In column D |
select seawater scale |
BT Value |
In column E |
select Lee et al., 2010 |
Step 4. Inputting data into CO2Sys and calculating carbon chemistry.
(a) Navigate to the DATA worksheet. In columns labelled INPUT CONDITIONS (A to E), you can insert data on up to 5 parameters: Salinity, Temperature (。C), Pressure (or depth, assumed 0 if left blank). Note that total Si (silicate) and Total P (phosphate) are optional (assumed 0 if left blank)
(b) In OUTPUT CONDITIONS (F and G), you can insert new temperature parameters. For example, if ocean temperature increases from 5 to 8。C, the INPUT condition would be 5。C and the OUTPUT condition would be 8。C. The results would be reported in the RESULTS (output conditions). However, you can change temperature and other properties in Input Conditions only (A to E).
(c) CO2Sys will only work if you input at least 2 out of the 4 carbon variables (H to L, fCO2 and pCO2 are the same). The carbon properties are as follows:
TA = total alkalinity
TCO2 = total inorganic carbon (equivalent to DIC)
pH = seawater pH
pCO2 or fCO2 = partial pressure of CO2 or fugacity of CO2
To clear data (yellow cell) and clear results (green cell), click in the relevant cell.
To run the macro, press the red START bar at the top of the spreadsheet. You will be asked if you have entered the data correctly. Say YES.
If you insert data into the INPUT CONDITIONS only, results will be reported in RESULTS (input conditions).
Now let’s check that you are using the spreadsheet correctly.
Check Test 1:
Insert the temperature, salinity, pH and pCO2 from Table 1 and task 2 into the INPUT and CONDITIONS columns on the DATA worksheet. Assume 0 total P and 0 total Si. The results for TA (total Alkalinity), TCO2 (or DIC), carbonate ion, and Revelle Buffer Factor are reported in Table 2.
If your answers are correct, move to check test 2.
If your answers are incorrect, please ask for help.
Table 1. Summary of parameters in the surface ocean (0m)
Temperature (。C) |
Salinity |
pH |
pCO2 |
10 |
34.5 |
8.05 |
398 |
20 |
34.5 |
8.05 |
398 |
20 |
27 |
8.05 |
398 |
5 |
27 |
8.05 |
398 |
|
|
|
|
Table 2. Summary of results using data from Table 1.
Temperature (。C) |
TA ( μmol kg- 1) |
TCO2 (μmol kg- 1) |
Carbonate ( μmol kg- 1) |
Revelle Buffer Factor |
Ωarg |
10 |
2334.3 |
2144.2 |
138.1 |
12.528 |
2.10 |
20 |
2333.2 |
2067.6 |
188.6 |
10.381 |
2.94 |
20 |
2183.1 |
1977.5 |
154.0 |
11.734 |
2.47 |
5 |
2243.1 |
2125.4 |
95.9 |
15.873 |
1.47 |
|
|
|
|
|
|
Check Test 2: Compare the difference in the Revelle buffer factor, pH and Ωarg in the low and high latitude ocean using properties in Table 3 when (a) pCO2 is 390 μatm and (b) when pCO2 increases to 420 μatm and (c) if pCO2 increases to 420 μatm and temperature increases by 5。C. Results are reported in Table 4.
Note: to change the temperature (part c), you can insert 25 and 8。C into the input column and 30 and 13。C into the output column, then run the macro. You then need to look at the OUTPUT results for the effect of increased temperature.
Table 3. Summary of input parameters in the surface of the low and high latitude ocean
Location |
Temperature (。C) |
Salinity |
Phosphate (μmol kg- 1) |
Silicate (μmol kg- 1) |
DIC (μmol kg- 1) |
Low latitude ocean |
25 |
35.2 |
0.01 |
1 |
2050 |
High latitude ocean |
8 |
32.5 |
1.8 |
8 |
2200 |
Table 4. Summary of results using information in Table 2.
Location |
Experiment conditions |
Revelle Buffer Factor |
pH |
Ωarg |
Low latitude ocean |
25。C, 390 μatm 25。C, 420 μatm 30。C, 420 μatm |
9.379 9.649 8.919 |
8.056 8.027 8.028 |
3.59 3.38 4.08 |
High latitude ocean |
8。C, 390 μatm 8。C, 420 μatm 13。C, 420 μatm |
13.293 13.752 12.381 |
8.068 8.038 8.048 |
2.00 1.87 2.31 |
If your answers are correct, you are ready to start the assignment
If your answers are incorrect, please ask for help.
2024-02-20