SCOPE 13 - The Global Carbon Cycle

12

Carbon in the Freshwater Cycle 

S. KEMPE
 
ABSTRACT
12.1 INTRODUCTION
12.2 THE WATER CYCLE
12.3 CARBON IN PRECIPITATION
12.4 VOLCANIC CO2
12.5 SOIL AIR
12.6 TYPES OF RUNOFF
12.7 GROUNDWATER IN KARST ROCKS
12.8 DEEP GROUNDWATER
12.9 RIVERS
12.10 REVISION OF LIVINGSTONE'S (1963) RIVER LOADS
12.11 POSSIBLE SINKS OF CARBON IN THE FRESHWATER CYCLE 
REFERENCES

ABSTRACT

Annually 1.08 x 1020 g of water precipitates on to the continents, of which 0.376 x 1020 g runs off in rivers. This water is the motor of erosion and terrigenous biological activity alike. The amount of CO2 in rain is small (0.065 x 1015 g C/year) compared with the amount which enters the-water in the soil (net flux reaching the oceans: 0.23 x 1015 g C/year). Some 0.8 % of all CO2 produced by root respiration and microbial activity reaches the ocean by water transport. The gross flux from soil to groundwater, however, is larger, because much of the CO2 is given off when groundwater reappears in springs. CO2 pressures are 100 times larger in soil air, and 10 times larger in shallow groundwater and in rivers, than in the atmosphere. Both ground and river waters show pronounced seasonal variability of CO2 pressure and alkalinity. An estimated 20% of carbonate rocks occurs in the crust, but only 4% of the continental area shows karst features. These areas play a distinct role in the carbonate dissolution rate; Europe, the continent with the highest karst percentage, has the largest carbonate concentrations in its rivers.

12.1 INTRODUCTION

The continental fresh waters contain the most diversified carbon pool on earth. Not only does fresh water act as a storage for various inorganic and organic compounds, but it is also a chief transport medium between some of the sizewise more important carbon storages in the biota, lithosphere, and oceans. Furthermore, water is a necessary reaction partner of carbon dioxide, both in the carbonate system and in the biological formation of organic matter. To carry out photosynthesis, land biota is totally dependent on the availability of water (see Chapter 8, this volume, Figure 8.2).

This chapter provides an introduction to the global water cycle, and a discussion of the different water compartments and their respective carbon contents.

12.2 THE WATER CYCLE

Fresh water is available as precipitation (rain, snow, or dew), ice, rivers, lakes, and ground water. The best values of the global water cycle currently available have been calculated by Baumgartner and Reichel (1975). Figure 12.1 illustrates the global water cycle: annually 4.25 x 1020 g of water evaporates from the oceans; this is equivalent to a water layer, 117.8 cm thick, spread over the total ocean surface. A total of 0.397 x 1020 g is transferred from the oceans to the continents (11 cm). The continents receive a precipitation of 1.081 x 1020 g annually, equivalent to a layer of 74.6 cm covering the total land area. Of this amount, two-thirds re-evaporate and one-third returns in rivers and glaciers to the sea. The global ratio between evaporation and runoff is 0.357.

 Figure 12.1 The global water cycle (drawn after data by Baumgartner and Reichel, 1975). P = precipitation, E = evaporation, r.t. = residence time

Table 12.1a lists some estimates of the sizes of various water reservoirs. It must be noted that there is no accepted value for the total volume of rivers and lakes available; estimates deviate by as much as a factor of 10. But it is this relatively tiny fraction of water which man most urgently needs, and which he most criminally misuses. Even less information is available for groundwater volumes. Estimates depend on assumptions of mean porosity, mean sediment depth, and mean depth of groundwater table.

Table 12.1b compiles the data for continental runoff and its carbon load. The most recent attempt to make a global calculation of the amount of river loads was made by Livingstone (1963). In Section 12.10 some corrections are made to the Livingstone calculations: these corrections include updated values for river discharge. Corrections on mean dissolved matter concentrations of the major rivers also seem to be necessary. For many rivers, long-term weighted means are available today, while Livingstone used only data from random sampling.

As rivers provide the bulk transportation from land to ocean (Table 12.2), their average loads are of primary interest for the carbon cycle. The second largest discharge rate of eroded material is due to glaciers, while minor amounts are transported by ground water, atmospheric dust, volcanic events, and shore erosion.

Table 12.1a Water volumes of the earth in solid, liquid, and gaseous forms. (Table 1, Baumgartner and Reichel, 1975. Reproduced by permission of the R. Oldenbourg-Verlag GmbH, München)


1015 g

%

Oceans
1 348000 000

97.39

Polar ice caps, icebergs, glaciers

27 820 000

2.01

Groundwater, soil moisture

8062000

0.58

Lakes and rivers

225 000

0.02

Atmosphere

13 000

0.001

Total
1 384 120 000

100.00

Fresh water

36 020 000

= 2.60

Fresh water as a percentage of total
Polar ice caps, icebergs, glaciers

77.23

Ground water to 800 m depth

9.86

Ground water from 800 to 4000 m depth

12.35

Soil moisture

0.17

Lakes (fresh water)

0.35

Rivers

0.003

Hydrated earth minerals

0.001

Plants, animals, humans

0.04

Total

100.00


Various estimates of masses of the hydrosphere (Units in 1015 g H2O)
 
Lvovich Flohn Garrels et al. %
(1970) (1973) (1975)

Oceans 1 370 000 000 1 370 000 000 1 370 000 000 80
Pore water in 60 000 000 320 000 000 18.8
rocks
Soil water 21 000
Groundwater 4 00 0000
(750 m)
Ice  24 400 000 27 000 000 16 500 000 1.2
Lakes and rivers 231 000 116 000 34 000 0.002
Atmosphere  12 400 10 500 0.0006

12.3 CARBON IN PRECIPITATION

The freshwater cycle starts with 1020 g of water falling annually as precipation on the continents (Fig. 12.1). The precipitation should chemically assume equilibrium with the atmosphere.

Table 12.1b Annual river discharge and carbon content


Livingstone (1963) (see Section 12.10)
Baumgartner and
Reichel (1975)

Discharge
Salinity  HCO3
C
Discharge
Continent
1015 g H2O
 ppm
ppm
ppm    1015 g H2 O

Europe
2 498.5
182
95
18.7
2 800
Asia
11 108.5
142
79
15.5
12 200
N. America
4 557.3
142
68
13.3
5 900
S. America
8008.4
69
31
6.1
11 100
Africa
5 901.3
121
43
8.5
3 400
Australia
316.3
59
31.6
6.2
2 400
(Antarctica as ice
2000)
Recalculated sums and
32 390.1
122.4
59.8
11.8
37 700
weighted means
Dissolved Organic Carbon (DOC):
3.28
(Garrels et al., 1975)
Particulate Organic Carbon (POC):
1.76
Discharge of carbon (in 1015 g/year) as
HCO3
DOC
POC
from Livingstone:
0.382
Baumgartner and Reichel:
0.445
0.123
0.066
Including new data from Amazon River:
0.454

With an atmospheric CO2 pressure (PCO2) of 0.0003 parts CO2 in 1 part of air, we can, according to Henry's law, expect the CO2 content of rainwater to be between I ppm at 0 °C and 0.5 ppm at 20 °C. Because continents and precipitation are unevenly distributed over the climatic zones of the globe (Figure 12.2 a, b), it is difficult to assess a global mean precipitation temperature. At 15 °C (0.6 ppm CO2), the annual flux with precipitation to the continent surface should amount to

0.6 x 10-6 CO2 x 1.08 x 1020 g H2 O = 0.065 x 1015 g CO2 (0.018 x 1015 g C/year or 0.0014 x 1015 moles CO2 )

Part of the CO2 is instantaneously lost again into the atmosphere, when the precipitation evaporates from rocks, soil, or vegetation. In this context, the ratio between physical evaporation and biological evapotranspiration should be of interest.

In the more densely populated humic zones of the earth, industrial fossil-fuel burning increases the CO2 content of air locally and rain may contain larger CO2 concentrations as compared with the global average.

Actual measurements of CO2 in rainwater (Table 12.3) show, in fact, a much larger content than expected from the equilibrium argument. The real flux of CO2 from atmosphere to land with precipitation may, therefore, be several times larger than previously anticipated. No global data are available on rain concentrations of organic carbon.

Table 12.2 Agents of material transport to oceans (After Garrels et al., 1975, Table 10. Reproduced by permission from Chemical Cycles and the Global Environment, by Garrels, Mackenzie and Hunt. Copyright © 1975 by William Kaufman, Inc., Los Altos, California. All rights reserved)


% of total 
Agent transport Remarks

Streams  89 Present dissolved load 17%, suspended 72%;
during geologic past more nearly equal.
 
Groundwater  2 Estimate poor, dissolved materials like those
of streams. Major area of ignorance with 
respect to possible contamination.
 
Dust  0 .2 Dust to ocean related to deserts and wind
patterns. Sahara major source for tropical
Atlantic. Composition similar to average
sedimentary rock; many dusts have high
(30%) organic content.
 
Shore erosion  1 Silts, muds, and sands eroded from shore-
lines by waves, tides, and currents. Composition
 like suspended load of streams.
Ice 7 Ground-up rock debris as well as material up
to sizes of boulders. Chiefly from Antarctica
and Greenland. Distributed in northern and
southern seas by icebergs. Composition
similar to average sediments.
Volcanic 0 .3 (?) Lavas and gases transported from earth's
interior. Amounts and compositions of
gases poorly known, but include CO2, CH4,
H2S, SO2, NH3, H2. Dusts from explosive
volcanoes may be important in climatic
control. No one knows how much material
from volcanoes is new to exogenic cycle.

CO2 in rain alone cannot account for the HCO3 found in rivers. According to 

H2O + CO2 + CaCO3 Ca(HCO3)2 Ca 2+ + 2HCO3

at least half of the stream HCO3 must originally be derived from the atmosphere, i.e. 0.23 x 1015  g C (Table 12.1b).

Figure 12.2 (a) Areas of oceans and continents versus latitude in 5° intervals. (b) Precipitation, evaporation, and runoff from continents versus latitude in 5 ° intervals (drawn after data by Baumgartner and Reichel, 1975)

Table 12.3 CO2 content of precipitation (Miotke 1968. Reproduced by permission of the author)


1. Rainwater
Mg CO2 /1
McLeod (1869)
2.7
Frankland (1874)
2.5
Baumert (1877) in Paul, B. H.
0.7
Peligot (1877) in Paul, B. H.
1.0
H. Lehmann (1956) Cuba
2.53.5
Bögli (1961) Switzerland
2.22.6
Miotke (1965) Picos, Spain
2.2
 
2. Meltwater from snow
Bögli (1951) fresh meltwater
0.851.76
Bögli (1961) from vegetated soil
2.643.63
Miotke (1965) Picos, fresh meltwater
0.4

12.4 VOLCANIC CO2

Whether or not volcanic CO2 or CO2 from hydrothermal sources play a quantitatively important role is not yet known. Recycled CO2, emitted by volcanic vents or moffettes, directly adds to the atmophere. The estimates of this process differ widely. In certain areas, however, there exists a diffuse CO2 flow through the crust. These areas are marked by rivers and lakes which carry more bicarbonate than alkaline earth metals. Soda lakes have developed by this process in Ethiopia, Kenya, Tanzania, and in Eastern Turkey (Lake Van; Kempe, 1977). Crater lakes often have sub-lacustrine CO2 sources, adding volcanic carbon (e.g. Laacher See, Federal Republic of Germany; Lake Kivu, East Africa).

12.5 SOIL AIR

Additional atmospheric CO2 dissolves in rainwater, only if rain comes into direct contact with alkaline earth metal carbonate rocks. Rain running down bare limestone creates karren by gradual dissolution of CaCO3 and slow uptake of atmospheric CO2 Bögli, 1975). The final carbonate concentration, however, which these runoff films reach, is far below that of ground or river water.

The major part of CO2 responsible for weathering of carbonates and silicates (see Chapter 13, this volume) enters precipitation water when it percolates through soil. Due to the immense total surface of soil particles, soil water is in equilibrium with soil air. At this point a substantial amount of CO2 enters the freshwater cycle.

In soils, bacterial oxidation decomposes the photosynthetically produced organic matter to CO2. Root respiration is an equally important source of CO2 in soil. CO2 is generated in the. upper few metres of soil, and diffuses upward to the atmosphere and downwards to groundwater. Concentrations and diffusion rates vary with seasons, climate, type of soil, and type of'vegetation cover. CO2 partial pressures (PCO2) of up to 0.2 have been measured in soils, and an average PCO2 of 0.03 is often quoted (Miotke, 1974). To monitor CO2 emission from soils, two methods have been followed: (i) the CO.2 increase in a plastic dome covering the soil surface is measured, and (ii) the apparent diffusion constants of the respective soil profile is determined in the laboratory, together with CO2 profiles under regular field conditions. Albertsen (1977) has obtained annual averages of CO2 emissions by this second method for four sites with different vegetation types in Northern Germany: coniferous forest 31 g C/m2 year; deciduous forest 86 g C/m2 year; forest clearing 196 g C/m2 year; and agricultural land 404 g C/m2 year. The bulk CO2 emission is directed towards the atmosphere; only in spring and early autumn was a concentration gradient towards the groundwater noticed.

The total CO2 content of soil air may be estimated by assuming an average groundwater table depth of 6 m (Garrels et al., 1975), an average PCO2 of 0.003 (which is certainly on the low side), an average gas effective porosity of 0.2, and a vegetation-covered total continental area of 100 x 106 km2 (see Chapter 5, this volume, Table 5.2, total continental area minus deserts, semideserts and ice). This calculation yields 36 x 109 m3 of CO2 of 9.7 x 1012 g C. As the annual net primary productivity amounts to about 60 x 1015 g C/year (Chapter 5, this volume), the residence time of CO2 in soil is approximately one hour. This result is obtained on the assumption that the total primary production is degraded in the soil. This is obviously not the case, as a residence time of one hour is in contradiction with the slow diffusion rate through soil. Therefore, large quantities of litter must degrade under direct contact with the free atmosphere.

The annual net flux of CO2 from soil to groundwater must at least equal the difference between half the HCO3 load of the streams and the CO2 in rain: 0.23 x 1015 g C 0.018 x 1015 g C = 0.21 x 1015 g C. Garrels et al. (1975) assume that this amount of CO2 used in the weathering of rocks is derived directly from the atmosphere. As discussed above, this does not seem to be the case. The CO2 in the process of weathering is derived from 26.5 x 1015 g C of CO2 in the model by Garrels et al. (1975) produced by respiration and decay, which then diminishes to about 26.3 x 1015 g C returning directly to the atmosphere. Consequently the flux of CO2 used in weathering from the atmosphere to land also diminishes.

The fractionation of degradation CO2 between atmosphere and groundwater is 122/1; 0.77% of the CO2 produced by degradation of terrestrial organic matter is transferred to groundwater.

With respect to the interaction of water and soil, one other aspect should be noted. If, indeed, most of the CO2 in groundwater is soil-derived, then this CO2 must be much older than that furnished by rain. The turnover of humus lapses over a long period, often several hundred years. Reiners (1973) has presented a good review of detritus sizes and turnover masses. Sinter chronology by 14C is limited by the inability to define the time of humus turnover. Double analyses of wood and its sinter cover show the flowstone cover to be much older than the wood, even after adjusting for the dead CO2 derived from the country rock of the cave roof. It should be possible, by this method, to estimate the fractionation of soil CO2 between percolated water and air.

12.6 TYPES OF RUNOFF

Three types of runoff have to be discerned: overland flow, interflow, and subsurface runoff by groundwater.

Surface runoff is mostly low in dissolved matter, but carries almost all of the suspended particles of the runoff. At times of maximal runoff, the concentration of dissolved substance decreases while suspended load increases.

Interflow occurs through the very upper part of the soil cover, through cracks in the soil, and small surface voids. Interflow behaves hydrologically much like surface runoff, but has geochemical similarities to groundwater.

Groundwater passes totally through the soil and can be clearly discerned geochemically from surface runoff. Its discharge peak after thunderstorms is much delayed or cannot be identified at all. Hydrologists have developed several methods to divide river stage curves into surface and groundwater runoff. These methods, however, are not widely used, and global estimates on the ratio of surface to groundwater discharge have not, as yet, become available.

Three types of groundwater aquifers can be discerned: cleft, porous, and conduit aquifers. The first two have diffuse flow patterns and occur in igneous rocks and shales (cleft), and in sandstones and unconsolidated alluvial rocks (porous). Conduit flow occurs mostly within carbonate karst rocks. Some 80% of the continental surface is covered by sediments, of which 60% are shales, 20% carbonate rocks, 15% sandstones, and 5% evaporite deposits. The average composition of the sediment cover is slightly different, with 65% shale and 15% carbonates (Garrels and Mackenzie, 1971; Garrels et al., 1975).

12.7 GROUNDWATER IN KARST ROCKS

The area governed by conduit flow is probably less than 20% of the continental surface, as many calcareous sandstones and shales do not develop conduits. Northern Germany, for example, is underlain by calcareous glacial till, and rivers show high alkalinity and alkaline earth contents similar to karst areas, though hydrologically and morphologically no criteria of karst have been developed.

Balázs (1977a) estimated the total area with morphological karst features to be 4% of the continents (Table 12.4). The bulk of the karst occurs in Europe and Asia between 60° N and 20° N latitude, with smaller areas in North America and even less in the southern hemisphere (Figure 12.3). Carbonates occur either in the orogenic belts, or as epicontinental sediments deposited from shallow seas. These areas can form vast karst plateaus once they are uplifted, as in the case of Southern China where the largest single karst area exists. The most pronounced karstification occurs where limestone and humid climate coincide. Figure 12.4 is a map which shows the main areas of karst: the AppalachianCaribbean, the European, the South-East Asian and Pacific Karst Provinces (Balázs, 1977b).

Table 12.4 Distribution of karst areas (Balázs, 1977a. Reproduced by permission of the 7th International Speleological Congress and by permission of the author)


Area of Area of karst/1000 km2 Percentage of karst/
continent continent =100%
Continent 106 km2 orogenic epeirogenic Total orogenic epeirogenic Total

Europe
10.5
528
890
1418
5.0
8.5
13.5
Asia
43.9
808
793
1601
1.8
1.8
3.6
Africa
30.3
67
923
990
0.2
3.0
3.2
North and Central 24.2
760
100
860
3.1
0.4
3.5
America
South America
17.9
90
90
0.5
0.5
Australia and the
8.5
132
250
382
1.6
2.9
4.5
Pacific
Total without
135.3
2385
2956
5341
1.8
2.2
4.0
Antarctica

Figure 12.3 Karst areas in comparison to orogen and platform areas in western and eastern hemispheres (top figure). Latitudinal distribution of karst areas in 5° intervals both for total and individual continents (Balász, 1977a. Reproduced by permission of the 7th International Speleological Congress and by permission of the author.)

Figure 12.4 Regions of largest geoclimatic potential for karstification (Balázs, 1977b. Reproduced by permission of the 7th International Speleological Congress and by permission of the author.)

These karst areas play an important role in the bicarbonate budget. Subsurface runoff is the main mode of discharge in karst. Here, high soil PCO2 (Miotke, 1974) can be neutralized to saturation with calcite (CaCO3) or dolomite (CaMg(CO3)2 ), binding large amounts of CO2. As limestone areas are more common in the humid areas of North America, Europe, and Asia than in the southern continents, we find a very low mean alkalinity in Australian and South American rivers (see Table 12.1) where hardly any limestone exists. In fact, Europe (with a karst area of 13.5%) has the largest average total dissolved solid (TDS) concentration (182 ppm) of all the continents. Half of this load is HCO3 (95 ppm).

The complex flow patterns within a carbonate conduit aquifer are illustrated by Figure 12.5. Very often, air-filled caves intersect the path of seepage water. Most caves exchange their air quite rapidly with the atmosphere. Cave air has, therefore, a lower PCO2 than soil water percolating from above. The water loses CO2, thus increasing the PCO2 of cave air. Air currents transport the CO2 back to the surface, where cave mouths are sources of CO2 . When the water loses CO2, sinter is formed. Where there is no soil above the cave, sinter cannot precipitate, a fact which is best illustrated by 14C flowstone analyses, which show an interruption of sintering in glacial times (Franke and Geyh, 1972).

Figure 12.5 Water types encountered in a carbonate aquifer and their interconnecting flowpatterns (Harmon et al., 1972. Reproduced by permission of the British Cave Assocation and by permission of the authors.)

Table 12.5 Karst water analyses (Harmon et al., 1972. Reproduced by permission of the British Cave Research Association and by permission of the authors)


Temp. Ca2+ Mg2+ HCO3- pH Hardness* Hardness† SIc PCO2 Number
log of
samples

GRAND AVERAGES OF KARST WATER TYPES
Surface
Soil water
15.5
14
6
100
5.91
60
83
-2.61 -0.98
11
 
Vadose Zone
Dripstones
13.5
58
23
300
7.64
239
246
+0.11 -2.31
18
Vertical shafts
11.4
31.3
4
91 7.47
94
75
-0.79 -2.62
31
 
Zones of fluctuation and phreatic storage
Cave streams
10.7
43
7
144
7.38
136
119
-0.60 -2.36
23
Standing cave pools
12.9
37
25
214
7.77
194
176
-0.08 -2.60
18
Wells
19.1
63
15
237
7.51
222
195
-0.02 -2.21
29
 
Groundwater discharge
Springs
12.4
62
10
200
7.30
196
164
-0.44 -2.14
66
 
REGIONAL AVERAGE OF ALL WATER TYPES
New York
10.8
76.2
8
222
7.57
224
182
+0.08 -2.31
8
Pennsylvania
10.0
43.5
17
183
7.39
178
150
-0.45 -2.24
37
West Virginia
12.9
13
1
44 7.23
37
37
-1.59 -2.65
5
Virginia
11.8
46
23
251
7.59
211
206
-0.14 -2.36
34
Kentucky
12.3
41
6
137
7.56
130
112
-0.39 -2.52
74
Alabama
10.2
18
2
65 7.20
52
53
-1.50 -2.59
13
Florida
24.6
50
7
153
7.78
154
125
+0.12 -2.61
10
Texas
22.6
73
20
316
6.96
264
259
-0.28 -1.49
5
Mexico
24.2
145
11
456
7.00
408
374
+0.17 -1.39 
26
Yucatan
26.2
118
50
382
7.13
503
314
+0.17 -1.57
11
REGIONAL AVERAGE OF ALL SPRINGS
New York
8.9
76
7
219
7.56
219
180
+0.04 -2.32
6
Pennsylvania
9.7
42
14
163
7.27
162
134
-0.61 -2.17
19
Virginia
10.6
37
4
116
7.44
108
95
-0.58 -2.46
7
Kentucky
11.9
46
6
151
7.45
140
124
-0.42 -2.38
21
Texas
22.3
83
16
319
6.81
274
261
-0.35 -1.33
4
Mexico
23.4
202
23
555
6.84
596
455
+0.23 -1.14
6

*Hardness as ppm Ca CO3 calculated from Ca2+ + Mg2+
†Hardness as ppm CaCO3 calculated from HCO3-.

With the availability of computer programs for calculating the carbonate equilibria in fresh waters, karst waters have been investigated intensively for their CO2 pressure and mineral saturation. Variation in PCO2, alkalinity, and saturation is due to the length of contact between water and rocks, to the kind of flow encountered, to the mean temperature (and hence to latitude and altitude), and to the seasons.

Table 12.5 compares soil water (top) with water from vadose and phreatic zones (Harmon et al., 1972). The PCO2 (expressed as a logarithm) of soil water is 0.10, while the saturation index of calcite (SIc) shows strong undersaturation (negative value). Deeper water still has a PCO2 of 0.05, while the SIc has reached values near saturation. This water, upon surfacing, loses CO2 until it reaches the atmospheric PCO2 of 0.000 33 (log PCO2 = -3.48) and will precipitate calcite.

The type of flow determines several parameters of the carbonate system. Parizek et al. (1971), who investigated karst springs with diffuse and conduit flow, found lower total hardness and lower calcite saturation in the case of rapid conduit flow. The variability of hardness is much larger in the conduit compared with the diffuse flow, but CO2 pressure stayed around 0.0032 (log = -2.5) in both groups.

Average composition of karst water is a function of mean temperature (Figure 12.6, Table 12.5). CO2 pressure and hardness increase at lower latitudes. This is due to a higher microbial CO2 production in soils of warmer climates, compared to those of cooler regions. For example, PCO2 increases from the northern United States (0.0032, log = -2.5) to Mexico (0.032, log = -1.5) by a factor of 10.

Figure 12.6 (a) Relation between average temperature of karst water samples in Table 12.5 and latitude (Harmon et al., 1972). (b) Variation of hardness and PCO2 with water temperature for averages of all samples within a region (compare Table 12.5) (Harmon et al., 1972). (c) Variation of hardness and PCO2 with temperature for spring samples only, averaged within the given regions (Harmon et al., 1972. All figures reproduced by permission of the British Cave Research Association.)

Seasonal variations in PCO2 of karstic waters have been monitored by Hess (1974) in a Pennsylvanian limestone karst and by the Arbeitsgemeinschaft für niedersächsische Höhlen (unpublished data) in a dolomite and gypsum karst in Germany. Figure 12.7 gives two log PCO2 curves from the German karst, depicting a spring from dolomitic source rock with a diffuse flow (annual mean log PCO2 1.7) and a gypsum karst spring with a more or less conduit type of flow (annual mean log PCO2 2.4). The main features of seasonal CO2 variation in Germany and Pennsylvania show a minimum in CO2 pressure in March-April, due to the beginning of the vegetation period, and a maximum in late summer and autumn, when respiration and decay of organic matter prevail. The amplitude is quite large, over log values of 1.0 in the German and 0.7 in the Pennsylvanian cases. These figures reveal a pronounced seasonal effect in the CO2 liberation rate from soil. It is, however, not yet possible to calculate from these figures the net flux of inorganic carbon and its seasonal variation. It is known from experience that alkalinity and saturation levels are not only a function of the PCO2, but also of the amount of water present.

Figure 12.7 Seasonal variation of PCO2 in springs of South Harz mountains, Germany. Top curve: Diffuse flow spring from a dolomite aquifer. Bottom curve: conduit flow spring from a gypsum aquifer (drawn from own data)

12.8 DEEP GROUNDWATER

So far, we have only been concerned with groundwater of the uppermost water-filled zone. Springs discharge water as long as the reservoirs behind them are filled. This gravitational resurging water is dependent on renewal by precipitation in the source area and by the transmissibility of the rock. Low transmissibility will cause slow outflow. The water discharged by springs is mainly recent percolation water, yet a certain amount of older deep groundwater is also involved.

Deep groundwater has very long residence times, up to 1000 or even 10 000 years. The chemistry of deep groundwater (White et al., 1963) is well known in certain areas because of exploration for oil or thermal water. Estimates on the mass of water within the crust differ widely (Table 12.1).

Many waters (conate waters) have been buried together with the sediments. In tectonically active regions a constant pore-water flux is maintained by compaction. By lowering fresh mud, with 80% porosity, 2 km down into the crust such compaction will occur, and only 10% of the porosity remains (Wunderlich, 1966). Under such pressures a variety of mineral exchange reactions occur, especially with carbonates. The resurgence of these pore waters is one of the factors which stimulate the circulation of deep groundwater.

12.9 RIVERS

In most cases, groundwater gives off CO2 at its reappearance. If the water stems from limestones, and is near saturation with respect to calcite, calcium carbonate precipitates downstream. Because of the influence of freshwater biota, most rivers never attain equilibrium of PCO2 with the atmosphere. A characteristic log value for PCO2 in rivers is 2.5. Figures 12.8 and 12.9 depict seasonal variations in PCO2 and alkalinity for the rivers Elbe (23-year mean) and the Amazon (for the year 1975). In the Elbe River, CO2 pressure and alkalinity decrease with the beginning of the phytoplankton bloom in early spring, with an annual minimum in April. Heterotrophic processes cause the CO2 pressure to rise quickly to a maximum in July, while alkalinity only slowly regains its former concentrations. A small minimum, caused by a late-summer plankton bloom, is encountered around August and September. The long-term mean of the log PCO2 in the Elbe River is 2.4.

The seasonal PCO2 and alkalinity variation in the Amazon is less regular, though large changes occur during the year. Alkalinity is high during times of spring floods, when more country rocks are exposed to erosion by water. PCO2 is low during periods of high waters, when nutrients are available to phytoplankton, while larger PCO2 values are encountered during periods of lower waters and the ceasing of photosynthesis.

Figure 12.8 Twenty-three-year mean of seasonal variation in PCO2 and alkalinity of the Elbe River at Hamburg (calculated from hydrochemical river records by the Hamburger Wasserwerke)

Figure 12.9 Seasonal variation of PCO2 and alkalinity of Amazon River at Obidos 1975 (drawn from data M.Sc. thesis M. L. Nossar Simöes de Dalgo, personal communication Dra. A. de Luca Rebello, Rio de Janeiro) 

Since the CO2 pressure of rivers is approximately 10 times larger than that of the atmosphere, there must be a considerable flux of CO2 from stream surfaces into the air. This flux could be estimated according to the laws of diffusion, if we knew the total surface area of streams and their average boundary layer thickness. However, neither figure is available at the moment.

Groundwater also discharges directly into the oceans. This is especially true for limestone areas bordering the sea; half of the drainage of the Yugoslavian Dinarides, for instance, is said to discharge below sea-level (personal communication, Dr. R. Gospodaric, Postojna). Water from caves, which developed during glacial times at low sea-level, continues to be transmitted directly to the sea. This water does not come into contact with air again, and, since it stems from carbonate rocks, it carries a substantially greater load of carbon than the global river mean. The size of this flux has not yet been estimated (see Table 12.2).

With respect to the various parts of the water cycle and the carbon flow interconnected with it, it becomes apparent that the gross flux of carbon from soil to fresh water is probably much larger than that of river loads discharged to the ocean (Table 12.1 b). How much more carbon that is remains obscure; it may be 50 or even close to 100 per cent of the final discharge.

An attempt is made in Figure 12.10 to show schematically the various fluxes found in the freshwater part of the carbon cycle. Instead of reservoir sizes and flux figures, however, only the PCO2 values and their range in the different compartments are given.

12.10 REVISION OF LIVINGSTONE'S (1963) RIVER LOADS

In 1963 Livingstone published the paper, `Chemical composition of rivers and lakes', which is still the fundamental work for estimating the erosion rates of the continents. By recalculating his world means, several errors, however, became apparent. Both the weighted mean of river salinities of the individual continents and of the world, as well as the HCO3 figure for world runoff are inaccurate.

The weighted mean is calculated by multiplying the volume of discharge with the salinity, adding these figures for continents (or total earth), and dividing this sum by the total discharge for all continents (or total earth):

 
(river discharge x salinity)
M (continents)=
river discharge

 Livingstone obtained a mean world river salinity of 120 ppm. By using his continental means, this figure should read 122.4 ppm. By using recalculated continent values, the salinity is 116.4 ppm and the HCO3 content becomes 59.8 ppm (instead of 58.4 ppm).

The runoff values quoted by Livingstone have been thoroughly revised by Baumgartner and Reichel (1975) (Table 12.1b). Using their continent discharge values, in combination with the revised continent salinities of Livingstone, we obtain a mean total salinity of only 111.06 ppm, i.e. 9 ppm lower (7.5%) than the original value of Livingstone. This is especially due to the present larger discharges of the Amazon and Orinoco in South America. On the other hand, the total discharged volume has increased (Livingstone, 32.4 x 103 km3 /year; Baumgartner andReichel, 37.8 x 103 km3 /year),which gives a total dissolved load of 3.77 x 1015 g/year (revised Livingstone value) and 4.20 x 1015 g/year (new value) respectively. For HCO3 the Livingstone value is 1.94 x 1015 g/year (= 0.381 x 1015 HCO3C g/year) and the new estimate becomes 2.26 x 1015 g/year (= 0.444 x 1015 HCO3C/year). There is new chemical data available for some of the larger rivers as well. Livingstone used a mean of 26.4 mg of HCO3 for the Amazon River. Data from the M.Sc. thesis of Miss Maria Lucia Nossar Simões de Dalgo (courtesy Dra. Angela de Luca Rebello, Rio de Janeiro, Brazil) indicates an annual weighted mean of 30.4 mg/l. As the Amazon supplies 15.1% of the world runoff, the new mean yields.

30.4

2.26x1015 +2.26 x1015 x0.151x


1   = 2.31x1015  g

26.4

(or 0.454 x 1015 g HCO3-C/year). This newer value is 19% higher than the revised Livingstone calculations.

Figure 12.10 -log PCO2 of different freshwater compartments. The CO2 pressure is expressed as its negative decadic logarithm: one unit less means a ten times larger pressure

12.11 POSSIBLE SINKS OF CARBON IN THE FRESHWATER CYCLE 

Industry accounts for approximately 5 x 1015 g of fossil carbon emitted per year into the atmosphere, and possibly the same amount is liberated by land management practices. Of this amount, only 2.4 x 1015 g stay in the atmosphere while 1.4 x 1015 g are believed to dissolve in the oceans. The rest enters the carbon cycle at points not yet known. Checking the freshwater cycle for such unnoticed `leaks', we have to keep in mind that the amount of carbon to be looked for surpasses the total flux rates discussed here. We can, therefore, expect only minor sinks in this part of the cycle. Nevertheless, these might exist at several points.

In the most heavily industrialized zones, SO2 emission has decreased the pH of rain below the value of 5.6, which is the normal value, due to the presence of CO2 and carbonic acid. At a lower pH, the solubility of CO2 increases (one pH unit causes a 100-fold solubility increase; Cole, 1975). CO2 could, therefore, be lost from the atmosphere faster in polluted areas than in non-polluted regions.

The atmospheric CO2 pressure has risen  from a preindustrial value of 0.000 29 to 0.000 33 today (see Chapter 3, this volume). The rise in CO2 pressure will certainly increase carbonate dissolution on the continents, though the increase due to carbonate equilibria is only very small. Nevertheless, there is some evidence of an increase in alkalinity and in CO2 pressure as in the case of the Elbe River (see Chapter 13, this volume), and other rivers can be checked for similarly rising trends. These increases, however, may also be induced by CO2 released from organic pollution.

In mature woods, CO2 resulting from respiration in soil is quickly recycled by plant uptake, thus keeping the bottom air low in CO2 . If a forest is cut, this balance is upset in two ways: (i) the total CO2 production increases (see examples in Section 12.4), and (ii) the CO2 content of bottom air rises due to lack of large plants to recycle the emitted CO2. Both processes should result in an enlarged CO2 diffusion to groundwater and hence to an enhanced weathering of rocks.

Also, groundwater usage causes a disturbance in the carbon cycle. Most pumping of groundwater involves the liberation of CO2. On the other hand, younger groundwater, with possibly larger amounts of CO2, is diverted from the upper levels of the subsurface water body into deeper parts, possibly carrying more CO2 downwards than the amount of CO2 liberated by the pumping of older ground water.

The eutrophication of fresh water is another process by which we can permanently fix CO2. Increasing amounts of organic matter are buried in lakes, marshes, and coastal seas.

Probably none of these fluxes will be larger than 1014 g C. Nevertheless, our attempt to investigate them shows how ignorant we are about the flow patterns and flux rates of carbon in connection with the freshwater cycle.

REFERENCES

Albertsen, M. (1977) Labor- and Felduntersuchungen zum Gasaustausch zwischen Grundwasser and Atmosphäre über natürlichen and verunreinigten Grundwässern. Ph.D. thesis, 1-145, Christen-Albrechts Universität Kiel.

Balázs, D. (1977a) The geographical distribution of karst areas. Proc. 7th Intern. Congr. Speleology, Sheffield 1977, 13-15, Sheffield.

Balázs, D. (1977b.) The optimal geo-climatic provinces of karstification. Proc. 7th Intern. Congr. Speleology, Sheffield 1977, 15-17, Sheffield.

Baumgartner, A. and Reichel, E. (1975) The World Water Balance, 1-179. R. Oldenbourg-Verlag, München, Wien.

Bögli, A. (1975) Solution of calcium carbonate and the formation of karren. Cave Geology 1(1), 1-28 (trans. of the original paper, 1960).

Cole, G. A. (1975) Textbook of Limnology, 1-283. The C. V. Mosby Company, St. Louis.

Flohn, H. (1973) Der Wasserhaushalt der Erde, Schwankungen and Eingriffe. Naturwissenschaften 60(7), 340-348.

Franke, H. W. and Geyh, M. A. (1972). Tropfsteinwachstum and Datierung. Mitt. Verb. Dt. Höhlen-u. Karstforscher 18(3), 59-60.

Garrels, R. M. and Mackenzie, F. T. (1971) Evolution of Sedimentary Rocks, 1-397. W. W. Norton and Company Inc., New York.

Garrels, R. M., Mackenzie, F. T. and Hunt, C. (1975) Chemical Cycles and the Global Environment, Assessing Human Influences, 3rd ed., 1-206. W. Kaufmann, Inc., Los Altos, Calif.

Harmon, R. S., Hess, J. W., Jacobson, R. W., Shuster, E. T., Haygood, C., and White, W. B. (1972) Chemistry of carbonate denudation in North America. Trans. Cave Res. Group of Great Britain 14(2), 96-103.

Hess, J. (1974) Hydrochemical investigations of the Central Kentucky Karst aquifer system. Ph.D. thesis, 1-219, Pennsylania State University.

Kempe, S. (1977) Hydrographie, Warven-Chronologie and organische Geochemie des Van Sees, Ost-Türkei. Mitt. Geol.-Paläont. Inst. Univ. Hamburg 47, 123-208.

Livingstone, D. A. (1963) Chemical composition of rivers and lakes. U.S. Geol. Surv. Prof. Pap. 440G, 1-64.

Lvovich, M. I. (1970) Publ. Assoc. Intern. Hydro. Scient. 93, 401. (Quoted in Flohn, H., 1973.)

Miotke, F. D. (1968) Karstmorphologische Studien in der glazial überformten Höhenstufe der `Picos de Europa' Nord-Spanien. Jb. Geogr. Ges. Hannover, Arb. Geograph. Inst. Techn. Univ. Hannover, Sonderheft 4, 1-161, Selbstverlag Geograph. Ges. Hannover, Hannover.

Miotke, F. D. (1974) Carbon dioxide and the soil atmosphere. Abh. Karst- u. Höhlenkunde, R. A 9, 1-49.

Parizek, R. R., White, W. B., and Langmuir, D. (1971) Hydrogeology and geochemistry of folded and faulted carbonate rocks of the Central Appalachian type and related land use problems. Min. Conser. Ser., Earth Miner. Sci. Exp. Station, Circular 82, 1-182, Pennsylvania State University.

Reiners, W. A. (1973) Terrestrial detritus and the carbon cycle. In: Woodwell, G. M. and Pecan, E. V. (eds), Carbon in the Biosphere. AEC Symp. Ser. 30, 303-327, NTIS U.S. Dept. of Commerce, Springfield, Virginia.

White, D. E., Hem, J. D., and Waring, G. A. (1963) Chemical composition of subsurface waters. U.S. Geol. Surv. Prof. Pap. 440F, 1-67.

Wunderlich, H. G. (1966) Wesen and Ursachen der Gebirgsbildung, 1-367 Hochschultaschenbücherverlag, Mannheim.

Back to Table of Contents |

 

The electronic version of this publication has been prepared at
the M S Swaminathan Research Foundation, Chennai, India.