The Rate of Air-Sea CO2 Exchange: Chemical Enhancement and Catalysis by Marine Microalgae.

## Chapter 6: Analysis of measured CO2 data from steady-state tank.

This chapter summarises how the data collected during three years of experiments in the steady-state gas exchange tank were analysed to produce the many CO2 transfer velocity and "carbon budget" plots, which will be discussed in the following two "results" chapter 7 and chapter 8. Many of the formulae and thermodynamic and kinetic constants listed here were also used for the global flux calculations both in chapter 3 and in chapter 9. Most of the key equations and concepts presented here have already been explained earlier in this thesis. Thus this is essentially a "reference" chapter which brings all the information together in a format which should make it clearer exactly how the calculations were made, thus enabling the following chapters to be dedicated to discussion of the significance of the results.

The data from the steady-state tank experiments was all processed using the computer program "Stanford Graphics". A multilayer spreadsheet and associated graphs were set up for processing all the tank data, with modifications for each experiment as appropriate. The section f this chapter follow approximately the order in which calculations were made in this spreadsheet. Firstly, Section 6.2 lists the formulae for the thermodynamic and kinetic constants based on temperature and salinity. Section 6.3 explains the carbonate speciation based on the pCO2 and alkalinity, and also discusses the relative contributions of the various species, and the internal consistency of the carbonate system. Section 6.4 gives the formulae of the predicted CO2 transfer velocities due to both diffusion and reaction, and Section 6.6 explains how the actual CO2 transfer velocities and their associated error ranges were calculated. Finally mass-balance estimates of total and biological carbon are explained

### 6.1 Measured Quantities

Table 6-1 shows the measured quantities and the variables representing them in calculations. Details of how the quantities were measured are given in the previous chapters.

Notes to table 6-1 :

TCO2, alkalinity, pH, chlorophyll, and SF6 and O2 transfer velocities were only measured occasionally, and were calculated as described in Section 5.9 , Section 5.10 , Section 5.10 , Section 5.4 , Section 5.11 and Section 5.12 respectively. The carbonate system calculations were therefore based on pCO2 and a constant total alkalinity (except when it was known to change). The measured TCO2 and pH have been plotted on the "Carbon budget" graphs in chapter 8 for comparison with the values calculated from the pCO2 measurements.

The baseline for predicted CO2 transfer velocities due to diffusion only was derived from the temperature and rpm calibration experiments (see Section 6.4 , Section 7.2 and Section 7.3 ). Therefore SF6 / O2 gas exchange measurements are plotted only for comparison.

Day and time were combined into one variable "time" = day + time / 24 which forms the x-axis of many graphs, the axis ticks representing midnight. Since many measurements were made in the evening, the data may be closer to the label for the following "day" than the date of measurement.

The typical values given are illustrative and not intended to be internally consistent. For more detail on measurement errors see the relevant section f chapter 4 and chapter 5. Factors which might increase the error of a particular CO2 transfer velocity measurement were noted and incorporated into the calculations as described in

### 6.2 Thermodynamic and Kinetic constants

Each temperature value was used to calculate several kinetic and thermodynamic constants, which are described in table 6-2 . Note that tk is the Kelvin temperature (t+273) The Salinity (s) was assumed to be 35 �, except for milliQ water blanks, in which case it was zero. The values of A, B, C etc. are different for each constant as shown in the fifth column, while the sixth column gives the literature reference of the data. The dissociation constants refer to the following reactions.:

k1 H2CO3 <==> HCO3- + H+

k2 HCO3- <==> CO32- + H+

kw H2O <==> H+ + OH-

kb B(OH)3 + H2O <==> H+ + B(OH)4-

kp1 H3PO4 <==> H+ + H2PO4-

kp2 H2PO4- <==> H+ + HPO42-

kp3 HPO42- <==> H+ + PO43-

Each constant k is defined as [H] [A] / [HA], where HA and A are the acid and base forms of each pair respectively. They are apparent equilibrium constants based on concentrations rather than activities, valid only for one specific temperature and salinity. In other words, the activity coefficients, which are determined predominantly by the bulk seawater medium rather than the species reacting, are built into the constant. The units are mol kg-1 (solution) rather than mol l-1. For more detail on this topic refer to papers by Hansson (1973) and Dickson (1993, 1994). All the thermodynamic constants used here can be found in Dickson (1994), where the chemistry of the carbonate system is discussed extensively.

Notes on table 6-2 :

1. The Oxygen and SF6 constants were not used routinely, but are included here for convenience

2. The rate constants for kCO2 and kOHKw are taken directly from the paper of Johnson (1982). -Emerson (1995) used a much lower value for kOHKw, which I believe to be incorrect, for reasons already explained in Section 1.5.2 . Moreover, the experimental results reported in Section 7.5 suggest that the value may need to be higher, rather than lower.

3. The Schmidt number formulae given by Wanninkhof (1992) are derived from Jahne (1987b) for CO2 and from Wilke and Chang (1955) for SF6 and Oxygen. King and Saltzman (1995) gave a slightly different formula for the SF6 Schmidt number.

### 6.3 Seawater chemistry

The measured pCO2, and constant values for total alkalinity, total phosphate and total borate were combined to calculate the speciation of the seawater carbonate system. The constant value for total alkalinity was based on the average values measured for the seawater, with various exceptions which are explained in the results (e.g. Section 8.3.3 ). The value for total borate was assumed to be 0.416 mmol kg-1, the average for surface seawater (Dickson 1994). Phosphate was only included when it had been added as a nutrient, in which case the concentration (typically 40 mmol l-1, approx as in f / 2 culture medium) was calculated from the weight added. Assuming that the contributions of HF, Si(OH)4 and HSO4- are too small to be significant for the purpose of these experiments, the parameters described above are sufficient to define the seawater carbonate system.

However, an iterative process is required to calculate this speciation exactly from these parameters. This is because the inclusion of borate, phosphate, and hydroxyl ions in the total alkalinity results in a set of seventh order differential equations, which cannot be solved algebraically. On the other hand, iteration is inconvenient in a spreadsheet with many data points.

A compromise was found, whereby the hydrogen ion concentration was estimated first assuming that the carbonate alkalinity was equal to the total alkalinity. The borate and phosphate were then speciated using this hydrogen ion concentration. The carbonate alkalinity was then recalculated by deducting the contributions of the borate, phosphate and hydroxyl ions from the total alkalinity. From this second carbonate alkalinity, a second hydrogen ion concentration and carbonate speciation were calculated. The method should become clearer from a study of table 6-3 .

It was found that a third iteration (i.e. recalculating the borate and phosphate speciation using the second hydrogen ion concentration, and from that recalculating a third carbonate alkalinity and hydrogen ion concentration) made no significant difference to the TCO2 and hydroxyl ion concentrations, which are the parameters of most interest.

Figure 6-1 shows the measured TCO2 as a function of the TCO2 calculated by this method. Generally, the internal consistency seems sufficient for the purposes of this investigation

It would be confusing to include the concentrations of every chemical species on all the carbon budget graphs, and as they are in equilibrium this would not convey independent information. However, it is interesting to see how the speciation changes as a function of TCO2. This is shown in figure 6-2 for a typical algal bloom culure.

The data in figure 6-2 is taken from the fourth Dunaliella bloom. Note that the phosphate concentration, similar to that in f/2 medium, is much higher than is typical of seawater, and therefore makes a significant contribution to the total alkalinity especially at low pCO2 (or TCO2). Before the phosphate was included in the calculations as described above, it was not possible to reconcile the predicted and measured TCO2, nor to explain the unusually high measured alkalinity.

### 6.4 Predicted CO2 transfer velocity, plus uncatalysed enhancement

The unenhanced transfer velocity for CO2, i.e. the transfer velocity expected if there were no chemical reaction with water, was predicted from the temperature and the motor speed. The dependence on temperature was based on an assumed Schmidt number dependence k µ Sc-1/2 (see Section 1.2.4 ). The dependence on motor speed was based on the calibration curve derived from the measurements of the CO2 transfer velocity, at a wide range of motor speeds, in acidified seawater (in which enhancement by reaction is negligible). This is discussed further in Section 7.2 . Note that the reduced diffusivity of CO2 in seawater due to its salinity is accounted for in the Schmidt number dependence.

Once the unenhanced transfer velocity has been calculated, the predicted enhancement due to reaction can be added by employing the formula of Hoover and Berkshire (1969). This uses the carbonate speciation as calculated in table 6-3 , and also the reaction rate constants from table 6-2 . The limitations of the Hoover and Berkshire equation have already been discussed in Sections 1.5.3-4 and will be considered further in Section 7.7.

The formula has been simplified by the use of a parameter "b" in the spreadsheet. "t" is a parameter used in Hoover and Berkshire's formulation. This should be clear from table 6-4

### 6.5 Measured CO2 transfer velocities and their error ranges

The actual CO2 transfer velocities were calculated from the measured pCO2 of air flowing in and out of the headspaces, and the flow rates, according to the principle outlined in Section 4.3 . The actual formulae used here are defined in table 6-5 .

As the CO2 transfer velocity is the key quantity of interest in this study, it is important to consider the error associated with its measurement. The errors in the various measured parameters have already been noted in the previous chapters. However, the total error lies not in any one measurement but chiefly in the ratio of the two pCO2 differences (pCO2 into headspace - pCO2 in headspace) / (pCO2 in headspace - pCO2 in water). It is very difficult to express this error as a standard deviation, even if the standard deviation of the individual measurements are known. This is because the top and bottom parts of the ratio are not statistically independent, and so it is not possible to use any simple statistical formula to calculate the error in this ratio.

On the other hand, it is relatively easy to express an error range, by combining the errors in the individual measurements to determine the most extreme cases. If the errors in the measurements are assumed to be standard deviations, the error range of these extreme cases will be considerably greater than the standard deviation. However, this quantity is useful to indicate the variation in the error between data points, i.e. showing which are the more reliable.

The errors in the individual measurements also vary. The error in measuring the flow rate was taken to be inversely proportional to the flow rate (it's easier to measure a slow flow). The error in the pCO2 entering the headspace was normally taken to be the same as the LiCOR error (1ppm), because the constant pressure flowing from the pump ensured that the LiCOR cell was well flushed. On the other hand, when the pure CO2 cylinder was used to deliver a high concentration of CO2 then this pCO2 varied more (as described in Section 5.7 ) and the error due to this variation was noted at the time. Generally, for such high-CO2 flow, the error was proportional to pCO2. The error in pCO2 leaving the headspace was taken to be higher (usually 3ppm) due to the lower pressure drop across the LiCOR and hence slower flushing. The error in water pCO2 was usually taken as 2ppm, but higher values were used when pCO2 was changing rapidly, as during intense algal growth (see plots in chapter 8).

An additional error might be caused by the headspace not having reached steady-state. An exponential formula describing the gradual approach towards steady-state has already been derived in Section 4.4 . Table 6-6 shows the calculation of non-steady-state error, based on this exponential formula. These non-steady-state errors were usually so small (<0.1ppm), that they were not incorporated into the error bars plotted on the graphs. However, these calculations gave a useful indication of whether the headspaces ought to be in steady state, and generally implied that discrepancies must be due to other factors.

### 6.6 Weighted average transfer velocities

In chapter 8, the measured CO2 transfer velocities from all the different blooms of algae must be compared with each other and with the measured transfer velocities at equivalent pCO2 from the no-algae "control" experiments (see Section 8.4 ). The composite plot of all the data contains so many datapoints with error bars that it is rather confusing. Moreover, when the water pCO2 was very low, at the peak of some of the algal blooms, the efflux transfer velocities have a very large error due to the small difference between pCO2 entering and leaving the headspace, and should not be considered of equal importance to the more reliable influx transfer velocities.

Therefore to make these composite plots less confusing, I also calculated "weighted average" transfer velocities for each datapoint. Whenever two headspaces were used simultaneously, usually one for "influx" and one for "efflux", the formula was:

weighted average transfer velocity=

(flow / area)1 (pCO2 in - pCO2 out)1 - (flow / area)2 (pCO2 in - pCO2 out)2
dimensionless solubility * (pCO2 out1 - pCO2 out2 )

Where the subscripts indicate the two different headspaces.

The weight is thus apportioned according to the difference between the pCO2 flowing into and out of each headspace, which is in most cases in inverse proportion to the maximum expected error. The water pCO2 is not used. This simpler formula was more convenient to calculate in the spreadsheet than one incorporating the full "expected error" calculation (as in table 6-5 ). Where only one headspace was used the transfer velocity was calculated in the normal way.

### 6.7 Calculation of Total and Biological Carbon

Since the pCO2 and flow rates of air entering and leaving the tank were known, and were the only inputs to an otherwise closed system, it was possible to estimate the total carbon in the tank from a cumulative mass-balance calculation. Subtracting the inorganic carbon (i.e. TCO2) from the total carbon in the tank water then gives the organic carbon, which is a useful indication of the algal biomass (hence I refer to it as "biological carbon"). Measured chlorophyll was plotted alongside this biological carbon for comparison, and gives an indication of the chlorophyll / carbon ratio (see figure 8-1 ).

The method of calculating the total carbon is simple but only approximate. Basically, the steady-state of each headspace at the time of each measurement is assumed to have persisted since the previous measurement. The CO2 flux leaving the tank is subtracted from the CO2 flux entering the tank, and this difference is multiplied by the time interval, and added to the previous value for total carbon. The formulas are shown in table 6-7 .

Notes on table 6-7 :

1. The time interval was normally calculated from the data automatically. However, on various occasions this was overridden following changes in the set-up of the gas flow system.

2. As this is only an approximate calculation, the air temperature was assumed constant at 20C. Note that the air temperature is not a factor in the transfer velocity calculation.

3. the factor of 10-12 comes from the conversions: ml --> m3 (in calculating gas volume) and ppm --> atm

4. The first value for totc was set equal to TCO2, i.e. bioc = 0

This method makes several simplifying assumptions, which are only reasonable for large blooms of algae in highly-buffered seawater. Chiefly, these are:

2. The initial total carbon is assumed to be all inorganic, i.e. equal to TCO2. This implies that there is negligible organic carbon in the initial filtered seawater, and in the starter culture. Again, the latter is obviously untrue, however if the initial organic carbon were measured, the only effect would be to add a fixed offset to the biological carbon figures.

3. The assumption that the total organic carbon is a good indicator of algal biomass, implies that there is not a large pool of dissolved organic carbon, bacteria, or dead cells etc. A comparison with the chlorophyll measurements (see Section 8.2 ) should help to indicate whether this is reasonable. Note that when the seawater was filtered, the bacteria should initially be low. At the end of each bloom the algae were observed to flocculate and settle out on the base of the tank, but before this occurred the base appeared clean and the water homogenous. Notes on such observations are included with the graphs.

With a dynamic computer model of the tank, combined with initial organic carbon measurements, it would be possible to calculate the total carbon more accurately. However, since we are only using the total carbon to give an indication of the algal biomass, such accuracy is not critical for the purposes of this study.

Continue to Chapter 7:
CO2 exchange without algae: Results and discussion

If you jumped into this page from elsewhere, you may find it more convenient to go to the "frames" version of this thesis, alternatively you can use the links below to jump straight to a particular chapter: Hits since 21st April 2000