Chemical Weathering and Riverine Carbonate System Driven by Human Activities in a Subtropical Karst Basin , South China

In the context of climate change, the input of acid substances into rivers, caused by human activities in the process of industrial and agricultural development, has significantly disrupted river systems and has had a profound impact on the carbon cycle. The hydrochemical composition and which main sources of the Lianjiang River (LR), a subtropical karst river in northern Guangdong Province, South China, were analyzed in January 2018. The objective was to explicate the influence on the deficit proportion of CO2 consumption, resulting from carbonate chemical weathering (CCW), driven by nitric acid (HNO3) and sulfuric acid (H2SO4), which is affected by exogenous acids from the industrial regions in north of the Nanling Mountains and the Pearl River Delta. The response of the riverine carbonate system to exogenous acid-related weathering was also discussed. HCO3 and Ca2+, respectively, accounted for 84.97% of the total anions and 78.71% of the total cations in the surface runoff of the LR, which was characterized as typical karst water. CCW was the most important material source of river dissolved loads in the LR, followed by human activities and silicate chemical weathering (SCW). Dissolved inorganic carbon (DIC), derived from CCW induced by carbonic acid (H2CO3), had the largest contribution to the total amount of DIC in the LR (76.79%), and those from CCW induced by anthropogenic acids (HNO3 and H2SO4) and SCW contributed 13.56% and 9.64% to the total DIC, respectively. The deficit proportion of CO2 consumption associated with CCW resulting from sulfuric acid and nitric acid (13.56%), was slightly lower than that of the Guizhou Plateau in rainy and pre-rainy seasons (15.67% and 14.17%, respectively). The deficit percentage of CO2 uptake associated with CCW induced by sulfuric acid and nitric acid, accounted for 38.44% of the total CO2 consumption related to natural CCW and 18.84% of the anthropogenic acids from external areas. DIC derived from CCW induced by human activities, had a significant positive correlation with the total alkalinity, SIc and pCO2 in river water, indicating that the carbonate system of the LR was also driven by exogenous acids, with the exception of carbonic acid. More attention should be paid to the effects of human activities on the chemical weathering and riverine carbonate system in the karst drainage basin.


Introduction
A river system is one of the most important transport channels of material and energy to the ocean in a drainage basin, which records the environment evolution processes of the river basin.Chemical weathering processes in a drainage basin can be reflected by the river hydrochemical characteristics [1,2], which is the main part of the mass exchange between the ocean and the continent.It is generally recognized that silicate chemical weathering (SCW) induced by carbonic acid is a net depletion process of the atmospheric/soil CO 2 .However, the CO 2 consumed by carbonate chemical weathering (CCW) associated with carbonic acid will return to the atmosphere with the precipitation of carbonate minerals on a geological time scale (>1 Ma).Thus, CCW only affects global climate change on a short time scale [3].Recently, some scholars have discovered that organic carbon can be produced by natural aquatic ecosystems through photosynthesis or directly by using dissolved inorganic carbon (DIC), which is derived from CCW and is rich in karst river (the fertilizing effect of DIC), that enters the lithosphere through sedimentation and burial [4][5][6].The CO 2 consumption associated with CCW have not only controlled climate change over a short time scale, but also dominated long-term climate change, since the emergence of aquatic photosynthetic organisms [7][8][9].
The karst basin ecosystem is subject to greater climate stress and intensifying human activities, which profoundly affect the river carbon cycle, and exhibits a higher sensitivity [10].With the development of urbanization, human activities (e.g., acid deposition, industrial and domestic sewage discharge, nitrogen fertilizer application, etc.) have gradually become an important factor influencing the chemical composition of river water.More and more studies have paid attention to the CCW influenced by nitric acid or sulfuric acid, which is derived from industrial and agricultural activities [11][12][13][14][15][16][17][18].Strong acid, such as sulfuric acid or nitric acid derived from anthropogenic sources, replaces carbonic acid agent and dissolves rocks in the process of rock chemical weathering.It also disturbs the rock weathering and the associated CO 2 consumption in natural conditions.Li et al. [19] analyzed the source of sulfate in the Jialingjiang River using stable isotopes (δ 34 S and δ 18 O) and found that domestic and industrial wastewater was an important source of sulfate.They also confirmed that sulfuric acid derived from the atmospheric input of high sulfur coal combustion was a weathering agent involved in rock weathering.
Chemical weathering in the drainage basin and fluvial chemical composition are more susceptible to human activities owing to the rapid dissolution kinetics of carbonate rocks and the special binary and three-dimensional geological structure in the karst area.A case study of two karst basins in Southern France showed that CCW induced by nitric acid, which was derived from nitrogen fertilizer, accounted for 30% [20] and 12-26% [21].The study on the glacial runoff by Anderson et al. [22] showed that the amount of carbonate rocks dissolved by sulfuric acid constituted 22% of the total dissolved carbonates.Chemical weathering of surface rocks is a basic geological process, which converts atmospheric/soil CO 2 to DIC [23].Generally, the DIC in the karst surface waters is mainly derived from carbonate rocks dissolved by carbonic acid.With the interference of human activities, the DIC produced by CCW induced by sulfuric acid and nitric acid may be an important part of the DIC flux in the karst surface water.Li et al. [24] reported that about 13% of DIC in the Yalongjiang River Basin in Southwest China was due to the involvement of sulfuric acid in CCW.
However, no atmospheric/soil CO 2 was consumed by CCW associated with strong acids.On the contrary, CCW with strong acids (especially sulfuric acid) has a net CO 2 release process on geological time scales [25], and Jiang [16] concluded that the atmospheric/soil CO 2 consumption of CCW should be strictly evaluated if sulfuric acid or nitric acid was involved by analyzing the main ions and δ 13 C DIC in the Nandong subterranean stream.Therefore, the carbon sink effect related to CCW needs to be reevaluated if sulfuric acid or nitric acid exists in a carbonate catchment.Perrin et al. [13] estimated that the CO 2 produced during CCW by nitric acid, which was derived from the nitrification of nitrogen fertilizers in the agricultural region of Southeastern France and around the world, were 2.2-5.3 × 10 5 t•a −1 and 1.2-2.9× 10 7 t•a −1 , respectively.Liu et al. [26] calculated the contribution of acid precipitation to the chemical composition of the Qiantangjiang River by the Galy Model and concluded that the estimated CO 2 consumption rate was 23% lower than that before deducting sulfuric acid weathering, if sulfuric acid involvement was considered.As a result, CCW is disturbed to varying degrees under the influence of human activities.This is of great significance for the accurate identification of the role of exogenous acids during the chemical weathering process in drainage basins in evaluating the interference of human activities with CCW and gaining an insight into the current global climate change and carbon cycle.However, the CCW processes involving strong anthropogenic acids have not yet received sufficient attention in the current global carbon cycle model [15].
Due to the slower dissolution kinetics and lower solubility of silicate than those of carbonate, the CO 2 consumption related to SCW in the karst watershed and clastic interlayers of carbonate rocks has been ignored, and all of the HCO 3 − was attributed to the contribution of CCW in previous studies [14,27,28].SCW consumes atmospheric/soil CO 2 to produce a net sink of atmospheric CO 2 over a geological time scale, and this cannot be neglected in the karst basins [29][30][31].In addition, the effects of anthropogenically sourced sulfuric acid and nitric acid are not considered in the common method of using hydrochemical compositions and runoff for the estimation of the karst carbon sink [6], and general studies only considered the effects of a single strong acid (sulfuric acid or nitric acid), which would affect the assessment accuracy of the CO 2 consumption related to CCW.
To accurately determine the compound influence of sulfuric acid and nitric acid resulting from human activities on the chemical weathering of carbonate basins and the riverine carbonate system, hydrogeochemical parameters of the surface runoff were investigated in the Lianjiang River Basin (LRB), which is in north of the Pearl River Delta (PRD) of China and under the dual influence of acidic gases transported from the PRD and the industrial regions in north of the Nanling Mountains.The main objectives were as follows: (1) to identify the hydrochemical characteristics of the Lianjiang River (LR) and the main sources of river dissolved loads; (2) to quantify the contribution proportion of SCW and CCW with different dissolution media to riverine DIC; (3) to evaluate the deficit proportion of CO 2 consumption caused by nitric acid and sulfuric acid weathering; and (4) to explicate the response characteristics of the riverine carbonate system to human activities.

Study Area
The LRB (112 • 10 E-113 • 18 E, 24 • 09 N-25 • 07 N) is located in the north of Guangdong Province, South China and at the southern edge of the important geographic and climatic dividing line, in the Nanling Mountains, East Asia.The LR belongs to the Pearl River Basin from the perspective of hydrology, and is located near the watershed between the Pearl River Basin and the Yangtze River Basin (Figure 1), with a mesoscale drainage area of 10,061 km 2 .The south end of the LRB is close to the Tropic of Cancer, and the climate type is a typical subtropical monsoon climate, with a humid climate, abundant rainfall between 1400 and 2200 mm (with an average of 1770 mm), and an average annual temperature of 19-20 • C. Most of the precipitation is observed from April to June, accounting for about 45-50% of rainfall throughout the whole year.After flowing from the headwater for 275 km, the LR flows into the Beijiang River at Jiangkouju, with an average gradient of 0.77%, and constitutes the largest tributary of the Beijiang River in the Pearl River System.The LR system belongs to the trellis drainage, extending from the northwest to the southeast of the drainage basin, and shows typical characteristics of a mountainous river.The main tributaries of the LR flow from the mountains at the edge of the drainage basin to the main stream (Table 1).The average annual runoff of the mainstream is 1.2 × 10 10 m 3 and mainly observed from March to August, accounting for about 70% of the total runoff.a µmol•L −1 and µeq•L −1 are expressed as µM and µEq, respectively; b Abbreviations of ML, UT, MT, DT and ST stand for the mainstream, upstream of tributaries, mid-stream of tributaries, downstream of tributaries and small tributaries of the LR, respectively.
The LRB is located in the east border of the boundary between the Yangtze Block and the Cathaysian Block, from the geological perspective [32], and is distributed in the top and west limb parts of the frontal arc in the Northern Guangdong ε-type Tectonic System (a complete example of the ε-type Tectonic System in China) [33].A large area of Carboniferous-Triassic limestone, dolomite and dolomitic limestone are outcropped in the central section of the LRB, which is characterized as a shallow marine carbonate formation and in which the limestone area accounts for more than 60% of the total basin area.There are Sinian metamorphic volcanic rocks, siliceous rocks and shales; Cambrian quartz sandstone and sandy shales; and Devonian purple and brown sandy shales distributed in the margin of upstream and mid-stream regions.In the former, the mineral composition is mainly potassium-sodium silicate, and there is a small amount of carbonate minerals, such as calcite and dolomite, in addition to silicate minerals in the latter.The Ordovician and Silurian strata in the LRB are missing, and the Jurassic and Cretaceous strata are scattered.A small fraction of Quaternary continental alluvial deposits and metamorphic rocks are observed in lower reaches (Figure 1).The terrains of the LRB are higher in the northwest and lower in the southeast with the main landform of karst plateau and basin.High mountains are distributed in the north and west, and karst basins and hilly areas are in the southeast.The karst plateau of the LRB is the eastern extension part of the karst area in Southwest China, which is one of the three largest concentrated karst areas in the world as well as the typical ecological vulnerable zone.The LRB is dominated by calcareous soil, which has developed as a consequence of the dissolution of carbonate rocks and was followed by red soil.There are numerous types of vegetation, including coniferous forests, evergreen broad-leaved forests, deciduous broad-leaved forests, mixed forests, bamboo forests and shrubs.Several cities, such as Lianzhou, Liannan, and Yangshan, are distributed in the upper and middle reaches of the mainstream (Figure 1).The upper reaches of the main tributaries are sparsely populated, yet the lower reaches are densely populated.Forest land accounts for the largest proportion of the drainage area (42%), and cultivated land accounts for about 10%.
The northern part of Guangdong Province has been seriously polluted by acid rain, and it has been a heavily acidic rain area (4.5 ≤ pH < 5.0 and acid rain frequency > 50%, Figure 1) in recent years [34].In addition to local pollutant gas emissions, acid rain in the LRB is also relevant to the long-distance transmission of pollutants, such as SO 2 and NO x from industrial activities of the economically developed PRD and the industrial districts in north of the Nanling Mountains [35].Rainfall events in the LRB are generally affected by the Quasi-stationary Front of South China, which is the main weather system regulating the rainfall in South China in winter and frequently stagnates between 22 • and 25 • N [36].It can easily form acid rain in rainfall events under the condition that acid substances are washed away by rainwater.The pH of acid rain in summer is higher than that in winter because of the large amount of rainfall from typhoons, which result in a high pH and low frequency of acid rain [37].Except for acid precipitation, local human activities, such as industrial and agricultural activities, will also affect the biogeochemical cycle in the drainage basin.It also deserves to be mentioned that biogeochemical processes in winter in the LRB proceed as usual because there are no hydrological freezing phenomena in the LR in January, which is annually the coldest month.

Sampling and Measurements
To fully understand the impact of human activities on the chemical weathering processes and riverine carbonate system in the LRB, 84 sampling sites in the mainstream and tributaries (including tertiary tributaries in headwater areas) were selected and sampled from 2 January to 8 January 2018.All of the water samples were collected in or near the middle channel of each river/stream, obviously avoiding pollution.Polyethylene sample bottles were fully rinsed using river water near the sampling sites more than three times before sampling and filled with river water ensuring that no air bubbles remain in the bottles.Samples were brought back to the laboratory as soon as possible and refrigerated at 4 • C. Alkalinity titration, ion content analysis and other measurements were completed within 48 h after sampling.The water temperature (T), pH and total dissolved solids (TDS) were tested in situ using an Ultrameter TM 6P portable conductivity meter (Myron L ® Company, Carlsbad, CA, USA), with a test accuracy of ±0.10 • C, ±0.01 and 0.01 mg•L −1 , respectively.The probe was calibrated before sampling and calibrated once during the sampling period according to the instructions.
Every water sample was divided into two parts.One part was filtered through a Whatman GF/F filter (GE Healthcare, Boston, MA, USA) (with a pore size of 0.45 µm), then used to determine the main cations and anions (K + , Na + , Ca 2+ , Mg 2+ , Cl − , NO 3 − , and SO 4 2− ) and dissolved silicon (expressed as DSi in this study).Samples for cations were acidified to pH < 2 with HNO 3 .The cations and DSi were measured using an IRIS (HR) type Inductively Coupled Plasma-Atomic Emission Spectrometric (ICP-AES, manufactured by TJA, Boston, MA, USA), and anions (Cl − , NO 3− , SO 4 2− ) were measured with ICS-900 Ion Chromatography, and the test accuracies were greater than 0.5%.The remaining unfiltered part of the sample was titrated to determine the total alkalinity (Alk.), according to Gran [38], with a Titralab-865 automatic titrator from the Hach Company (Loveland, CO, USA), the experimental error of which was less than 5%.Before titration, the pH meter was calibrated with the particular buffer solutions (with pH values of 4.005, 7.000 and 10.012).Using the titrant (HCl, 0.05 mol•L −1 ), the filtrated sample (100 mL) was then taken and titrated to pH = 4.5, 4.2, 3.8 and 3.5, and the corresponding HCl volumes were recorded.The final titration point was controlled by a glass electrode in the titration process [38].Each sample was titrated 3 times, and the sample alkalinity was calculated by the Graphics of Gran.The total alkalinity was approximately equal to the concentration of HCO 3 − , when the pH ranged from 6 to 10 in the water sample.

Data Processing
The DEM used in this study came from the digital elevation data of GDEMDEM 30 M resolution from the Geospatial Data Cloud website [39].The geological map used was based on the "Geological Map of Guangdong Province, Hong Kong and Macao Special Administrative Region" (1:1,750,000), and the drainage system and rock distribution maps were extracted using ArcGIS 10.0 (Environmental Systems Research Institute, Inc., Redlands, CA, USA).IBM SPSS Statistics 19.0 (IBM, Armonk, NY, USA) was used for the data statistics and the Pearson correlation analysis.Graphs were drawn using Grapher 12.0 software (Golden Software, LLC., Golden, CO, USA) and Origin Pro 2018 (OriginLab Corporation, Northampton, MA, USA).The dissolved inorganic carbon (DIC) and partial pressure of carbon dioxide (pCO 2 ) were calculated through the CO2SYS based on the pH, water temperature and alkalinity of the water obtained from the sampling and measuring period [2].The calcite saturation index (SIc) of the river water was calculated through the hydrogeochemical program of WATSPEC15 [28].

Physical and Chemical Parameters
The hydrochemical compositions of the mainstream and main tributaries in the LR are listed in Table 1.The temperature of the river water samples varies from 9.3 to 20.5 • C, with an average of 14.3 • C in January 2018.The minimum value appears at the sample 4 in the upper reach of the mainstream, and the maximum value appears in the Zhutianhe River, a tributary in the lower reach of the LR, showing the coupled effects of the southern subtropical monsoon climate and elevation on the water temperature.The range of pH in the LR varies from 6.79 to 9.02, with an average of 7.82, and most of the water samples are characterized as weakly alkaline, which indicates the dominant effect of the karst water on the pH of the LR.The variation range of the total dissolved solids (TDS) are from 22.54 to 357.5 mg•L −1 , with an average of 131.8 mg•L −1 , which is higher than the average of 100 mg•L −1 of most rivers worldwide [3], reflecting that the rock chemical weathering is relatively intense in the LRB.
Equivalent concentrations of total cations (TZ + = 2 Ca 2+ + 2 Mg 2+ + K + + Na + ) of all samples from the LR are between 326 µeq•L −1 and 5906 µeq•L −1 , with a mean of 2210 µeq•L −1 , which is higher than the average value of 61 rivers in the world [40].Equivalent concentrations of total anions (TZ − = HCO 3 − + 2SO  77) are not in equilibrium but vary from 10.25% to 18.29%, which is probably due to the presence of undetected anions, such as organic acids and phosphates, resulting from domestic sewage.

Cations and Anions
To conveniently clarify the hydrochemical type of the LR water and determine the main material sources of the dissolved loads, 84 water samples were categorized into five groups: the mainstream of the LR (ML), upstream of tributaries (UT), mid-stream of tributaries (MT), downstream of tributaries (DT) and small tributaries (ST) (Table 1).Ternary graphs of the main cations and anions in the LR were plotted according to the measured chemical compositions (Figure 2).Cations of the LR in January 2018 are dominated by calcium, accounting for 78.71% of the total cations, followed by magnesium, sodium and potassium, accounting for 9.35%, 6.92%, and 5.03%, respectively.The anions are dominated by bicarbonate, accounting for 84.97% of the total anions.Sulfate ion is the second largest anion, accounting for 8.98% of the total cations, followed by nitrate and chloride ions, accounting for 3.78% and 2.27%, respectively.Therefore, the hydrochemical type of the mainstream and tributaries in the LR belongs to the HCO 3 -Ca type, indicating that the typical karst water characterized by "higher calcium and carbon" dominates the river hydrochemical composition, which is consistent with the geological background, dominated by carbonate rocks in the LRB (Figure 1).equilibrium state of river water can be reflected by the inorganic charge balance index [NICB = (TZ + − TZ − )/TZ + × 100].The NICB of 72 out of the 84 samples are in the range from −10% to 10%, which could be explained by the fact that the anions and cations in the LR are generally in equilibrium, and those of the remaining 12 samples (12, 13, 17, 39, 40, 46, 54, 56, 65, 67, 71, and 77) are not in equilibrium but vary from 10.25% to 18.29%, which is probably due to the presence of undetected anions, such as organic acids and phosphates, resulting from domestic sewage.

Cations and Anions
To conveniently clarify the hydrochemical type of the LR water and determine the main material sources of the dissolved loads, 84 water samples were categorized into five groups: the mainstream of the LR (ML), upstream of tributaries (UT), mid-stream of tributaries (MT), downstream of tributaries (DT) and small tributaries (ST) (Table 1).Ternary graphs of the main cations and anions in the LR were plotted according to the measured chemical compositions (Figure 2).Cations of the LR in January 2018 are dominated by calcium, accounting for 78.71% of the total cations, followed by magnesium, sodium and potassium, accounting for 9.35%, 6.92%, and 5.03%, respectively.The anions are dominated by bicarbonate, accounting for 84.97% of the total anions.Sulfate ion is the second largest anion, accounting for 8.98% of the total cations, followed by nitrate and chloride ions, accounting for 3.78% and 2.27%, respectively.Therefore, the hydrochemical type of the mainstream and tributaries in the LR belongs to the HCO3-Ca type, indicating that the typical karst water characterized by "higher calcium and carbon" dominates the river hydrochemical composition, which is consistent with the geological background, dominated by carbonate rocks in the LRB (Figure 1).It is worth mentioning that the SO4 2− concentrations of Samples 28 and 82 are higher (1122 and 1114 μmol•L −1 ), accounting for 37.13% and 47.41% of the total anions, respectively.In addition, some proportions of (K + + Na + ) in the upstream of tributaries to the total cations are greater than 40%, and the DSi concentrations vary from 55.71 to 333.57μmol•L −1 , with an average of 159.22 μmol•L −1 (Table 1), which is slightly higher than the average concentration of DSi in the global rivers (158 μmol•L −1 ) [41].The DSi concentrations of Samples 25, 26, 48, 51, and 52 (303.93,317.86, 309.64, 333.57and 301.43 μmol•L −1 , respectively) are close to that of the Changhuajiang River (376.39 μmol•L −1 ), a silicatedominated basin on Hainan Island, China [42], indicating that the chemical compositions in the upstream of tributaries were characterized by the chemical weathering of silicate rocks.It is worth mentioning that the SO 4 2− concentrations of Samples 28 and 82 are higher (1122 and 1114 µmol•L −1 ), accounting for 37.13% and 47.41% of the total anions, respectively.In addition, some proportions of (K + + Na + ) in the upstream of tributaries to the total cations are greater than 40%, and the DSi concentrations vary from 55. 71 [42], indicating that the chemical compositions in the upstream of tributaries were characterized by the chemical weathering of silicate rocks.

Spatial Patterns of pCO 2 , SIc and DIC
The pCO 2 in the mainstream of the LR varies from 712 to 1624 µatm, with an average of 1025 ± 313 µatm, and those of the main tributaries are all supersaturated, relative to those of the atmosphere (407.54 µatm) [43], increasing from upstream to downstream by 1177, 1242 and 1394 µatm, respectively, which are lower than the average value of global rivers (2900 µatm) [44].The SIc spatial distribution characteristics of tributaries are similar to those of pCO 2 , and the SIc values from upstream to downstream in tributaries are −0.74,−0.51 and −0.03, respectively, showing an increasing trend from upstream to downstream, which reflects that the calcite saturation state of the river water varies from an unsaturated to a saturated state from upstream to downstream in the main tributaries of the LR.The SIc of the LR mainstream varies from −0.59 to 0.35, with an average of −0.11 ± 0.29, which is generally considered as a calcite unsaturation state.
The calculated DIC in the mainstream and tributaries vary from 263 to 4982 µmol•L −1 , with an average of 1740 µmol•L −1 .HCO 3 − is the most abundant dissolved carbonate in river water, and its proportion to DIC varies from 70.24% to 98.06%, with an average of 94.6% (Table 2).
The DIC contents in tributaries gradually increase from upstream to downstream (1461, 1645 and 1932 µmol•L −1 , respectively).37 07" N) is only 158 km.Therefore, the LRB is obviously affected by the warm and humid air mass from the South China Sea and the Western Pacific Ocean, as well as that from the Indian Ocean.Cl − is the main component of rainwater, with a conservative behavior during the water cycle process, and is gradually used as a reference element to estimate the contribution of sea salt precipitation to the river hydrochemical composition [45].It is generally believed that the cyclic chloride in surface water is mainly derived from sea salt precipitation, evaporite dissolution and human activities (i.e., the common salt in domestic sewage) [46,47].
There are hardly any saline formations but distributes several cities and towns at the intersections of rivers in the LRB.Sample 51, which was collected in a small tributary of the Qinglianshui River flowing through a single lithology (granite) region that is not influenced by salt formation and human activities, has the lowest value of Cl − concentrations (11.27 µmol•L −1 ) of all samples, and it can be assumed that the Cl − of Sample 51 completely originates from sea salt precipitation.Therefore, the Cl − concentration of Sample 51 is used to represent the Cl − concentration in the rainfall of the LRB.Sea salt correction is performed with reference to the concentration ratios of other ions in standard seawater, where Cl − /Na + = 1.15,Ca 2+ /Na + = 0.02, Mg 2+ /Na + = 0.11, and HCO 3 − /Na + = 0.004 [3].
The calculated proportions of sea salt precipitation and human activities to Cl − concentrations are 23.23% and 76.77%, respectively.Correlation analysis between ions shows that there is a significant positive correlation between K + , Na + , and Cl − (Table 3), indicating that there is a consistent source for K + , Na + and Cl − , which may result from sea salt precipitation or domestic sewage discharge.Furthermore, the correlation coefficient R between Na + and Cl − is 0.472, and that of K + and Cl − is 0.813 at a 0.01 test level, which indicates that the sources of K + and Cl − are more consistent and largely due to the sea salt precipitation.However, the fluctuation ranges of Na + /Cl − molar ratios are relatively great, from 0.55 to 16.86, and the mean value of 2.42 is greater than that of standard seawater (0.86) [48].The excessive Na + indicates that there are other sources with regard to Na + , such as the chemical weathering of silicate minerals (e.g., sodium feldspar) and common salt in domestic wastewater.

Carbonate Chemical Weathering
Ca 2+ and HCO 3 − are dominant in the river water ions of the LR (Figure 2).It is found that there are significant positive correlations between Ca 2+ and Mg 2+ , HCO 3 − , TDS (Table 3), with correlation coefficients of 0.764, 0.961 and 0.981 (p < 0.01), which indicates that there is a consistent geological source of CCW for Ca 2+ , Mg 2+ and HCO 3 − .Therefore, it has an important controlling effect for CCW on the geochemistry of the LR, which is consistent with the geological background of carbonate rocks, which is the main type of the basin bedrock (the distribution area of carbonate terrains accounts for more than 60% of the total basin area).(Ca 2+ + Mg 2+ ):HCO 3 − molar ratios, induced by calcite and dolomite dissolution, can be expressed as Equations ( 1) and (2) [49]: The relative contribution of calcite or dolomite dissolution to river solutes can be judged by the molar ratios of Mg/Ca and Mg/HCO 3 due to the fact that Mg 2+ is mainly derived from the dissolution of dolomite minerals in drainage basins, and carbonates are predominant and hardly affected by the dissolution of gypsum.It can be observed that the molar ratio of Mg/Ca, produced by calcite dissolution, is less than 0.1, while that of dolomite dissolution is equal to 1 (Equations ( 1) and ( 2)).The equal dissolution of calcite and dolomite yields a Mg/Ca ratio of 0.33 and Mg/HCO 3 ratio of 0.125, and the resultant Mg/HCO 3 ratios are 0.08 and 0.16, which indicates that the amount of dolomite dissolution is half or twice that of calcite based on mass balance, respectively [50,51].The calculated average values of Mg/Ca and Mg/HCO 3 in the LR are 0.18 (0.1-0.2) and 0.09 (0.05-0.11), respectively, indicating that the hydrochemical compositions of the LR, induced by CCW, should be attributed to the compound contributions of calcite and dolomite dissolution.Furthermore, calcite dissolution is the dominant source, which is also consistent with the large-scale exposure of limestone in the basin, especially in the area where the mainstream of the LR drains.

Chemical Weathering of Silicates
It has shown that the contribution of silicate weathering to CO 2 consumption has often been ignored in previous studies of karst carbon sink [14,27,28,52,53], which may lead to an overestimation of the CCW rate and CO 2 consumption associated with CCW in karst watersheds.The contribution of SCW to riverine DIC in karst watersheds should not be ignored.The exogenous water from the non-carbonate terrains around the karst watershed plays an important role in the karst process.For instance, the exogenous water with a relatively low pH and SIc will result in the severe dissolution of carbonate rocks in karst areas, which will greatly strengthen the potential of karst carbon sink [54][55][56].Ford and William [57] found that the amount of carbonates dissolved in the CCW processes will be enhanced by 50% when the recharge area of exogenous water accounts for 50%.On the other hand, clastic interlayers are often observed in carbonate formations and in which, the export of HCO 3 − induced by silicate weathering cannot simply be attributed to the contribution of CCW.
Although the rate of SCW is relatively slow, it consumes atmospheric/soil CO 2 on a geological time scale (>1 Ma), and which process constitutes the natural thermostat of the Earth [58].The proportion of non-carbonate formations in the LRB is more than 30%, according to the computational geometry module in ArcMap (Environmental Systems Research Institute, Inc., Redlands, CA, USA).The silicate formations are mainly distributed at the edge of the river basin, in other words, in the upstream of the mainstream and main tributaries, as well as in the lower reaches of south bank of the LR, including igneous rocks, metamorphic rocks, clastic rocks and Quaternary alluvial sediments (Figure 1), and the former dominates the non-carbonate formations.
It is generally known that SCW is one of the important natural sources of K + and Na + , for example, the chemical weathering of silicate rocks, such as granites and metamorphic rocks, can produce Na + , and that of K-feldspar and mica can produce K + .Being less polluted by human activities in the upper reaches of the LR's main tributaries, the average concentrations of Na + and K + are 113.86 and 53.12 µmol•L −1 , respectively, and their proportions to the total cations, respectively, account for 19.6% and 7.07%.With the dilution effect of the karst water input (poor in Na + and K + ), this shows a decreasing tendency in the proportions of Na + and K + to the total cations in the midstream of tributaries to the downstream, which are 14.91% and 5.32% in the middle reaches, respectively, and 13.62% and 5.06% in the lower reaches.
The DSi in the runoff mainly originates from SCW.The average concentration of DSi in the tributaries of the LR is 138.15 µmol•L −1 , which is slightly lower than the average value of DSi in the global rivers (158 µmol•L −1 ) [41] and the upper reaches of the Sorocaba River in Southeastern Brazil, where igneous and metamorphic rocks predominate (240.75 µmol•L −1 ) [59].It is close to that of the Ganges (134 µmol•L −1 ) [60], but higher than that of the Yellow River (95.29 µmol•L −1 ) [61], the Sanchahe River (80-110 µmol•L −1 , upstream of the Wujiang River) [62] and the Kochechum River (87.5 µmol•L −1 ) in the permafrost region of Central Siberia [63].This indicates that the ratio of silicate formation to the drainage area dominates the output of DSi in the river basin, and it also shows that the SCW process in the LRB cannot be ignored from the perspective of the DSi concentration.
It was reported above that sea salt precipitation is the main source of K + and Na + in the LR.In addition, SCW in non-carbonate formations and carbonate interlayers also contribute significantly to K + and Na + .The ratio of (Na + + K + ) and (Na + + K + + Cl − ) is used to determine the relative contributions of SCW to the riverine hydrochemical composition [64].Calculations show that (Na + + K + )/(Na + + K + + Cl − ) varies from 0.47 to 0.95, with an average of 0.71.It can be assumed that the products of SCW relatively contribute 71% of K + and Na + in the LR, without considering other sources of K + and Na + .
There is a positive correlation between DSi and (Na + + K + )/(Na + + K + + Cl − ) for the LR (R 2 = 0.41, p < 0.01) (Figure 3), which shows that the DSi content controlled by SCW increases simultaneously with the ratio of (Na + + K + )/(Na + + K + + Cl − ), indicating that the significant contribution of SCW in non-carbonate formations to K + and Na + is further confirmed in the karst basin of the LR.The lower correlation coefficient of DSi and (Na + + K + )/(Na + + K + + Cl − ) may be attributed to the DSi uptake from diatoms in the case of slower velocities, as in January [65].It consumes atmospheric/soil CO 2 during the SCW process and produces HCO 3 − while releasing Na + and K + .The contents of HCO 3 − , Na + and K + should simultaneously increase.It can be seen in Table 3 that there is a significant positive correlation between HCO 3 − and K + at 0.01 levels, with a correlation coefficient of 0.342.However, there is no correlation between HCO 3 − and Na + , and the correlation coefficient is only 0.059, which indicates that Na + also has other sources, except for sea salt precipitation and SCW, which is probably related to the discharge of domestic sewage from the residents of cities and towns near the riverbank.
the ratio of (Na + + K + )/(Na + + K + + Cl − ), indicating that the significant contribution of SCW in noncarbonate formations to K + and Na + is further confirmed in the karst basin of the LR.The lower correlation coefficient of DSi and (Na + + K + )/(Na + + K + + Cl − ) may be attributed to the DSi uptake from diatoms in the case of slower velocities, as in January [65].It consumes atmospheric/soil CO2 during the SCW process and produces HCO 3 − while releasing Na + and K + .The contents of HCO 3 − , Na + and K + should simultaneously increase.It can be seen in Table 3 that there is a significant positive correlation between HCO 3 − and K + at 0.01 levels, with a correlation coefficient of 0.342.However, there is no correlation between HCO 3 − and Na + , and the correlation coefficient is only 0.059, which indicates that Na + also has other sources, except for sea salt precipitation and SCW, which is probably related to the discharge of domestic sewage from the residents of cities and towns near the riverbank.

Human Activities
Human activities, including atmospheric acid precipitation, urban and industrial wastewater, agricultural fertilization and other sources of pollution will affect river hydrochemical composition [66].Acid precipitation and industrial wastewater mainly impact SO4 2− and NO3 − in surface runoff.Na + , Cl − and NO3 − will increase if the domestic sewage drains into the river system.The application of potassium and nitrogen fertilizer in agriculture activities results in non-point source pollution and increases the concentrations of K + and NO3 − in the river water, respectively [66,67].
The LRB is influenced by a large amount of SO2 and NOx from the industrial districts in north of the Nanling Mountains and the PRD, which is one of the most developed areas in East Asia.The middle and long distance transmission of acid substances has a profound impact on the northern part of Guangdong Province.The contribution rate of acid substances, transported from external lands to local acid rain, can reach 49%, according to previous research on rainfall by stable isotope techniques [68].Therefore, the urbanization/industrialization of the PRD and industrial districts in north of the Nanling Mountains is a non-negligible factor that results in a low pH of rainfall in the LRB.Several cities, such as Lianzhou, Liannan, and Yangshan, are distributed in the upper and middle reaches of the LR, where the population density is relatively high.In addition, it is inevitable that agricultural activities, scattered near the main tributaries, will also influence the chemical composition of the LR.
The concentrations of SO4 2− and NO3 − in the mainstream of the LR (Sample 3) are 165.63 and 91.94 μmol•L −1 , respectively, which is close to those of the Qiantangjiang River (223 and 146 μmol•L −1 ), which is affected by severe acid precipitation [26], but much less than those of the upper reaches in

Human Activities
Human activities, including atmospheric acid precipitation, urban and industrial wastewater, agricultural fertilization and other sources of pollution will affect river hydrochemical composition [66].Acid precipitation and industrial wastewater mainly impact SO 4 2− and NO 3 − in surface runoff.Na + , Cl − and NO 3 − will increase if the domestic sewage drains into the river system.The application of potassium and nitrogen fertilizer in agriculture activities results in non-point source pollution and increases the concentrations of K + and NO 3 − in the river water, respectively [66,67].
The LRB is influenced by a large amount of SO 2 and NO x from the industrial districts in north of the Nanling Mountains and the PRD, which is one of the most developed areas in East Asia.The middle and long distance transmission of acid substances has a profound impact on the northern part of Guangdong Province.The contribution rate of acid substances, transported from external lands to local acid rain, can reach 49%, according to previous research on rainfall by stable isotope techniques [68].Therefore, the urbanization/industrialization of the PRD and industrial districts in north of the Nanling Mountains is a non-negligible factor that results in a low pH of rainfall in the LRB.Several cities, such as Lianzhou, Liannan, and Yangshan, are distributed in the upper and middle reaches of the LR, where the population density is relatively high.In addition, it is inevitable that agricultural activities, scattered near the main tributaries, will also influence the chemical composition of the LR.
The concentrations of SO 4 2− and NO 3 − in the mainstream of the LR (Sample 3) are 165.63 and 91.94 µmol•L −1 , respectively, which is close to those of the Qiantangjiang River (223 and 146 µmol•L −1 ), which is affected by severe acid precipitation [26], but much less than those of the upper reaches in the Wujiang River (590 and 200 µmol•L −1 ), which is drained in the Guizhou karst plateau and is significantly influenced by industrial activities and acid rain [53].The ratio between SO 4 2− /Ca 2+ and NO 3 − /Ca 2+ is generally used to determine the anthropogenic sources of SO 4 2− and NO 3 − .NO 3 − /Ca 2+ is not correlated with SO 4 2− /Ca 2+ in the mainstream and major tributaries of the LR (Figure 4), indicating that there are different sources for SO 4 2− and NO 3 − .Several samples have a relatively high SO 4 2− /Ca 2+ (>0.3) and low NO 3 − /Ca 2+ (<0.5), indicating that these samples are more subject to wastewater, which is relevant to industry and mining in the upstream of some tributaries, while another small group of samples from the upstream and downstream of some tributaries have a lower SO 4 2− /Ca 2+ (<0.3) and higher NO 3 − /Ca 2+ (>1.5) and are connected with domestic sewage discharge, which can be demonstrated by the human activities of town residents in karst depressions and karst valleys within the LRB.As the sampling was conducted during a period with hardly any agricultural activities (January), the effect of nitrogen application on NO 3 − was neglected in this study.Most of the samples have a relatively low SO 4 2− /Ca 2+ (average value is 0.16) and NO 3 − /Ca 2+ (average value is 0.68) (Figure 4), therefore, SO 4 2− and NO 3 − in the river water are largely due to acid precipitation, and this phenomenon is consistent with the fact that the LRB is located in a severe acid rain area in East Asia (Figure 1).
the Wujiang River (590 and 200 μmol•L −1 ), which is drained in the Guizhou karst plateau and is significantly influenced by industrial activities and acid rain [53].The ratio between SO 4 2− /Ca 2+ and NO 3 − /Ca 2+ is generally used to determine the anthropogenic sources of SO 4 2− and NO 3 − .NO 3 − /Ca 2+ is not correlated with SO 4 2− /Ca 2+ in the mainstream and major tributaries of the LR (Figure 4), indicating that there are different sources for SO /Ca 2+ (average value is 0.16) and NO 3 − /Ca 2+ (average value is 0.68) (Figure 4), therefore, SO 4 2− and NO 3 − in the river water are largely due to acid precipitation, and this phenomenon is consistent with the fact that the LRB is located in a severe acid rain area in East Asia (Figure 1).] are distributed between the endmember of carbonic acid-involved CCW and nitric acid and sulfuric acid-involved weathering in Figure 6, which demonstrates that nitric acid and sulfuric acid play important roles in CCW, and once again proves that carbonic acid is not the only weathering medium in the LRB from this distribution pattern.produced by CCW associated with sulfuric acid and nitric acid, is 2, and that of distributed between the endmember of carbonic acid-involved CCW and nitric acid and sulfuric acid-involved weathering in Figure 6, which demonstrates that nitric acid and sulfuric acid play important roles in CCW, and once again proves that carbonic acid is not the only weathering medium in the LRB from this distribution pattern.] are distributed between the endmember of carbonic acid-involved CCW and nitric acid and sulfuric acid-involved weathering in Figure 6, which demonstrates that nitric acid and sulfuric acid play important roles in CCW, and once again proves that carbonic acid is not the only weathering medium in the LRB from this distribution pattern.Based on the above analysis, the contribution of sea salt precipitation and SCW to DIC must be subtracted from the total riverine DIC to calculate the effect of acid precipitation on DIC.Sea salt correction (see Section 4.1.1)was used to evaluate the contribution of a sea salt precipitation endmember to riverine DIC in this work.The calculated contribution proportion of sea salt precipitation to DIC concentration is very small, only accounting for 0.01%, which is because the pH of rainwater in the LRB is below 5.6 (Figure 1), and the pCO 2 in the river water is generally higher than that in the atmosphere.
The hydrochemical characteristics of the LR are mainly controlled by the geological background of the watershed.The chemical weathering of carbonate has a dominant role in riverine DIC.Additionally, SCW cannot be ignored.The reactions of SCW followed Mortatti and Probst [72], which are shown as follows: Since SCW associated with sulfuric acid and nitric acid cannot produce HCO 3 − , the carbonic acid-involved SCW was only considered to estimate the CO 2 consumption of SCW in this paper.
According to Equations ( 4) and ( 5), the HCO 3 − generated from SCW can be calculated as half of the DSi concentration, and the ions of Ca 2+ and Mg 2+ , derived from SCW, can be calculated as 1/5 of the DSi concentration [73].The calculated contribution proportions of SCW to riverine DIC range from 0.59% to 54.89%, with a mean value of 9.64%.The average DIC produced by SCW is 49.97 µmol•L −1 , indicating that the SCW plays a relatively important role in the DIC concentration of the LR.

CO 2 Consumption
Deficit Related with Anthropogenic Acid-Involved CCW Carbonic acid, nitric acid and sulfuric acid were involved in CCW in the LRB.To estimate the relative contribution to riverine DIC of CCW induced by nitric and sulfuric acid, we assumed that k 1 (mol) carbonic acid, k 2 (mol) nitric acid and k 3 (mol) sulfuric acid were involved in the CCW process, rewriting Equation (3) as [53,74]: Afterwards, the molar concentrations of (Ca 2+ + Mg 2+ ) and HCO 3 − produced by CCW in the LRB are: Thus: In the above equations, the subscripts of "riv", "carb", "sil" and "cyc" represent the river water, carbonate weathering, silicate weathering and cyclic salt, respectively.
According to the above equations, the concentration of carbonic acid (k 1 ) involved in CCW can be calculated by subtracting the corresponding values of cyclic salt (sea salt precipitation) and SCW from the measured HCO 3 − , Ca 2+ , and Mg 2+ concentration data at each sampling point (Equations ( 7)-( 11)).
The difference between the residual amount of HCO 3 − after deducting those from the sea salt precipitation and SCW in river water and 2k 1 (the molar concentration of HCO 3 − derived from CCW by carbonic acid is double that of carbonic acid) is the sum of k 2 and 2k 3 (the values of k 2 and k 3 cannot be calculated separately).Therefore, the k 1 value can be used to calculate the molar concentrations of carbonic acid that participates in the CCW and the DIC concentrations from each endmember, respectively, and then the deficit proportion of CO 2 consumption induced by nitric acid and sulfuric acid can also be calculated according to the above results.As mentioned in Section 3.2, SO 4 2− is abnormally high for Samples 28 and 82.To prevent the calculated results from significantly deviating from the actual ones, these two samples were not considered here.The calculation results show that the relative proportions of carbonic acid and anthropogenic acids (nitric acid and sulfuric acid) in CCW are 76.97% and 23.03%, respectively.DIC from carbonic acid weathering accounts for 76.79% of the total DIC in rivers and streams, of which 38.4% is derived from atmospheric/soil CO 2 , and the deficit proportion of CO 2 consumption, associated with CCW resulting from anthropogenic acids, is approximately 13.56%.Thus, the proportion attributed to acid substances from other areas outside the LRB is 6.64% (see Section 4.1.4),and the above results show that the exogenous acid sourced from human activities has a non-negligible interference with the water-rock interaction and the geochemical cycle of carbon across the basin.The calculation of the deficit proportion is slightly lower than the loss of CO 2 consumption caused by nitric acid and sulfuric acid, in both the rainy and pre-rainy seasons in the Guizhou Plateau (15.67% and 14.17%, respectively) [53], which is mainly attributed to the severe acid rain pollution, frequent industrial activity, and high percentage of carbonate rock distribution.The output of DIC, caused by anthropogenic acids in the LR, represents 38.44% of the CO 2 consumption associated with natural carbonate weathering, and the proportion will reach 18.84% for the external acids, both of which are higher than that of SCW associated with nitric acid in France (6-15%) [13].Thus, chemical weathering influenced by human activities in the carbonate basin, especially in which affected by the heavy acid rain, need to be paid more attention.

Response of Riverine Carbonate System to Human Activities
It consumes CO 2 from the atmosphere/soil in the chemical weathering process and directly releases the carbon sequestered in carbonate rocks during CCW and then, the yielded HCO 3 − are transported to the river system through surface runoff.Therefore, the chemical weathering of rocks will directly influence the riverine carbonate system [75].A correlation analysis between DIC riv and the parameters related with the riverine carbonate system in the LR shows that riverine TDS, pCO 2 , alkalinity and SIc all increase with the DIC riv concentration, and there are extremely significant positive correlations among those parameters, with correlation coefficients of 0.975, 0.633, 0.999 and 0.8, respectively, p < 0.01 (Table 4).Carbon dioxide, as an important driving force of the occurrence of the karst process, will promote the dissolution of carbonate rocks after entering the karst dynamic system, which represents an increase of Ca 2+ , Mg 2+ and HCO 3 − , and an enhancement of pCO 2 .On the contrary, the escape of CO 2 from the karst system may result in the carbonate precipitation, and the HCO 3 − concentration and pCO 2 in river water will reduce [76].Moreover, a constantly increasing HCO 3 − promotes the combination of Ca 2+ and Mg 2+ with HCO 3 − in the river water, which will lead the saturated state of calcite to be over-saturated and provide favorable conditions for calcite precipitation, thus causing a rise in pCO 2 .The correlation between the DIC riv and pH in the LR water is relatively weak, with a correlation coefficient of 0.259 (p < 0.05), and there is no correlation between the water temperature and DIC riv .However, there is a weak positive correlation between the water temperature and DIC C , with a correlation coefficient of 0.249 (p < 0.05, Table 4), showing that temperature constrains the chemical weathering of carbonate rocks and have an impact on the output of DIC.Within a certain range of temperature conditions, higher temperatures give rise to faster dissolution and yields more HCO 3 − [77].
In addition, it can also reflect that there are different sources of DIC in the LR by the weak correlation between the water temperature and DIC C and the uncorrelated relationship between water temperature and DIC N−S , and the CCW associated with sulfuric acid and nitric acid derived from human activities, which is discussed in this paper, is one of the important sources of riverine DIC.The participation of nitric acid and sulfuric acid in the CCW results in a weakening of carbon sequestration capacity induced by chemical weathering.In addition, it promotes the activation of carbon sequestered in carbonate rocks [78,79].HCO 3 − , which derived from the product of chemical weathering, contributes significantly to the riverine DIC, which constitutes a non-ignorable part of the river carbon flux that transports to the ocean.There is a significant positive correlation between the DIC produced in CCW induced by carbonic acid in the LR and the total riverine DIC, with a correlation coefficient of 0.984, and the DIC produced by anthropogenic acid-involved CCW is also positively correlated with the total DIC in river waters (R = 0.787, p < 0.01) (Table 4 and Figure 7a,b), which shows that the riverine DIC is mainly constrained by the mechanism of the CCW process induced by carbonic acid, followed by nitric acid and sulfuric acid.The DIC derived from the latter is also important for the riverine DIC in the LR.The scattered data points of DIC N−S /DIC riv show that carbonate systems in different regions of the LRB are driven by varying degrees of weathering processes associated with nitric acid and sulfuric acid (Figure 7b).DICN−S is also significantly positively correlated with the major parameters related with the riverine carbonate system, such as TDS, pCO2, alkalinity and SIc, with the correlation coefficients of 0.853, 0.517, 0.785 and 0.615 (p < 0.01) (Table 4).The logarithmically increasing relationship between DICN−S and SIc (R = 0.615, p < 0.01) also indicates that the carbonate system of the LR responds positively to CCW induced by anthropogenic nitric acid and sulfuric acid (Figure 7c,d).As CCW influenced by anthropogenic acids continues to occur, carbon sequestered in carbonate rocks will continuously release into the river system, and DIC derived from CCW will also increase.The saturation state of calcite will also be gradually transformed from being unsaturated to saturated, or even oversaturated accompanied with the increase of DICN−S.The increase in the saturation of calcite in river water will drive CO2 release to the atmosphere through the water-air interface, prompting the transition of the river to a carbon source [80], and, in turn, affecting the source and sink effects of drainage basins.The associated event should be considered in the estimation of the CO2 consumption induced by chemical weathering.Therefore, a systematic view of the river carbon transport process is crucial for the accurate estimation of the CO2 consumption in the drainage basin.
There is a weak exponentially increasing relationship between DICN−S and pCO2 (R = 0.517, p < 0.01), which reflects that pCO2 may be affected by other factors, except for the effects of the CCW, for instance, the terrestrial CO2 produced by root respiration will enter the river system through surface runoff, subsurface flow and groundwater [81,82]; in addition, pCO2 can also increase due to the precipitation of calcite in supersaturated water and the decomposition of organic carbon in river sediments and river waters [83,84].In brief, chemical weathering caused by anthropogenic sulfuric acid and nitric acid is also the critical factor driving the increase of SIc and pCO2 in the LR and constrains the carbonate system of the karst river, together with the chemical weathering involved with carbonic acid.

Conclusions
Anions and cations in the mainstream and main tributaries of the LR were dominated by HCO3 − and Ca 2+ in January 2018, accounting for 84.97% of the total anions and 78.71% of the total cations, DIC N−S is also significantly positively correlated with the major parameters related with the riverine carbonate system, such as TDS, pCO 2 , alkalinity and SIc, with the correlation coefficients of 0.853, 0.517, 0.785 and 0.615 (p < 0.01) (Table 4).The logarithmically increasing relationship between DIC N−S and SIc (R = 0.615, p < 0.01) also indicates that the carbonate system of the LR responds positively to CCW induced by anthropogenic nitric acid and sulfuric acid (Figure 7c,d).As CCW influenced by anthropogenic acids continues to occur, carbon sequestered in carbonate rocks will continuously release into the river system, and DIC derived from CCW will also increase.The saturation state of calcite will also be gradually transformed from being unsaturated to saturated, or even oversaturated accompanied with the increase of DIC N−S .The increase in the saturation of calcite in river water will drive CO 2 release to the atmosphere through the water-air interface, prompting the transition of the river to a carbon source [80], and, in turn, affecting the source and sink effects of drainage basins.The associated event should be considered in the estimation of the CO 2 consumption induced by chemical weathering.Therefore, a systematic view of the river carbon transport process is crucial for the accurate estimation of the CO 2 consumption in the drainage basin.
There is a weak exponentially increasing relationship between DIC N−S and pCO 2 (R = 0.517, p < 0.01), which reflects that pCO 2 may be affected by other factors, except for the effects of the CCW, for instance, the terrestrial CO 2 produced by root respiration will enter the river system through surface runoff, subsurface flow and groundwater [81,82]; in addition, pCO 2 can also increase due to the precipitation of calcite in supersaturated water and the decomposition of organic carbon in river sediments and river waters [83,84].In brief, chemical weathering caused by anthropogenic sulfuric acid and nitric acid is also the critical factor driving the increase of SIc and pCO 2 in the LR and constrains the carbonate system of the karst river, together with the chemical weathering involved with carbonic acid.

Conclusions
Anions and cations in the mainstream and main tributaries of the LR were dominated by HCO 3 − and Ca 2+ in January 2018, accounting for 84.97% of the total anions and 78.71% of the total cations, respectively.The average concentration of DSi in the drainage basin (159.22 µmol•L −1 ) was slightly higher than that of global rivers (158 µmol•L −1 ).Therefore, SCW in karst basins should not be ignored.HCO 3 − was the most abundant dissolved carbonate species in the LR, accounting for an average of 94.6% of DIC.Most of the dissolved loads originated from mineral weathering, dominated by calcite dissolution.SO 4 2− and NO 3 − were mainly derived from acid precipitation, induced by acid substances from external transmission (from the PRD and industrial districts in north of the Nanling Mountains) and local supply, but the contribution from local industrial wastewater and domestic sewage was less.
The SCW contributed an average of 9.64% to riverine DIC.Sulfuric acid and nitric acid from human activities were the weathering agents that participated in the CCW.The relative proportions of carbonic acid and anthropogenic acids, participating in CCW, were 76.97% and 23.03%, respectively.The deficit proportion of CO 2 consumption, associated with CCW resulting from sulfuric acid and nitric acid in winter, was approximately 13.56%, which was slightly lower than those of the Guizhou Plateau in the rainy and pre-rainy seasons (15.67% and 14.17%, respectively).The deficit percentage of CO 2 consumption by nitric acid and sulfuric acid to the CO 2 consumption associated with natural CCW was found to be 38.44% and 18.84% for anthropogenic acids from external areas, and these values may be lower than the actual values because the utilization of DIC by aquatic photosynthetic organisms was not considered in this study.In addition, the relative proportion of sulfuric acid and nitric acid, which is involved in CCW, needs to be further identified.
The DIC produced by CCW associated with nitric acid and sulfuric acid showed a significant positive correlation with riverine DIC, and there was a logarithmically increasing relationship with SIc and an exponentially increasing relationship with pCO 2 , indicating that the carbonate system of the LR was driven by human activities.The exogenous acids derived from human activities had a non-negligible interference with the water-rock interaction and the geochemical cycle of carbon across the basin.Therefore, the chemical weathering and riverine carbonate system, influenced by anthropogenic acids in the carbonate basin, should be paid more attention.

Figure 1 .
Figure 1.Characteristics of location, sampling sites, river system, geological background, terrains and rainfall pH in the LRB.Nos.1-84 stand for the serial number of 84 sampling sites distributed in the Lianjiang River System.A-N stand for the mainstream and main tributaries of the LR: mainstream of the LR (A), the Xingzihe River (B, XZHR), the Baoanhe River (C, BAHR), the Chaotianhe River (D, CTHR), the Dongbeihe River (E, DBHR), the Sanjiangshui River (F, SJSR), the Dongguanshui River (G, DGSR), the Guanbeihe River (H, GBHR), the Qigonghe River (I, QGHR), the Qinglianshui River (J, QLSR), the Shuibianhe River (K, SBHR), the Datanhe River (L, DTHR), the Jintanhe River (M, JTHR), and the Zhutianhe River (N, ZTHR).The Pearl River Delta and the industrial districts in north of the Nanling Mountains are represented as lowercase letters a and b, given in the upper left figure.

Figure 1 .
Figure 1.Characteristics of location, sampling sites, river system, geological background, terrains and rainfall pH in the LRB.Nos.1-84 stand for the serial number of 84 sampling sites distributed in the Lianjiang River System.A-N stand for the mainstream and main tributaries of the LR: mainstream of the LR (A), the Xingzihe River (B, XZHR), the Baoanhe River (C, BAHR), the Chaotianhe River (D, CTHR), the Dongbeihe River (E, DBHR), the Sanjiangshui River (F, SJSR), the Dongguanshui River (G, DGSR), the Guanbeihe River (H, GBHR), the Qigonghe River (I, QGHR), the Qinglianshui River (J, QLSR), the Shuibianhe River (K, SBHR), the Datanhe River (L, DTHR), the Jintanhe River (M, JTHR), and the Zhutianhe River (N, ZTHR).The Pearl River Delta and the industrial districts in north of the Nanling Mountains are represented as lowercase letters a and b, given in the upper left figure.

Figure 2 .
Figure 2. Ternary graphs of cations and anions of the LR in January 2018.

Figure 2 .
Figure 2. Ternary graphs of cations and anions of the LR in January 2018.
(<0.5), indicating that these samples are more subject to wastewater, which is relevant to industry and mining in the upstream of some tributaries, while another small group of samples from the upstream and downstream of some tributaries have a lower SO 2+ (>1.5) and are connected with domestic sewage discharge, which can be demonstrated by the human activities of town residents in karst depressions and karst valleys within the LRB.As the sampling was conducted during a period with hardly any agricultural activities (January), the effect of nitrogen application on NO 3 − was neglected in this study.Most of the samples have a relatively low SO 4 2−

4. 2 .
Evidence of Chemical Weathering Associated with Sulfuric and Nitric Acids Equivalent data points of [Ca 2+ + Mg 2+ ] and [HCO 3 − ] in the LR deviate from the 1:1 line (Figure 5a), with most situated above the 1:1 line, indicating that extra anions are necessary to balance the excess Ca 2+ and Mg 2+ .Nevertheless, those of [Ca 2+ + Mg 2+ ] and [HCO 3 − + SO 4 2− + NO 3 − ] areclose to the 1:1 line (Figure5b), of which average ratio is 0.96.Therefore, the equilibrium between the anions and cations in the LR requires SO 4 2− and NO 3 − to compensate.

Figure 7 .
Figure 7. Relationship between DIC from different sources and the relationship among DICN−S, SIc and pCO2.

Figure 7 .
Figure 7. Relationship between DIC from different sources and the relationship among DIC N−S , SIc and pCO 2 .

Table 1 .
Main physical and chemical parameters of the LR water in January 2018.

Table 2 .
Variations and average values for pCO 2 , SIc, DIC and other parameters in the mainstream and tributaries of the LR.Atmospheric precipitation is the dominant recharge source of the surface runoff in the LRB in the context of a humid subtropical monsoon climate.The distance between the southernmost border of the LRB (113 • 05 04" E, 23 • 55 33" N) and the Pearl River Estuary (113 • 42 06" E, 22

Table 3 .
Correlations of the main ion concentrations in the LR in January 2018.
** Correlation is significant at the 0.01 level (two-tailed); * Correlation is significant at the 0.05 level (two-tailed).
4.3.Chemical Weathering and the Related CO 2 Consumption Driven by Anthropogenic Acids 4.3.1.Influence from Sea Salt Precipitation and SCW

Table 4 .
Correlation analysis of parameters related with carbonate system in the LR.* Correlation is significant at the 0.01 level (two-tailed); * Correlation is significant at the 0.05 level (two-tailed); DIC riv , DIC C and DIC N−S represent the total amount of DIC in water samples, DIC derived from carbonic acid-involved CCW and DIC derived from nitric acid and sulfuric acid-involved CCW, respectively. *