Increases in terrestrially derived carbon stimulate organic carbon processing and CO2 emissions in boreal aquatic ecosystems

Abstract

The concentrations of terrestrially derived dissolved organic carbon have been increasing throughout northern aquatic ecosystems in recent decades, but whether these shifts have an impact on aquatic carbon emissions at the continental scale depends on the potential for this terrestrial carbon to be converted into carbon dioxide. Here, via the analysis of hundreds of boreal lakes, rivers and wetlands in Canada, we show that, contrary to conventional assumptions, the proportion of biologically degradable dissolved organic carbon remains constant and the photochemical degradability increases with terrestrial influence. Thus, degradation potential increases with increasing amounts of terrestrial carbon. Our results provide empirical evidence of a strong causal link between dissolved organic carbon concentrations and aquatic fluxes of carbon dioxide, mediated by the degradation of land-derived organic carbon in aquatic ecosystems. Future shifts in the patterns of terrestrial dissolved organic carbon in inland waters thus have the potential to significantly increase aquatic carbon emissions across northern landscapes.

ARTICLE

Received 3 Sep 2013 | Accepted 18 Nov 2013 | Published 13 Dec 2013

The concentrations of terrestrially derived dissolved organic carbon have been increasing throughout northern aquatic ecosystems in recent decades, but whether these shifts have an impact on aquatic carbon emissions at the continental scale depends on the potential for this terrestrial carbon to be converted into carbon dioxide. Here, via the analysis of hundreds of boreal lakes, rivers and wetlands in Canada, we show that, contrary to conventional assumptions, the proportion of biologically degradable dissolved organic carbon remains constant and the photochemical degradability increases with terrestrial inuence. Thus, degradation potential increases with increasing amounts of terrestrial carbon. Our results provide empirical evidence of a strong causal link between dissolved organic carbon concentrations and aquatic uxes of carbon dioxide, mediated by the degradation of land-derived organic carbon in aquatic ecosystems. Future shifts in the patterns of terrestrial dissolved organic carbon in inland waters thus have the potential to signicantly increase aquatic carbon emissions across northern landscapes.

DOI: 10.1038/ncomms3972

Increases in terrestrially derived carbon stimulate organic carbon processing and CO2 emissions in boreal aquatic ecosystems

Jean-Franois Lapierre1, Franois Guillemette1,2, Martin Berggren1,3 & Paul A. del Giorgio1

1 Groupe de Recherche Interuniversitaire en Limnologie et en Environnement Aquatique (GRIL), Dpartement des Sciences Biologiques, Universit du Qubec Montral, Case Postale 8888, succursale Centre-Ville, Montral, Quebec, Canada H3C 3P8. 2 Department of Ecology and Genetics/Limnology, Uppsala University, SE-752 36 Uppsala, Sweden. 3 Department of Physical Geography and Ecosystem Science, Lund University, SE-223 62 Lund, Sweden. Correspondence and requests for materials should be addressed to J-F.L. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 4:2972 | DOI: 10.1038/ncomms3972 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

& 2013 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3972

Emissions of carbon dioxide (CO2) from inland waters to the atmosphere are now recognized as signicant relative to other components of the Earth carbon system14. These

emissions are fuelled by land-derived carbon1,5,6, and one of the pathways leading to widespread aquatic CO2 supersaturation involves the transport of dissolved organic carbon (DOC) from land to continental waters and its subsequent degradation within the aquatic network. There is now mounting evidence that shifts in climate and other environmental factors are driving increasing aquatic concentrations of terrestrially derived DOC in many landscapes of the world79, in what has been termed the browning10 of inland waters. Whereas the magnitude and regulation of the browning phenomenon are the focus of intense research and debate7,1117, the potential implications of such a shift in terrestrial carbon (C) concentrations for the functioning of the receiving aquatic systems have been less explored, particularly in terms of the impacts on aquatic CO2 emissions.

In this context, understanding the mechanistic underpinnings of the links between terrestrial DOC and aquatic CO2 becomes critical to our ability to predict the potential impact of current and projected shifts in the delivery of this terrestrial carbon on global aquatic CO2 emissions18. DOC and CO2 concentrations are indeed positively related across temperate and boreal lakes1922 but the mechanisms underlying this relationship are still debated19. Should the positive relationship between CO2 and DOC be mostly driven by the covariation between the loading of terrestrial DOC and that of soil-derived CO2, then the use of this relationship to understand and predict the effects of increasing DOC on surface water CO2 dynamics is questionable.

If, in contrast, the relationship between DOC and CO2 is to a large extent driven by the actual degradation of terrestrially derived DOC into CO2 within aquatic systems, then this relationship has a causal underpinning and can be used to predict future scenarios involving changes in the concentrations of terrestrial DOC in continental waters.

Understanding the potential for the degradation of terrestrially derived DOC is therefore essential to establish the link between DOC and CO2 dynamics. It has long been recognized that terrestrial DOC is susceptible to photochemical degradation in the water23, but it has traditionally been assumed that this DOC is generally recalcitrant to biological degradation2426, and therefore that the contribution of this latter pathway to aquatic CO2 production is modest. This view is rapidly changing24,25, and a growing volume of evidence indicates that terrestrial DOC may be readily degraded by biological processes in aquatic environments, such as the Amazon basin27 and the Arctic8,28, temperate29,30 and boreal31,32 biomes. Boreal biomes contain the highest density of freshwater on Earth, yet this evidence for biologically degradable DOC is scattered and based on isolated components of the continental aquatic network. As a result, the overall susceptibility of DOC to photochemical and biological degradation, and how this susceptibility varies over gradients of increasing terrestrial inuence and at scales relevant for continental C budgets, has not been well established19,3336.

Here we explicitly address the mechanistic underpinnings of the relationship between aquatic CO2 and DOC and, in particular, we explore the concentrations of photochemically degradable DOC and Bd-DOC, and surface water CO2 dynamics in boreal aquatic networks. We base this work on simultaneous measurements of Biologically degradable DOC and photochemically degradable DOC, as well as ambient CO2 concentrations and uxes, in hundreds of lakes, wetlands and rivers across the northern boreal landscapes of Qubec (Canada), spanning a very large gradient of both DOC and terrestrial inuence. We rst assess how these two facets of DOC degradability vary along a gradient of increasing terrestrial inuence (using chromophoric

dissolved organic matter (CDOM) as a proxy (see Supplementary Note 1) and then explore to what extent DOC degradability mediates the relationship between terrestrial DOC and aquatic CO2 concentrations and uxes.

ResultsEnvironmental gradients covered. The conditions covered in this study encompass much of the range in environmental conditions reported in the literature for northern freshwaters (Supplementary Table S1). In particular, surface water concentrations of DOC ranged from 2.1 to 123.9 mg l 1 (median 9.4 mg l 1) across all sites, and CDOM ranged from

0.13 to 30 m 1 (median 4.4 m 1, absorption coefcient mea

sured at 440 nm in Naperian units); both were strongly, positively correlated (r2 0.60, Supplementary Fig. S1), suggesting that

much of the variability in ambient DOC concentrations is driven by differences in the loading of terrestrial DOC. The vast majority of ecosystems was supersaturated in CO2 (average partial pressure of CO2 (pCO2) was 754, 3,396 and 5,896 matm for lakes, rivers and wetlands, respectively), and thus emitted CO2 to the atmosphere (average CO2 ux was 450, 1,459 and 960 mg C m 2 d 1 for lakes, rivers and wetlands, respectively).

DOC degradability patterns. To reconstruct the links between terrestrial DOC and aquatic CO2 uxes, we rst assessed the patterns in DOC degradability based on measurements of the amount of DOC lost to biological and photochemical processes under standardized conditions (see Methods). There were positive relationships between the amounts of both Biologically degradable DOC (Bd-DOC) and photochemically degradable DOC (Pd-DOC) and CDOM concentrations across lakes, wet-lands and rivers (Fig. 1). Although it could be a priori expected that higher amounts of degradable DOC would result from increasing overall amounts of DOC, it is the nature of the relationships reported here that is of particular interest. The slope of the type II regression (see Methods) between Bd-DOC and CDOM (Fig. 1a; log-slope 1.15, 95% condence interval

0.971.35, Po0.001) indicates that across the whole gradient of systems studied, the DOC pool in the brownest systems with the highest terrestrial inuence was as biologically degradable (in terms of the proportion that can be removed) as the DOC pool in clearer waters. In systems with higher terrestrial inuence, DOC was in fact proportionately more photochemically degradable (Fig. 1b; log-slope 1.54, 95% CI 1.401.68, Po0.001). As a

result, the percentage of Bd-DOC shows no signicant relation

ship with CDOM (P 0.82) but the percentage of Pd-DOC

increases signicantly with CDOM (r2 0.30, Po0.0001,

Supplementary Fig. S2).

Linking DOC processing to CO2 dynamics. Stable or increasing DOC degradability over gradients of terrestrial inuence suggest that increasing CDOM should lead to enhanced conversion of DOC to CO2 in aquatic ecosystems, and thus to higher concentrations and uxes of CO2. There was indeed a positive relationship between ambient pCO2 and CDOM across all our systems ( Supplementary Fig. S3), but whether or not this pattern may be directly attributed to the conversion of DOC to CO2 within those systems cannot be determined by this simple relationship. Alternative mechanisms unrelated to direct DOC degradation, including the covariation in the loading of soil-derived DOC and CO2 to aquatic systems, could also explain this relationship. To explicitly address these different scenarios, we further explored the causal links that may exist between terrestrially derived DOC, DOC degradability, and excess aquatic CO2 concentrations and ux (fCO2, in mg C m 2 d 1) using

2 NATURE COMMUNICATIONS | 4:2972 | DOI: 10.1038/ncomms3972 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2013 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3972 ARTICLE

structural equation modelling (SEM)37. In addition to empirically determining the existence of a link between CDOM and fCO2, this approach allows to determine the pathways by which CDOM is most probably inuencing CO2. We tested several alternative structures linking CDOM to CO2 for all the sampled sites (see Supplementary Fig. S4 for key examples), and the only statistically signicant structure is one where CDOM-related increases in both Bd-DOC and Pd-DOC positively inuence aquatic pCO2, and

subsequently aquatic fCO2 (Fig. 2a). The two paths leading from CDOM to fCO2 were highly signicant, suggesting that both biological and photochemical processing of DOC contribute signicantly to the observed CO2 dynamics. Note that this structure did not contain a signicant direct link between CDOM and pCO2, the latter reecting variables that inuence pCO2 and that may co-vary with CDOM, such as soil-derived CO2.

Interestingly, although there was a signicant positive relationship between pCO2 and DOC across all systems (r2 0.28, Po0.001), there were no signicant structures linking

DOC, DOC degradability and CO2 (Supplementary Fig. S4), thus suggesting that it is the degradation of the terrestrial component of DOC (CDOM) that is more directly inuencing aquatic CO2 dynamics.

The effect of physical forcing. While Fig. 2 suggests an important role of DOC processing for aquatic CO2 dynamics, it also shows that at this scale, the regulation of CO2 dynamics is not homogenous across all systems. Roughly half of the variation in pCO2 was explained by factors other than DOC degradation, and, considering that uxes are a direct function of concentration, the relationship between excess pCO2 and fCO2 was unexpectedly weak (r2 0.22). Physical forcing, however, strongly mediates the

pCO2 versus fCO2 relationship, with low-turbulence (expressed by the gas exchange coefcient, kCO2 in m d 1) sites sustaining higher average CO2 concentrations for a given CO2 production or input rate than systems that typically have higher turbulence38. In this regard, there was a clear bimodal distribution of the empirically determined kCO2 around a median of 1 m d 1 across all sampled sites (Fig. 3a), and the regression of fCO2 against pCO2 drastically improved when sites with kCO2 over and below 1 m d 1 were considered separately (Fig. 3b). This suggests that, at the expense of the generality of the relationship presented in Fig. 2, separating sites into more homogenous groups with regard to physical structure provides a more robust quantication of the impact of changing terrestrial DOC concentrations on aquatic CO2 uxes.

Quantifying the effect of CDOM on CO2 uxes. On the basis of above results, we repeated the SEM analysis separately for these two subgroups of sites that were split based on their kCO2 (sites

below and above 1 m d 1, Fig. 4). The structures that emerged as signicant for each separate group generally resembled that obtained for the ensemble of data shown in Fig. 2, but there were some differences in the resulting links and parameters with implications for the prediction of aquatic CO2 emissions. As mentioned above, there was no signicant direct link between CDOM and pCO2 when all the systems are considered (Fig. 2), suggesting that the effect of CDOM on pCO2 at the whole-network scale is predominantly degradation mediated. This direct link became signicant, however, when high and low kCO2

a

Wetlands

LakesRivers and streams

Biologically degradable DOC

(mg Cl1 )

1

0.1

1 10

1 10

CDOM (a440 nm ; m1)

CDOM (a440 nm ; m1)

b

100

NS* 0.77

Photochemically degradable DOC

(mg Cl1 )

10

10

r 2=0.35

0.59

0.53 0.30

pCO2

fCO2

0.47

Pd-DOC

Bd-DOC r 2=0.59

r 2=0.51

r 2=0.22

1

0.1

0.01

Figure 1 | Concentrations of degradable DOC across gradients of terrestrial inuence. Concentrations of Bd-DOC and Pd-DOC increase as a function of ambient CDOM for natural water samples from a wide range of aquatic systems in the boreal and temperate landscapes of Qubec, Canada. (a) Log(Bd-DOC) 1.15 Log(CDOM) (95% condence interval

(CI) 0.971.35); 1.00, Po0.0001, n 345, r2 0.30. (b) Log(Pd-

DOC) 1.54 Log(CDOM)(95% CI 1.401.68); 0.90, Po0.0001, n 359,

r2 0.59. The dashed lines represent a slope of 1, that is, a proportional

increase in degradable DOC as CDOM increases.

CDOM

Figure 2 | Linking terrestrially derived DOC to aquatic uxes of CO2.

Structural equation model linking changes in aquatic CDOM to changes in CO2 uxes to the atmosphere. Arrows denote directional links between two variables, and the coefcients represent the rate at which a response variable varies as a function of a given change in its predictor, based on centered and standardized data (see Methods). All paths are signicant at Po0.001. The dashed arrow denotes a non-signicant (NS*) path that was not included in the nal model presented here; it is included only for presentation purposes to emphasize the absence of a signicant direct link between CDOM and pCO2 or CO2 ux. pCO2 represents excess pCO2.

Degradation-mediated CDOM indirect effect coefcient on fCO2 is 0.25. n 216. w2 5.20, P 0.39, df 5.

NATURE COMMUNICATIONS | 4:2972 | DOI: 10.1038/ncomms3972 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

& 2013 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3972

systems were analysed separately, implying that a fraction of the CDOM effect was not due to DOC degradation when more homogenous groups of sites were considered. Moreover, the path connecting Pd-DOC and pCO2 was non-signicant in low-turbulence systems (Fig. 4a). As a result, most of the overall CDOM effect on fCO2 (0.39, Fig. 4b) was mediated by Bd-DOC and Pd-

DOC (combined coefcient 0.30, Fig. 4b) for high kCO2 systems

(mostly lakes and fast-owing rivers), while roughly half (coefcient 0.20, Fig. 4a) of the overall CDOM effect on fCO2

(coefcient 0.37) was not mediated by either Bd-DOC or Pd-

DOC for low kCO2 systems (mostly wetlands and slow-owing rivers). These coefcients can then be used to estimate how CO2 uxes respond to increases in CDOM and, together, the different SEM analyses thus provide insights on the links that exist between terrestrial inuence and CO2 emissions in aquatic ecosystems, in particular, on the role that DOC processing plays in that framework.

DiscussionCollectively, our results highlight the importance of terrestrially derived DOC in fuelling both Bd-DOC and Pd-DOC pools across a wide range of aquatic environments. The variability shown in Fig. 1 suggests that such patterns may have been masked in previous studies conducted at smaller scales or across narrower environmental gradients where autochthonous DOC may have emerged as a more important driver of overall DOC degradability. At the scale of the present study, however, increasing terrestrial inuence (as reected in CDOM) clearly implies higher absolute and relative potentials for DOC to be degraded across lakes, rivers and wetlands. Although those patterns are based on potential degradation rates under standardized conditions, they unambiguously show that both the total amount and the proportion of DOC that can potentially be converted to CO2 tend to increase with terrestrial inuence and, as a corollary, it would be expected that surface water CO2 concentrations should increase as a function of CDOM.

The SEM analyses (Figs 2 and 4) demonstrate that there is a link from CDOM to CO2 uxes that is mediated by the

a

100

80

60

Number of sites

40

20

0

0.1 10 100

1

kCO2 (m d1)

KCO2 < 1

KCO2 > 1

b

CO 2fluxes (mg m2 d1 )

10,000

1,000

100

1010 100 1,000Excess pCO2 (atm)

10,000

Figure 3 | Relationship between concentrations and uxes of CO2.

(a) Frequency distribution of the measured diffusive kCO2 across the sampled sites. Note the bimodal distribution around kCO2 1 m d 1.

(b) Regression of CO2 uxes to the atmosphere as a function of surface water CO2 excess partial pressure for systems with kCO2 over and below 1.

Analysis of covariance reveals different intercepts (Po0.001) but similar slopes (P 0.08) for both groups. kCO241: Log(fCO2) 0.73 (0.03)

Log(excess CO2) 0.77 (0.08), n 247, Po0.0001, r2 0.74.kCO2o1:

Log(fCO2) 0.73 (0.03)Log(excess CO2) 0.26 (0.13), n 172,

Po0.0001, r2 0.68. Dashed line illustrates the inadequacy of an overall t

in the absence of grouping: Log(fCO2) 0.36 (0.03)Log(excess

CO2) 1.67 (0.09), n 419, Po0.001, r2 0.27. Note that sample size

differs here compared with Figs 2 and 4, hence the different coefcients of determination between pCO2 and fCO2.

CDOM

a b

CDOM

r 2=0.26

0.47

Pd-DOC

Bd-DOC r 2=0.40

r 2=0.36

r 2=0.58

r 2=0.30

0.50

0.35 0.24

pCO2

fCO2

0.82

0.80

0.11

Bd-DOC r 2=0.64

r 2=0.36

r 2=0.67

Pd-DOC

0.46 NS*

pCO2

fCO2

0.76

0.56

0.26

Figure 4 | The pathways linking CDOM to CO2 uxes in high- and low-turbulence systems. Structural equation models linking changes in aquatic CDOM to changes in CO2 uxes to the atmosphere for the two subsets of sites separated on the basis of kCO and therefore of average turbulence. Arrows denote directional links between two variables, and the coefcients represent the rate at which a response variable varies as a function of a given change in its predictor, based on centered and standardized data (see Methods). All paths are signicant at Po0.001. The dashed arrow denotes a non-signicant (NS*) path that was not included in the nal model presented here; it is included only for presentation purposes to emphasize the absence of a signicant direct link between Pd-DOC and pCO2 in low-turbulence systems. pCO2 represents excess pCO2. (a) Low-turbulence systems (kCO2o1 m d 1). CDOM indirect effect coefcient on fCO2 is 0.37, degradation-mediated CDOM indirect effect on fCO2 is 0.17. n 102. w2 5.07, P 0.41, df 5. (b) High-

turbulence systems (kCO241 m d 1). CDOM indirect effect coefcient on fCO2 is 0.39, degradation-mediated CDOM indirect effect on fCO2 is 0.30. n 114. w2 5.45, P 0.24, df 4.

4 NATURE COMMUNICATIONS | 4:2972 | DOI: 10.1038/ncomms3972 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2013 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3972 ARTICLE

degradation of DOC when very large gradients are considered, such that increasing concentrations and subsequent degradation of terrestrial DOC in aquatic ecosystems should indeed result in rising CO2 emissions. These results, however, do not imply that coloured DOC itself is the only pool being degraded, but rather that systems with high CDOM concentrations are associated with a higher absolute potential for conversion of DOC, both coloured and uncoloured, into CO2 in surface waters. They do not imply either that DOC degradation is the only factor inuencing aquatic CO2 dynamics in these northern regions. For example, recent inter-regional assessments have highlighted the importance of terrestrial inorganic C inputs in shaping aquatic CO2 dynamics6,33,39. Similarly, studies that have targeted specic components of the aquatic network, such as low DOC lakes40 or headwater streams41, have also highlighted the importance of soil CO2 inputs. There is no actual contradiction between these various studies and the results we present here: the direct inuence of terrestrial DOC is superimposed onto a background of CO2 dynamics driven by external inputs of CO2. This background is always present but varies both regionally and also locally between different types of aquatic systems19,39, determined by variations in the relative proportion of DOC versus CO2 delivered from land to the inland water network.

The large heterogeneity that exists among the studied systems was at least partly captured by the SEM analyses performed on subsets of sites sharing comparable patterns in gas exchange coefcient (Fig. 4).When more homogenous sets of systems, in terms of physical structure, were analysed, factors other than DOC degradation emerged as modulating the effect of CDOM on CO2 concentrations and uxes. We can only speculate on the causes underlying the differences in pathways linking CDOM to CO2 in high kCO2 versus low kCO2 systems. It is possible that physical structure integrates other environmental factors that may modulate the impact of CDOM on CO2 dynamics, such as area and depth42, which may in turn be related to residence time, exposure to sunlight or the relative importance of groundwater inputs in regulating CO2. Regardless of the actual processes, the different SEM analyses converge to essentially the same conclusion, which remains adequate whether the entire data set or more homogenous subsets of sites are considered. The effect of terrestrially derived DOC on CO2 evasion is not entirely but signicantly driven by DOC processing in continental watersheds.

Our results provide a mechanistic underpinning to previous studies that had arrived to similar conclusions concerning the role of DOC but based on regional C mass balances32, and also provide a process-based framework underlying the regional relationships that have been observed between DOC concentrations and aquatic CO2 uxes19,21,22. We show that

DOC degradation strongly mediates this relationship, across rivers ranging from stream order 0 to 6, lakes ranging in size and morphometry from a few hectares to several thousands of square kilometres, and wetlands and ponds that also vary greatly in position and features, within a heterogeneous boreal landscape covering roughly 900,000 km2. Our results may not necessarily apply to some specic aquatic components or portions of the landscape but rather provide an integrated perspective of the role of DOC processing in sustaining aquatic CO2 concentrations and evasion in continental watersheds.

The SEM analyses further allow a rst-order quantication of the potential effect of shifts in terrestrial DOC concentrations on aquatic fCO2. Although the long-term DOC trends are highly variable in time and across landscapes11,12,15,16, DOC concentrations have consistently increased across temperate and boreal regions over the past 20 years7. More specically, in our own study area, DOC concentrations have increased between 20

and 50% since 1990 (ref. 43), whereas CDOM has increased by up to 155% within the same period in some European landscapes15. The coefcients provided by the SEMs (Fig. 4) suggest that for every percent increase in CDOM there should be a 0.17 and0.30% increase in CO2 evasion from low and high kCO2 systems, respectively, owing to enhanced DOC processing in those systems. Our results imply that as a consequence, aquatic CO2 uxes may have been increasing in northern Qubec and elsewhere across northern temperate and boreal landscapes over the past decades, conceivably from 10% to up to 50% where terrestrial DOC concentrations have increased the most. There are, to our knowledge, no long-term data sets on aquatic CO2 that can be used to test this hypothesis, which we suggest should be incorporated into future long-term research programmes.

Our results further suggest that aquatic CO2 emissions will keep increasing as a direct function of increasing terrestrially derived DOC concentrations in aquatic ecosystems. Considering the major contribution of the aquatic networks in these northern, water-rich landscapes to the global continental aquatic C emissions1,2,4 (ca. 1.5 Pg y 1), these potential shifts could signicantly alter the ne balance that exists between freshwater

CO2 emissions and the estimated continental (2.8 Pg C y 1) and oceanic carbon sinks (2.2 Pg C y 1)4. Although more focused measurements should be undertaken to adequately quantify the magnitude of the enhanced aquatic CO2 emissions at continental or global scales, the scenario that we propose here undoubtedly involves the mobilization and degradation of terrestrial organic C, which, contrary to historical assumptions, appears to be highly reactive in freshwaters but would otherwise experience much slower degradation in more stable environments such as soil or oceanic refuges1,44,45. The absence of evidence for generalized long-term declines in discharge, and the positive relationships that have been reported between long-term trends in discharge and DOC concentrations in landscapes where browning is happening11,12,43, suggest that greater amounts of DOC are circulating in inland waters today compared with a few decades ago, yet no increase in the export of DOC from continents to the oceans have been observed so far in regions where browning is happening46. These observations are consistent with our ndings that this additional carbon is readily processed and converted to CO2 as it enters the aquatic networks, which may further enhance the erosion of the terrestrial C (refs 4749) sink beyond a direct effect of increasing global temperatures.

Methods

Study area and sampling. Over the course of 8 years, we sampled lakes, streams, rivers and wetlands (mostly beaver ponds) located in seven temperate and boreal regions of Qubec, eastern Canada (4453N, 7080W, Supplementary Fig. S5). Sites were sampled from late May to early October and a subset of sites was sampled up to three times in a given year.

A more detailed description of the regional features can be found in Lapierre and del Giorgio19. In short, these regions cover a total area of B900,000 km2 and are characterized by a wide range of geography, climate and dominant plant cover: mean annual temperature ranges from 4.9 C in the Laurentians to 2.0 C in

the Eastmain region, and mean annual precipitations ranges from 695 mm y 1 in the Eastmain to 1,150 mm y 1 in the Eastern Townships regions. Dominant vegetation ranges from deciduous forest in the Eastern Townships and Laurentians regions to mixed deciduous and conifer forest in the Abitibi and Chicoutimi regions; the James Bay, Chibougamau and Eastmain regions are characterized by spruce-moss forest. Topography is generally at and lakes and rivers are typically shallow and turbid in the low-altitude (typically 200300 m elevation) Abitibi and James Bay regions, which are located in the Abitibi clay belt. Water discharge was very low in rivers and streams in the latter region, whereas it was the highest in the most elevated and steepest (200900 m elevation) Laurentians and Chicoutimi regions.

Sampling took place from 2003 to 2005 in Eastern Townships, and from 2009 to 2011 in the remaining regions. The Eastern Townships sites were sampled in June, July and August each year, the Eastmain sites were sampled twice from early June to early July 2009, the Laurentians sites were sampled once between mid-June and mid-August 2010, and three sampling campaigns were performed during the ice-

NATURE COMMUNICATIONS | 4:2972 | DOI: 10.1038/ncomms3972 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

& 2013 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3972

free period for the remaining regions in 2010 and 2011. For these latter regions, the rst campaign occurred from late May to early June, the second campaign from mid-July to mid-August and the third campaign from late September to early October. A subset of sites were sampled at each campaign and the remaining ones were sampled only once during mid-summer. For this analysis, we have used the indiviudal sampling points, without any averaging. Because of various logistical constraints, not all variables could be measured at every site all the time, which explains the different number of observations presented for the main variables. In addition, DOC photochemical degradability was not measured for the Eastern Townships sites.

Water was sampled at 0.5 m from the surface from the deepest measured point of lakes and near the shore of streams, rivers and wetlands. Ambient pCO2 was measured in situ using an EGM-4 infrared gas analyser (PP-systems, Boston, MA, USA) after equilibration by pumping water through a Liqui-Cel Mini Module (Charlotte, NC, USA) contactor membrane (for lakes) or by headspace equilibrium (rivers, streams and wetlands). Excess pCO2 was calculated by subtracting ambient atmospheric pCO2 at the corresponding sites. CO2 ux, and the associated diffuse gas exchange coefcients for CO2 (kCO2), were measured for every site on the basis of changes in pCO2 with time in oating chambers38. The chamber was connected in closed loop to an EGM-4 infrared gas analyser (PP-systems), pCO2 in the chamber was measured every minute for 10 min and accumulation rate was calculated by linear regression.

Chemical and optical analyses. Coloured dissolved organic matter (CDOM) was determined spectrophotometrically in 0.45-mm ltered natural water samples. We measured absorbance at 440 nm with a ultravioletvisible Ultrospec 2100 spectrometer (Biochrom, Cambridge, UK) using a 2-cm quartz cuvette and corrected by subtracting nanopure water blanks. The absorption coefcient (a440, in m 1)

was reported in Naperian units by dividing the optical absorbance at 440 nm by the path length in metres and multiplying by 2.303 (ref. 50).

Concentrations of DOC from ltered (0.45 mm) water samples were measured in duplicate (taken from a unique vial) on an OI-1010 TIC-TOC Analyzer (OI Analytical, College Station, TX, USA). Total phosphorous ( Supplementary Table S1) was analysed spectrophotometrically using the standard molybdenum blue method after persulfate digestion51, and Chl a samples ( Supplementary Table S1) were analysed spectrophotometrically following ltration on Whatman (Maidstone, UK) GF/F lters and hot ethanol (90%) extraction52. Filters were sonicated before extraction. The surface and catchment areas of the lakes we sampled were determined using the ArcView 3.2 and ArcGIS v.9 software applied on the DEM derived from (1:50,000) maps.

DOC degradability. The DOC degradation data presented here reect the potential for DOC to be mineralized under standardized conditions, and not to actual degradation rates under ambient conditions. Our approach aimed at removing or minimizing site-specic features such as local temperature and light exposure and thus to ensure that the resulting DOC degradability was comparable across the entire spectrum of sites.

Biological degradation experiments were carried out in the dark at a xed incubation time (14 days) and temperature (20 C), and the data presented in Fig. 1a thus correspond to the amount of Bd-DOC over a period of 14 days. We ltered 1.5 l of ambient water through 2.7 mm pore-size GF/F lters (Whatman)

for biological degradation experiments to preserve most bacteria and avoid re-inoculation53. We used the same pool of water for photochemical degradation experiments to minimize the differences in manipulations and thus maximize the comparability of both properties over the large range of samples. Two 500-ml subsamples were used as replicates for biological degradation experiments and dispensed in acid-washed glass bottles (500 ml). Figure 1a presents the average of both replicates. The remainder of the water was used for photochemical degradation experiments, which were carried out in a solar simulator (Qsun XE1-BC, Qlab, FL, USA) under a standard light dose (0.68 W m 2 at 340 nm, spectrum representative of natural sunlight) over a xed incubation time (24 h) and temperature (24 C). The data presented in Fig. 1b thus correspond to the amount of photochemically degradable DOC in 1 day at the maximum sustainable light intensity from the solar simulator; considering diel light cycles and the natural light intensity at the studied latitudes, this would roughly correspond to 6 days of natural light exposure in the sampled regions.

Irradiated samples were contained in 40 ml glass tubes, and thus the data presented would correspond to loss rates in the rst centimetres of the water column. Roughly 80% of the light between 360 and 700 nm passes through the glass, decreasing linearly to reach 10% at 315 nm; no light below 300 nm went through. The measured loss rates are thus conservative compared with what would have been obtained with quartz tubes, but still represent a level of exposure representative of the rst centimetres of the water column. Given the very strong light intensity and the small size of the vials, the light dose within the vials was high and the loss rates because of photochemical degradation were high, thus a 24-h period was sufcient to obtain reliable estimates of DOC losses. We assumed that all bacteria were killed by the very high ultraviolet dose. Even if that were not the case, the DOC loss rates owing to photochemical degradation were systematically at least an order of magnitude higher than those measured in the biological degradation experiments (when compared on a per day basis), implying that if

bacteria did survive, they could not contribute importantly to the measured loss rates. DOC concentrations were measured within 1 month after incubation.

These standardized biological and photochemical degradation rates could potentially be extrapolated to ambient waters by adjusting them for in situ temperature, light exposure and water retention time, among other factors, but this is not the focus of this study. Here we rather concentrate on the comparative analysis of this standardized DOC degradability across sites.

Statistical analyses. Data were log transformed to attain normality and homo-scedasticity. Single regressions were performed using JMP software version 7.0 (SAS Institute, NC, USA). We used major axis (type II) regressions to obtain a representative slope for the regression between two log-transformed variables containing measurement errors in Fig. 1. We performed SEM to identify the sequences of presumed causal relationships that may exist37 between CDOM, the amounts of Bd-DOC and Pd-DOC, ambient DOC, excess pCO2 and fCO2, using

EQS 6.1 (Multivariate Sofware, CA, USA). This analysis tests the signicance of individual paths between variables as well as that of the structure as a whole. Similar to most analyses, a path with Po0.05 is typically considered signicant.

Regarding the overall model, however, a value of P40.05 means that the proposed structure effectively captures the variability in the modelled variables and the directional links between them, whereas Po0.05 would imply a signicant difference between the observed and the modelled data37.

The coefcients provided by the SEM for the links between two variables (noted beside the arrows in Figs 2 and 4) are conceptually analogous to the slope of the linear regression between those variables, but the computation differs from major axis regressions, thus the coefcients provided by the different types of analyses cannot be directly compared37. Furthermore, the log-transformed data have been centred and scaled for the SEM analyses, such that different coefcients within a given structure can be directly compared. A higher coefcient between two variables denotes a stronger link than between two variables linked by a lower coefcient. Further details on the SEM analyses are provided in the Supplementary Note 2.

References

1. Cole, J. J. et al. Plumbing the global carbon cycle: Integrating inland waters into the terrestrial carbon budget. Ecosystems 10, 171184 (2007).

2. Tranvik, L. J. et al. Lakes and reservoirs as regulators of carbon cycling and climate. Limnol. Oceanogr. 54, 22982314 (2009).

3. Battin, T. J. et al. The boundless carbon cycle. Nat. Geosci. 2, 598600 (2009).4. Canadell, J. G. et al. Contributions to accelerating atmospheric CO(2) growth from economic activity, carbon intensity, and efciency of natural sinks. Proc. Natl Acad. Sci. USA 104, 1886618870 (2007).

5. Battin, T. J. et al. Biophysical controls on organic carbon uxes in uvial networks. Nat. Geosci. 1, 95100 (2008).

6. Butman, D. & Raymond, P. A. Signicant efux of carbon dioxide from streams anr rivers in the United States. Nat. Geosci. 4, 839842 (2011).

7. Monteith, D. T. et al. Dissolved organic carbon trends resulting from changes in atmospheric deposition chemistry. Nature 450, 537539 (2007).

8. Schuur, E. A. G. et al. The effect of permafrost thaw on old carbon release and net carbon exchange from tundra. Nature 459, 556559 (2009).

9. Larsen, S., Andersen, T. O. M. & Hessen, D. O. Climate change predicted to cause severe increase of organic carbon in lakes. Glob. Change Biol. 17, 11861192 (2011).

10. Roulet, N. & Moore, T. R. Environmental chemistry - browning the waters. Nature 444, 283284 (2006).

11. Erlandsson, M. et al. Thirty-ve years of synchrony in the organic matter concentrations of Swedish rivers explained by variation in ow and sulphate. Glob. Change Biol. 14, 11911198 (2008).

12. Clark, J. M. et al. The importance of the relationship between scale and process in understanding long-term DOC dynamics. Sci. Total Environ. 408, 27682775 (2010).

13. Evans, C. D., Chapman, P. J., Clark, J. M., Monteith, D. T. & Cresser, M. S. Alternative explanations for rising dissolved organic carbon export from organic soils. Glob. Change Biol. 12, 20442053 (2006).

14. Freeman, C., Evans, C. D., Monteith, D. T., Reynolds, B. & Fenner, N. Export of organic carbon from peat soils. Nature 412, 785785 (2001).

15. Chapman, P. J. et al. Changes in water colour between 1986 and 2006 in the headwaters of the River Nidd, Yorkshire, UK. Biogeochemistry 101, 281294 (2010).

16. Zhang, J. et al. Long-term patterns of dissolved organic carbon in lakes across eastern Canada: evidence of a pronounced climate effect. Limnol. Oceanogr. 55, 3042 (2010).

17. Schindler, D. W. et al. The effects of climatic warming on the properties of boreal lakes and streams at the Experimental Lakes Area, northwestern Ontario. Limnol. Oceanogr. 41, 10041017 (1996).

18. Billings, S. A., Ziegler, S. E., Schlesinger, W. H., Benner, R. & Richter, D. de B. Predicting carbon cycle feedbacks to climate: Integrating the right tools for the job. Eos Trans. AGU 93, 188188 (2012).

6 NATURE COMMUNICATIONS | 4:2972 | DOI: 10.1038/ncomms3972 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2013 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms3972 ARTICLE

19. Lapierre, J. F. & del Giorgio, P. A. Geographic and environmental drivers of regional differences in the lake pCO2 versus DOC relationship across northern landscapes. J. Geophys. Res.-Biogeosci. 117, G03015 (2012).

20. Sobek, S., Tranvik, L. J. & Cole, J. J. Temperature independence of carbon dioxide supersaturation in global lakes. Glob. Biogeochem. Cycle 19, GB2003 (2005).

21. Roehm, C. L., Prairie, Y. T. & del Giorgio, P. A. The pCO(2) dynamics in lakes in the boreal region of northern Quebec, Canada. Glob. Biogeochem. Cycle 23, GB3013 (2009).

22. Jonsson, A., Karlsson, J. & Jansson, M. Sources of carbon dioxide supersaturation in clearwater and humic lakes in northern Sweden. Ecosystems 6, 224235 (2003).

23. Graneli, W., Lindell, M. & Tranvik, L. Photo-oxidative production of dissolved inorganic carbon in lakes of different humic content. Limnol. Oceanogr. 41, 698706 (1996).

24. Holmes, R. M. et al. Lability of DOC transported by Alaskan rivers to the arctic ocean. Geophys. Res. Lett. 35, L03402 (2008).

25. Bianchi, T. S. The role of terrestrially derived organic carbon in the coastal ocean: a changing paradigm and the priming effect. Proc. Natl Acad. Sci. USA 108, 1947319481 (2011).

26. Hedges, J. I., Clark, W. A. & Cowie, G. L. Fluxes and reactivities of organic-matter in a coastal marine bay. Limnol. Oceanogr. 33, 11371152 (1988).

27. Ward, N. D. et al. Degradation of terrestrially derived macromolecules in the Amazon River. Nat. Geosci. 6, 530533 (2013).

28. Cory, R. M., Crump, B. C., Dobkowski, J. A. & Kling, G. W. Surface exposure to sunlight stimulates CO2 release from permafrost soil carbon in the Arctic. Proc. Natl Acad. Sci. USA 110, 34293434 (2013).

29. McCallister, S. L. & del Giorgio, P. A. Evidence for the respiration of ancient terrestrial organic C in northern temperate lakes and streams. Proc. Natl Acad. Sci. USA 109, 1696316968 (2012).

30. Guillemette, F., McCallister, S. L. & del Giorgio, P. A. Differentiating the degradation dynamics of algal and terrestrial carbon within complex natural dissolved organic carbon in temperate lakes. J. Geophys. Res. Biogeosci. 118, 963973 (2013).

31. Algesten, G. et al. Role of lakes for organic carbon cycling in the boreal zone. Glob. Change Biol. 10, 141147 (2004).

32. Spencer, R. G. M., Aiken, G. R., Wickland, K. P., Striegl, R. G. & Hernes, P. J. Seasonal and spatial variability in dissolved organic matter quantity and composition from the Yukon River basin, Alaska. Global Biogeochem. Cycles 22, GB4002 (2008).

33. Humborg, C. et al. CO2 supersaturation along the aquatic conduit in Swedish watersheds as constrained by terrestrial respiration, aquatic respiration and weathering. Glob. Change Biol. 16, 19661978 (2010).

34. Dubois, K. D., Lee, D. & Veizer, J. Isotopic constraints on alkalinity, dissolved organic carbon, and atmospheric carbon dioxide uxes in the Mississippi River.J. Geophys. Res. Biogeosci. 115, G02018 (2010).35. Berggren, M., Lapierre, J.-F. & del Giorgio, P. A. Magnitude and regulation of bacterioplankton respiratory quotient across freshwater environmental gradients. ISME J. 6, 984993 (2012).

36. Lauerwald, R., Hartmann, J., Moosdorf, N., Kempe, S. & Raymond, P. A. What controls the spatial patterns of the riverine carbonate system? A case study for North America. Chem. Geol. 337338, 114127 (2013).

37. Shipley, B. Cause and Correlation in Biology: A Users Guide to Path Analysis, Structural Equations and Causal Inference 2nd edn (Cambridge University Press, 2002).

38. Vachon, D., Prairie, Y. T. & Cole, J. J. The relationship between near-surface turbulence and gas transfer velocity in freshwater systems and its implications for oating chamber measurements of gas exchange. Limnol. Oceanogr. 55, 17231732 (2010).

39. McDonald, C. P., Stets, E. G., Striegl, R. G. & Butman, D. Inorganic carbon loading as a primary driver of dissolved carbon dioxideconcentrations in the lakes and reservoirs of the contiguous United States. Glob. Biogeochem. Cycle 27, 285295 (2013).

40. Maberly, S. C., Barker, P. A., Stott, A. W. & De Ville, M. M. Catchment productivity controls CO2 emissions from lakes. Nat. Clim. Change 3, 391394 (2012).

41. Wallin, M. B. et al. Evasion of CO2 from streams The dominant component of the carbon export through the aquatic conduit in a boreal landscape. Glob. Change Biol. 19, 785797 (2013).

42. Vachon, D. & Prairie, Y. T. The ecosystem size and shape dependence of gas transfer velocity versus wind speed relationships in lakes. Can. J. Fish. Aquat. Sci 70, 17571764 (2013).

43. Couture, S., Houle, D. & Gagnon, C. Increases of dissolved organic carbon in temperate and boreal lakes in Quebec, Canada. Environ. Sci. Pollut. Res. 19, 361371 (2012).

44. Ekschmitt, K., Liu, M. Q., Vetter, S., Fox, O. & Wolters, V. Strategies used by soil biota to overcome soil organic matter stability - why is dead organic matter left over in the soil? Geoderma 128, 167176 (2005).

45. Bianchi, T. S. et al. Enhanced transfer of terrestrially derived carbon to the atmosphere in a ooding event. Geophys. Res. Lett. 40, 116122 (2013).46. Raike, A., Kortelainen, P., Mattsson, T. & Thomas, D. N. 36-year trends in dissolved organic carbon export from Finnish rivers to the Baltic Sea. Sci. Total Environ. 435-436, 188201 (2012).

47. Davidson, E. A. & Janssens, I. A. Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature 440, 165173 (2006).

48. Dorrepaal, E. et al. Carbon respiration from subsurface peat accelerated by climate warming in the subarctic. Nature 460, 616679 (2009).

49. Gudasz, C. et al. Temperature-controlled organic carbon mineralization in lake sediments. Nature 466, 478473 (2010).

50. Cuthbert, I. D. & del Giorgio, P. Toward a strandard method of measuring color in fresh-water. Limnol. Oceanogr. 37, 13191326 (1992).

51. Cattaneo, A. & Prairie, Y. T. Temporal variability in the chemical characteristics along the Riviere De Lachigan - how many samples are necessary to describe stream chemistry. Can. J. Fish. Aquat. Sci. 52, 828835 (1995).

52. Marker, A., Nusch, E., Rai, H. & Riemann, B. The measurement of photosynthetic pigments in freshwaters and standardization of methods: conclusion and recommendation. Arch. Hydrobiol. Beih. Ergebn. Limnol. 14, 91106 (1980).

53. Guillemette, F. & del Giorgio, P. A. Reconstructing the various facets of dissolved organic carbon bioavailability in freshwater ecosystems. Limnol. Oceanogr. 56, 734748 (2011).

Acknowledgements

This work is part of the Carbon Biogeochemistry in Boreal Aquatic Systems (CarBBAS) research chair, co-funded by the Natural Science and Engineering Research Council of Canada (NSERC) and Hydro-Qubec. We thank A. Campeau, A. St-Pierre, A. Parkes and the CarBBAS team for the contribution to the eld and laboratory components of this study, and we also thank M. Pace, S. Ziegler, M. Jansson, A. Stubbins, C. Ruiz Gonzalez and A. Lindroth for comments on the early versions of the manuscript, as well as J. Comte, Y.T. Prairie, D. Vachon and N. Fortin St-Gelais for advice on analyses.

Author contributions

J.-F.L. and P.d.G. designed the study; all authors contributed to the manuscript structure;J.-F.L. and F.G. processed and analysed the data; J.-F.L. and P.d.G. wrote the manuscript; all authors commented on the manuscript.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Lapierre, J.-F. et al. Increases in terrestrially derived carbon stimulate organic carbon processing and CO2 emissions in boreal aquatic ecosystems.

Nat. Commun. 4:2972 doi: 10.1038/ncomms3972 (2013).

NATURE COMMUNICATIONS | 4:2972 | DOI: 10.1038/ncomms3972 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

& 2013 Macmillan Publishers Limited. All rights reserved.