A lake classification concept for a more accurate global estimate of the dissolved inorganic carbon export from terrestrial ecosystems to inland waters

The magnitude of lateral dissolved inorganic carbon (DIC) export from terrestrial ecosystems to inland waters strongly influences the estimate of the global terrestrial carbon dioxide (CO2) sink. At present, no reliable number of this export is available, and the few studies estimating the lateral DIC export assume that all lakes on Earth function similarly. However, lakes can function along a continuum from passive carbon transporters (passive open channels) to highly active carbon transformers with efficient in-lake CO2 production and loss. We developed and applied a conceptual model to demonstrate how the assumed function of lakes in carbon cycling can affect calculations of the global lateral DIC export from terrestrial ecosystems to inland waters. Using global data on in-lake CO2 production by mineralization as well as CO2 loss by emission, primary production, and carbonate precipitation in lakes, we estimated that the global lateral DIC export can lie within the range of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {0.70}_{-0.31}^{+0.27} $$\end{document}0.70−0.31+0.27 to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {1.52}_{-0.90}^{+1.09} $$\end{document}1.52−0.90+1.09 Pg C yr−1 depending on the assumed function of lakes. Thus, the considered lake function has a large effect on the calculated lateral DIC export from terrestrial ecosystems to inland waters. We conclude that more robust estimates of CO2 sinks and sources will require the classification of lakes into their predominant function. This functional lake classification concept becomes particularly important for the estimation of future CO2 sinks and sources, since in-lake carbon transformation is predicted to be altered with climate change.

Integrating inland waters into Earth system models Earth system models (ESMs) simulate the interactions between global climate and biogeochemical cycles based on the physical, chemical, and biological properties of the three main components of the Earth system: land, atmosphere, and ocean. Connecting atmospheric transport, ocean circulation, and terrestrial biosphere models (TBMs; land) allows simulation of carbon stores and fluxes at the global scale (Falkowski et al. 2000;IPCC 2013). TBMs simulate biogeochemical and physical processes of terrestrial ecosystems, including up to 25 key processes (Fisher et al. 2014). While global terrestrial gross primary production (GPP land ) can be measured using satellite remote sensing data, autotrophic and heterotrophic respiration of land ecosystems (R a and R h , respectively) need to be simulated in current TBMs. In TBMs, the biomass production on land is quantified as net primary production (NPP land ; NPP land = GPP land − R a ), while the amount of carbon stored in or released from terrestrial ecosystems, terrestrial net ecosystem production (NEP land ), is obtained by subtracting the total terrestrial ecosystem respiration (R a + R h ) from GPP land (Fisher et al. 2014). Inland waters connect the Earth system components land and ocean (Cole et al. 2007;Drake et al. 2017;Tranvik et al. 2009). However, current ESMs do not simulate carbon fluxes in inland waters (Bauer et al. 2013). Instead, inland waters were for a long time regarded as Bpassive pipes^between land and ocean. Recently, aquatic carbon fluxes were integrated into TBMs on a regional scale (Langerwisch et al. 2016).
In recent years, the view of inland waters being passive carbon transporters between land and ocean has changed. Instead, lakes have been identified as important regulators of carbon processing along the land to ocean aquatic continuum (LOAC) (Battin et al. 2009;Biddanda 2017;Cole et al. 2007;Tranvik et al. 2009). The functioning of lakes in carbon transport and transformation controls both lateral (i.e., hydrologic transport) and vertical (i.e., emission and burial) carbon fluxes along the LOAC and thus influences the global carbon balance (Battin et al. 2009;Biddanda 2017;Cole et al. 2007;Tranvik et al. 2009). These findings have important implications for the calculation of the NEP land in TBMs. When the proportion of R h from terrestrial biomass that leaves terrestrial ecosystems through lateral hydrologic export to streams and lakes (Oquist et al. 2014) is not accounted for when simulating R h , NEP land is overestimated. The recognition of the importance of lateral aquatic carbon transport for continental carbon budgets has led to the realization that terrestrial ecosystems are less efficient in sequestering carbon than previously assumed (Butman et al. 2016;Ciais et al. 2008). A recent study showed that the NEP land of the conterminous USA might have been overestimated by more than 25%, as lateral aquatic carbon fluxes were not accounted for in present TBMs (Butman et al. 2016). Thus, realistic estimates of the terrestrial carbon sink/source require accurate quantification of the lateral carbon export from terrestrial ecosystems to inland waters. These estimates are currently not available and are unrealistic to measure over large geographic regions. Thus, the inclusion of lateral inland water carbon fluxes into global ESMs remains difficult, but essential if we seek to reconcile global carbon budgets (Battin et al. 2009;Butman et al. 2016;Weyhenmeyer et al. 2015).
Lake functioning along the aquatic continuum Rivers, floodplains, and lakes control carbon transport as well as transformation along the LOAC (Cole et al. 2007;Raymond et al. 2013;Tranvik et al. 2009). The integration of lakes into global carbon dioxide (CO 2 ) budgets is difficult, since lakes function differently depending on their characteristics and location ). Nutrient conditions, hydrology, catchment characteristics, lake morphology, and regional climate are important factors determining the functioning of lakes in the global carbon cycle (Lewis Jr. 2011;Tranvik et al. 2009;Weyhenmeyer et al. 2015). The role of lakes in dissolved inorganic carbon (DIC) transport and transformation ( Fig. 1) depends on the characteristics of each lake. In-lake CO 2 consumption and production might, for example, be the most important drivers of lake carbon dynamics in warm eutrophic lakes (Almeida et al. 2016), while lateral CO 2 transport can be highly significant in boreal lakes (Weyhenmeyer et al. 2015). Although decomposition rates of organic carbon are higher in lakes with short water residence times (Catalán et al. 2016), short water residence times generally result in lower in-lake CO 2 production and consumption, if the majority of carbon is transported downstream before being processed .
At present, different assumptions regarding the functioning of lakes in the global carbon cycle are found in the literature. Many global carbon estimates, including those currently used for policy decisions (IPCC 2013), assume that in lakes, CO 2 is efficiently produced by mineralization of terrestrial dissolved organic carbon (DOC) (Aufdenkampe et al. 2011;Battin et al. 2009;Cole et al. 2007;IPCC 2013;Raymond et al. 2013;Tranvik et al. 2009). However, assuming that all lakes function similarly can be problematic since the transformation of DOC to CO 2 in many lakes is less efficient than previously thought Stets et al. 2009;Weyhenmeyer et al. 2015). In numerous lakes, a large proportion of the emitted CO 2 originates from terrestrial ecosystem respiration (R h ) and has been transported to inland waters via discharge (Weyhenmeyer et al. 2015). Assuming this CO 2 to be produced by in-lake DOC mineralization results in an overestimation of NEP land in TBMs and ESMs. Thus, assumptions about lake functioning have a large impact on the calculated lateral DIC export from terrestrial ecosystems to inland waters.
The aim of this study was to develop a conceptual model to quantify variations in the estimate of the global lateral DIC export from terrestrial ecosystems to inland waters, depending on the assumed predominant function of lakes.
A lake classification concept and effects on the global terrestrial DIC export To integrate lakes into ESMs, we classified lakes into three functional categories: (1) lakes as active carbon transformers, (2) lakes as intermediate active carbon transformers, and (3) lakes as passive open channels, where classes 1 and 3 represent the ends of a continuum of possible lake functioning depending on lake characteristics (Fig. 2).
To demonstrate the effect of different assumptions about lake functioning on calculated estimates of DIC export from terrestrial ecosystems to inland waters (streams, rivers, lakes, and reservoirs), we established the following mass-balance equation by accounting for all main DIC fluxes along the LOAC (Fig. 1): where DIC export is the DIC exported from terrestrial ecosystems to inland waters, DIC ocean is the DIC entering the oceans via lateral surface and groundwater transport, CO 2 _emission lake is the net CO 2 emission from lakes and reservoirs (hereafter lake), CO 2 _emission lotic is the net CO 2 emission from streams and rivers (lotic systems), GPP lake is the CO 2 consumption by lake primary production, GPP lotic is the CO 2 consumption by primary production in streams and rivers, CCP lake is the in-lake calcium carbonate precipitation, CCP lotic is the calcium carbonate precipitation in streams and rivers, MIN lake is the amount of CO 2 produced by lake mineralization, and MIN lotic is the amount of CO 2 produced by stream and river mineralization. Organic carbon sedimentation and burial is not included in our conceptual model, since we restricted the analysis to DIC fluxes along the LOAC. Our model ( Fig. 2) is meant to provide a conceptual framework that can be applied at different spatial scales with any of the most comprehensive available estimates describing the fluxes stated in Eq. 1.
In lakes that predominantly function as active carbon transformers (e.g., warm eutrophic lake ecosystems (Almeida et al. 2016)), GPP lake is substantial for the lake carbon budget. The terrestrial DIC export for landscapes in which lakes predominantly function as active carbon transformers can be estimated using Eq. 1.
When assuming lakes to function predominantly as intermediate active carbon transformers, CO 2 emissions from lakes would mainly be sustained by in-lake DOC mineralization, and GPP lake as well as CCP lake would be close to zero. In that case, DIC export can be estimated as: This intermediate lake function reflects that some lakes, often referred to as heterotrophic lakes, have high bacterial and photochemical CO 2 production, and DOC-derived CO 2 is the main source of lake CO 2 emissions. These lakes have little phytoplankton and zooplankton production, which is typical for some boreal lakes (Jonsson et al. 2001) and probably even for nutrient-poor and deep lakes. Assuming CO 2 _emission lake and MIN lake approach zero in Eq. 2 does not indicate that these lakes have no in-lake carbon transformation, but that the CO 2 production in these lakes is mainly sustained by mineralization of allochthonous organic carbon. Thus, for this lake type, CO 2 _emission lake cannot be included as flux term when calculating the DIC export from terrestrial ecosystems to inland waters. This lake function reflects the way global lakes are currently accounted for in ESMs, where it is assumed that CO 2 in inland waters originates mainly from in-lake mineralization of allochthonous organic carbon (Aufdenkampe et al. 2011;Battin et al. 2009;Cole et al. 2007;IPCC 2013;Raymond et al. 2013;Tranvik et al. 2009).
When instead assuming that lakes predominantly function as passive open channels (e.g., small boreal and temperate lakes (Jonsson et al. 2003;Stets et al. 2009)), GPP lake , CCP lake , and MIN lake are minor, and CO 2 emissions are mainly sustained by hydrologic DIC inflow to lakes that is derived from terrestrial inorganic carbon export. In this case, we assumed GPP lake , CCP lake , and MIN lake to approach zero in Eq. 1. Thus, with lakes as passive open channels, the DIC export can be estimated as: Sensitivity of the global terrestrial DIC export to lake functioning To demonstrate the sensitivity of the global DIC export to the assumed functioning of global lakes using Eqs. 1, 2, and 3, we collected published data on global DIC ocean , CO 2 _emission lake , CO 2 _emission lotic , GPP lake , MIN lake , and MIN lotic from the literature (Table 1). For our calculations, we chose the most recent estimate of each respective flux, since these were the most accurate available flux estimates on the global scale. While the estimate for GPP lake from Lewis Jr. (2011) used in Eq. 1 already entails a scaling of lake GPP based on the latitudinal distribution of lakes and prevailing nutrient conditions, Raymond et al. (2013) simulated CO 2 _emission lake from non-tropical lakes using DOC and lake area and used a median value for tropical lakes in their estimate. Thus, the Raymond et al. (2013) model does not account for CO 2 emissions derived from lateral DIC inputs to lakes; hence, lakes that function as passive open channels are excluded. Consequently, scaling global CO 2 _emission lake and MIN lake according to the assumed predominant lake functions (Eqs. 1, 2, and 3) is reasonable.
Our sensitivity analysis illustrates the likely range of the calculated DIC export from terrestrial ecosystems to inland waters under the assumption of different lake functions. When assuming that all lakes on Earth function as active carbon transformers, DIC export according to Eq. 1 equaled 1:52 þ1:09 −0:90 P g C y r − 1 . DIC e x p o r t b e c a m e s m a l l e s t , i.e., 0:70 þ0:27 −0:31 Pg C yr −1 , when we considered lakes as intermediate active carbon transformers (Eq. 2). When we considered lakes as passive open channels (Eq. 3), the DIC export turned to 1:02 þ0:79 −0:57 Pg C yr −1 . Thus, we found that calculations of DIC export can vary between 0:70 þ0:27 −0:31 and 1:52 þ1:09 −0:90 Pg C yr −1 , depending on the assumed predominant function of lakes (Fig. 3). These numbers for the lateral DIC export are about 25 and 50% of the total carbon export from terrestrial ecosystems to inland waters estimated by previous studies (Aufdenkampe et al. 2011;Battin et al. 2009;Tranvik et al. 2009). We suggest that the variability in calculated lateral DIC fluxes from terrestrial ecosystems to inland waters has a strong influence on estimates of the terrestrial CO 2 sink and might explain a share of the residual terrestrial CO 2 sink of approximately 2 Pg C  (Houghton 2003;Nakayama 2017;Schimel 1995). It should be noted that DIC export comprises DIC originating from rock weathering as well as soil-derived CO 2 .

Lakes as active transformers
Only about 70% of DIC export might be soil derived (i.e., of atmospheric origin), while the other 30% is derived from rock weathering and is therefore part of the slow carbon cycle and has no atmospheric origin (Ciais et al. 2008;Einsele et al. 2001).
In our analysis, we used the flux estimate of Lauerwald et al. (2015) for CO 2 _emission lotic , since this study simulates CO 2 emissions from running waters at a much higher resolution than earlier approaches. To our knowledge, no global estimates for GPP lotic and CCP lotic exist. We set GPP lotic and CCP lotic to zero when testing the sensitivity of DIC export to lake functioning, since the fluxes in streams and rivers were kept constant for all three functional lake classes and thus did not influence the result of the sensitivity analysis. Except for the estimates for GPP lake and MIN lake (Pace and Prairie 2005) as well as MIN lotic (Caraco and Cole 1999), which were indicated in moles C yr −1 , the values were directly used from the cited publications. The fluxes given in moles C yr −1 were converted into fluxes in g C yr −1 by multiplying the given values by the atomic mass of carbon. For flux estimates reported as a range in the original study, we used the mid-range for our calculations. The uncertainties presented here indicate the range of the estimates reported in the respective studies ( Table 1). The uncertainties of the respective flux estimates used in Eqs. 1, 2, and 3 were summed, resulting in the reported uncertainties of DIC export . In our calculations, GPP lake comprises CO 2 uptake by phytoplankton from all available sources (Lewis Jr. 2011), and MIN lake is a lumped value for the total pelagic and sediment mineralization of autochthonous and allochthonous organic carbon (Pace and Prairie 2005). For a more detailed description of the values used to calculate DIC export , consult the respective studies cited in Table 1.

Sensitivity of terrestrial DIC export to uncertainties in global flux estimates
The lateral DIC export from terrestrial ecosystem to inland waters is currently highly uncertain (Drake et al. 2017), complicating the calculation of global terrestrial carbon sinks/sources (Butman et al. 2016). The different global estimates used to calculate DIC export (Table 1) each entail an uncertainty that sums to the uncertainty of our calculated DIC export . Uncertainties in DIC export have previously been estimated as~±1.0 Pg C yr −1 , which is close to the uncertainty of our calculated DIC export of 1:52 þ1:09 −0:90 (Eq. 1). Although DIC export lies most likely between 0.70 and 1.52 Pg C yr −1 , the full range of the calculated DIC export when considering the upper and lower boundary value for active carbon transformers and intermediate active carbon transformers is 0.39 to 2.61 Pg C yr −1 (Fig. 3).
To demonstrate the relevance of our functional lake classification concept in comparison to the uncertainties of the global flux estimates used in Eqs. 1, 2, and 3, we performed an analysis on the sensitivity of DIC export to uncertainties in the global flux estimates used (Table 1). We varied each respective value used in Eqs. 1, 2, and 3 by ± 25% and compared our calculated DIC export to the DIC export calculated with an error of ± 25% (DIC export_error ). While the differences in lake functioning resulted in a variation of DIC export of 0.82 Pg C yr −1 , DIC export_error differed from DIC export by 0.76, 0.35, and 0.51 Pg C yr −1 for active carbon transformers, intermediate active carbon transformers, and passive open channels, respectively. Thus, when assuming an error of ± 25% for each respective CO 2 production by mineralization in rivers MIN lotic 0.40 (Caraco and Cole 1999) flux term in Eqs. 1, 2, and 3, the variability in DIC export resulting from different assumptions on lake functioning was larger than the variability between DIC e x p o r t and DIC export_error . These results further highlight the importance of considering lake functioning when calculating DIC export . When testing the sensitivity of DIC export to lake functioning, we set GPP lotic and CCP lotic to zero, since to our knowledge, no global estimates on GPP lotic and CCP lotic exist. To examine the effect of setting GPP lotic (that is definitely > 0) to zero when testing the sensitivity of DIC export to lake functioning, we assumed GPP lotic as 25% of GPP lake and re-calculated DIC export . Adding GPP lotic resulted in an increase of DIC export for active carbon transformers (Eq. 1) of 21%. We did not test for the effect of non-zero GPP lotic for the other lake types, since we had assumed GPP lake to approach zero for intermediate active and passive lakes, and GPP lotic is usually considerably smaller than GPP lake .

Spatial variations in lake functioning
Our calculated DIC export estimates were based on the assumption that all lakes on Earth are either active carbon transformers, intermediate active carbon transformers, or passive open channels. This assumption does not depict reality since in some regions on Earth, lakes will rarely function as passive open channels or active carbon transformers. Although it is beyond the scope of this study to allocate a lake function to each of the 117 million lakes on Earth, we made a rough estimate of how many lakes potentially can function as active carbon transformers or as passive open channels. Based on Lewis Jr. (2011), we assumed that all lakes located between 39°N and 39°S have the potential to function as active carbon transformers. We chose this latitude as a borderline for potentially active lakes, since the modeled median lake gross primary production of global lakes increases sharply between 42.5°and 37.5°latitude from~400 to~800 g C m −2 yr −1 (Lewis Jr. 2011). Accordingly, we assumed that all lakes between 54°N and 84°N may function as passive open channels as these lakes are located in the boreal and subarctic zone, are usually small and shallow, and often function as passive open channels (Weyhenmeyer et al. 2015). Even if not taken into consideration in this estimate, it must be noted that a substantial number of humic-rich boreal lakes most probably function as intermediate active carbon transformers. We did not assign lakes located between 39°and 54°latitude to any of our three categories, since the functioning of lakes within these latitudinal bands might vary strongly depending on nutrient conditions, hydrology, catchment characteristics, lake morphology, and regional climate.
Using the abundance and total area of lakes for 3°latitudinal bands from the Global Water Body database (GLOWABO) (Verpoorter et al. 2014), we found that about 25% of lakes on Earth, corresponding to 35% of the global lake area, potentially function as passive open channels, while 60% of lakes, corresponding to 45% of the global lake area, might act predominantly as active carbon transformers. Since the functioning of different lakes within an ecoregion varies   Raymond et al. (2013), and Lauerwald et al. (2015) widely Weyhenmeyer et al. 2015), this estimate is a first-order approximation of the global distribution of lake functions according to our conceptual model and can serve as an exemplifying application of our classification concept. This classification according to latitudinal distribution only accounts for the control of climate on lake functioning. Additional factors including hydrologic regime, land-use, and regional geography can exert strong influences on lake functioning and should be considered in future, more accurate estimates of the functioning of lakes in the global carbon cycle.

Refining terrestrial DIC export estimates
We focused our analysis on the role of lakes as modulators of DIC transport and transformation. It has to be noted that the functioning of streams, rivers, floodplains, and wetlands is also of high importance for DIC transport and transformation along the LOAC (Aufdenkampe et al. 2011;Raymond et al. 2013). At the current stage, the use of our conceptual model for calculating the global lateral DIC export from terrestrial ecosystems to inland waters has limitations, which future studies should address, when more whole-lake carbon budgets and estimates of global carbon fluxes are available. At present, our lake classification concept (Figs. 1 and 2) applies only to open lake systems along the LOAC. While this is the prevalent lake type globally, closed basins can play an important role in arid and semi-arid regions (Einsele et al. 2001;Li et al. 2017). The function of these lakes in the global carbon cycle is presently unknown; however, a recent study estimated that in global closed drainage basins, about 0.15 Pg C yr −1 of the incoming DIC is stored (Li et al. 2017). Our estimates are also limited by the lack of data for global GPP lotic and CCP lotic . We assumed these fluxes to be zero in Eqs. 1, 2, and 3, but performed a sensitivity analysis to demonstrate the effect of GPP lotic on our DIC export estimate according to Eq. 1. As GPP lotic and CCP lotic are larger than zero, including estimates for GPP lotic and CCP lotic would increase the calculated DIC export rates. Further, the global estimate of CCP lake (Table 1) used in Eq. 1 is relatively low, considering that some lakes, e.g., Attersee and Lake Constance, precipitate 4-17% of the incoming carbonate (Einsele et al. 2001). Thus, we suggest that our estimates for DIC export are conservative, particularly since we used the most recent global estimate for CO 2 _emission lotic from Lauerwald et al. (2015) that is substantially lower than earlier estimates of CO 2 _emission lotic (Table 1).
A significant amount of uncertainty in our calculated DIC export arises from the variability in lake primary production that complicates quantification of global lake CO 2 consumption. Lake primary production along a global latitudinal gradient can vary by a factor of 1000 (Jonsson et al. 2003;Melack and Kilham 1974) and significantly differs for individual lakes on decadal scale in relation to climate variation (Pettersson et al. 2003). Depending on the estimates for global lake abundance used, Lewis Jr. (2011) reports a variation of global lake GPP between 1.05 and 1.51 Pg C yr −1 . This uncertainty of 0.46 Pg C yr −1 accounts for a variation in our calculated DIC export for active carbon transformers (1.52 Pg C yr −1 acc. to Eq. 1) of about 30%. Neglecting GPP lake in the global DIC budget that considers lakes as active carbon transformers (Eq. 1) would result in a reduction of the calculated DIC export by about 85%. In our calculations, a large share of GPP lake is from GPP in warm, nutrient-rich tropical lakes (Lewis Jr. 2011), which we classified as active carbon transformers. We assumed that all photosynthetically fixed CO 2 in lakes is of terrestrial origin. This assumption is supported by the fact that a large majority of inland waters is supersaturated with CO 2 and CO 2 uptake from the atmosphere is minor in global lakes. However, the relative contribution of terrestrial vs. atmospheric-derived CO 2 to photosynthetic CO 2 fixation in inland waters is presently unknown (Drake et al. 2017).
Assuming that a share of the photosynthetically fixed CO 2 is taken up by lakes directly from the atmosphere would lower our estimates of DIC export slightly. In lakes with high hydrologic DIC inputs, primary production can be influenced by terrestrially derived DIC, and CO 2 emissions in net autotrophic lakes can mainly arise from hydrologic DIC inputs (Stets et al. 2009). Consequently, CO 2 consumption by primary production in lakes along the LOAC needs to be considered when calculating terrestrial DIC export rates.
The single published estimates for the same global inland water DIC flux differ significantly (Table 1). The accuracy of the global estimates for inland water DIC fluxes partly depends on the resolution at which spatial variabilities are accounted for (Lauerwald et al. 2015;McDonald et al. 2013). The spatial resolution currently used to calculate carbon fluxes along the LOAC of maximum 0.5°is too coarse to account for the diversity and regional distribution of soil types, wetlands, streams, rivers, and lakes (Lauerwald et al. 2015;Regnier et al. 2013). However, attempts to resolve, for example, CO 2 emissions from the global river network at a higher resolution have progressed constantly during the past 10 years (Lauerwald et al. 2015). With ongoing refinements of global estimates of inland water DIC fluxes, our conceptual model, accounting for the predominant lake functions, can be a valuable tool for more robust estimates of the global DIC export from terrestrial ecosystems to inland waters.
Overall, we show here that the consideration of lake functioning is necessary to estimate the magnitude of the global DIC export from terrestrial ecosystems to inland waters. Similar to our estimations for DIC, we expect that our lake classification concept (Fig. 2) can also be applied to other carbon forms (e.g., DOC) as its transport and transformation processes vary strongly among lakes depending on lake characteristics. Since lake functioning differs widely across the globe (Lewis Jr. 2011;Tranvik et al. 2009), accurate estimates of the lateral DIC export will require the prediction of lake functioning for each of the 117 million lakes on Earth Verpoorter et al. 2014). With constant changes in anthropogenic carbon outputs altering the global carbon cycle , and disturbances of natural conditions causing eutrophication or a global temperature rise, lake functioning is likely to change (Gudasz et al. 2010;Lewis Jr. 2011 ; Fig. 2). The continuously high activity in dam construction on global scale (Zarfl et al. 2015) will probably shift numerous riverine systems towards a state at which they predominantly act as active carbon transformers. Global climate change will likely increase the activity of already existing lakes and reservoirs in carbon transformation, especially in the temperate and boreal region (Flanagan et al. 2003;Gudasz et al. 2010;Tranvik et al. 2009). Thus, our functional lake classification concept becomes particularly important for the calculation of the future lateral DIC export from soils to inland waters and future estimations of terrestrial CO 2 sinks and sources.
Author information Correspondence and requests for materials should be addressed to Fabian.Engel@ebc.uu.se Funding information Financial support was received from the Swedish Research Council (Grant No. 2016-04153), the European Union's Horizon 2020 research and innovation program under the Marie Sklodowska-Curie grant agreement No 643052 (C-CASCADES project), and from the Knut and Alice Wallenberg Foundation (KAW project). This work profited from the Global Lake Ecological Observatory Network (GLEON). The Estonian partners were supported by institutional research funding IUT 21-02 of the Estonian Ministry of Education and Research, and American partners were supported, in part, through NSF EF-1137327.

Compliance with ethical standards
Conflict of interest The authors declare that they have no conflict of interest.
Open Access This article is distributed under the terms of the Creative Comm ons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.