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

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 (δCDIC) also showed strong spatial variability (−13.0 ‰ to −1.4 ‰) in the Indian monsoonal rivers with relatively depleted δCDIC 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 δCDIC 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.


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 (CO 2 ), while the latter releases CO 2 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 CO 2 contribute DIC to rivers.On the other hand, autotrophic production by aquatic plants (photosynthesis) and evasion of CO 2 to atmosphere withdraw DIC from rivers.All these pro-Published by Copernicus Publications on behalf of the European Geosciences Union.
cesses in the catchments and within the rivers are strongly coupled to atmospheric CO 2 , because they act as either sinks or sources of atmospheric CO 2 (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 = 10 12 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 CO 2 at the airwater 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 δ 13 C DIC as a potential tracer, and to estimate the riverine export flux of DIC to the northern Indian Ocean from the Indian subcontinent.
The stable isotopic composition of DIC (δ 13 C DIC ) 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 CO 2 is about 0 ‰ (Coplen et al., 2002), whereas it is about −27 ‰ to −26 ‰ if the DIC is derived from the oxidation of or-ganic matter produced by C 3 plants (O'Leary, 1988).The δ 13 C of DIC generated by the carbonic-acid (formed by soil CO 2 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 δ 13 C DIC in the range of −8 ‰ to −7 ‰ and −4 ‰ to −3 ‰, respectively, if the carbonic acid is formed by the dissolution of atmospheric CO 2 .Though the δ 13 C of DIC derived from different sources is well separable (Deines et al., 1974), the isotopic fractionation by instream physical and biological processes alters the δ 13 C of the DIC source (Fig. 2).For example, photosynthesis and degassing enriches (O'Leary, 1988;Finlay, 2004;Parker et al., 2005Parker et al., , 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 δ 13 C 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 CO 2 (Raymond et al., 2013) and emit CO 2 to atmosphere due to oversaturation (Oquist et al., 2009;Campeau et al., 2017), the isotopic equilibration between the DIC and CO 2 in the atmosphere significantly influences the δ 13 C DIC in rivers (Abongwa and Atekwana, 2014;Deirmendjian and Abril, 2018) due to selective fluxes of 12 CO 2 and 13 CO 2 at the air-water interface.Hence, the influence of biogeochem- ical processes within the rivers must be considered while interpreting the δ 13 C DIC 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

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 × 10 12 m 3 yr −1 ) than into the Arabian Sea (0.3 × 10 12 m 3 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; https://www.gsi.gov.in/, last access: 5 November 2018).Based on discharge, the monsoonal rivers in this study were divided into two types, namely the major (> 150 m 3 s −1 ) and medium-sized (< 150 m 3 s −1 ) rivers.

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 GF/F (nominal pore size of 0.7 µm) under moderate vacuum and were stored in liquid nitrogen for Chl a analysis.

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 highpurity 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 CO 2 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 13 C, and X stands for 13 C / 12 C.The highpurity tank of CO 2 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 δ 13 C. 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 andRagu (1995, 1996), the Central Water Commission of New Delhi (2006Delhi ( , 2012)), and Kumar et al. (2005).Catchment area of the rivers was obtained from Water Resources Information System of India (WRIS, http://www.india-wris.nrsc.gov.in/, 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.

Hydrographic characteristics
Surface water temperatures were higher in rivers from the NE and SE regions (mean 30.9 ± 1.2 • C) than the rivers from the SW and NW regions (27.3 ± 1.5 • C) 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 mediumsized (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 δ 13 C DIC 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 δ 13 C DIC varied from −13.0 to −1.4 ‰, with a significant spatial variability (Fig. 3d; Table 1) in the rivers sampled.Relatively depleted δ 13 C DIC values were observed in rivers of the NW region (−11.1 ± 2.3 ‰), while enriched δ 13 C DIC was found in rivers of the SE region (−3.5 ± 2.3 ‰) of India (Fig. 3d).The δ 13 C DIC 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).

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 www.biogeosciences.net/16/505/2019/Biogeosciences, 16, 505-519, 2019 Table 1.Catchment 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 (δ 13 C DIC ) of Indian rivers were given.Measured pH and calculated δ 13 C values of CO 2 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.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).

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 (http://www.imd.gov.in, 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).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 m 3 m −2 ) than the rivers from the SE (0.17 m 3 m −2 ), NE (0.6 m 3 m −2 ) and NW (0.32 m 3 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 (r 2 = 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 ×10 3 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 (r 2 = 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 CO 2 with the atmosphere, autotrophic removal and heterotrophic addition of DIC.Since the Indian monsoonal estuaries have been reported to be a source of CO 2 to the atmosphere during the discharge period due to heterotrophic decomposition of organic matter (Sarma et al., 2001(Sarma et al., , 2011(Sarma et al., , 2012;;Gupta et al., 2008Gupta et al., , 2009;;Bhavya et al., 2018), the DIC input from the dissolution of atmospheric CO 2 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(Sarma et al., , 2012)).In contrast, significant negative correlation between chlorophyll a and DIC (r 2 = −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 δ 13 C DIC and Chl a (r 2 = 0.49; p<0.01;Fig. 4d), supporting this argument, because preferential uptake of 12 C over 13 C during photosynthesis leaves the residual DIC enriched in 13 C. On the other hand, δ 13 C DIC showed significant positive correlation with DO saturation (r 2 = 0.50, p<0.01;Fig. 4e; depleted δ 13 C DIC values with more undersaturation of DO) and DOC concentrations (r 2 = 0.43, p<0.01;Fig. 4f), as was observed in the Xi river (Zou et al., 2016).Altogether, enriched δ 13 C DIC values are associated with higher DOC, less undersaturation of DO and higher phytoplankton biomass (Chl a), while the depleted δ 13 C DIC 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 δ 13 C DIC data set, as the δ 13 C DIC in rivers is also influenced by pollution, catchment lithology and outgassing of CO 2 (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 (r 2 = 0.41, p<0.01;Fig. 4g), suggesting a considerable influence of pollution from the megacities and industries on DIC in the Indian 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 eastflowing 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 be- tween total alkalinity (TA) and δ 13 C DIC (r 2 = 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.

Major sources of DIC in the Indian monsoonal rivers
Though the δ 13 C DIC 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 δ 13 C of DIC source.The influence of major in-stream processes on the δ 13 C DIC 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 δ 13 C 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 δ 13 C of the source can be interpreted as the influence of the in-stream processes.Further, we approximated the δ 13 C of CO 2 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 δ 13 C DIC 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 (δ 13 C DIC as a function of 1/DIC; Fig. 5a) and the Miller-Tans plot (δ 13 C DIC × DIC as a function of DIC; Fig. 5b; r 2 = 0.61, p<0.01, and r 2 = 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 δ 13 C 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 CO 2 .Calculated δ 13 C of CO 2 ranged from −21.5 ‰ to −9.6 ‰ in the Indian rivers, with a mean value of −13.0 ± 2.7 ‰. Calculated δ 13 C of CO 2 is linearly correlated with the measured δ 13 C DIC , but correlation coefficient (r 2 ) is only 0.51 (Fig. 5c), suggesting significant spatial variability in the influence of in-stream prowww.biogeosciences.net/16/505/2019/Biogeosciences, 16, 505-519, 2019 cesses on the δ 13 C DIC .The Miller-Tans plot of CO 2 (δ 13 C-CO 2 × CO 2 as a function of CO 2 ) showed a highly significant linear regression model with a slope of −10.7 ‰ (r 2 = 0.97; p<0.001;Fig. 5d).These results indicated that chemical weathering of carbonate and silicate minerals by soil CO 2 (−10 ‰ to −9 ‰) is the major source of DIC in the Indian rivers.Deviations of the measured δ 13 C DIC (−13.0 ‰ to −1.4 ‰) from that of the approximated δ 13 C of DIC source (−3.0 ‰ to −2.0 ‰) and δ 13 C of CO 2 (−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 δ 13 C of the DIC source is towards the negative side (depletion; δ 13 C DIC <−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 δ 13 C DIC : −3.1 ‰), this could be due to the competition between autotrophy, degassing and heterotrophy, as these processes influence the δ 13 C DIC in opposite directions (Fig. 2); the former two processes cause enrichment, while the latter depletes δ 13 C DIC .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 δ 13 C DIC values of these rivers.

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 km 3 yr −1 ) from the catchment area of about 970 km ×10 3 km than the Arabian Sea (122 km 3 yr −1 from the catchment area of 244 km ×10 3 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 km 3 yr −1 ) compared to the American (11 799 km 3 yr −1 ) and African (3786 km 3 yr −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 river-ine 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.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 ×10 3 km to as high as 313 km ×10 3 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.

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 km 3 yr −1 ) and relatively smaller catchment (20 km ×10 3 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 (r 2 = 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 (r 2 = −0.49,p<0.001 -Fig.6b; r 2 = −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 (r 2 = 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 ×10 3 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 CO 2 that increases CO 2 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 (r 2 = 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.

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 δ 13 C of DIC source from the Keeling and Miller-Tans plots (−2.0 ‰ and −3.0 ‰, respectively), and the calculated δ 13 C of CO 2 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 www.biogeosciences.net/16/505/2019/Biogeosciences, 16, 505-519, 2019 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 km 3 yr −1 while the latter receives only 122 km 3 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.

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

Figure 3 .
Figure 3. Spatial 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 (δ 13 C DIC ; 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.

Figure 4 .
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 δ 13 C of DIC, respectively.Significant linear relationships of δ 13 C 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 δ 13 C 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).

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

Figure 6 .
Figure 6.Significant 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.