The impact of disturbed peatlands on river outgassing in Southeast Asia

River outgassing has proven to be an integral part of the carbon cycle. In Southeast Asia, river outgassing quantities are uncertain due to lack of measured data. Here we investigate six rivers in Indonesia and Malaysia, during five expeditions. CO2 fluxes from Southeast Asian rivers amount to 66.9±15.7 Tg C per year, of which Indonesia releases 53.9±12.4 Tg C per year. Malaysian rivers emit 6.2±1.6 Tg C per year. These moderate values show that Southeast Asia is not the river outgassing hotspot as would be expected from the carbon-enriched peat soils. This is due to the relatively short residence time of dissolved organic carbon (DOC) in the river, as the peatlands, being the primary source of DOC, are located near the coast. Limitation of bacterial production, due to low pH, oxygen depletion or the refractory nature of DOC, potentially also contributes to moderate CO2 fluxes as this decelerates decomposition.

T he importance of inland waters in the global carbon cycle has gained more awareness since the last decade through studies which have revealed that inland waters (rivers, streams, lakes, reservoirs and estuaries) are not a passive conduit but play an integral role both for carbon storage and greenhouse gas emissions to the atmosphere [1][2][3][4] . These studies estimate, in line with the 5th Assessment IPCC Report 5 , that on a global scale B45-60% (0.9-1.4 Pg C per year) of carbon entering the freshwater system is decomposed and emitted back to the atmosphere as CO 2 . Another 0.2-0.6 Pg C per year is buried in freshwater sediments and about 0.9 Pg C 1 reaches the coastal ocean. However, an estimate of inland water outgassing by Raymond et al. 6 revealed an emission of 2.1 Pg C per year, of which 1.8 Pg C per year from streams and rivers, which is significantly larger than previous estimates. On the other hand, a more recent study by Lauerwald et al. 7 estimates a global river outgassing of 0.65 Pg C per year, however, they have excluded stream orders o2. These variable findings challenge our current understanding of the global carbon cycle and in particular that of the terrestrial biosphere as a sink for anthropogenic CO 2 . Still, large uncertainties remain in outgassing fluxes due to scarcity of data, which we aim to resolve for Southeast (SE) Asia. Indonesia and Malaysia are areas of particular interest due to their peatlands, which together store 66.5 Pg C (ref. 8). It has been shown recently that the fluvial organic carbon flux increases once these tropical peatlands are disturbed 9 . However, it remains unclear to what extent this influences the CO 2 emissions from these aquatic systems.
In this study, river outgassing fluxes are quantified for SE-Asia by using measurements from four rivers in Sumatra, Indonesia, and two rivers in Sarawak, Malaysia (Fig. 1). CO 2 fluxes and yields for all rivers are determined and related to peat coverage, which uncovers a positive relationship. Based on this correlation, CO 2 fluxes for SE-Asia are calculated, which reveal moderate fluxes and show that SE-Asia is not a hotspot for river outgassing.

Results
River processes and variability. In estuaries, salinities higher than zero can be measured 18-50 km downstream with increasing salinities towards the coastal ocean (Fig. 2a), whereas pCO 2 concentrations show an opposite pattern with increasing concentrations going upstream (Fig. 2c). Indeed, in the estuaries, salinity is inversely correlated with CO 2 concentrations (Fig. 2b) as a consequence of mixing between river and oceanic waters and the associated decrease of pH in rivers. To calculate mean river parameter values, estuaries were excluded to avoid the influence from ocean waters and only data points with salinity 0-0.1 in the respective rivers were taken into consideration (Fig. 2d, Table 1). Figure 2d shows the variability of the pCO 2 data in the Sumatran rivers in 2009 and 2013. The highest pCO 2 concentrations were measured in the Siak river, which also reveals the highest dissolved organic carbon (DOC) and the lowest O 2 concentrations (Table 1, Fig. 3). This agrees with results derived from numerical model and DOC decomposition experiments, which show that DOC leaching from peat soils and the decomposition of its labile fraction are the main factors controlling DOC and O 2 concentrations in the Siak 10 . The resulting input of dissolved inorganic carbon and the already low pH, caused by the organic acids from peat, shift the carbonate system (CO 2 3HCO À 3 3CO 2 À 3 ) towards CO 2 and explain the high CO 2 concentrations in the Siak. Increasing pCO 2 concentrations associated with increasing DOC concentrations denote that DOC leaching and decomposition are prime factors controlling pCO 2 concentrations also in all other studied rivers (Fig. 3).
Contrary to the Siak, which is a typical black-water river where its dark brown colour reduces light penetration to a few centimeters, non-black-water rivers have a higher light availability. Accordingly, photosynthesis plays a role and therefore these rivers have a pronounced day and night cycle, as seen in the Musi river during our expedition in 2013. As we entered this river during the day the pCO 2 concentrations were lower than during the departure, which took place during the night when photosynthesis and CO 2 uptake could not take place, but decomposition and CO 2 production prevailed. In total, the observed variability of pCO 2 in the Musi was ± 21.5%, which was reduced to ± 6.2% in the Siak due to a reduced impact of photosynthesis in peat-draining rivers.   Figure 2 | Patterns of salinity and pCO 2 underway measurements on the way in and out of the estuaries and rivers in 2009 and 2013. Estuaries start at 0 km, which is the border with the river at salinity zero, with increasing (negative) distance towards the coastal ocean. Measurements in the rivers start at 0 km at salinity o0.1 with increasing distance upstream. (a) Decrease of salinity across estuary towards river. (b) Inverse relationship of pCO 2 versus salinity (40.1). (c) Increase of pCO 2 across estuary towards river. (d) pCO 2 data points at zero salinity (0-0.1) against river distance, with number of data points in the lower right corner.  The Musi also shows an interannual variability with lower pCO 2 concentrations in 2009 as compared with 2013 due to lower precipitation rates and lower DOC leaching.
Although this supports our former finding that DOC leaching is a prime factor controlling pCO 2 concentrations, there are also exceptions, which point to other processes, as indicated by the high pCO 2 but low DOC concentrations in the Siak in 2013.
Contrary to non-peat-draining rivers, enhanced discharge lowers DOC concentrations in rivers draining undisturbed peat. This reverse behaviour is caused by the high water saturation of peat due to which enhanced precipitation leads to flooding and the resulting surface runoff increases discharge, which dilutes the DOC concentrations in the peat-draining rivers [11][12][13] . However, the peatlands in the Siak catchment as well as elsewhere in Indonesia are largely drained. When precipitation is relatively low, the ground water table could fall below the peat and the DOC concentrations in the river would be reduced 12 . This might explain the exceptionally low DOC concentration in April 2013, which can be associated with low precipitation rates (Table 1). However, despite the low DOC concentration, pCO 2 concentrations remained high pointing to an additional carbon source. This could have been a decaying plankton bloom favoured by the enhanced light availability due to the lower DOC concentrations. A similar situation might also explain the Lupar in 2013, where low DOC concentrations are associated with low precipitation rates, but relatively high pCO 2 concentrations still occur. This indicates that the dampening effects of reduced carbon input from soil leaching on the pCO 2 concentrations could be counterbalanced by input of CO 2 from presumably plankton blooms, which reduces the temporal variability of pCO 2 in the rivers during time characterized by extremely low DOC leaching.
River outgassing and DOC export. Quantification of both river outgassing and riverine DOC export sheds light on their flux ratio and consequently the share of river outgassing with respect to the carbon entering the freshwater system. CO 2 emissions were calculated for all investigated rivers using equation 1 (Methods section) with their respective piston velocities and respective DpCO 2 , which were based on their averaged pCO 2 values (Table 1) and an atmospheric CO 2 concentration of 390 matm. Subsequently, these fluxes were multiplied by their respective river surface area (Methods section and Table 2) to calculate the flux per river. Riverine DOC exports were calculated by multiplying the DOC concentrations with the discharges, which were based on the averaged monthly precipitations (Table 1) and assuming an evapotranspiration of 38% (ref. 9). The resultant net water-air CO 2 fluxes and riverine DOC exports are presented in Table 3. Note that both CO 2 and DOC fluxes are shown in grams 10 9 mth À 1 . These results suggest that on average 53.3±6.5% of the organic carbon leached from soils is decomposed in the river and emitted as CO 2 into the atmosphere. This percentage is based on monthly averages from six rivers, has a relatively small s.e. and is similar to the 45-60% that was reported as the global average in  *For the CO 2 calculations of Lupar and Saribas, only the catchment area upstream from the measurement point was considered, as the measurement data is representative for the upstream river characteristics and not the peatland area further downstream. The river coverage of the entire catchment area of the Lupar is 2.85% and that of the Saribas 3.64%. wRiver coverages for Malaysia, Indonesia and Southeast Asia are derived from the average pCO 2 , piston velocity and CO 2 yield found for these locations (Table 6).
the IPCC report, as mentioned before. Therefore we may assume that this percentage is also a representative value for river outgassing in SE-Asia.
CO 2 fluxes and peat coverage in SE-Asia. As the data of this study is limited to a peat coverage of 21.9%, the lateral DOC exports of 95.7 g C m À 2 per year for disturbed peat and 61.7 g C m À 2 per year for pristine peat, as determined in rivers on Borneo by Moore et al. 9 , were used to extend the data set to 100% peat coverage. Based on the fact that carbon leaching equals the sum of CO 2 outgassing and riverine DOC export, this DOC export yield for 100% peat coverage was then converted to CO 2 yield by assuming a river outgassing of 53.3 ± 6.5% as found in the SE-Asian rivers. This results in CO 2 yield for 100% peat coverage of 109.2 g C m À 2 per year (disturbed) and 70.5 g C m À 2 per year (pristine). Knowing that 90% of SE-Asian peatlands are disturbed and 10% pristine 14 , the CO 2 yield for 100% peat coverage amounts to 105.3 ± 27.6 g C m À 2 per year. This ratio for disturbed and pristine is naturally accounted for in the in situ river measurements, assuming an even distribution.
The lateral DOC export of 95.7 g C m À 2 per year and the net ecosystem C loss of disturbed peatlands of 433 g C m À 2 per year derived from eddy covariance measurements 15 together amount to 529 g C m À 2 per year. Including the river outgassing flux of 109.2 g C m À 2 per year, which accounts for a 20% increase, raises the net ecosystem C loss to a total of 638 g C m À 2 per year (net C source).
The extent of peat coverage in the investigated rivers differs in each catchment and correlates with the DOC and CO 2 yield, emphasizing the importance of peat as DOC and CO 2 source (Fig. 4). The correlation found between peat coverage and CO 2 yield shows that the correlation initially appears linear, but levels off after a peat coverage of 25%, which indicates that the rate of pCO 2 production declines with increasing peat coverage. This may be attributed to a limitation of bacterial production as a consequence of the low pH caused by the acidic organic environment 16 or oxygen depletion (this study, refs 10,16). Based on the peat coverages and respective CO 2 yields, the annual CO 2 fluxes were interpolated for Indonesia, Malaysia and SE-Asia (Table 4) using the regression. As a whole, SE-Asia releases 66.9 ± 15.7 Tg C per year through river outgassing, of which the majority is emitted by Indonesia with 53.9±12.4 Tg C per year. This is due to the fact that Indonesia holds 83% of SE-Asian peatlands in addition to a large land surface area. River outgassing in Malaysia amounts to 6.2 ± 1.6 Tg C per year. Although no CH 4 measurements were taken in the river, CH 4 fluxes are not considered to play a significant role in river carbon emissions. This is supported by CH 4 fluxes from the Saribas and Lupar estuaries, which emit 27 ± 24 t CH 4 -C per year and 84 ± 24 t CH 4 -C per year, respectively, and are a minute fraction of the estuary CO 2 fluxes 0.09±0.08 and 0.31±0.09 Tg C per year (ref. 17). CH 4 fluxes from peat soils are with 0.02 g C m À 2 per year also much lower than those of peat soil CO 2 fluxes 250 g C m À 2 per year (ref. 18). The error of the DpCO2 is the s.e., the Kco2 has the largest error from the mean, the discharge has the s.d. and the DOC flux, CO2 flux and outgassing errors are based on the best/worst case deviation, with the average having a s.e. Averaged and used values are in bold.  Considering that the data is collected during the transitional stages during the year with respect to precipitation, the mean values represent a suitable yearly average. Although extreme events such as droughts or heavy rainfall have not been measured, their influence will not alter the estimated CO 2 fluxes significantly. As discussed earlier, the effect of increased discharge due to high water saturated peatlands dilutes the DOC concentration, which dampens the effect of enhanced discharge 11,19,20 . Therefore a decrease is expected in DOC and hence pCO 2 concentrations during extreme events in peatland areas. In nonpeat areas, enhanced precipitation will increase DOC leaching and hence, DOC and pCO 2 concentrations in the river waters. However, as seen in Table 4, the impact of peat on the CO 2 fluxes, wherein the impact is defined as the share of CO 2 emissions from peatlands expressed in percentage, is much larger than that of non-peat areas. The regression indicates that non-peat areas have a CO 2 yield of 2.0 g C m À 2 per year and only contribute 7.1% (4.8 Tg per year) to the CO 2 fluxes in SE-Asia, as opposed to 92.8% (62.1 Tg per year) from peatlands (Table 3). Assuming a case of heavy rainfall, wherein the CO 2 yield would double to 4.0 g C m À 2 per year, would increase the CO 2 flux from non-peat areas to 9.6 Tg per year, but would not significantly increase the total annual CO 2 fluxes for SE-Asia. Additionally, the effect of the large peat impact and that CO 2 fluxes from peatlands will decrease due to dilution from increased discharge, will outweigh the increased CO 2 flux of non-peat areas, thereby maintaining a steady or even decreased CO 2 flux during extreme rainfall. Therefore, our annual CO 2 flux estimates can be assumed to be representative for SE-Asia, wherein fluctuations due to extreme events are accounted for in the error range.
Study comparison. For data comparison, the CO 2 flux and piston velocity estimates of Raymond et al. 6 were used, who have calculated global inland water-CO 2 effluxes by means of a global set of calculated piston velocities, pCO 2 values 21 and COSCAT areas (Coastal Segmentation and related CATchments). These COSCAT areas are characterized by their coastal segment limits and length and by catchment characteristics, such as runoff direction and physiographic units 22 . Raymond et al. 6 assigned an average piston velocity and an average pCO 2 value to each of these COSCATs from which they estimated the CO 2 yield per m 2 land per year. Several of these COSCAT areas, 10 of which are relevant for our study area, overlap SE-Asia and are summarized in Table 5. Based on the COSCATs in Table 5 the pCO 2 , K CO2 and CO 2 yield were averaged for Malaysia (row 1-3), Indonesia (row 2-8) and SE-Asia (row 1-10), after which the corresponding CO 2 flux was calculated and compared with our data (Table 6).
In addition, an indirect comparison was made based on recent estimates by Lauerwald et al. 7 (Table 6). Although no data is available for SE-Asia, they do provide estimates of pCO 2 concentrations and piston velocities for the tropical zone (o25°), which we used to predict their CO 2 fluxes in SE-Asia (Table 6), assuming a river surface coverage as determined in our study ( Table 2).
The comparison of Raymond et al. 6 with our data reveals that the CO 2 flux estimates for Indonesia and SE-Asia by Raymond et al. 6 are almost three times higher than those found in this study with 144.7 Tg C per year and 181.5 Tg C per year, respectively. The Malaysian estimate peaks almost eightfold with 48.8 Tg C per year. These overestimations were already anticipated by Raymond et al. 6 and are explained by the small number of calculated CO 2 values they had available in these COSCATs. Moreover, compared with our findings, their piston velocities appear to be largely overestimated (Table 5) and are mainly responsible for the resulting CO 2 flux overestimations of Raymond et al. 6 . As for Lauerwald et al. 7 , their data for the tropical zone (o25°) include piston velocities similar to ours, but show pCO 2 concentrations that are 21, 40 and 35% lower in Malaysia, Indonesia and SE-Asia, respectively, and result in lower CO 2 emissions.
Explanatory arguments for moderate fluxes. The overall statement of this study conveys that SE-Asia is in fact not such a river CO 2 outgassing hotspot as one could assume due to the The errors of the CO2 yield and flux represent the best/worst case deviation. *The Lupar and Saribas stations were located upstream from the peatland area; their data are therefore typical for a peat coverage of 0%. With the Lupar pCO2 measurements at 0% peat coverage, the correlation between annual CO2 yield and peat coverage was derived (Fig. 4b). The Saribas, having only one (seasonal) data point, was excluded from this correlation, but is shown in grey in Fig. 4b.   24), and show that CO 2 outgassing from SE-Asian rivers is rather moderate. The main reason for these moderate fluxes is the relatively short residence time 10 of DOC in the river water due to the location of the peatlands near the coast, which are the main source of DOC. Limitation of bacterial production, as a consequence of the low pH caused by the acidic organic environment 16 , oxygen depletion (this study, refs 10,16), or the refractory nature of DOC (refs 10,25,26), potentially also contributes to the moderate CO 2 fluxes as this decelerates decomposition, especially with increasing peat coverage.
In conclusion, this study shows that river outgassing fluxes in SE-Asia are in fact moderate and, in line with the 5th Assessment IPCC report, suggest that B53.3±6.5% of carbon entering the freshwater system is decomposed and emitted back to the atmosphere as CO 2 . Globally there are three tropical regions, of which Africa 16 and Amazonia 24 are significantly important contributors to the river outgassing budget. However, due to the fact that the main source of DOC is located near the coast, which further shortens the residence time of DOC in rivers, SE-Asia is a moderate emitter and global river outgassing estimates can therefore be scaled down. In future assessments, it needs to be considered that rivers function in a different way with respect to discharge, depending on residence time of DOC in the river and the soil type as it strongly affects the DOC leaching and consequently respiration and CO 2 emission.

Methods
Study area. SE-Asian peatlands cover 27.1 million hectares and are located in the coastal plains of the islands of Sumatera, Borneo and Ian Jaya 27 . Tropical peatlands are particularly vulnerable to anthropogenic stressors, such as deforestation and drainage, which are used to convert peat swamp forests into cropland. Today, a large fraction of this peat is already found under oil palm plantations 28 and only 10% remains undisturbed 14 . The investigated rivers in Sumatra, Indonesia, are the Musi, Batanghari, Indragiri and Siak. All four rivers originate from the Barisan Mountains with a short steep descent and continue to flow towards the Malacca Strait. Towards the river mouths, the rivers cut through peatland areas which scatter the low-lying areas along the northern coast and are subject to leaching of organic matter 29 . The Siak is a typical black-water river 30 ; its brown colour is derived from dissolved organic matter leached from adjacent disturbed peatlands, which cover 21.9% of its catchment area. The Indragiri, Batanghari and Musi have a peat coverage of 11.9, 5.0 and 3.5%, respectively 31 Fig. 1 (ref. 38).
Sampling methods. In Sumatra, pH, salinity, temperature, dissolved oxygen and pCO 2 were measured continuously by means of underway instruments. All sensors were arranged in a flow through system and supplied with surface water from an approximate depth of 1 m. Salinity was measured using a Seabird SBE 45 Micro TSG sensor. Temperature and pH were measured with a Meinsberg EGA 140 SMEK with integrated temperature sensor. Oxygen measurements were conducted with an Aanderaa Optode 3835. pCO 2 was measured with two devices: the Li-Cor 7,000 pCO 2 analyzer (October 2009 and October 2012) and the Contros HydroC CO 2 Flow Through sensor (October 2012 and April 2013). Before the expeditions both devices were calibrated, of which the Contros at 100, 448 and 800 p.p.m. The Li-Cor 7,000 analyzer was calibrated with certificated NOAA reference gases (#CB08923 with 359.83 p.p.m., #CA06265 with 1,021.94 p.p.m. and another certified calibration gas with 8,000 p.p.m. CO 2 ).
In addition to the continuous measurements, water samples were collected at each station using a Niskin bottle at an approximate depth of 1.5 m. DOC samples were filtered (0.45 mm), stored in 60 ml high-density polyethylene (HDPC) bottles and acidified with phosphoric acid (20%) to a pH-value of 2. After a total storage time during and after the expeditions of maximum 3 weeks, DOC samples were analysed upon arrival in the laboratory in Bremen, Germany, with a Dohrmann DC-190 Total Organic Carbon Analyzer. The samples were combusted at 680°C within a quartz column, filled with Al 2 O 3 -balls covered with platinum. The released CO 2 was purified, dried and measured by a non-dispersive infrared detection system. The relative s.d. for the method was ± 2%.
In Sarawak, pCO 2 was measured with an in situ FTIR analyzer 39 , using a Weiss equilibrator 40 . In the headwater region, which was not accessible by boat, a Li-820 CO 2 analyzer was used together with headspace equilibration in a 10-l water bottle (June 2013) and a 600-ml conical flask (March 2014). The FTIR and the Li-820 were calibrated with the same set of secondary standards, ranging from 380 to 10,000 p.p.m. CO 2 . Samples were taken for DOC following the same procedure as described above. Dissolved oxygen, pH and conductivity were measured using a WTW Multi 3420 with an FDO 925 oxygen sensor, a SenTix 940 IDS pH sensor and a TetraCon 925 conductivity sensor. Salinity was calculated from conductivity using the equations by Bennet 41 . Floating chamber measurements were performed to determine the CO 2 flux. The chamber had a volume of 8.7 l and a surface area of ARTICLE 0.05 m 2 . The CO 2 flux was determined from the increase of the CO 2 mixing ratio over time, which was monitored with the Li-820.
Calibration experiment. To calibrate the Contros river measurements, a CO 2 calibration experiment was conducted during which different concentrations of CO 2 gas were delivered using a gas mixing system. Before the calibration with water measurements, the gas concentrations delivered by the gas mixing system were controlled by the mixing system regulator, the Li-Cor 7000, the Li-820 and the cavity ring-down spectrometer (Picarro G2201-i) in a range from circa 500 to 6,000 p.p.m. (Supplementary Fig. 1a). The gas was then used to calibrate freshwater in a range of 1,500-5,500 p.p.m. that was pumped into the Li-Cor 7,000 equilibrator and the Contros sensor. The measured pCO 2 concentrations were correlated and the regression equation (Supplementary Fig. 1b) was used to calibrate the Contros river data measured during the expedition in 2013.
CO 2 flux calculation. CO 2 fluxes (F) were calculated from pCO 2 using: where K CO2 is the CO 2 piston velocity (cm h À 1 ), K 0 the solubility coefficient of CO 2 in seawater 42 and DpCO 2 (matm) is the sea-air pCO 2 difference. Usually, the piston velocity is poorly constrained and spatially and temporally highly variable. Raymond and Cole 43 have pointed out that flux estimates can easily be altered by considerable amounts depending on the choice of the piston velocity. Meanwhile, the ways to determine the piston velocity are multifarious. Relationships with wind speed, hydraulic characteristics 44 , in situ measurements using floating chambers 45,46 or dual tracer techniques 47 are commonly used.
Empirical models neglect small-scale fluctuations and allow for estimates on larger scales. These models relate environmental parameters to the piston velocity. In coastal systems and the ocean, this is mostly wind speed 48 . In rivers, stream velocity, stream slope, depth, discharge and bedrock roughness have been identified as the main drivers of in-stream turbulence and consequently the piston velocity 44 . In this study, floating chambers were used to determine piston velocities. The performance of floating chambers is a matter of debate, where both over-and underestimations occur due to artificial turbulence and shielding of wind, respectively, are suggested. On the other hand, sometimes, a relatively good agreement between floating chamber measurements and other techniques is reported 49 . The floating chamber measurements conducted on the Lupar and Saribas river are described in detail in Müller et al. 50 . There, we also discuss potential biases. The piston velocities of the Lupar and Saribas rivers were determined with nine floating chamber measurements and are averaged at 26.5 ± 9.3 cm h À 1 and 17.0 ± 13.6 cm h À 1 , respectively, after normalization to a CO 2 Schmidt number of 360 (30°C) 44 . For the Siak river, Rixen et al. 10 derived a similar piston velocity of 22.0 cm h À 1 using an oxygen balance model. Given that these rivers have piston velocities of a similar range, despite having different catchment sizes 51 and slopes, we believe that they provide a good representation of the piston velocities in the Indragiri, Batanghari and Musi, Therefore, their average of 21.8 ± 4.7 cm h À 1 is used to calculate the fluxes in these three rivers, whereas the river-specific piston velocities were applied accordingly in the other rivers.
Peat coverages. Catchment area was calculated using a relief model of the Earth's surface in ArcGIS 9.3 with the ArcHydro extension. Digital Elevation Data (DEM) such as the SRTM90mDEM can be obtained from the web site of the Consortium for Spatial Information (CGIAR-CSI) of the Consultative Group for International Agricultural Research (CGIAR, http://www.cgiar-csi.org) 52 . The SRTM90mDEM has a spatial resolution of 90 m at the equator and was originally produced in the framework of the NASA Shuttle Radar Topographic Mission (SRTM). Peat coverage (%) for each catchment was estimated using a combination of the FAO Soil Map of the World 31 and the catchment area derived from the ArcGIS DEM model.
The peat coverage for Malaysia, Indonesia and SE-Asia was derived from Hooijer et al. 27 , who based their peat coverage percentages for these areas on field surveys provided by Wetlands International and the FAO Digital Soil Map of the World 31 , as well as Miettinen et al. 28 , who used satellite images from as recent as 2010. However, data of Miettinen et al. 28 only covers Malaysia, Sumatra and Kalimantan, whereas that of Hooijer et al. 27 covers entire SE-Asia. Therefore, peat coverage for Malaysia, Indonesia and SE-Asia is based on the combination of Hooijer et al. 27 and Miettinen et al. 28 , wherein the most recent data was integrated where possible.
River surface coverage. For all six rivers, surface area estimation has been conducted based on the length and width of their primary course and main tributaries. Length was estimated using HydroSheds stream lines and Google Earth was used to estimate width along multiple sections in the primary course and main tributaries for each river. River area (%) was calculated based on the share of river surface area with respect to the catchment area. River coverages for Malaysia, Indonesia and SE-Asia are derived from the average pCO 2 , piston velocity and CO 2 yield found for these locations ( Table 6). The results are summarized in Table 2.
Uncertainty estimates. The errors 53 associated with the averaged parameter per expedition are presented as the s.d. The errors of the averaged parameter per river are calculated as the s.e. where possible; otherwise the error is presented as the deviation of the two averages from the mean. Throughout the calculations of the CO 2 yield and fluxes, the s.e. as derived from the averages were integrated. Therefore, the errors of the CO 2 yields and fluxes are representative for the best and worst case scenarios, named as the best/worst case deviation.