Temporary and net sinks of atmospheric CO2 due to chemical weathering in subtropical catchment with mixing carbonate and silicate lithology

The study provided the major ion chemistry, chemical weathering rates and temporary and net CO2 sinks in the Bei Jiang, which was characterized as a hyperactive region with high chemical weathering rates, carbonate and silicate mixing lithology, and abundant sulfuric acid chemical weathering agent of acid deposition and acid mining drainage (AMD) origins. The total chemical weathering rate of 85.46 t km−2 a−1 was comparable to that of other rivers in the hyperactive zones between the latitudes 0 and 30. A carbonate weathering rate of 61.15 t km−2 a−1 contributed to about 70 % of the total. The lithology, runoff, and geomorphology had a significant influence on the chemical weathering rate. The proportion of carbonate outcrops had a significant positive correlation with the chemical weathering rate. Due to the interaction between dilution and compensation effect, a significant positive linear relationship was detected between runoff and total carbonate and silicate weathering rates. The geomorphology factors such as catchment area, average slope, and hypsometric integral value (HI) had nonlinear correlation with chemical weathering rate and showed significant scale effect, which revealed the complexity in chemical weathering processes. Dissolved inorganic carbon (DIC) apportionment showed that CCW (carbonate weathering by CO2) was the dominant origin of DIC (35 %–87 %). SCW (carbonate weathering by H2SO4) (3 %– 15 %) and CSW (silicate weathering by CO2) (7 %–59 %) were non-negligible processes. The temporary CO2 sink was 823.41× 103 mol km−2 a−1. Compared with the temporary sink, the net sink of CO2 for the Bei Jiang was approximately 23.18× 103 mol km−2 a−1 of CO2 and was about 2.82 % of the “temporary” CO2 sink. Human activities (sulfur acid deposition and AMD) dramatically decreased the CO2 net sink, even making chemical weathering a CO2 source to the atmosphere.

Abstract. The study provided the major ion chemistry, chemical weathering rates and temporary and net CO 2 sinks in the Bei Jiang, which was characterized as a hyperactive region with high chemical weathering rates, carbonate and silicate mixing lithology, and abundant sulfuric acid chemical weathering agent of acid deposition and acid mining drainage (AMD) origins. The total chemical weathering rate of 85.46 t km −2 a −1 was comparable to that of other rivers in the hyperactive zones between the latitudes 0 and 30 • . A carbonate weathering rate of 61.15 t km −2 a −1 contributed to about 70 % of the total. The lithology, runoff, and geomorphology had a significant influence on the chemical weathering rate. The proportion of carbonate outcrops had a significant positive correlation with the chemical weathering rate. Due to the interaction between dilution and compensation effect, a significant positive linear relationship was detected between runoff and total carbonate and silicate weathering rates. The geomorphology factors such as catchment area, average slope, and hypsometric integral value (HI) had nonlinear correlation with chemical weathering rate and showed significant scale effect, which revealed the complexity in chemical weathering processes. Dissolved inorganic carbon (DIC) apportionment showed that CCW (carbonate weathering by CO 2 ) was the dominant origin of DIC (35 %-87 %). SCW (carbonate weathering by H 2 SO 4 ) (3 %-15 %) and CSW (silicate weathering by CO 2 ) (7 %-59 %) were non-negligible processes. The temporary CO 2 sink was 823.41 × 10 3 mol km −2 a −1 . Compared with the temporary sink, the net sink of CO 2 for the Bei Jiang was approximately 23.18 × 10 3 mol km −2 a −1 of CO 2 and was about 2.82 % of the "temporary" CO 2 sink. Human activities (sulfur acid deposition and AMD) dramatically decreased the CO 2 net sink, even making chemical weathering a CO 2 source to the atmosphere.

Introduction
Chemical weathering driven by weak carbonic acid (H 2 CO 3 ) that originates from atmosphere CO 2 or soil respiration under natural conditions is a fundamental geochemical process regulating the atmosphere-land-ocean carbon fluxes and Earth's climate (Guo et al., 2015). Carbonate and silicate weathering define the two typical categories of chemical weathering. From the view of the global carbon cycle, the CO 2 consumption due to carbonate weathering is recognized as the "temporary" sink because the flux of CO 2 consumed by carbonate dissolution on the continents is balanced by the flux of CO 2 released into the atmosphere from the oceans by carbonate precipitation on the geological timescale (Cao et al., 2015;Garrels, 1983). The consumption of CO 2 during chemical weathering of silicate rocks has been regarded as the net sink Published by Copernicus Publications on behalf of the European Geosciences Union.
of CO 2 and regulates the global carbon cycle (Hartmann et al., 2009(Hartmann et al., , 2014bKempe and Degens, 1985;Lenton and Britton, 2006). Thus in a carbonate-silicate mixing catchment, it is essential to distinguish proportions of the two most important lithological groups, i.e., carbonates and silicates, and evaluate the net CO 2 sink due to chemical weathering of silicate (Hartmann et al., 2009).
In addition to the chemical weathering induced by H 2 CO 3 , sulfuric acid (H 2 SO 4 ) of anthropogenic origins produced by sulfide oxidation such as acid deposition caused by fossil fuel burning and acid mining discharge (AMD) also becomes an important chemical weathering agent in the catchment scale. Many studies have shown the importance of sulfide oxidation and subsequent dissolution of other minerals by the resulting sulfuric acid at catchment scale (Hercod et al., 1998;Spence and Telmer, 2005). Depending on the fate of sulfate in the oceans, sulfide oxidation coupled with carbonate dissolution could facilitate a release of CO 2 to the atmosphere (Spence and Telmer, 2005). The carbonate weathering by H 2 SO 4 plays a very important role in quantifying and validating the ultimate CO 2 consumption rate. Thus, under the influence of human activities, the combination of silicate weathering by H 2 CO 3 and carbonate weathering by H 2 SO 4 controlled the net sink of atmospheric CO 2 .
Numerous studies on chemical weathering of larger rivers have been carried out to examine hydrochemical characteristics, chemical erosion and CO 2 consumption rates, and long-term climatic evolution of the Earth, for example for the Yangtze River (Chen et al., 2002;Ran et al., 2010), the Huang He (Zhang et al., 1995), the Pearl River (Gao et al., 2009;Xu and Liu, 2010;Zhang et al., 2007), the Huai River (Zhang et al., 2011), the rivers of the Qinghai-Tibetan Plateau Li et al., 2011;Wu et al., 2008), the Mekong River , the rivers of the Alpine region (Donnini et al., 2016), the Sorocaba River (Fernandes et al., 2016), the rivers of Baltic Sea catchment (Sun et al., 2017), the Amazon River (Gibbs, 1972;Mortatti and Probst, 2003;Stallard and Edmond, 1981;Edmond, 1983, 1987), the Lena River (Huh and Edmond, 1999), and the Orinoco River (Mora et al., 2010). For simplicity of calculation, most of the research has ignored sulfuric-acidinduced chemical weathering and resulted in an overestimation of CO 2 sink. To overcome this shortcoming of the traditional mass balance method, we applied a dissolved inorganic carbon (DIC) source apportionment procedure to discriminate the contribution of sulfuric-acid-induced chemical weathering to validate the temporary and net sink of CO 2 in a typical hyperactive region with carbonate-silicate mixing lithology to give a further understanding of basin-scale chemical weathering estimation.
About half of the global CO 2 sequestration due to chemical weathering occurs in warm and high-runoff regions (Ludwig et al., 1998), called the hyperactive regions and hot spots (Meybeck et al., 2006). The Pearl River located in the subtropical area in south China includes three principal rivers: the Xi Jiang, Bei Jiang, and Dong Jiang. The warm and wet climatic conditions make the Pearl River a hyperactive region in China. The three river basins have distinct geological conditions. The Xi Jiang is characterized as the carbonatedominated area, and the Dong Jiang has silicate as the main rock type. The Bei Jiang, which is the second largest tributary of the Pearl River, is characterized as a typical carbonatesilicate mixing basin. In addition, as the severe acid deposition (Larssen et al., 2006) and active mining area (Li et al., 2019), chemical weathering induced by sulfuric acid make it necessary to reevaluate the temporary and net sinks of atmospheric CO 2 . Thus, the Bei Jiang in southeast China with a typical subtropical monsoon climate and carbonate-silicate mixing geologic settings was selected as the study area. The three main objectives are summarized as follows: (1) reveal spatial-temporal variations in major element chemistry of the river water, (2) calculate the chemical weathering rate and unravel the controlling factors on chemical weathering processes, and (3) determine the temporary sink of CO 2 and evaluate the influence of sulfide oxidation on the net sink of CO 2 by DIC apportionment.

Study area
The Bei Jiang basin, which is the second largest tributary of the Pearl River basin, is located in the southeast of China (Fig. 1). It covers an area of 52 068 km 2 and has a total length of 573 km. The river basin is located in subtropical monsoon climate zone, with the mean annual temperature across the drainage basin ranging from 14 to 22 • C and the mean annual precipitation ranging from 1390 to 2475 mm. The average annual runoff is 51 × 10 9 m 3 , with 70 %-80 % of the flux occurring from April to September. This can be attributed to the fact that more than 70 % of the annual precipitation (about 1800 mm yr −1 ) is concentrated in the wet season (April to September).
Lithology in the river basin is composed of limestone, sandstone, gneiss, and glutenite. In the upper basin, carbonate rock (mainly of limestone) outcrops in the west and center, while sandstone of the Devonian era and mudstone of the Paleogene era outcrop in the east of the upper stream. In the middle of basin, limestone and sandstone cover most of the area, and Cretaceous volcanic rocks are found in the tributary (Lian Jiang), mainly granite. In the lower basin, Archean metamorphic rocks outcrop in the west and are composed of gneiss and schist; sandstone covers the rest of area of the lower basin. Quaternary sediments are scattered along the main stream of the river. The carbonate and silicate rock outcrops in the Bei Jiang basin are 10 737 km 2 (28 %) and 24 687 km 2 (65 %), respectively.

Sampling procedure and laboratory analysis
Water samples were collected monthly at 15 hydrologic stations from January to December in 2015 (Fig. 1). The river waters were sampled by a portable organic class water sampler along the middle thread of the channel on the first day of each month. In addition, to discriminate the contribution of rain inputs, the daily rainwater was also sampled at five stations (SJs, FLXs, YDs, XSs, and XGLs) along the main stream. The rainwater collector is consisted of a funnel with a diameter of 20 cm and a 5 L plastic bottle. A rubber ball is set up in the funnel to prevent evaporation. All the river and rainwater were filtered through a 0.45 µm glass fiber filter and stored in 100 mL tubes and stored below 4 • C until analysis.
Electric conductivity (EC), pH, and temperature (T ) were measured using a multi-parameter water quality meter (HACH-HQ40Q), and alkalinity (HCO − 3 ) was measured in filtered water samples by titration in situ. The dissolved SiO 2 was measured with the molybdenum yellow method and was analyzed using an ultraviolet spectrophotometer (Shimadzu UV-2600). The cations (Na + , K + , Ca 2+ , Mg 2+ ) and anions (Cl − , SO 2− 4 ) were analyzed with ion chromatography (Thermo Fisher ICS-900) with a limit of detection (L.O.D) of 0.01 mg L −1 . Reference, blank, and replicate samples were × 100 of the water samples was less than 5 %.

Chemical weathering rates
The mass balance equation for element X in the dissolved load can be expressed as (Galy and France-Lanord, 1999) (1) Here [X] denotes the elements of Ca 2+ , Mg 2+ , Na + , K + , Cl − , SO 2− 4 , and HCO − 3 in millimoles per liter. The subscripts riv, pre, eva, sil, car, and anth denote the river, precipitation source, evaporite source, silicate source, carbonate source, and anthropogenic source.
In the study area, the anthropogenic source of major ions except for SO 2− 4 was ignored due to the following two reasons. (1) Two main characteristics of much-polluted rivers are that total dissolved solid (TDS) is greater than 500 mg L −1 and the Cl − /Na + molar ratio is greater than that of sea salts (about 1.16) (Cao et al., 2016a;Gaillardet et al., 1999). The TDS in the study area ranged from 73.79 to 230.16 mg L −1 , and the low TDS implied that the anthropogenic origins of major ions could be ignored in the study. However, the Bei Jiang is characterized as a typical region suffering from severe acid deposition (Larssen et al., 2006) and are active mining area (Li et al., 2019). The acid deposition and acid mining discharge contribute to the highest concentration of SO 2− 4 .
(2) Natural origin of SO 2− 4 is the dissolution of evaporite, such as gypsum, while no evaporite was found in the study area. If SO 2− 4 comes from the gypsum dissolution, the ratios of Ca 2+ and SO 2− 4 should be close to 1 : 1. The stoichiometric analysis (Fig. 2) showed that the ratio of Ca 2+ and SO 2− 4 deviated from 1 : 1 and also proved this point.
Therefore, on the basis of the theory of rock chemical weathering and ignoring the anthropogenic origins of major ions (except for SO 2− 4 ), the major elements of river water can be simplified as follows.
Firstly, the measured ion concentrations of the rainwater are rectified by evaporation coefficient α = 0.63 = P /R (with P the precipitation and R the runoff), and the contributions of atmospheric precipitation are calculated. Secondly, the molar ratios of Ca 2+ /Na + (0.4) and Mg 2+ /Na + (0.2) for the silicate end-member (Zhang et al., 2007) are used to calculate the contribution of Ca 2+ and Mg 2+ from silicate weathering, and then, residual Ca 2+ and Mg 2+ were attributed to carbonate weathering. For monthly data, the contributions of different sources can be calculated as follows.
Here R denotes the proportions of dissolved cations from different sources. S denotes the total concentrations of cations for river water in millimoles per liter. The total, carbonate and silicate chemical weathering rates (TWR, CWR, and SWR) of a year can be estimated as follows.
Here TWR, CWR, and SWR in metric tons per square kilometer per year, Q i denotes discharge in cubic meters per month, and A denotes the catchment area in square kilometers.

DIC apportionments
In the Bei Jiang, the pH values of water samples ranged from 7.5 to 8.5 with an average of 8.05. Under these pH conditions, the major species of DIC is HCO − 3 . In addition, HCO − 3 accounted for more than 95 % at all sampling sites based on calculation; thus the concentration of HCO − 3 (mmol L −1 ) was used to represent the DIC concentration in this study. The riverine DIC originates from several sources including carbonate minerals, respired soil CO 2 , and atmospheric CO 2 , and it could be affected by processes occurring along the water pathways (Khadka et al., 2014;Li et al., 2008). Four dominant weathering processes, including (1) carbonate weathering by carbonic acid (CCW), (2) carbonate weathering by sulfuric acid (SCW), (3) silicate weathering by carbonic acid (CSW), (4) and silicate weathering by sulfuric acid (SSW), can be described by the following reaction equations.
Here CaSiO 3 represents an arbitrary silicate. According to the study of Galy and France-Lanord (1999) and Spence and Telmer (2005), carbonate and silicate weathering by carbonic acid have the same ratio as carbonate and silicate weathering by sulfuric acid, and for monthly data the mass balance equations are as follows: where the subscripts CCW, SCW, CSW, and SSW denote the four end-members defined by carbonate weathering by carbonic acid, carbonate weathering by sulfuric acid, silicate weathering by carbonic acid, and silicate weathering by sulfuric acid, respectively. The parameter α denotes the proportion of DIC derived from each end-member process. According to the above description, the ion balance equations are as follows.
Combining the above equations, the proportions of HCO − 3 derived from the three end-members (CCW, SCW, and CSW) can be calculated, and the DIC (equivalent to HCO − 3 ) fluxes by different chemical weathering processes are calculated by the following equations.

CO 2 consumption rate and CO 2 net sink
According to the Reactions (R1)-(R4), only the processes of CCW and CSW can consume the CO 2 from the atmosphere or soil and only half of the HCO − 3 in the water due to carbonate weathering by carbonic acid from atmospheric CO 2 . Thus, the CO 2 consumption rates (CCRs) for CCW and CSW can be calculated as follows (Zeng et al., 2016).
Here Q is discharge in cubic meters per year, [HCO − 3 ] is concentration of HCO − 3 in millimoles per liter, and A is catchment area in square kilometers, so that the CCR is in 10 3 mol km −2 a −1 .
According to the classical view of the global carbon cycling (Berner and Kothavala, 2001), the CCW is not a mechanism that can participate in the amount of CO 2 in the atmosphere because all of the atmospheric CO 2 fixed through CCW is returned to the atmosphere during carbonate precipitation in the ocean. However, when sulfuric acid is involved as a proton donor in carbonate weathering, half of the dissolved carbon is re-released to the atmosphere during carbonate precipitation. Thus, SCW leads to a net release of CO 2 in the ocean-atmosphere system over a timescale typical of residence times of HCO − 3 in the ocean (10 5 years). Meanwhile, in the case of CSW, followed by carbonate deposition, 1 of the 2 mol of CO 2 involved is transferred from the atmosphere to the lithosphere in the form of carbonate rocks, while the other 1 mol returns to the atmosphere, resulting in a net sink of CO 2 . Therefore, the net CO 2 consumption rate (CCR Net ) due to chemical weathering can be concluded as follows.

Spatial and statistical analysis
The hypsometric integral value (HI) (Pike and Wilson, 1971) was employed in this study to evaluate the influence of terrain on the chemical weathering. HI is an important index to reveal the relationship between morphology and development of landforms and can be used to establish the quantitative relationship between the stage of geomorphological development and the material migration in the basin (Pike and Wilson, 1971;Singh et al., 2008;Strahler, 1952). The HI value of each watershed is calculated by the elevation / relief ratio method and can be obtained by the following equation (Pike and Wilson, 1971): Here HI is the hypsometric integral, Mean.elevation is the mean elevation of the watershed, Min.elevation is the minimum elevation within the watershed, and Max.elevation is the maximum elevation within the watershed. According to the hypsometric integral value (HI), the geomorphological development can be divided into three stages: nonequilibrium or young stage (HI > 0.6), equilibrium or mature stage (0.35 < HI ≤ 0.6), and monadnock or old age (HI ≤ 0.35), which can reflect the erodible degree and erosion trend of the geomorphology (Xiong et al., 2014). The watershed of the study area was divided by using the hydrological analysis module of ArcGIS. The average slope and HI were calculated by the spatial analysis module of Ar-cGIS. The area of silicate-carbonate outcrops was calculated by the hydrological module of ArcGIS based on a geology map provided by the China Geological Survey. The data of river water discharge were provided by the local hydrology bureau.
All statistical tests were conducted using SPSS version 22.0. One-way analysis of variance (ANOVA) was performed to check the differences of monthly major ion concentrations and dissolved inorganic carbonate isotopes with significance at p < 0.05. Principal component analysis (PCA) was employed to unravel the underlying dataset through the reduced new variables, analyzing the significant factors affecting the characteristics of water chemistry.

Chemical compositions
The major physical-chemical parameters of river water samples are presented in Table 1. In Table 1, the chemical parameters of river water are the flow-weighted average over 12 months. For every sampling station, the flow-weighted average of ion concentration can be expressed followed the elements of Ca 2+ , Mg 2+ , Na + , K + , Cl − , SO 2− 4 , and HCO − 3 in millimoles liters. Q denotes average monthly discharge in cubic meters per second. The subscripts i denotes 12 months from January to December. For all the monthly samples, the pH values ranged from 7.5 to 8.5 with an average of 8.05. Average EC was 213 µs cm −1 , ranging from 81 to 330 µs cm −1 . The TDS of river water samples varied from 73.8 to 230.2 mg L −1 , with an average of 157.3 mg L −1 , which was comparable with the global average of 100 mg L −1 (Gaillardet et al., 1999). Compared with the major rivers in China, the average TDS was significantly lower than the Yangtze (Chen et al., 2002), the Huang He (He et al., 2017) the Zhu Jiang (Zhang et al., 2007), the Huai He (Zhang et al., 2011), and the Liao He (Ding et al., 2017). However, the average TDS was higher than the rivers draining silicate-rock-dominated areas, e.g., the Dong Jiang (59.9 mg L −1 ) in southern China (Xie et al., 2013), the North Han River (75.5 mg L −1 ) in South Korea, (Ryu et al., 2008), the Amazon River (41 mg L −1 ), and the Orinoco River (82 mg L −1 ) draining the Andes (Dosseto et al., 2006;Edmond et al., 1996).
Major ion compositions are shown in the Piper plot (Fig. 3). Ca 2+ was the dominant cation with concentration ranging from 199 to 1107 µmol L −1 , accounting for approximately 49 % to 81 %, with an average of 66 % (µEq) of the total cation composition in the river water samples. HCO − 3 was the dominant anion, with concentration ranging from 640 to 2289 µmol L −1 . On average, it comprised 77 % (59 %-92 %) of total anions, followed by SO 2− 4 (16 %) and Cl − (6 %). The major ionic composition indicated that the water chemistry of the Bei Jiang basin was controlled by both carbonate and silicate weathering.
The PCA was used to extract the factors controlling the chemical compositions. The varimax rotation was used to reduce the number of variables to two principal components (PCs), which together explain 76.88 % of the total variance in the data. The first PC (PC1) explained approximately 50.02 % of the total variations and was considered to represent the carbonate weathering factor because of the high loadings of EC, TDS, Ca 2+ , Mg 2+ , and HCO − 3 concentrations. The second PC (PC2) explained 26.85 % of the total variance and presented high loadings for Na + and K + con- centrations. Thus, the PC2 represented a silicate weathering factor. These two PCs were considered to be two important sources of major ions in the Bei Jiang basin.
The hydrochemical compositions of rainwater are presented in Table S1 in the Supplement. Ca 2+ was the dominant cation with concentration ranging from 6.9 to 282.6 µmol L −1 , accounting for approximately 65 % of the total cation composition in the rainwater samples. SO 2− 4 was the dominant anion, with concentration ranging from 21.9 to 1462 µmol L −1 , accounting for approximately 67 % of the total anion composition in the rainwater samples.

Seasonal and spatial variations
There were significant seasonal variations in the major ion concentrations (Fig. 4). Two basic patterns of temporal variations could be observed. The first one was related to the carbonate-weathering-derived ions such as Ca 2+ and HCO − 3 , which showed high values in November and low values in June. The second one was for the silicate-weathering-derived ions such as Na + and K + , which showed high values in February and low values in June. The minimums occurred in June, and all the ions showed a significant dilution effect during the high-flow periods.
It was clear that the Ca 2+ and HCO − 3 concentrations had a decreasing trend from upstream to downstream (Fig. 5); this characteristic agrees with the trends observed in the Yangtze River and the Huai River, where the major elements or TDS concentrations of the main channel showed a general decreasing trend, and the tributaries display the dilution effect to the main channel. For other silicate-weathering-derived ions such as Na + , there was a slight increasing trend implying the chemical inputs from the tributaries. These trends were in accordance with the lithology in the study area. The Table 1. The major physical-chemical parameters of river water samples at 15 hydrological stations in the Bei Jiang (mean ± SD). The total dissolved solid (TDS, mg L −1 ) expressed as the sum of major inorganic species concentrations (Na  carbonate dominates in the upper stream basin; when the river drains this area, carbonate weathering contributes to the elevation of Ca 2+ and HCO − 3 . As the river enters into the downstream dominated by silicate, the relatively low ion concentrations due to silicate weathering contributed to diluting the Ca 2+ and introducing extra Na + to the main channel. Atmospheric precipitation inputs, anthropogenic inputs (here referring to the acid deposition and AMD), and chemical weathering of rocks and minerals are the major sources contributing to the hydrochemistry in the river basin. Previous studies have shown that rock weathering contributions to major element compositions of the river can be interpreted in terms of mixing among three main end-members: the weathering products of carbonates, silicates, and evaporites (Cao et al., 2016b;Négrel et al., 1993;Ollivier et al., 2010). The river water samples in the Bei Jiang basin are displayed on the plots of Na-normalized molar ratios (Fig. 6). In these plots, the contributions from carbonate weathering correspond to the trend toward the high-Ca 2+ end-member close to the top right corner, while silicate weathering corresponds to the trend toward the high-Na + end-member close to the lower left corner. It was clear that the samples with high ratio of carbonate outcrop had the highest molar ratios of Ca 2+ /Na + , Mg 2+ /Na + , and HCO − 3 /Na + , which were the samples located toward the carbonate weathering end-member. However, the samples with low Ca 2+ /Na + , Mg 2+ /Na + , and HCO − 3 /Na + ratios showed the influence of silicate weathering. In addition, major ion compositions of the Bei Jiang were mainly contributed by the weathering of carbonates and silicates and showed little contribution of evaporite weathering.
Based on the chemical balance method, the calculated contributions of different sources to the total cationic loads are presented in Fig. 7. The results show that carbonate weathering was the most important mechanism controlling the local hydrochemistry and contributed approximately 50.06 % (10.96 %-79.96 %) of the total cationic loads. Silicate weathering and atmospheric precipitation inputs accounted for 25.71 % (5.55 %-70.38 %) and 17.92 % (0 %-46.95 %), respectively. Evaporite weathering had the minimum contribution with an average of 6.31 % (0 %-24.36 %) to the total cationic loads.
The result of chemical weathering rates is listed in Table 2. The carbonate weathering contributes about 70 % of the total chemical weathering, and the averages of carbonate and silicate weathering rate in the Bei Jiang basin were 61.15 and 25.31 t km −2 a −1 , respectively. In addition, chemical weathering rates showed significant seasonal variations, with the highest carbonate and silicate weathering rates in May (16.75 and 5.50 t km −2 per month, respectively) and the lowest carbonate and silicate weathering rates in February (0.95 and 0.39 t km −2 per month, respectively). Gaillardet et al. (1999) reported the chemical weathering rate of major rivers all over the world and found that the hyperactive zone with a high chemical weathering rate is generally located between the latitude 0 and 30 • and our study belongs to this area (Fig. 8). The factors influencing the balance between CWR and SWR will be further discussed in the following sections.

Factors affecting chemical weathering
Many factors control the chemical weathering rates, including terrain, geotectonic properties, lithology, land cover, climatic conditions (temperature, precipitation, etc.), and hydrological characteristics (Ding et al., 2017;Gislason et al., 2009;Hagedorn and Cartwright, 2009). For this study, the lithology, hydrological characteristics, and geomorphology were selected as the major factors to be discussed.

Lithology
Among all the factors controlling the chemical weathering rates, lithology is one of the most important factors because different types of rocks have different weathering abilities (Viers et al., 2014). The TWR had a significant positive correlation (p < 0.01) with the ratios of the proportion of carbonate and a nonsignificant positive correlation with that of silicate outcrops (Fig. 9a, b). Furthermore, a significant correlation (p < 0.01) was found between the CWR and proportion of carbonate outcrops (Fig. 9c), but the correlation between the SWR and the proportion of silicate outcrops was low and not statistically significant (p > 0.05, Fig. 9d). The correlation analysis confirmed that carbonate outcrop ratios were the sensitive factor controlling the chemical weathering rates, and the rapid kinetics of carbonate dissolution played an important role in weathering rates in the Bei Jiang basin.

Runoff
Chemical weathering is a combination of two processes, including dissolution of primary minerals and precipitation of secondary mineral growth (Eiriksdottir et al., 2011;Hartmann et al., 2014a;Liu et al., 2013). The dissolution process is related to precipitation and runoff. In general, river water chemistry is usually diluted by river runoff (Q), and this dilution effect is variable in different basins (Rao et al., 2019). The dilution effects of major elements caused by increasing water flow can be expressed by a log linear equation, the standard rating relationship Walling, 1986;Zhang et al., 2007): Figure 6. Mixing diagrams using Na-normalized molar ratios: (a) Mg 2+ /Na + vs. Ca 2+ /Na + and (b) HCO − 3 /Na + vs. Ca 2+ /Na + for the Bei Jiang basin. The color ramp showed the percentage of carbonate outcrops. where C i is the concentration of element i (mmol L −1 ), Q is the water discharge (m 3 s −1 ), a is the regression constant, and b is the regression exponent. The linear fitting result is shown by Fig. 10, and the parameters b for major elements obtained from the dataset were 0.08 (Na + ), 0.05 (K + ), 0.08 (Ca 2+ ), 0.02 (Mg 2+ ), 0.06 (HCO − 3 ), 0.12 (Cl − ), 0.11 (SO 2− 4 ), and −0.005 (SiO 2 ). In many cases, b ranges from −1 to 0 due to the chemical variables that are influenced in various ways and various extents. However, in our study area, the values of b were positive and not comparable to the observations in major Asian rivers such as the Yangtze (Chen et al., 2002), the Yellow (Chen et al., 2005), the Pearl (Zhang et al., 2007), and the Mekong . This suggested additional and significant solute sources in the river basin that might contribute to and compensate for the effect of dilution by precipitation. The difference of slope for individual dis-solved components at different stations reflected the different sources and the solubility of source materials.
Due to the compensation effect of chemical weathering, significant positive linear relationship was detected between Q and TWR, CWR, and SWR. Thus, the linear regression analyses between Q and TWR, CWR, and SWR were conducted to further reveal the effect of runoff on chemical weathering rate. The slope of the liner regression equations for all 15 hydrological station watersheds in the Bei Jiang basin are summarized in Table 3. The linear relations indicated that the increase in runoff could accelerate the chemical weathering rates, but the variations in K values revealed that the degrees of influences were different due to multiple factors, such as the influence of geomorphology.

Geomorphology
The geomorphology factors including catchment area, average slope, and HI, which could influence the runoff generation process and physical and chemical weathering, were selected to give a further explanation of the variation in K values. As shown in Fig. 11a, the K values were found to have a nonlinear relationship with the areas of subcatchment and could be fitted by an exponential decay model, which showed that the K values decreased dramatically with the initial increase in area and quickly become stable after reaching the threshold. The threshold value for K 1 , K 2 , and K 3 was about 5000 km 2 . This indicated that the compensation effect was more significant in small catchments. The average topographic slope of each subcatchment ranged from 37 to 63 • . With the increasing average slope, the residence time of both surface water and groundwater decreases. Kinetics of carbonate and silicate reactions was determined by the reaction time, which could be related by the residence time of water. In our study area, the K values showed nonlinear negative correlation with average slope (Fig. 11e, f, g). When the average slope increases, the resulting small residence time (time of water-rock reactions) makes the compensation effect weak in the study area.  Hypsometric analysis showed that the HI ranged from 0.18 to 0.34. According to the empirical classification by HI (HI > 0.6, nonequilibrium or young stage; 0.35 < HI ≤ 0.6, equilibrium or mature stage; HI ≤ 0.35, monadnock or old age), the geomorphological development in the Bei Jiang was recognized as old age, which reflects that the erodible degree and erosion trend of the geomorphology were high. Furthermore, the nonlinear positive correlations between HI and K values (Fig. 11g, h, i) also addressed the fact that geomorphology development has a significant influence on chemical weathering and related CO 2 consumption processes.

Sulfate origin and DIC apportionment
The successful application of DIC apportionment calculation mentioned in Sect. 3.2.2 is dependent on the origins of sulfate (SO 2− 4 ). Three origins of SO 2− 4 should be discriminated, including atmospheric acid deposition (Larssen and Carmichael, 2000), acid mining discharge (AMD) (Li et al., , 2019, and chemical weathering of evaporite such as the dissolution of gypsum (Appelo and Postma, 2005). Acid rain events occurred frequently in south and east China after 1980 (Larssen et al., 2006). The pH isolines based on data from 86 monitoring stations (Larssen et al., 2006) showed that in the Bei Jiang the rain pH was lower than 4.5, and our monitoring dataset also proved this result. Sulfur wet deposition estimated based on the observed bulk wet sulfur deposition data and the RAINS-Asia model (Larssen and Carmichael, 2000) ranged from 2000 to 5000 eq ha −1 a −1 , which showed that the acid sulfur deposition was one of the most important sources of river sulfate. In addition, considering the abundant ore resources in the Bei Jiang, the second possible source of SO 2− 4 is sulfide oxidation due to mining. In our previous study, the SO 2− 4 with AMD origin mainly came from the tributary Wen Jiang (Wen et al., 2018). These two sources could offer a sufficient amount of the chemical weathering agent H 2 SO 4 and be actively involved in the chemical weathering due to the following reaction mechanism (take carbonate for example) (Taylor et al., 1984;van Everdingen and Krouse, 1985).
The third source came from dissolution of gypsum and could not offer active H 2 SO 4 to induce carbonate and silicate dissolution. Two pieces of evidence were summarized to indicate the absence of gypsum in the study area.
(1) Lithology in the river basin is composed of limestone, sandstone, gneiss, and glutenite. HI showed that geomorphology development has entered into the old age, and the evaporite such as halite and gypsum has been consumed by dissolution.
(2) The stoichiometric relationship between Ca 2+ and SO 2− 4 ( Fig. 2) showed that all of the samples in the study area were located below the 1 : 1 gypsum dissolution line, and they are also below the 1 : 2 carbonate weathering line induced by sulfuric acid (SCW). These two points combined gave the evidence to prove the absence of contribution of gypsum dissolution to river SO 2− 4 . Thus, the DIC apportionment could be calculated according to Eqs. (23) to (25), and the results of the three main processes (CCW, CSW, and SCW) contributing to the DIC origin in the Bei Jiang water are shown in Table 4. It was found that CCW was the dominant origin of DIC (35 %-87 %) and that SCW (3 %-15 %) and CSW (7 %-59 %) were non-negligible weathering processes.

3886
Y. Cao et al.: Temporary and net sinks of atmospheric CO 2 due to chemical weathering Figure 11. The relationships between K values and catchment area (a, b, c), average slope (d, e, f), and HI (g, h, i) for the Bei Jiang.

Temporary and net CO 2 sink
According to the classical view of the global carbon cycling (Berner and Kothavala, 2001), the CO 2 sink induced by chemical weathering varies for different timescales. At a short-term timescale, carbonic-acid-based carbonate and silicate weathering (CCW and CSW) and transport of the HCO − 3 to oceans through rivers is an important temporary carbon sink (Khadka et al., 2014) and can be calculated by the sum of CCR CCW and CCR CSW . Thus, it was significant to estimate the CCR of CCW and CSW (Liu and Dreybrodt, 2015;Liu et al., 2011). However, at the geological timescale (> 10 6 years), over the timescale typical of the residence time of HCO − 3 in the ocean (10 5 years), the CCW is not a mechanism that can participate in the net sink of CO 2 in the atmosphere because all of the atmospheric CO 2 fixed through CCW is returned to the atmosphere during carbonate precipitation in the ocean. Meanwhile, in the case of CSW, followed by carbonate deposition, 1 of the 2 mol of CO 2 involved is transferred from the atmosphere to the lithosphere in the form of carbonate rocks, while the other 1 mol returns to the atmosphere. The CSW is recognized as the net sink of atmosphere CO 2 . In addition, when sulfuric acid is involved as a proton donor in carbonate weathering, half of the carbon is dissolved to the atmosphere during carbonate precipitation. Thus, SCW leads to a net release of CO 2 in the ocean-atmosphere system. Thus the net CO 2 sink (expressed by CCR Net in this study) is controlled by the DIC apportionment according to Eq. (28). The results of CCR Total , CCR CCW , CCR CSW , and CCR Net were summarized in Table 4. The CCR Total was 823.41 × 10 3 mol km −2 a −1 . Comparing with other Chinese rivers, such as the Songhua River (189 × 10 3 mol km −2 a −1 ) (Cao et al., 2015) and other rivers calculated by Gaillardet et al. (1999) including the Heilong River (53 × 10 3 mol km −2 a −1 ), the Yangtze River (609 × 10 3 mol km −2 a −1 ), the Huang He (360×10 3 mol km −2 a −1 ), the Xi Jiang (960 × 10 3 mol km −2 a −1 ), the Jinsha Jiang (420 × 10 3 mol km −2 a −1 ), the Langcang Jiang (980 × 10 3 mol km −2 a −1 ), the Nu Jiang (1240×10 3 mol km −2 a −1 ), the Yalong Jiang (870 × 10 3 mol km −2 a −1 ), the Dadu He (1280 × 10 3 mol km −2 a −1 ), and Min Jiang (660 × 10 3 mol km −2 a −1 ), our study area showed relatively high CCR due to a high chemical weathering rate. In addition, the CCR CCW and CCR CSW were 536.59 × 10 3 (65 %) and 286.82 × 10 3 (35 %) mol km −2 a −1 , respectively. Compared with the temporary sink, the net sink of CO 2 for the Bei Jiang was approximately 23.18×10 3 mol km −2 a −1 of CO 2 sinking in the global carbon cycling. It was about 3 % of the temporary CO 2 sink. In addition, the CO 2 net sink of each subbasin was also different and showed large spatial variations due to heterogeneity of geology and human activities. The geology showed weak correlation with the CO 2 net sink (Fig. 12a), while the [SO 2− 4 ] SCW and [SO 2− 4 ] SSW have weak negative correlation with the CO 2 net sink (Fig. 12b). This proved that human activities (sulfur acid deposition and AMD) decreased the CO 2 net sink and even made chemical weathering a CO 2 source to the atmosphere.

Conclusions
This study revealed the temporary and net sinks of atmospheric CO 2 due to chemical weathering in a subtropical hyperactive catchment with mixing carbonate and silicate lithology under the stress of chemical weathering induced by anthropogenic sulfuric acid agent. During the sampling period, the pH values ranged from 7.5 to 8.5 and TDS varied from 73.8 to 230.2 mg L −1 . Ca 2+ and HCO − 3 were the dominant cation and anion. Water chemical patterns and PCA showed that carbonate and silicate weathering were the most important processes controlling the local hydrochemistry. On average, carbonate and silicate weathering contributed approximately 50.06 % and 25.71 % of the total cationic loads, respectively.
The averages of carbonate and silicate weathering rate in the Bei Jiang basin were 61.15 and 25.31 t km −2 a −1 , respectively. The high rate was comparable to other rivers located in the hyperactive zone between the latitude 0 and 30 • . The lithology, runoff, and geomorphology had significant influences on the chemical weathering rate. (1) Due to the difference between kinetics of carbonate and silicate dissolution processes, the proportion of carbonate outcrops had significant positive correlation with the chemical weathering rate and confirmed that carbonate outcrop ratios were the sensitive factor controlling the chemical weathering rates, and the rapid kinetics of carbonate dissolution played an important role in weathering rates. (2) Runoff mainly controlled the season variations, and the dilution effect was weak in the study area. Due to the compensation effect of chemical weathering, a significant positive linear relationship was detected between Q and TWR, CWR, and SWR. (3) The geomorphology factors such as slope and HI had a nonlinear correlation with chemical weathering rate and showed a significant scale effect, which revealed the complexity in chemical weathering processes.
Compared with the temporary sink, the net sink of CO 2 for the Bei Jiang was approximately 23.18 × 10 3 mol km −2 a −1 of CO 2 sinking in global carbon cycling. It was about 2.82 % of the temporary CO 2 sink. Human activities such as sulfur acid deposition and AMD have significantly altered the CO 2 sinks.
Data availability. The data in this study have been presented in the tables of the article and can also be requested from the corresponding author.