Introduction

Chemical weathering of silicate minerals consumes atmospheric CO2 and stoichiometrically produces dissolved silica (DSi, referred to silicic acid, H4SiO4) and HCO3 1,2. These solutes are transported to the coastal zone and lakes by rivers where the carbon can be synthesized into organic matter by phytoplankton and the silica can be used to produce diatom frustules. Ultimately, some of these organisms sink to the bottom sediments where carbon (C) and silicon (Si) can be sequestered. Thus, silicate weathering has a net-sink effect on atmospheric CO21,3. However, this is not the case for carbonate mineral weathering by carbonic acid, because CO2 is re-released to the atmosphere when carbonate mineral is re-precipitated on geological time scales (>10 to 100 ka). Besides, strong acids (i.e. sulfuric and nitric acid) from anthropogenic sources have enhanced carbonate weathering and CO2 evasion4,5,6. Therefore, C and Si biogeochemical cycles are important in regulating atmospheric CO2 concentrations and hence climate change.

Since the industrial revolution, anthropogenic perturbation such as nutrient enrichment and construction of dams have altered chemical weathering and C and Si fluxes through the river systems7,8, and the related mechanisms are suggested to further research9. C sequestration coupled with the Si cycle in agricultural ecosystems has recently been summarized10, however, coupling of C and Si geochemical cycles have not yet been highlighted in freshwaters. This may be because (1) riverine inorganic C and dissolved Si (DSi) concentrations vary greatly on a global scale11 and there is not a conspicuous relationship between chemical species such as calcium and bicarbonate12, (2) C and Si are traditionally considered to be less important than nitrogen (N) and phosphorus (P) in limiting the growth of photosynthetic organism.

Nutrient biogeochemical cycles in aquatic ecosystem linked through the metabolic activity of living organisms13,14. For example, diatoms require DSi to build their frustules15 and C for their organic molecules. To test the extent of coupling between C and Si in inland waters, we compiled data from global rivers, and lakes from the United States of America and the English Lake District. The specific aims of this study were to (1) characterize the relationship between dissolved inorganic C and Si in rivers and (2) determine how phytoplankton couples the C and Si biogeochemical cycles in lakes.

Results

Variation of carbon and silicon in rivers

Major ion concentrations ranged over three orders of magnitude in the 175 global rivers (Supplementary Tables 4). The average concentrations of DSi and HCO3 were 165 and 1296 μmol l−1, respectively. Sodium (Na)-normalized DSi and HCO3 concentrations showed a significant positive correlation (Fig. 1a–c), however, this was not the case for their absolute concentrations in the 175 global rivers (Supplementary Table 7). As expected, Na-normalized DSi and HCO3 ratios increased from watersheds dominated by silicate weathering to those dominated by carbonate weathering, while DSi:HCO3 ratios showed the reverse tendency, suggesting that they can be used to discern the source of mineral weathering in global rivers (Figs 1a–c and 2a,b). DSi:HCO3 ratios were positively correlated with 87Sr/86Sr (Fig. 2a,b) and negatively correlated with total dissolved solid (TDS) (Fig. 2c). DSi and HCO3 in the Changjiang catchment showed a similar geochemical behaviour to those in the 175 global rivers (Figs 1a,b and 2a,b).

Figure 1
figure 1

Relationships between Na-normalized DSi and HCO3.

(a) From the 175 rivers; (b) From the Changjiang River; (c) From the GLORICH database (Supplementary Methods); (d) From the lakes in United States of America (AL) and in the English Lake District (EL). The data from the American Lakes are the mean of each dataset (Supplementary Table 2), whereas the data from the English Lake District are the mean for each lake every five years between 1991 and 2010 (Supplementary Methods).

Figure 2
figure 2

Relationship between DSi:HCO3 and 87Sr/86Sr, TDS (total dissolved solid), and chlorophyll a, respectively.

Data in (a,c) are from the 175 rivers; data in (b) are from the Changjiang River, data in (d) are from the two lake regions (referred to Fig. 1).

Variation of carbon and silicon in lakes

Average DSi and HCO3 concentrations in the American lakes were of a similar magnitude to those in the 175 global rivers, and an order of magnitude greater than the concentrations in the English Lake District (Supplementary Tables 4–6) because their catchments were small and mainly comprised slowly weathering igneous and silicate rocks. Average contribution of HCO3 to dissolved inorganic carbon (DIC) in the English Lake District (78.8%) was less than that in the American lakes (91%) because of the lower HCO3 concentrations. For both lake regions, both absolute and Na-normalized DSi and HCO3 concentrations were correlated significantly, but the sign of the correlation (i.e. positive or negative) was different (Fig. 1d and Supplementary Tables 8 and 9). The DSi:HCO3 ratio was negatively correlated to chlorophyll a, and this relationship was more significant in the English Lake District than in the American lakes (Fig. 2d), explaining the different sign of the correlation. In the English Lake District, the relationship between the DSi:DIC and the algal Si:C ratio varied throughout the year and took the form of a bell-shaped curve (Fig. 3).

Figure 3
figure 3

Relationship between the DSi:DIC and the algal Si:C ratios in the English Lake District.

The data represent the means from January, April, July and October for 20 lakes every 5 years (1991, 1995, 2000, 2005, and 2010). The cross represents the seasonal means of the 20 lakes over 5 years.

Discussion

Riverine DSi and HCO3 mainly originate from chemical weathering. The variation in their concentration results from heterogeneity in factors such as land cover, runoff, and evaporation in different rivers at the global scale11, but Na-normalized concentrations can eliminate the influences of these factors as Na is mostly originating from weathering processes16,17. This is the reason that the Na-normalized DSi and HCO3 ratios are linked and can reflect their origins. The similar tracing function of the DSi:HCO3 ratio is supported by its significant relation to 87Sr/86Sr, which is known to identify the nature of watershed rock and local weathering processes12,16. 87Sr/86Sr is a better tracer of the end member of carbonate weathering than that of silicate weathering (Fig. 2a), as this ratio will depend on the age of the geological province18,19. The DSi:HCO3 ratio is influenced by the contribution of silicate weathering, which in turn is strongly controlled by runoff12,20,21. Data from the 175 global rivers are consistent with this: the DSi:HCO3 ratios were positively correlated with runoff (r = 0.373, n = 167, p < 0.001). The DSi:HCO3 ratio was also negatively correlated with TDS (Fig. 2c) since the rate of silicate chemical weathering is a function of runoff and TDS12.

Lacustrine DSi:HCO3 ratios do not reflect the characteristics of chemical weathering determined by the catchment because algal activity differentially takes up Si vs carbon and the ratio therefore decreases with an increase in phytoplankton biomass (i.e. chlorophyll a). Diatoms can use both CO2 and HCO3 as an inorganic C source22. Thus, DSi:DIC ratio responds to algal activity more significantly than DSi:HCO3 ratio, although HCO3 is a dominant species in DIC. The average Si:C atomic ratio of freshwater diatoms, 0.7923, is usually larger than that of inflowing rivers. Therefore, with an increase in diatom biomass, residual solutes show decreasing Si:C atomic ratios24. This conclusion is also strengthened by the data from the English Lake District where the inflowing streams had average DSi and DIC concentrations of 43.8 and 370 μmol l−1 respectively (i.e. DSi:DIC = 0.12)25,26, whereas the average DSi and DIC concentrations in the lakes themselves were 22.3 and 243.9 μmol l−1 respectively (i.e. DSi:DIC = 0.09; Supplementary Tables 6).

The long-term data from the English Lake District (1991–2010) provide further insights into the interactions between the DSi:DIC and the algal Si:C ratios (Figs 3 and 4). In winter, low temperature and light limits algal growth. In January, the algal assemblages have their lowest biovolume, of which diatoms represent 37%, and have an Si:C ratio similar to the water DSi:DIC ratio. In spring, diatoms bloom and dominate the phytoplankton (49% of total biovolume), resulting in an increase of the algal Si:C ratio and a decrease in the water DSi:DIC ratio. In summer, Si depletion, stratification and nutrient shifts favour non-diatom phytoplankton27. Although diatoms decrease, the total algal biovolume is higher than in spring and the algal Si:C and DSi:DIC ratios were at a minimum. In the autumn, continued inputs of Si from the catchment and a weakening of stratification, start to permit diatom growth again. This results in an increase in the algal Si:C ratio, bringing them into equilibrium with the water DSi:DIC ratio.

Figure 4
figure 4

Relationships between selected variables of lakes in the English Lake District.

The data represent the average and standard deviation of Jan, Apr, Jul and Oct for 20 lakes over 5 years. The standard deviation was calculated from the average of each year.

Overall, this bell-shaped relation through the year (Fig. 3) has an average DSi:DIC ratio of 0.14 that is approximately equal to the average algal Si:C ratio of 0.15, values very close to the well-known Redfield ratio, 106C:16N:16Si:1P28. This bell-shaped curve suggests that there is a negative feedback regulation on Si and C stoichiometry by lacustrine phytoplankton. When the algal Si:C value deviates from the bioavailable Si:C ratio, the relative contribution from diatoms changes to adjust for the deviation. This may bring the system into an equilibrium state, which means that the Si:C atomic ratio of the algal assemblage is similar to that in the water column (Fig. 3). It is likely that this stoichiometric homeostasis may apply to other nutrients in lacustrine system, such as N and P. Similar stoichiometric homeostasis of N and P ([NO3]:[PO43−]≈15:1) in the oceans is also thought to be regulated by algae29.

The coupling of Si and C in rivers reflects their origins and controls by chemical weathering, whereas in lakes this signature is modified by algal activity. Many human-induced environmental problems, such as climate change and eutrophication, are the result of disruptions of natural nutrient stoichiometry14. Our study supports this. Once lacustrine systems deviate from the equilibrium state, they are prone to algal blooms. Furthermore, our homeostatic theory derived from the English Lake District also suggests that in the equilibrium state, C fixation by algae in lakes is quantitative and depends on the DSi concentration, and this hypothesis maybe requires to be tested more widely with other data. These findings provide some new insights into global C and Si biogeochemical cycles.

Methods

Data from global rivers were collected from the GEMS-GLORI world river discharge database30, from which 175 rivers with SiO2 (i.e. DSi) concentrations were selected (Supplementary Table 1). Data from small to medium-sized river catchments were collected from GLORICH database, and the samples having DSi, HCO3 and Na concentrations were selected. Data from American lakes were derived from the 2007 National Lakes Assessment, and divided into 19 datasets according to chlorophyll a concentration (Supplementary Table 2). Surveys of the lakes in the English Lake District were conducted seasonally in 1991, 1995, 2000, 2005 and 2010, and chemical variables and phytoplankton species composition were investigated. Algae were identified and enumerated using a microscope31. The biovolume of each species was calculated geometrically32, and the C and Si contents per cell were calculated on the basis of the relationship between cell volume and cellular C and Si contents33,34, respectively. Concentrations of CO2, HCO3 and CO32− were calculated based on alkalinity, pH and temperature, with equilibrium constants corrected for temperature and ionic strength35. Further details about the Methods are found in the Supplementary Methods.

Additional Information

How to cite this article: Wang, B. et al. Coupling of carbon and silicon geochemical cycles in rivers and lakes. Sci. Rep. 6, 35832; doi: 10.1038/srep35832 (2016).