Articles | Volume 16, issue 2
Biogeosciences, 16, 505–519, 2019

Special issue: Human impacts on carbon fluxes in Asian river systems

Biogeosciences, 16, 505–519, 2019

Research article 30 Jan 2019

Research article | 30 Jan 2019

Export fluxes of dissolved inorganic carbon to the northern Indian Ocean from the Indian monsoonal rivers

Export fluxes of dissolved inorganic carbon to the northern Indian Ocean from the Indian monsoonal rivers
Moturi S. Krishna, Rongali Viswanadham, Mamidala H. K. Prasad, Vuravakonda R. Kumari, and Vedula V. S. S. Sarma Moturi S. Krishna et al.
  • CSIR – National Institute of Oceanography, Regional Centre, Visakhapatnam, 530017, India

Correspondence: Moturi S. Krishna (


Rivers are an important source of dissolved inorganic carbon (DIC) to the adjacent coastal waters. In order to examine the spatial variability in the distribution and major sources of DIC in the Indian monsoonal rivers and to quantify their export flux to the northern Indian Ocean, 27 major and medium-sized rivers were sampled during the discharge period. Significant spatial variability in concentrations of DIC (3.4–73.6 mg L−1) was observed, and it is attributed to spatial variations in the precipitation pattern, the size of rivers, pollution and lithology of the catchments. The stable isotopic composition of DIC (δ13CDIC) also showed strong spatial variability (−13.0 ‰ to −1.4 ‰) in the Indian monsoonal rivers with relatively depleted δ13CDIC values in rivers of the northwest of India (-11.1±2.3 ‰) and enriched values in the southeast of India (-3.5±2.3 ‰). Results of the linear least-squares regression models of Keeling and Miller–Tan's plots indicated that the chemical weathering of carbonate and silicate minerals by soil CO2 is the major source of DIC in the Indian monsoonal rivers. Spatial variability in the deviation of δ13CDIC from the approximated δ13C of the source may probably be due to dominant autotrophic production in rivers of the southeastern region, whereas heterotrophic decomposition of organic matter largely influences the other Indian monsoonal rivers. It is estimated that the Indian monsoonal rivers annually export ∼10.3 Tg of DIC to the northern Indian Ocean, of which the major fraction (75 %) enters into the Bay of Bengal, and the remaining fraction reaches to the Arabian Sea. This is consistent with the freshwater flux, which is 3 times higher for the Bay of Bengal (∼378 km3 yr−1) than for the Arabian Sea (122 km3 yr−1). Despite discharge from the Indian monsoonal rivers accounting for only 1.3 % of the global freshwater discharge, they disproportionately export 2.5 % of the total DIC exported by the world's major rivers. Despite rivers from the region in the southwest (SW) of India exporting DIC that is an order of magnitude lower (0.3 Tg yr−1) than the rivers from other regions of India, the highest yield of DIC was found in the rivers of the SW region of India. It is attributed to intense precipitation (∼3000 mm), favorable natural vegetation of tropical moist deciduous and tropical wet evergreen and semi-evergreen forests, tropical wet climate, high soil organic carbon, and the dominance of red loamy soils in catchments of the rivers of the SW region.

1 Introduction

Dissolved inorganic carbon (DIC) is one of the major constituent of carbon species in rivers. DIC in rivers mainly originates from the geogenic (weathering of carbonate and silicate rocks) and biogenic (decomposition of organic matter in soils) sources (Meybeck, 1987; Mook and Tan, 1991; Gaillardet et al., 1999, Dessert et al., 2001; Viers et al., 2007; Raymond et al., 2008; Tamooh et al., 2013). The former consumes atmospheric carbon dioxide (CO2), while the latter releases CO2 fixed by the terrestrial plants. In addition to these major sources in the catchment, DIC is also contributed by various physical and biological processes within the rivers. For instance, heterotrophic decomposition of organic matter, photo-oxidation of dissolved organic carbon (DOC), autotrophic respiration and dissolution of atmospheric CO2 contribute DIC to rivers. On the other hand, autotrophic production by aquatic plants (photosynthesis) and evasion of CO2 to atmosphere withdraw DIC from rivers. All these processes in the catchments and within the rivers are strongly coupled to atmospheric CO2, because they act as either sinks or sources of atmospheric CO2 (e.g. Berner et al., 1983; Mook and Tan, 1991; Gaillardet et al., 1999; Richey et al., 2002). The DIC in rivers and its export to the coastal oceans is thus intimately linked to the global carbon cycle (Campeau et al., 2017).

Riverine export fluxes of DIC to coastal regions of the world's oceans have been estimated on the global (Gaillardet et al., 1999) and regional scales (Richey et al., 2002; Wallin et al., 2013; Crawford et al., 2014; Campeau et al., 2014; Kokic et al., 2015) to understand the component of DIC in the global carbon budget. Annual export flux of DIC from the world's major river systems to the global ocean has been estimated as ∼327–385 Tg (1 Tg =1012 g; Ludwig et al., 1998; Meybeck and Vorosmarty, 1999). However, many of the regional studies on DIC export fluxes were limited only to the major river systems (e.g. Gaillardet et al., 1999; Joesoef et al., 2017), for example, the Mississippi River (Raymond and Cole, 2003; Raymond et al., 2008; Cai et al., 2008), Yangtze and Pearl rivers (Cai et al., 2008), Congo (Wang et al., 2013) River, etc. Regional studies on the riverine export fluxes of DIC are very important for the global carbon cycle and budget, as the export fluxes are largely dependent on the hydrological, lithological and environmental conditions, which are highly variable on the regional scales. However, DIC measurements are still lacking in several medium-sized rivers from different regions of the world in general, in particular from Asia.

Studies on the sources and export fluxes of DIC from the Indian rivers are very limited. Though DIC measurements were conducted in some Indian estuaries, for example, Mandovi and Zuari (Sarma et al., 2001), Godavari (Sarma et al., 2011), Cochin (Gupta et al., 2009; Bhavya et al., 2018), Hooghly (Mukhopadhyay et al., 2002; Samanta et al., 2015), Mahanadi (Pattanaik et al., 2017) and Chilka (Gupta et al., 2008; Muduli et al., 2013), they were confined only to the internal cycling of DIC and exchange of CO2 at the air–water interface but were not focused on the sources and export fluxes of DIC. The major sources of DIC in the Indian rivers remain unclear, except only a couple of rivers, Krishna (Das et al., 2005; Laskar et al., 2014) and Ganges (Samanta et al., 2015). Further, the quantity of annual DIC export by the Indian rivers to the coastal regions is unknown. Here, we made an attempt to understand the major sources of DIC in the Indian monsoonal rivers (Fig. 1), using δ13CDIC as a potential tracer, and to estimate the riverine export flux of DIC to the northern Indian Ocean from the Indian subcontinent.

Figure 1Map showing the study region. Rivers sampled in this study were indicated by solid black line. Distribution of soils in catchments of the Indian monsoonal rivers sampled was also shown. Rivers draining the four regions, i.e., northwest (NW), southwest (SW), southeast (SE) and northeast (NE), were shown by solid black arrows. The source is found at (last access: 12 December 2018).


Figure 2Schematic diagram showing the typical range of δ13C of different sources of dissolved inorganic carbon (DIC) in rivers. Various major processes influencing the δ13C of DIC (δ13CDIC) within the rivers were also shown. Black arrows indicate the direction of change in δ13CDIC due to the influences of different in-stream process, mentioned next to the arrows.


The stable isotopic composition of DIC (δ13CDIC) is widely used to identify the major sources of DIC in the aquatic systems (e.g. Singh et al., 2005; Tamooh et al., 2013; Samanta et al., 2015; Zou, 2016). The isotopic composition of DIC originated by dissolution of atmospheric CO2 is about 0 ‰ (Coplen et al., 2002), whereas it is about −27 ‰ to −26 ‰ if the DIC is derived from the oxidation of organic matter produced by C3 plants (O'Leary, 1988). The δ13C of DIC generated by the carbonic-acid (formed by soil CO2 dissolution) weathering of silicates is about −21 ‰ to −17 ‰ (Solomon and Cerling, 1987), while it is in the range of −10 ‰ to −9 ‰ for carbonate rocks, because half of the carbon comes from carbonate rocks (0 ‰, Land, 1980) during weathering. The weathering of silicate and carbonate minerals yield δ13CDIC in the range of −8 ‰ to −7 ‰ and −4 ‰ to −3 ‰, respectively, if the carbonic acid is formed by the dissolution of atmospheric CO2. Though the δ13C of DIC derived from different sources is well separable (Deines et al., 1974), the isotopic fractionation by in-stream physical and biological processes alters the δ13C of the DIC source (Fig. 2). For example, photosynthesis and degassing enriches (O'Leary, 1988; Finlay, 2004; Parker et al., 2005, 2010; Polsenaere and Abril, 2012; Venkiteswaran et al., 2014), while the heterotrophic decomposition of organic matter and photo-oxidation of dissolved organic carbon depletes the δ13C of DIC (Opsahl and Zepp, 2001; Finlay, 2003; Waldron et al., 2007; Vähätalo and Wetzel, 2008; Fig. 2). Though rivers are generally in disequilibrium with atmospheric CO2 (Raymond et al., 2013) and emit CO2 to atmosphere due to oversaturation (Oquist et al., 2009; Campeau et al., 2017), the isotopic equilibration between the DIC and CO2 in the atmosphere significantly influences the δ13CDIC in rivers (Abongwa and Atekwana, 2014; Deirmendjian and Abril, 2018) due to selective fluxes of 12CO2 and 13CO2 at the air–water interface. Hence, the influence of biogeochemical processes within the rivers must be considered while interpreting the δ13CDIC results for identification of DIC sources. The main objectives of this study are to (i) identify the major sources of DIC in the Indian monsoonal rivers, (ii) estimate the export flux and yield of DIC to the northern Indian Ocean, and (iii) examine the major processes in the catchments and within the rivers controlling DIC in the Indian monsoonal rivers.

2 Study region, sampling and methods

2.1 Study area

The Indian peninsula bifurcates the northern Indian Ocean into the Bay of Bengal and the Arabian Sea. Although these two basins occupy the same latitudinal belt, their oceanographic processes were reported to be remarkably different due to higher freshwater flux into the Bay of Bengal (1.63×1012 m3 yr−1) than into the Arabian Sea (0.3×1012 m3 yr−1; Subramanian, 1993; Gauns et al., 2005). The large freshwater influx leads to the formation of a strong vertical salinity stratification in the Bay of Bengal (Varkey et al., 1996) that prevents vertical mixing of nutrient rich subsurface water with the surface (Prasanna Kumar et al., 2004). As a result, the Bay of Bengal is considered to be relatively less productive (Prasanna Kumar et al., 2002) than the adjacent Arabian Sea, which is one of the highly productive zones in the world (Madhupratap et al., 1996; Barber et al., 2001; Smith, 2001) due to injection of nutrients into the surface through the seasonal upwelling and convective mixing (Shetye et al., 1994; Madhupratap et al., 1996; Muraleedharan and Prasanna Kumar, 1996).

Discharge from the Indian monsoonal rivers is largely fed by the monsoon-induced precipitation over the Indian subcontinent, which receives > 80 % of its annual rainfall during the southwest (SW) monsoon period (June–September; Soman and Kumar, 1990). Though some amount of rainfall occurs during the northeast (NE) monsoon (December–March), it does not generate discharge, as it will be stored within the dam reservoirs for domestic, industrial and irrigation purposes. Discharge from the Indian monsoonal rivers mainly occurs during the SW monsoon season (Vijith et al., 2009; Sridevi et al., 2015); hence, these rivers are called monsoonal rivers. Since the major portion of the annual freshwater discharge occurs only during the SW monsoon, the entire estuary is filled with freshwater (Vijith et al., 2009; Sridevi et al., 2015) during this period. As discharge is small during the rest of the year, the discharge during the SW monsoon (wet period) is considered to be equivalent to the annual discharge of the monsoonal rivers. Based on rainfall intensity, forest cover, vegetation and soil type in the catchment, rivers sampled in the present study were categorized into four groups, namely the northwest (NW), SW, southeast (SE) and NE rivers of India (Fig. 1). The SW region of India is characterized by the intense rainfall during SW monsoon (∼3000 mm) following the NE (1000–2500 mm), SE (300–500 mm) and NW (200–500 mm) regions of India (Soman and Kumar, 1990). The SW rivers drain red loamy soils, while the NW rivers drain black soils. Except the major rivers Godavari and Krishna, all the rivers reaching the Bay of Bengal (NE and SE rivers) drain red loamy and alluvial soils in their upper and lower catchments, respectively. The Godavari and Krishna rivers drain black soils in their upper catchment, whereas they drain red loamy and alluvial soils in their middle and lower catchments, respectively (Geological Survey of India;, last access: 5 November 2018). Based on discharge, the monsoonal rivers in this study were divided into two types, namely the major (> 150 m3 s−1) and medium-sized (< 150 m3 s−1) rivers.

2.2 Sample collection

Water samples were collected from the freshwater regions of the estuaries to obtain reliable export fluxes of DIC to the coastal ocean. Samples were collected at two to three locations to minimize the spatial variability within the freshwater zone of the estuary. Further, to minimize the interannual variability in DIC concentrations, sampling was conducted in 2 different years, and the mean was used for export flux estimations. Further, samples were collected in midstream in the river, using a local mechanized boat to avoid the contamination from river banks.

In situ measurements and sample collection were conducted in 27 rivers of the Indian subcontinent (Fig. 1) during the SW monsoon season of the years 2011 and 2014. Surface water samples at each location were collected for phytoplankton biomass (Chl a), DIC and dissolved oxygen (DO). Samples for DIC were collected in airtight crimp-top glass bottles, and poison (mercuric chloride) was added to arrest the biological activity. DO analysis was carried out at a temporary shore laboratory setup for sample processing after the completion of sampling on each day. Water samples were filtered through GFF (nominal pore size of 0.7 µm) under moderate vacuum and were stored in liquid nitrogen for Chl a analysis.

2.3 Methods

Temperature and salinity at the sampling locations were measured using a conductivity–temperature–density (CTD) profiling system (Sea Bird Electronics, SBE 19plus, United States of America). Concentration of DO was determined by Winkler's method (Carritt and Carpenter, 1966) using an auto titrator (Metrohm, Switzerland) with potentiometric endpoint detection. The analytical precision of the method was ±0.07 % (RSD). Dissolved oxygen saturation is computed following formulations given by Garcia and Gordon (1992). DIC concentrations in water samples were measured at our institute's laboratory using a coulometer (UIC Inc., USA) connected to an automatic subsampling system. Based on the repeated analysis of samples and standards, the precision of the method was ±0.02 mg L−1. The certified reference materials (CRMs) supplied by Andrew G. Dickson from Scripps Institute of Oceanography, USA, and internal standards were used to test the accuracy of our DIC measurements; they were found to be within ±0.2 % to 0.3 %. The potentiometric Gran titration method (Metrohm, Switzerland) was used for determination of pH and total alkalinity and followed the standard operating procedures given by the Department of Energy (DOE; 1998).

The stable carbon isotopic composition of DIC in the water was measured on a gas-bench–isotope-ratio mass spectrometer (EA-IRMS-Delta V, Finnigan, Germany). 50 mL airtight bottles with rubble septa were filled with 0.5 mL of high-purity orthophosphoric acid and purged with high-purity helium. About 1 mL of a water sample is injected into the bottle and incubated at constant temperature of 50 C for 12 h. The CO2 extracted into the head space is injected into the IRMS through gas bench. The results are expressed relative to conventional standards, that is, with Pee Dee Belemnite (PDB) limestone for carbon (Coplen, 1996) as δ values, defined as


where R refers to 13C, and X stands for 13C12C. The high-purity tank of CO2 was used as a working standard for carbon. These gases were calibrated with standards of the International Atomic Energy Agency (IAEA). Standard deviation on 20 aliquots of the same sample was lower than 0.05 ‰ for δ13C. Chlorophyll a (Chl a) on the filter was extracted into dimethylformamide (DMF) and measured the extract fluorometrically using a spectrofluorophotometer (Varian Eclipse, Varian Electronics, UK) following Suzuki and Ishimaru (1990). Annual mean discharge data of the rivers were taken from Meybeck and Ragu (1995, 1996), the Central Water Commission of New Delhi (2006, 2012), and Kumar et al. (2005). Catchment area of the rivers was obtained from Water Resources Information System of India (WRIS,, last access: 10 October 2018). Soil organic carbon data were taken from Kishwan et al. (2009) and Sreenivas et al. (2016), and the rainfall data were obtained from Soman and Kumar (1990). Dissolved organic carbon (DOC) data for the Indian rivers were taken from Krishna et al. (2015)

Total export flux of DIC from each river was estimated by multiplying the mean concentrations of DIC at near-zero salinity (river endmember) with the annual discharge. Spatial variability of DIC concentrations within the river was minimized to a large extent by collecting samples from two to three locations in each river, while the interannual variability was minimized by collecting samples during discharge periods of 2 years. However, variability in DIC concentrations within the discharge period results in some uncertainties in our estimations of DIC export fluxes. Time series measurements in the Godavari estuary (our unpublished results) revealed that the variability in DIC concentrations within the discharge period is up to 10 %. Therefore, the error associated with our DIC flux estimates may be about 10 %. The DIC flux normalized by catchment area (yield) was calculated by dividing the total DIC export flux of the river by its catchment area.

Figure 3Spatial variability in the concentration (mg L−1; a), export flux (Tg yr−1; b) and yield (g m−2 yr−1; c) of dissolved inorganic carbon (DIC) and its stable isotopes (δ13CDIC; d) in the Indian monsoonal rivers studied. Rivers geographically located in the northwest (NW), southwest (SW), southeast (SE) and northeast (NE) regions of India were also shown. Rivers draining into the Bay of Bengal (east-flowing rivers) were shown with gray shading, while rivers draining into the Arabian Sea (west-flowing) were shown with no shading.


3 Results

3.1 Hydrographic characteristics

Surface water temperatures were higher in rivers from the NE and SE regions (mean 30.9±1.2C) than the rivers from the SW and NW regions (27.3±1.5C) of India. Dissolved oxygen saturation varied, from as low as 63 % to as high as 105 %, with a mean saturation of 90±11 %. The rivers from the SW region of India recorded more undersaturation of DO (82±7 %) than the rivers located in the NE (89±15 %), NW (93±3 %) and SE (96±11 %) regions of India. Chlorophyll a (Chl a) concentrations varied broadly, from 0.8 to 7.5 mg m−3, with relatively higher mean concentrations in rivers of the SE region (4.7±2.5 mg m−3) followed by the SW (2.8±0.7 mg m−3) regions of India. On the other hand, relatively low Chl a was observed in the medium-sized (2.6±1.3 mg m−3) when compared to the major estuaries (3.2±2.1 mg m−3).

3.2 DIC concentrations and δ13CDIC

DIC concentrations in the Indian monsoonal rivers widely varied, from 3.4 (Bharathapuzha) to 73.6 mg L−1 (Vellar), with a significant spatial variability (Fig. 3a; Table 1). The highest mean DIC concentration was observed in rivers of the SE region (37.4±6.3 mg L−1), while the lowest DIC was found in the SW region (5.2±2.1 mg L−1) of India. Intermediate values were found in rivers of the NW (28.4±8.9 mg L−1) and NE (17.1±6.2 mg L−1) regions of India. DIC concentrations were found to be similar in the major (22.7±13.6 mg L−1) and medium (21.1±13.2 mg L−1) rivers (homoscedastic Student's t test; p=0.76). The mean DIC concentration found in this study (21.4±16.3 mg L−1) is similar to those observed earlier in the major river systems of India (Brahmaputra; Singh et al., 2005) and elsewhere in the world, for example, in British rivers (Jarvie et al., 2017) and Swedish rivers (Campeau et al., 2017). However, DIC concentrations in the present study are higher than the global mean DIC (10.3 mg L−1, Meybeck and Vorosmarty, 1999; Table 1) but are lower than those reported in the rivers draining into the Gulf of Trieste (northern Adriatic; 37–66 mg L−1; Tamse et al., 2014).

The δ13CDIC varied from −13.0 to −1.4 ‰, with a significant spatial variability (Fig. 3d; Table 1) in the rivers sampled. Relatively depleted δ13CDIC values were observed in rivers of the NW region (-11.1±2.3 ‰), while enriched δ13CDIC was found in rivers of the SE region (-3.5±2.3 ‰) of India (Fig. 3d). The δ13CDIC values found in this study are well within the range of values reported earlier in rivers of India (Das et al., 2005) and elsewhere in the world, for example, in Swedish streams (−27.6 ‰ to −0.6 ‰; Campeau et al., 2017) and rivers from Italy and Slovenia (−12.8 ‰ to −7.7 ‰, Tamse et al., 2014).

Table 1Catchment area, discharge, length and elevation of each river, and annual mean rainfall, number of dams and population density in each watershed of the Indian monsoonal rivers sampled. Concentrations, export fluxes and yields of dissolved inorganic carbon (DIC) and its stable isotopes (δ13CDIC) of Indian rivers were given. Measured pH and calculated δ13C values of CO2 from isotopic fractionation factors across DIC speciation were also provided. Rivers located in the northern region (north of 16 N) of India were shown in italics. Of these rivers, Mahisagar, Sabarmati, Tapti and Narmada are located in the northwestern (NW) region, whereas rivers from Godavari to Hyadri are located in the northeastern (NE) region of India.

Download Print Version | Download XLSX

3.3 Export fluxes and yield of DIC

Annual export flux of DIC to the coastal ocean from the individual rivers varied broadly from 0.01 Tg (Chalakundi) to as high as 2.33 Tg (Krishna; Fig. 3b; Table 1). Among the rivers sampled, rivers of the NE region of India export higher DIC (6.52 Tg yr−1) while the lowest was found from rivers of the SW region (0.24 Tg yr−1; Table 1). The Indian monsoonal rivers together export about 10.32 Tg yr−1 of DIC to the northern Indian Ocean, of which 7.81 Tg (75 %) enters into the Bay of Bengal, and the remaining enters into the Arabian Sea (2.51 Tg). The yield of DIC ranged from 2.8 (Bharathapuzha) to 20.7 g m−2 yr−1 (Baitarani; 3c; Table 1), excluding the exceptionally high yield of 119 g m−2 yr−1 from the Haldia River. The mean yield was found to be more or less similar in rivers from all the four regions of India, i.e, the NW (8.4 g m−2 yr−1), SW (8.8 g m−2 yr−1), SE (6.6 g m−2 yr−1) and NE (7.7 g m−2 yr−1) regions. Despite the export flux of DIC being lowest from rivers of the SW region (0.24 Tg yr−1), interestingly, the yield from rivers of this region is on par with (even slightly higher than) the other Indian monsoonal rivers (Table 1; Fig. 3b, c). Yields of DIC found in this study are similar to those found earlier in rivers elsewhere in the world (Huang et al., 2012).

4 Discussion

4.1 Distribution of DIC in the Indian monsoonal rivers

Distribution of DIC in the Indian monsoonal rivers showed large spatial variability, with the lowest values in rivers from the SW region of India (Fig. 3a). DIC concentrations in rivers are known to be influenced by the intensity of precipitation over the catchment, basin lithology (Giesler et al., 2013; Löfgren et al., 2014), length of the fluvial network (Hotchkiss et al., 2015), and in-stream physical and biological processes (Mook and Tan, 1991; Raymond et al., 2008). The spatial distribution of rainfall over the Indian subcontinent (, last access: 8 October 2018) shows that the SW region receives the highest annual rainfall (∼3000 mm) when compared to the rest of India (Soman and Kumar, 1990).

Figure 4(a) Exponential decrease and (b) linear increase of dissolved inorganic carbon (DIC) concentrations with increasing rainfall over the catchment and the length of the river, respectively. (c) Inverse and (d) linear relationships of chlorophyll a (Chl a) with concentrations and δ13C of DIC, respectively. Significant linear relationships of δ13C of DIC with (e) dissolved oxygen (DO) saturation and (f) dissolved organic carbon (DOC) concentration. Linear relationships observed between (g) DIC concentrations and population density in the catchment, and (h) total alkalinity and δ13C of DIC in the Indian monsoonal rivers during the study period. Ovals with dashed lines indicate the outliers which were not included in the regression equations. Rivers of the northwestern region of India showed linear relationships, as shown by the other Indian rivers but with a different slope (f–h).


The intense precipitation over the SW region is expected to cause higher weathering rates and thus higher DIC in rivers (e.g., Gupta et al., 2011), but lower DIC concentrations were found in rivers of this region. It could be due to the influence of dilution, because the dense precipitation over the small catchment area (Table 1) might have diluted DIC concentrations in rivers of this region. In order to understand the influence of the density of rainfall on DIC in rivers, we normalized the volume of discharge from the river with its catchment area. The catchment area-normalized volume of discharge was found to be much higher in rivers from the SW region (1.71 m3 m−2) than the rivers from the SE (0.17 m3 m−2), NE (0.6 m3 m−2) and NW (0.32 m3 m−2) regions of India. Higher catchment area-normalized discharge (about 3 times) might have diluted DIC concentrations in the rivers of the former region. A strong exponential decrease in DIC concentrations with increasing rainfall over the catchment (r2=0.72, p<0.001; Fig. 4a) also suggests that DIC concentration in the Indian rivers is strongly influenced by density of precipitation over the catchment. Rivers of the SW region are relatively small in size, both in terms of catchment area (total catchment area: 20 km ×103 km) and the length of the river (mean length: 126 km), when compared to the rivers from other regions (SE, NE and NW) of India (Table 1). Since the contribution of DIC from in-stream processes, such as decomposition of organic matter, has been demonstrated to increase along the course of the fluvial network (Hotchkiss et al., 2015), possibly due to an increase in the residence time of water (Catalan et al., 2016), the lowest DIC concentrations found in rivers from the SW region may also, at least partly, be due to their small size. Fairly good positive correlation between DIC concentrations and length of the rivers (r2=0.38, p<0.01; Fig. 4b) also supports this argument.

The major physical and biological processes controlling DIC concentrations in rivers are the exchange of CO2 with the atmosphere, autotrophic removal and heterotrophic addition of DIC. Since the Indian monsoonal estuaries have been reported to be a source of CO2 to the atmosphere during the discharge period due to heterotrophic decomposition of organic matter (Sarma et al., 2001, 2011, 2012; Gupta et al., 2008, 2009; Bhavya et al., 2018), the DIC input from the dissolution of atmospheric CO2 may be unlikely. On the other hand, organic matter decomposition is expected to add a significant amount of DIC, as enhanced bacterial respiration rates were reported during this period (Sarma et al., 2011, 2012). In contrast, significant negative correlation between chlorophyll a and DIC (r2=-0.44, p<0.01; Fig. 4c), except few SE rivers where elevated phytoplankton biomass (Chl a > 5 mg m−3) was recorded, suggests that autotrophic removal of DIC is also significant in the Indian monsoonal rivers during the study period. A significant positive relationship was observed between the δ13CDIC and Chl a (r2=0.49; p<0.01; Fig. 4d), supporting this argument, because preferential uptake of 12C over 13C during photosynthesis leaves the residual DIC enriched in 13C. On the other hand, δ13CDIC showed significant positive correlation with DO saturation (r2=0.50, p<0.01; Fig. 4e; depleted δ13CDIC values with more undersaturation of DO) and DOC concentrations (r2=0.43, p<0.01; Fig. 4f), as was observed in the Xi river (Zou et al., 2016). Altogether, enriched δ13CDIC values are associated with higher DOC, less undersaturation of DO and higher phytoplankton biomass (Chl a), while the depleted δ13CDIC values are associated with more undersaturation of DO and less DOC. This suggests that both autotrophic removal and heterotrophic addition control DIC in the Indian rivers during the discharge period, with a considerable spatial variability. However, the influence of these processes on DIC concentrations is difficult to separate with this bulk δ13CDIC data set, as the δ13CDIC in rivers is also influenced by pollution, catchment lithology and outgassing of CO2 (Brunet et al., 2005; Bouillon et al., 2009; Shin et al., 2011; Zeng et al., 2011; Tamooh et al., 2013). Excluding Sabarmati and Mahisagar rivers, DIC concentrations showed a fairly good linear relationship with population density over the catchment of the river (r2=0.41, p<0.01; Fig. 4g), suggesting a considerable influence of pollution from the megacities and industries on DIC in the Indian rivers.

Figure 5Linear least-squares regression models of (a) δ13CDIC as a function of 1/DIC (Keeling plot) and (b) δ13CDICxDIC as a function of DIC concentrations (Miller–Tans plot) in the Indian monsoonal rivers. (c) Linear relationship between calculated δ13C of CO2 and the measured δ13CDIC values, and (d) Miller–Tans linear regression model of δ13C-CO2×CO2 as a function of CO2 concentration in the Indian monsoonal rivers.


Spatial distribution of soils shows that rivers of the NW region of India and upper reaches of Krishna and Godavari rivers drain the lime-rich black soils (Fig. 1), while rivers from the SW region drain red loamy soils, whereas the east-flowing rivers drain the lime-poor red sandy soils in the upper but lime-rich alluvial soils in the lower reaches (Fig. 1). Lateritic soils, which are poor in lime and silicate, occupied the catchment of the rivers in the SW region of India. Relatively lower chemical weathering rates of the lateritic than the non-lateritic soils could be one of the reasons for the observed lower DIC concentration in the rivers from the SW region of India. A significant positive correlation was found between total alkalinity (TA) and δ13CDIC (r2=0.52; p<0.01; Fig. 4h), suggesting that a significant contribution of DIC is from weathering of carbonate minerals in the catchment. Though the higher chemical weathering rates were reported for the Deccan Trap basalts (Das et al., 2005; Singh et al., 2005), which occupied the catchments of rivers of the NW region of India and upper reaches of Godavari and Krishna, higher DIC concentrations were also observed in rivers draining over the metamorphic rocks. This suggests that the influences of factors other than bedrock are also significant in the concentrations of DIC in the Indian rivers.

4.2 Major sources of DIC in the Indian monsoonal rivers

Though the δ13CDIC is a promising tool for deciphering the sources of DIC, its interpretation for source material identification in rivers is still challenging, because multiple physical and biological processes within the rivers significantly alter the δ13C of DIC source. The influence of major in-stream processes on the δ13CDIC must be separated before interpreting the results for major sources of DIC, and failing to do this leads to erroneous conclusions. In order to identify and separate DIC sources, here we used two different graphical mixing model techniques, the Keeling (Keeling, 1958; Pataki et al., 2003) and Miller–Tans plots (Miller and Tans, 2003). These models approximate the hypothetical δ13C of source material as an intercept (in Keeling plot) and slope (in Miller–Tans plot) of the linear least-squares regression equations (Pataki et al., 2003; Campeau et al., 2017). The deviations from the approximated δ13C of the source can be interpreted as the influence of the in-stream processes. Further, we approximated the δ13C of CO2 using a set of enrichment factors of isotopic fractionation across the carbonate species (Zhang et al., 1995) in order to filter the impact of DIC speciation and pH on the bulk δ13CDIC values. This approach has already been used by Quay et al. (1992), Mayorga et al. (2005) and recently by Campeau et al. (2017).

Significant negative relationships were observed in both the Keeling plot (δ13CDIC as a function of 1/DIC; Fig. 5a) and the Miller–Tans plot (δ13CDIC× DIC as a function of DIC; Fig. 5b; r2=0.61, p<0.01, and r2=0.72, p<0.01, respectively) of DIC in the Indian rivers, except in the rivers draining the Deccan Trap basalts. Both graphical mixing models, the Keeling and Miller–Tan's plot, approximated the similar δ13C of source material (−3.0 ‰ and −2.0 ‰, respectively; Fig. 5a, b), suggesting that weathering of carbonate minerals is the predominant source of DIC in the Indian monsoonal rivers rather than biogenic soil CO2. Calculated δ13C of CO2 ranged from −21.5 ‰ to −9.6 ‰ in the Indian rivers, with a mean value of -13.0±2.7 ‰. Calculated δ13C of CO2 is linearly correlated with the measured δ13CDIC, but correlation coefficient (r2) is only 0.51 (Fig. 5c), suggesting significant spatial variability in the influence of in-stream processes on the δ13CDIC. The Miller–Tans plot of CO2 (δ13C-CO2×CO2 as a function of CO2) showed a highly significant linear regression model with a slope of −10.7 ‰ (r2=0.97; p<0.001; Fig. 5d). These results indicated that chemical weathering of carbonate and silicate minerals by soil CO2 (−10 ‰ to −9 ‰) is the major source of DIC in the Indian rivers. Deviations of the measured δ13CDIC (−13.0 ‰ to −1.4 ‰) from that of the approximated δ13C of DIC source (−3.0 ‰ to −2.0 ‰) and δ13C of CO2 (−10.7 ‰) could be due to the influence of in-stream process. In more than 75 % of the Indian rivers sampled, the deviation from the δ13C of the DIC source is towards the negative side (depletion; δ13CDIC<−3.0 ‰), suggesting that heterotrophic decomposition of organic matter is the dominant process controlling DIC in these rivers. While a lack of (or very little) deviation was observed only in rivers from the SE region of India (mean δ13CDIC: −3.1 ‰), this could be due to the competition between autotrophy, degassing and heterotrophy, as these processes influence the δ13CDIC in opposite directions (Fig. 2); the former two processes cause enrichment, while the latter depletes δ13CDIC. Relatively higher phytoplankton biomass (mean Chl a: 4.6 mg m−3) and less undersaturation of DO (98.7 %) was observed in these rivers compared to the mean of the rest of the Indian rivers (2.4 mg m−3 and 87.5 %, respectively), suggesting that autotrophy is one of the dominant processes controlling DIC in rivers from the SE region of India. The total number of dams on the rivers from this (SE) region (mean 155, Table 1) is not significantly higher from that of the mean of total number of dams on the Indian rivers sampled (mean 135), suggesting that degassing due to storage of water may not be the dominant process responsible for enrichment in δ13CDIC values of these rivers.

4.3 Total DIC export by the Indian monsoonal rivers to the northern Indian Ocean

Indian monsoonal rivers annually export ∼10.3 Tg of DIC to the northern Indian Ocean. Nearly three-fourths of this amount (7.8 Tg) reaches to the Bay of Bengal, while the Arabian Sea receives only one-fourth (2.5 Tg). This is consistent with the higher magnitude of freshwater discharge to the Bay of Bengal (378 km3 yr−1) from the catchment area of about 970 km ×103 km than the Arabian Sea (122 km3 yr−1 from the catchment area of 244 km ×103 km). The total DIC exported by the Indian monsoonal rivers (10.3 Tg yr−1) is lower than the DIC exported by the American (61.4 Tg yr−1) and African (17.7 Tg yr−1) rivers and major rivers draining to the tropical Atlantic from South America and Africa (53 Tg yr−1; Araujo et al., 2014). It is mainly due to the fact that freshwater discharge from the Indian monsoonal rivers is very low (∼500 km3 yr−1) compared to the American (11 799 km3 yr−1) and African (3786 km3yr−1) rivers. However, the Indian monsoonal rivers are exporting disproportionately higher DIC, because they account for only 1.3 % of the global river discharge but export 2.5 % of the global riverine DIC to the oceans (400 Tg yr−1). Though American and African rivers account for 30 % and 10 % of the global river discharge, they export only 15 % and 4.4 % of global riverine DIC to oceans, respectively. A disproportionately higher DIC flux from the Indian rivers could be due to relatively higher weathering rates of silicate and carbonate minerals in their drainage basins (Das et al., 2005; Gurumurty et al., 2012; Pattanaik et al., 2013). Higher DIC fluxes from the tropical regions are mainly attributed to the favorable climatic conditions, lithology and land use change (Huang et al., 2012) of this region for higher dissolution.

Figure 6Significant relationships of dissolved inorganic carbon (DIC) yield with (a) river discharge in medium estuaries, (b) catchment areas of the medium-sized rivers, (c) catchment areas of the major rivers, (d) rainfall over the catchment of all the rivers sampled and (e) soil organic carbon (OC) content in catchments of the Indian monsoonal rivers studied. Since the data on soil OC are not available for each watershed, (e) was plotted using the available soil OC data on regional scale (NW, SW, SE and NE regions of India). Hence, it contains only four points.


Krishna et al. (2015) reported that Indian monsoonal rivers export 2.32 Tg yr−1 of dissolved organic carbon (DOC) to the northern Indian Ocean. The total fluvial dissolved carbon flux (DIC + DOC) would be 12.6 Tg yr−1, to which DIC flux contributed up to ∼81 %. The predominance of DIC has also been found in rivers elsewhere in the world, for example, in British rivers (Jarvie et al., 2017) and high-altitude Swedish rivers (Campeau et al., 2017). Since the catchment area of the Indian monsoonal rivers ranged widely, from as low as 1 km ×103 km to as high as 313 km ×103 km, the export fluxes of DIC were normalized with the catchment area of the river (yield) in order to examine various factors controlling the DIC export to the northern Indian Ocean.

4.4 Yield of DIC from the Indian monsoonal rivers

The yield of DIC found in this study (mean 8.7±5.2 g m−2 yr−1) is close to those found in rivers from the tropical region of Asia but is significantly higher than those reported from tropical region of the American and African continents (Huang et al., 2012). The yield was highest (8.8±5.6 g m−2 yr−1) in rivers from the SW region of India, despite the fact that they export relatively lower DIC (0.3 Tg yr−1) due to their low volume of discharge (46 km3 yr−1) and relatively smaller catchment (20 km ×103 km) than the rivers from SE, NE and NW regions of India (Table 1). DIC yield showed a significant positive correlation with the volume of discharge (r2=0.67, p<0.001; Fig. 6a) in medium-sized rivers, and no such relationship was found in the major rivers. Significant negative relationship was observed between DIC yield and catchment area of river (r2=-0.49, p<0.001 – Fig. 6b; r2=-0.43, p<0.001 – Fig. 6c; for medium-sized and major rivers, respectively), suggesting that the smaller rivers export more DIC per unit area of catchment compared to the major river systems, and thus inclusion of DIC data from medium-sized rivers in the world significantly alters the global estimations of DIC. A fairly good linear relationship between the yield of DIC and the intensity of precipitation (r2=0.43, p<0.01 – Fig. 6d) was observed only in the rivers which receives > 2000 mm of annual mean precipitation. Higher precipitation over the catchment increases the yield of DIC, because the dense precipitation enhances the extraction of DIC from soils and rocks in their catchment. Therefore, high precipitation (∼3000 mm) over the small catchment (20 km ×103 km) could have increased DIC yield from the rivers of the SW region of India.

Sreenivas et al. (2016) and Krishwan et al. (2009) found that the soil organic and inorganic carbon contents in the surface (100 cm) soils in the catchment of rivers in the SW region were higher and lower, respectively, than the catchments of the rivers from the SE, SW and NE regions of India. Decomposition of soil organic matter releases excess CO2 that increases CO2 in soils, leading to the formation of acidic conditions in soils. This would increase the DIC yield by more dissolution of soil carbonates and chemical weathering of carbonate and silicate rocks (Zou et al., 2016). A significant linear correlation was found between soil organic carbon content and DIC yield in this study (r2=0.65, p<0.001; Fig. 6e), suggesting that higher soil organic carbon in the catchment of the rivers from the SW region could have elevated the yield of DIC from rivers of this region. The basin-scale studies are, however, required for a comprehensive understanding of the influence of environmental and anthropogenic factors on export fluxes and yield of DIC from the Indian monsoonal rivers.

5 Summary

In order to examine the spatial variability in the sources and distribution of dissolved inorganic carbon (DIC) in the Indian monsoonal rivers, and to estimate their export fluxes of DIC to the northern Indian Ocean, we sampled a total of 27 major and medium-sized rivers during the wet period. An order of magnitude variability was found in DIC concentrations among the rivers sampled (3.4–73.6 mg L−1), with a lower mean concentration of 6.6±2.1 mg L−1 in rivers located in the SW region of India. It is attributed to significant spatial variability in the size of rivers, precipitation pattern, pollution and lithology in their catchments. The approximated δ13C of DIC source from the Keeling and Miller–Tans plots (−2.0 ‰ and −3.0 ‰, respectively), and the calculated δ13C of CO2 suggested that DIC in the Indian rivers mainly originates from chemical weathering of carbonate minerals but is largely affected by autotrophic production in rivers from the southeastern region of India and heterotrophic decomposition of organic matter in rivers from other regions of India. Indian monsoonal rivers together export ∼10.3 Tg yr−1 of DIC to the northern Indian Ocean, of which 7.8 Tg yr−1 enters into to the Bay of Bengal, while the Arabian Sea receives only 2.5 Tg yr−1. It is mainly attributed to the volume of river discharge, as the former receives ∼378 km3 yr−1 while the latter receives only 122 km3 yr−1 from the Indian monsoonal rivers. Dense rainfall and higher soil organic carbon content in the catchment of rivers from the SW region than in the catchment of the other Indian rivers resulted in the highest yield of DIC from the former than the latter.

Data availability

The data set used in the current study can be obtained from the corresponding author by an e-mail request.

Author contributions

MSK and VVSSS planned, executed and prepared the paper. RV, MHKP and URK were involved in the sample collection, data analysis, graphical presentation and scientific discussions.

Competing interests

The authors declare that they have no conflict of interest.

Special issue statement

This article is part of the special issue “Human impacts on carbon fluxes in Asian river systems”. It is not associated with a conference.


We would like to thank the director of the CSIR – National Institute of Oceanography (NIO), Goa, and the scientist in charge of the NIO Regional Centre in Visakhapatnam for their kind support and encouragement. We also acknowledge Maripi Dileep Kumar, NIO Goa, for his guidance and encouragement. The work is part of the Council of Scientific and Industrial Research (CSIR) funded research project. This publication has the NIO contribution number 6339.

Edited by: Gwenaël Abril
Reviewed by: Garuda V. M. Gupta and two anonymous referees


Abongwa, P. T. and Atekwana, E. A.: Assessing the temporal evolution of dissolved inorganic carbon in waters exposed to atmospheric CO2: a laboratory approach, J. Hydrol., 505, 250–265, 2014. 

Araujo, M., Noriega, C., and Lefevre, N.: Nutrients and carbon fluxes in the estuaries of major rivers flowing into the tropical Altantic, Front. Mar. Sci., 1, 1–16, 2014. 

Barber, R. T., Marra, J., Bidigare, R. C., Codispoti, L. A., Halpern, D., Johnson, Z., Latasa, M., Goericke, R., and Smith, S. L.: Primary productivity and its regulation in the Arabian Sea during 1995, Deep-Sea Res. Pt. II, 48, 1127–1172, 2001. 

Berner, R. A., Lasaga, A. C., and Garrels, R. M.: The carbonate–silicate geochemical cycle and its effect on atmospheric carbon dioxide over the past 100 million years, Am. J. Sci., 283, 641–683, 1983. 

Bhavya, P. S., Sanjeev Kumar, Gupta, G. V. M., Sudharma, K. V., and Sudheesh, V.: Spatio-temporal variation in δ13CDIC of a tropical eutrophic estuary (Cochin estuary, India) and adjacent Arabian Sea, Cont. Shelf Res., 153, 75–85,, 2018. 

Bouillon, S., Abril, G., Borges, A. V., Dehairs, F., Govers, G., Hughes, H. J., Merckx, R., Meysman, F. J. R., Nyunja, J., Osburn, C., and Middelburg, J. J.: Distribution, origin and cycling of carbon in the Tana River (Kenya): a dry season basin-scale survey from headwaters to the delta, Biogeosciences, 6, 2475–2493,, 2009. 

Brunet, F., Gaiero, D. M., Probst, J. L., Depetris, P. J., Gauthier-Lafaye, F., and Stille, P.: δ13C tracing of dissolved inorganic carbon sources in Patagonian rivers (Argentina), Hydrol. Process., 19, 3321–3344,, 2005. 

Cai, W. J., Guo, X. H., Chen, C. T. A., Dai, M. H., Zhang, L. J., Zhai, W. D., Lohrenz, S. E., Yin, K. D., Harrison, P. J., and Wang, Y. C.: A comparative overview of weathering intensity and HCO3 flux in the world's major rivers with emphasis on the Changjiang, Huanghe, Zhujiang (Pearl) and Mississippi Rivers, Cont. Shelf Res., 28, 1538–1549, 2008. 

Campeau, A., Wallin, M. B., Giesler, R., Löfgren, S., Mörth, C.-M., Schiff, S., Venkiteswaran, J. J., and Bishop, K.: Multiple sources and sinks of dissolved inorganic carbon across Swedish streams, refocusing the lens of stable C isotopes, Nature Sci. Rep., 7, 9158,, 2017. 

Carritt, D. E. and Carpenter, J. H.: Comparison and evaluation of currently employed modifications of the Winkler method for determining dissolved oxygen in seawater: A NASCO report, J. Mar. Res., 24, 286–318, 1966. 

Catalan, N., Marce, R., Kothawala, D. N., and Tranvik, L. J. Organic carbon decomposition rates controlled by water retention time across inland waters, Nat. Geosci., 9, 501–504,, 2016. 

Central Water Commission: Integrated Hydrological Data Book, New Delhi, India, 383 pp., 2006. 

Central Water Commission: Integrated Hydrological Data Book, New Delhi, India, 680 pp., 2012. 

Coplen, T. B.: New guidelines for reporting stable hydrogen, carbon and oxygen isotope-ratio data, Geochim. Cosmochim. Ac., 60, 3359–3360, 1996. 

Coplen, T. B., Hopple, J. A., Böhlke, J. K., Peiser, H. S., Rieder, S. E., Krouse, H. R., Rosman, K. J. R., Ding, T., Vocke Jr., R. D., Révész, K. M., Lamberty, A., Taylor, P., and De Bièvre, P.: Compilation of minimum and maximum isotope ratios of selected elements in naturally occurring terrestrial materials and reagents, U.S. Department of the Interior and U.S. Geological Survey, USA, 2002. 

Crawford, J. T., Lottig, N. R., Stanely, E. H. Walker, J. F., Hanson, P. C., Finlay, J. C., and Strieg, R. G.: CO2 and CH4 emissions from streams in a lake-rich landscape: Patterns, controls, and regional significance, Global Biogeochem. Cy., 28, 197–210,, 2014. 

Das, A., Krishnaswami, S., and Bhattacharya, S. K.: Carbon isotope ratio of dissolved inorganic carbon (DIC) in rivers draining the Deccan Traps, India: Sources of DIC and their magnitudes, Earth Planet. Sc. Lett., 236, 419–429,, 2005. 

Deines, P., Langmuir, D., and Harmon, R. S.: Stable carbon isotope ratios and the existence of a gas phase in the evolution of carbonate ground waters, Geochim. Cosmochim. Ac., 38, 1147–1164,, 1974. 

Deirmendjian, L. and Abril, G.: Carbon dioxide degassing at the groundwater-stream atmosphere interface: isotopic equilibration and hydrological mass balance in a sandy watershed, J. Hydrol., 558, 129–143,, 2018. 

Dessert, C., Dupre, B., Francois, L. M., Schott, J., Gaillardet, J., Chakrapani, G., and Bajpai, S.: Erosion of Deccan Traps determined by river geochemistry: impact on the global climate and the 87Sr86Sr ratio of seawater, Earth Planet. Sc. Lett., 188, 459–474, 2001. 

Finlay, J. C.: Controls of stream water dissolved inorganic carbon dynamics in a forested watershed, Biogeochemistry, 62, 231–252, 2003. 

Finlay, J. C.: Patterns and controls of lotic algal stable carbon isotope ratios, Limnol. Oceanogr., 49, 850–861, 2004. 

Gaillardet, J., Dupre, B., Louvat, P., and Allegre, C. J.: Global silicate weathering and CO2 consumption rates deduced from the chemistry of large rivers, Chem. Geol., 159, 3–30, 1999. 

Garcia, E. H. and Gordon, L. I.: Oxygen solubility in seawater better fitting equations, Limnol. Oceanogr., 37, 1307–1312,, 1992. 

Gauns, M., Madhupratap, M., Ramaiah, N., Jyothibabu, R., Fernandes, V., Paul, J. T., and Kumar, S. P.: Comparative accounts of biological productivity characteristics and estimates of carbon fluxes in the Arabian Sea and the Bay of Bengal, Deep-Sea Res. Pt. II, 52, 2003–2017, 2005. 

Giesler, R., Morth, C.-M., Karlsson, J., Lundin, E. J., Lyon, S. W., and Humborg, C.: Spatiotemporal variations of pCO2 and δ13C-DIC in subarctic streams in northern Sweden, Global Biogeochem. Cy., 27, 176–186,, 2013. 

Gupta, G. V. M., Sarma, V. V. S. S., Robin, R. S., Raman, A. V., Jai Kumar, M., Rakesh, M., and Subramanian, B. R.: Influence of net ecosystem metabolism in transferring riverine organic carbon to atmospheric CO2 in a tropical coastal lagoon (Chilka Lake, India), Biogeochemistry, 87, 265–285,, 2008. 

Gupta, G. V. M., Thottathil, S. D.,Balachandran, K. K., Madhu, N. V., Madeswaran, P., and Nair, S.: CO2 supersaturation and net heterotrophy in a tropical estuary (Cochin, India): influence of anthropogenic effect, Ecosystems, 12, 1145–1157,, 2009. 

Gupta, H., Chakrapani, G. J., Selvaraj, K., and Kao, S.-J.: The fluvial geochemistry, contributions of silicate, carbonate and saline–alkaline components to chemical weathering flux and controlling parameters: Narmada River (Deccan Traps), India, Geochim. Cosmochim. Ac., 75, 800–824, 2011. 

Gurumurthy, G. P., Balakrishna, K., Riotte, J., Braun, J.-J., Audry, S., Shankar, H. N. U., and Manjunatha, B. R.: Controls on intense silicate weathering in a tropical river, southwestern India, Chem. Geol., 300–301, 61–69, 2012. 

Hotchkiss, E. R., Hall Jr., R. O., Sponseller, R. A., Butman, D., Klaminder, J., Laudon, H., Rosvall, M., and Karlsson, J.: Sources of and processes controlling CO2 emissions change with the size of streams and rivers, Nat. Geosci., 8, 696–699,, 2015. 

Huang, T.-H., Fu, Y.-H., Pan, P.-Y., and Arthur, C. T.: Fluvial carbon fluxes in tropical rivers, Curr. Opin. Environ. Sust., 4, 162–169, 2012. 

Jarvie, H. P., King, S. M., and Neal, C.: Inorganic carbon dominates total dissolved carbon concentrations and fluxes in British rivers: Application of the THINCARB model – Thermodynamic modeling of inorganic carbon in freshwaters, Sci. Total Environ., 575, 496–512, 2017. 

Joesoef, A., Kirchman, D. L., Sommerfield, C. K., and Cai, W.-J.: Seasonal variability of the inorganic carbon system in a large coastal plain estuary, Biogeosciences, 14, 4949–4963,, 2017. 

Keeling, C. D.: The concentration and isotopic abundances of atmospheric carbon dioxide in rural areas, Geochim. Cosmochim. Ac., 13, 322–334,, 1958. 

Kishwan, J., Pandey, R., and Dhadwal, V. K.: India's forest and tree cover: Contribution as a carbon sink, Tech. Pap. 130, ICFRE BL-23, 2009. 

Kokic, J., Wallin, M. B., Chmiel, H. E., Denfeld, B. A., and Sobek, S.: Carbon dioxide evasion from headwater systems strongly contributes to the total export of carbon from a small boreal lake catchment, J. Geophys. Res. Biogeosciences, 120, 13–28,, 2015. 

Krishna, M. S., Prasad, V. R., Sarma, V. V. S. S., Reddy, N. P. C., Hemalatha, K. P. J., and Rao, Y. V.: Fluxes of dissolved organic carbon and nitrogen to the northern Indian Ocean from the Indian monsoonal rivers, J. Geophys. Res.-Biogeo., 120, 2067–2080, 2015. 

Kumar, R., Singh, R. D., and Sharma, K. D.: Water resources of India, Curr. Sci. India, 89, 794–811, 2005. 

Land, L. S.: The isotopic and trace element geochemistry of dolomite, in: the state of the art, Concepts and Models of Dolomitization, edited by: Zenger, D. H., Dunham, J. B., and Ethington, R. A., SEPM Special Publication No. 28, 63, 485,, 1980. 

Laskar, A. H., Gandhi, N., Thirumalai, K., Yadava, M. G., Ramesh, R., Mahajan, R. R., and Kumar, D.: Stable carbon isotopes in dissolved inorganic carbon: extraction and implications for quantifying the contributions from silicate and carbonate weathering in the Krishna River system during peak discharge, Isot. Environ. Healt. S., 50, 156–168, 2014. 

Löfgren, S., Froberg, M., Yu, J., Nisell, J., and Ranneby, B.: Water chemistry in 179 randomly selected Swedish headwater streams relatedto forest production, clear-felling and climate, Environ Monit Assess., 186, 8907–8928,, 2014. 

Ludwig, W., Amiotte-Suchet, P., Munhoven, G., and Probst, J. L.: Atmospheric CO2 consumption by continental erosion: present-day controls and implications for the last glacial maximum, Global Planet. Change, 17, 107–120, 1998. 

Madhupratap, M., Prasanna Kumar, S., Bhattathiri, P. M. A., Kumar, M. D., Raghukumar, S., Nair, K. K. C., and Ramaiah, N.: Mechanism of the biological response to winter cooling in the northeastern Arabian Sea, Nature, 384, 549–552, 1996. 

Mayorga, E., Aufdenkampe, A., Masiello, C. A., Krusche, A. V., Hedge, J. I., Quay, P. D., Richey, J. E., and Thomas, A. B.: Young organic matter as a source of carbon dioxide outgassing from Amazonian rivers, Nature, 436, 538–541,, 2005. 

Meybeck, M.: Global chemical weathering of surficial rocks estimated from river dissolved loads, Am. J. Sci., 287, 401–428, 1987. 

Meybeck, M. and Ragu, A.: GEMS/water contribution to the Global Register of River Inputs (GLORI), Provisional Final Rep., UNEP/WHO/UNESCO, Geneva, Switzerland, 245 pp., 1995. 

Meybeck, M. and Ragu, A.: River discharges to the oceans. An assessment of suspended solids, major ions, and nutrients, Environ. Inf. and Assess. Rep., 240,, 1996. 

Meybeck, M. and Vorosmarty, C. J.: Global transfer of carbon by rivers, Global Change News Lett., 37, 18–19, 1999. 

Miller, J. B. and Tans, P. P.: Calculating isotopic fractionation from atmospheric measurements at various scales, Tellus B, 55, 207–214,, 2003. 

Mook, W. and Tan, F.: Stable carbon isotopes in rivers and estuaries, Biogeochemistry of major world rivers, Wiley, Chichester, UK, 42, 245–264, 1991. 

Muduli, P., Vardhan, K. V., Robin, S. RCharankumar, B., Sivaji, P., Akkur, R., Gollapalli, R., Subramanian, B. R.: Distribution of dissolved inorganic carbon and net ecosystem production in a tropical brackish water lagoon, India. Cont. Shelf. Res., 64, 75–87,, 2013 

Mukhopadhyay, S. K., Biswas, H., De, T. K., Sen, S., and Jana, T. K.: Seasonal effects on the air–water carbon dioxide exchange in the Hooghly estuary, NE coast of Bay of Bengal, India, J. Environ. Monitor., 4, 549–552, 2002. 

Muraleedharan, P. M. and Prasanna Kumar, S.: Arabian Sea upwelling – A comparison between coastal and open ocean regions, Curr. Sci. India, 71, 842–846, 1996. 

O'Leary, M. H.: Carbon Isotopes in Photosynthesis, BioScience, 38, 328–336,, 1988. 

Opsahl, S. P. and Zepp, R. G.: Photochemically-induced alteration of stable carbon isotope ratios (δ13C) in terrigenous dissolved organic carbon, Geophys. Res. Lett., 28, 2417–2420,, 2001. 

Öquist, M. G., Wallin, M., Seibert, J., Bishop, K., and Laudon, H.: Dissolved inorganic carbon export across the soil/stream interface and its fate in a boreal headwater stream. Environ. Sci. Tech., 43, 7364–7369,, 2009. 

Parker, S. R., Poulson, S. R., Gammons, C. H., and DeGrandpre, M. D.: Biogeochemical controls on diel cycling of stable isotopes of dissolved O2 and dissolved inorganic carbon in the Big Hole River, Montana, Environ. Sci. Technol., 39, 7134–7140,, 2005. 

Parker, S. R., Poulson, S. R., Smith, M. G., Weyer, C. L., and Bates, K. M.: Temporal variability in the concentration and stable carbon isotope composition of dissolved inorganic and organic carbon in two Montana, USA Rivers, Aquat Geochem., 16, 61–84,, 2010. 

Pataki, D. E., Ehleringer, J. R., Flanagan, L. B., Yakir, D., Bowling, D. R., Still, C. J., Buchmann, N., Kaplan, J. O., and Berry, J. A.: The application and interpretation of Keeling plots in terrestrial carbon cycle research, Global Biogeochem. Cy., 17, 1022,, 2003. 

Pattanaik, J. K., Balakrishnan, S., Bhutani, R., and Singh, P.: Estimation of weathering rates and CO2 drawdown based on solute load: Significance of granulites and gneisses dominated weathering in the Kaveri River basin, Southern India, Geochim. Cosmochim. Ac., 121, 611–636, 2013. 

Pattanaik, S., Sahoo, R. K., Satapathy, D. R., Panda, C. R., Choudhury, S. B., and Mohapatra, P.: Intra-annual Variability of CO2 Flux in the Mahanadi Estuary – A Tropical Estuarine System, India. Ann. Mar. Sci., 1, 005–012, 2017. 

Polsenaere, P. and Abril, G.: Modelling CO2 degassing from small acidic rivers using water pCO2, DIC and δ13C-DIC data. Geochim. Cosmochim. Ac., 91, 220–239,, 2012. 

Prasanna Kumar, S., Muraleedharan, P. M., Prasad, T. G., Gauns, M., Ramaiah, N., de Souza, S. N., Sardesai, S., and Madhupratap, M.: Why is the Bay of Bengal less productive during summer monsoon compared to the Arabian Sea?, Geophys. Res. Lett., 29, 2235,, 2002. 

Prasanna Kumar, S., Nuncio, M., Narvekar, J., Kumar, A., Sardessai, S., DeSousa, S. N., Gauns, M., Ramaiah, N., and Madhupratap, M.: Are eddies nature's trigger to enhance biological productivity in the Bay of Bengal?, Geophys. Res. Lett., 31, L07309,, 2004. 

Quay, P. D., Wilbur, D. O., Richey, J. E., Hedges, J. I., Devol, A. H., and Victoria, R.: Carbon cycling in the Amazon River: Implications from the 13C compositions of particles and solutes, Limnol. Oceanogr., 37, 857–871,, 1992. 

Raymond, P. A. and Cole, J. J.: Increase in the export of alkalinity from North America's largest river, Science, 301, 88–91,, 2003. 

Raymond, P. A., Oh, N. H., Turner, R. E., and Broussard, W.: Anthropogenically enhanced fluxes of water and carbon from the Mississippi River, Nature, 451, 449–452,, 2008. 

Raymond, P. A., Hartmann, J., Lauerwald, R., Sobek, S., McDonald, C., Hoover, M., Butman, D., Striegl, R., Mayorga, E., Humborg, Ch., Kortelainen, P., Durr, H., Meybeck, M., Ciais, P., and Cuth, P.: Global carbon dioxide emissions from inland waters, Nature, 503, 355–359,, 2013. 

Richey, J. E., Melack, J. M., Aufdenkampe, A. K., Ballester, V. M., and Hess, L. L.: Outgassing from Amazonian rivers and wetlands as a large tropical source of atmospheric CO2, Nature, 416, 617–620, 2002. 

Samanta, S., Dalai, T. K., Pattanai K. J. K., Rai, S. K., and Mazumdar, A.: Dissolved inorganic carbon (DIC) and its δ13C in the Ganga (Hooghly) River estuary, India: Evidence of DIC generation via organic carbon degradation and carbonate dissolution, Geochim. Cosmochim. Ac., 165, 226–248, 2015. 

Sarma, V. V. S. S., Kumar, M. D., and Manerikar, M.: Emission of carbon dioxide from a tropical estuarine system, Goa, India, Geophys. Res. Lett., 28, 1239–1242,, 2001. 

Sarma, V. V. S. S., Kumar, N. A. , Prasad, V. R., Venkataramana, V., Appalanaidu, S., Sridevi, B., Kumar, B. S. K., Bharati, M. D., Subbaiah, Ch. V., Acharyya, T., Rao, G. D., Viswanadham, R., Gawade, L., Manjary, D. T., Kumar, P. P., Rajeev, K., Reddy, N. P. C., Sarma, V. V., Kumar, M. D., Sadhuram, Y., and Murty, T. V. R.: High CO2 emissions from the tropical Godavari estuary (India) associated with monsoon river discharges, Geophys. Res. Lett., 38, L08601,, 2011. 

Sarma, V. V. S. S., Viswanadham, R., Rao, G. D., Prasad, V. R., Kumar, B. S. K., Naidu, S. A. Kumar, N. A., Rao, D. B., Sridevi, T., Krishna, M. S., Reddy, N. P. C., Sadhuram, Y., and Murty, T. V. R.: Emission of carbon dioxide from the Indian monsoonal estuaries, Geophys. Res. Lett., 39, L03602,, 2012. 

Shetye, S. R., Gouveia, A. D., and Shenoi, S. S. C.: Circulation and water masses of the Arabian Sea, P. Indian A. S.-Earth, 103, 107–123, 1994. 

Shin, W. J., Chung, G. S., Lee, D., and Lee, K. S.: Dissolved inorganic carbon export from carbonate and silicate catchments estimated from carbonate chemistry and δ13CDIC, Hydrol. Earth Syst. Sci., 15, 2551–2560,, 2011. 

Singh, S. K., Sarin, M. M., and France-Lanord, C.: Chemical erosion in the eastern Himalaya: Major ion composition of the Brahmaputra and 13C of dissolved inorganic carbon, Geochim. Cosmochim. Ac., 69, 3573–3588, 2005. 

Smith S. L.: Understanding the Arabian Sea: reflections on the 1994–1996 Arabian Sea expedition, Deep-Sea Res. Pt. II, 48, 1385–1402, 2001. 

Solomon, D. K. and Cerling, T. E.: The annual carbon dioxide cycle in a montane soil: observations, modeling and implications for weathering, Water Resour. Res., 23, 2257–2265, 1987. 

Soman, M. K. and Kumar, K. K.: Some aspects of daily rainfall distribution over India during the southwest monsoon season, Int. J. Climatol., 19, 299–311, 1990. 

Sreenivas, K., Dadhwal, V. K., Suresh, K., Sri Harsha, G., Tarik, M., Sujatha, G, Suresh, J. R. G., Fyzee, M., and Ravisankar, T.: Digital mapping of soil organic and inorganic carbon status in India, Geoderma, 269, 160–173,, 2016. 

Sridevi, B., Sarma, V. V. S. S., Murty, T. V. R., Sadhuram, Y., Reddy, N. P. C., Vijayakumar, K., Raju, N. S. N., Jawahar Kumar, C., Raju, Y. S. N., Luis, R., Kumar, M. D., and Prasad, K. V. S. R.: Variability in stratification and flushing times of the Gautami–Godavari estuary, India, J. Earth. Syst. Sci., 124, 993–1003, 2015. 

Subramanian, V.: Sediment load of Indian rivers, Curr. Sci. India, 64, 928–930, 1993. 

Suzuki, R. and Ishimaru, T.: An improved method for the determination of phytoplankton chlorophyll using N,N-dimethyl formamide, J. Oceanogr., 46, 190–194, 1990. 

Tamooh, F., Borges, A. V., Meysman, F. J. R., Van Den Meersche, K., Dehairs, F., Merckx, R., and Bouillon, S.: Dynamics of dissolved inorganic carbon and aquatic metabolism in the Tana River basin, Kenya, Biogeosciences, 10, 6911–6928,, 2013.  

Tamse, S., Ogrinc, N., Walter, L. M, Turk, D., and Faganeli, J.: River Sources of Dissolved Inorganic Carbon in the Gulf of Trieste (N Adriatic): Stable Carbon Isotope Evidence, Estuar. Coast., 38, 151–164,, 2014. 

Vähätalo, A. V. and Wetzel, R. G.: Long-term photochemical and microbial decomposition of wetland-derived dissolved organic matter with alteration of 13C:12C mass ratio, Limnol. Oceanogr., 53, 1387–1392,, 2008. 

Varkey, M. J., Murty, V. S. N., and Suryanarayana, A.: Physical oceanography of the Bay of Bengal and Andaman Sea, Oceanogr. Mar. Biol., 34, 1–70, 1996. 

Venkiteswaran, J. J., Schiff, S. L., and Wallin, M. B.: Large Carbon Dioxide Fluxes from Headwater Boreal and Sub-Boreal Streams, PLoS ONE, 9, e101756,, 2014. 

Viers, J., Oliva, P., Dandurand, J. L., Dupré, B., Gaillardet, J., Heinrich, D. H., and Karl, K. T.: Chemical weathering rates, CO2 consumption, and control parameters deduced from the chemical composition of Rivers, Treatise on Geochemistry Pergamon, Oxford, UK, 2007. 

Vijith, V., Sundar, D., and Shetye, S. R.: Time-dependence of salinity in monsoonal estuaries, Estuar. Coast. Shelf S., 85, 601–608,, 2009. 

Waldron, S., Scott, E. M., and Soulsby, C.: Stable isotope analysis reveals lower-order river dissolved inorganic carbon pools are highly dynamic, Environ. Sci. Technol., 41, 6156–6162,, 2007. 

Wallin, M. B., Grabs, T., Buffam, I., Laudon, H., Agren, A., Oquist, M. G., and Bishop, K.: Evasion of CO2 from streams – The dominant component of the carbon export through the aquatic conduit in a boreal landscape, Glob. Change Biol., 19, 785–797,, 2013. 

Wang, Z. A., Bienvenu, D. J., Mann, P. J., Hoering, K. A., Poulsen, J. R., Spencer, R. G. M., and Holmes, R. M.: Inorganic carbon speciation and fluxes in the Congo River, Geophys. Res. Lett., 40, 511–516, 2013. 

Zeng, F.-W., Masiello, C. A., and Hockaday, W. C.: Controls on the origin and cycling of riverine dissolved inorganic carbon in the Brazos River, Texas, Biogeochemistry, 104, 275–291,, 2011. 

Zhang, J.: Carbon isotope fractionation during gas-water exchange and dissolution of CO2, Geochim. Cosmochim. Ac., 59, 107–114,, 1995. 

Zou, J.: Sources and Dynamics of Inorganic Carbon within the Upper Reaches of the Xi River Basin, Southwest China, PLoS One, 11, e0160964,, 2016. 

Short summary
An order-of-magnitude variability in DIC was found within the Indian estuaries due to significant variability in size of rivers, precipitation pattern and lithology in the catchments. Indian monsoonal estuaries annually export ∼ 10.3 Tg of DIC to the northern Indian Ocean, of which 75 % enters into the Bay of Bengal. Our results indicated that chemical weathering of carbonate and silicate minerals by soil CO2 is the major source of DIC in the Indian monsoonal rivers.
Final-revised paper