Spatiotemporal Dynamics of Dissolved Organic Carbon and Freshwater Browning in the Zoige Alpine Wetland, Northeastern Qinghai-Tibetan Plateau

: The concentrations of dissolved organic carbon (DOC) and its light-absorbing fraction (chromophoric dissolved organic matter; CDOM) in surface waters, particularly those draining organic-rich peatlands, have dramatically increased over the past decade due to climate change and human disturbance. To explore the spatiotemporal dynamics of DOC and CDOM in surface waters of the northeastern Qinghai-Tibetan Plateau, we collected water samples from two rivers in the Zoige alpine wetland and from two rivers in its adjacent alpine-gorge region, during wet and dry seasons. DOC concentration ranged from 4.82 mg · L − 1 to 47.83 mg · L − 1 , with a mean value of 15.04 mg · L − 1 , 2.84 times higher than the global average. The Zoige rivers had higher DOC concentration and highly terrigenous CDOM. Signiﬁcantly higher DOC concentration was observed for the Zoige rivers in the wet season compared to the dry season. In contrast, the alpine-gorge rivers had higher DOC levels in the dry season. No signiﬁcant correlations were observed between DOC and CDOM at all rivers due to the inﬂuence of autochthonous sources on the alpine-gorge rivers and intensive photochemical degradation of terrigenous DOM in the Zoige rivers. Signiﬁcant relationships between CDOM and speciﬁc ultraviolet absorbance at 254 nm (SUVA 254 ) and between CDOM / DOC and SUVA 254 were observed, indicating that the aromaticity of DOM in the rivers was mainly determined by CDOM. Moreover, the DOC / CDOM properties of the Hei River indicate critical human-induced water quality degradation. High DOC level and high browning degree were found in rivers in the Zoige alpine wetland, indicating that large amounts of terrigenous DOC were released to the aquatic systems of the region.


Introduction
Dissolved organic matter (DOM), operationally measured as dissolved organic carbon (DOC) and defined as the organic fraction in solution that passes through a 0.45 µm filter [1], is a ubiquitous component in soils and natural waters which establishes a connection between terrestrial and aquatic carbon pools [2,3]. DOC is related to the energy supply, light penetration, the biogeochemical cycles of carbon (C) and elements coupled to C, and the transportation of metals and nonpolar organic contaminants in waters [4][5][6], and thus plays an important role in environmental processes in aquatic ecosystems. Moreover, DOC acts as the precursor of the carcinogenic by-products (e.g., trihalomethanes) produced in the water disinfection process, exerting adverse effects on the quality of drinking water absorption) and identify the characteristics and sources of DOM; (iii) evaluate the relationships between DOC and CDOM absorption parameters, and further investigate DOM and CDOM quality, sources, and their roles in environmental processes.

Study Area
Two rivers in the Zoige wetland, the Bai River and the Hei River (major tributaries of the Yellow River), and two rivers in the alpine-gorge area of the northeastern QTP, the Zagunao River and the Min River, were examined ( Figure 1). With an average altitude over 3500 m, the Zoige wetland has a cold, humid continental monsoon climate. The mean annual temperature is around 1 °C, with the highest monthly mean temperature recorded in July (10.7 °C) and the lowest in January (−10.1 °C). The mean annual precipitation is around 650mm, 86% of which falls during April and October (data were obtained from 1971-2000 climate data of the Zoige meteorological station). In the Zoige region, the growing season is very short (May to September), and so is the non-frozen period (16-25 days) [25]. The main land cover types of the Zoige region are alpine meadow and swamp meadow [31]. The Min River is the largest tributary in the upper reaches of the Yangtse River. In this study, we focused on the upstream of the Min River (tracing to the headwater) and the Zagunao River (the largest tributary of the Min River). Compared to the climate conditions of the Zoige wetland, this watershed region also has distinct dry and wet seasons, but has slightly higher temperature and precipitation levels. In addition, forests and scrubs are the main land cover types in the alpine-gorge river basins. Land vegetation data were obtained from Resource and Environment Data Cloud Platform of Institute of Geographical Sciences and Natural Resources Research, Chinese Academy of Sciences.

Water Sample Collection
Water samples were collected from 35 sites across the four rivers in June and November (representing wet and dry seasons) of 2018 (Supplementary Materials Table S1). Teflon sampling bottles were washed once with soap water and then several times with ultra-pure water before being dried in an oven at 60 • C [32]. After pre-rinsing the Teflon sampling bottles three times with the surface water, we collected three parallel samples and preserved them in a dark refrigerator at 4 • C. After the collection, the samples were immediately transported to the laboratory, and were filtered through 0.45 µm filters (PVDF, Millipore, MA, USA) before analysis.

DOC and CDOM Absorption Determination
A Milti N/C 2100S Total Organic Carbon (TOC) analyzer (Analytik Jena, Germany) was used to measure the DOC and total dissolved nitrogen (TDN) concentrations (mg·L −1 ) [33]. Absorption spectra of the samples were measured between 240 and 550 nm at 3 nm intervals using an Aqualog ® spectrophotometer (Horiba, HORIBA Instruments Incorporated, NJ, USA) with a 1 cm quartz cuvette, with Milli-Q water used as reference [33]. The CDOM absorption coefficient a(λ 0 ) (m −1 ) was calculated from: where A(λ) is the measured absorbance at wavelength λ and l is the cuvette path length 1 cm [12].
The absorption coefficient at 355 nm was selected to represent the relative concentration of CDOM [33][34][35]. SUVA 254 (L·mg C −1 ·m −1 ) was defined as dividing the UV absorbance at λ = 254 nm by the DOC concentration (mg·L −1 ), and higher values indicate greater aromaticity [14]. The spectral slope parameter (S) was calculated from: where a(λ) is the CDOM absorption coefficient at wavelength λ, a(λ 0 ) is the absorption coefficient at a reference wavelength λ 0 [13]. The spectral slope ratio (S R ) was calculated as the ratio of S 275-295 to S 350-400 , with higher values indicative of lower aromaticity and lower molecular weight [12].

Fluorescence Determination
Excitation-emission matrices (EEMs) fluorescence spectra were also measured on the Aqualog ® spectrophotometer (Horiba, HORIBA Instruments Incorporated, NJ, USA) with Milli-Q water as reference [33]. The excitation wavelength (E x ) and emission wavelength (E m ) ranged from 240 to 550 nm and 214 to 619 nm, respectively. The intervals and slits were at 3.0 nm and 2.5 nm, respectively. The Raman scattering and Rayleigh scattering were automatically eliminated from analysis by the system [33]. The fluorescence index (FI) was calculated by dividing the fluorescence intensity at 450 nm by that at 500 nm with E x = 370, which can be used to distinguish between terrestrial DOM (FI < 1.4) and microbial DOM (FI > 1.9) [36]. The humification index (HIX) was defined as dividing the cumulative fluorescence intensity at 435-480 nm by that at 300-345 nm with E x = 254, with higher values indicative of higher humification degree of DOM [37]. Data of FI and HIX are presented in the Supplementary Materials Figure S1 online.
Statistical analysis was carried out using the Python 3.7 software package and the IBM SPSS Statistics 21 software (IBM, Armonk, NY, USA). The Shapiro-Wilk (SW) test was applied to determine if the response variable of each group of data would have a normal distribution. Based on the distribution pattern of each variable, ANOVA (for a normal distribution) or Mann-Whitney (for non-normal distributions) were used to test for significance differences among groups, and Pearson's correlation (normal) or Spearman's correlation (non-normal) was used for correlation analysis. Assuming that the terrestrial DOM input and the microbial activity related with autochthonous processes were both high in summer, data from June for all sampling sites were used to conduct the hierarchical clustering to determine the DOM types in the four rivers. Hierarchical clustering analysis was performed through Ward's method with squared Euclidian distance as the measurement of similarity, and the variables involved in the process were DOC concentration, a (355), SUVA 254 , S R , TDN, FI, and HIX.

DOC Concentration
DOC concentration in rivers is affected by climate, land cover, soil characteristics and structure, runoff, and anthropogenic activities [38][39][40][41]. In our study area, DOC concentration ranged from 4.82 mg·L −1 in the Zagunao River to 47.83 mg·L −1 in the Bai River, with a mean value of 15.04 mg·L −1 . (Figure 2). The mean DOC concentration in the rivers of our study area was not only higher than that in other rivers of the Qinghai-Tibetan Plateau (1.16-6.57 mg·L −1 ) [42][43][44] and the global average (5.3 mg·L −1 ) [45], but also higher than that in African (8.0 mg·L −1 ) [45] and Arctic rivers (12 mg·L −1 ) [46]. Although differences may exist in the research time scale between our study and others, it can be concluded that the DOC levels in these examined rivers were relatively high in our study periods.
Water 2020, 12, x FOR PEER REVIEW 5 of 16 similarity, and the variables involved in the process were DOC concentration, a (355), SUVA254, SR, TDN, FI, and HIX.

DOC Concentration
DOC concentration in rivers is affected by climate, land cover, soil characteristics and structure, runoff, and anthropogenic activities [38][39][40][41]. In our study area, DOC concentration ranged from 4.82 mg·L −1 in the Zagunao River to 47.83 mg·L −1 in the Bai River, with a mean value of 15.04 mg·L −1 . (Figure 2). The mean DOC concentration in the rivers of our study area was not only higher than that in other rivers of the Qinghai-Tibetan Plateau (1.16-6.57 mg·L −1 ) [42][43][44] and the global average (5.3 mg·L −1 ) [45], but also higher than that in African (8.0 mg·L −1 ) [45] and Arctic rivers (12 mg·L −1 ) [46]. Although differences may exist in the research time scale between our study and others, it can be concluded that the DOC levels in these examined rivers were relatively high in our study periods.

Spatial Variation
The Hei River had the highest mean DOC concentration (19.90 ± 8.85 mg·L −1 ), followed by the Min (16.58 ±5 .32 mg·L −1 ), Bai (15.20 ± 9.38 mg·L −1 ), and the Zagunao rivers (8.20 ± 2.77 mg·L −1 ) ( Figure  2). The varying DOC concentrations could be partly attributed to the varying catchment land cover types, which is further related to the stock of organic carbon in the adjacent terrestrial ecosystems. It has been shown that in the Qinghai-Tibetan Plateau, rivers surrounded by wet meadows and alpine meadows had a higher DOC concentration than those surrounded by alpine steppes and alpine deserts [42]. Similar DOC patterns have been reported for Irish streams as follows: peatland and forest dominated watershed > grassland dominated watershed > arable land dominated watershed

Spatial Variation
The Hei River had the highest mean DOC concentration (19.90 ± 8.85 mg·L −1 ), followed by the Min (16.58 ± 5.32 mg·L −1 ), Bai (15.20 ± 9.38 mg·L −1 ), and the Zagunao rivers (8.20 ± 2.77 mg·L −1 ) ( Figure 2). The varying DOC concentrations could be partly attributed to the varying catchment land cover types, which is further related to the stock of organic carbon in the adjacent terrestrial ecosystems. It has been shown that in the Qinghai-Tibetan Plateau, rivers surrounded by wet meadows and alpine meadows had a higher DOC concentration than those surrounded by alpine steppes and alpine deserts [42]. Similar DOC patterns have been reported for Irish streams as follows: peatland and forest dominated watershed > grassland dominated watershed > arable land dominated watershed [40].
Our results further support these findings. The Hei and Bai rivers, flowing across the Zoige Peatland with catchments dominated by wet meadows and alpine meadows, both had higher DOC levels than the Zagunao Rivers, which is mainly surrounded by forests and scrubs. Besides, exhibiting relatively high values in November, the Min River, which is mainly surrounded by forests and scrubs, had similar mean DOC concentrations to the Bai River. This can be attributed to the high Lignin and cellulose decomposition rates of forest litter during the freeze-thaw period (October to April) around the catchments of the Min River [47], which results in abundant soil DOC being produced and exported to surface streams in November.
As for the spatial fluctuation across sites, since the sites of a specific river have similar climatic and hydrological conditions in a specific season, their variation of DOC concentration is mostly influenced by land cover types and pollution sources (e.g., domestic sewage and other organic pollutants) [48].

Seasonal Variation
All examined rivers showed significant seasonal variation in DOC concentration (p < 0.05) ( Figure 2). Higher DOC concentrations were found in June and lower values were found in November in the Zoige rivers (Figure 2c,d). This is consistent with the results reported by research on rivers in a permafrost wetland catchment in northeastern China [49], which found that maximum DOC concentration in rivers coincided with the maximum flood peaks (generally during spring and summer). A study of DOM in several large Arctic Rivers also reported higher DOC concentration in June than compared to November [50]. According to the study of Guo et al. [49], in wet and warm seasons, more DOC in the upper organic soil layer is released to surface streams with the increase of rainfall and fast runoff, and increased thawing depth enables deeper soil layers to release DOC into streams, while in the catchments of alpine-gorge rivers where the carbon stock of the organic soil layer is lower and limited, the terrestrial input of DOC is no longer the crucial factor in determining DOC levels in waters. Thus, the dilution effect of the increased runoff in the alpine-gorge rivers, which contributes to lower DOC concentrations in June, cannot be neglected anymore. In the meantime, the presence of lower DOC concentrations in June than November in the alpine-gorge rivers may also result from the greater influences of photochemical and microbial degradation, considering that the increased light level and longer water residence time in June could both facilitate the degradation of DOM [35].

CDOM and CDOM/DOC
Quantified as a(355), CDOM concentrations ranged from 0 to 15.895 m −1 in the examined rivers (Figure 3a-d). The alpine-gorge rivers had an average a(355) level of 1.26 ± 1.08 m −1 , which was close to aquatic systems dominated by autochthonous sources (usually < 1.0 m −1 ) [51], while the Zoige rivers had an average a(355) level of 8.37 ± 4.71 m −1 , belonging to the terrigenous-dominated type (usually > 5.0 m −1 ) [35]. The Bai River stands out as having the highest average a(355)/DOC, SUVA 254 , and lowest average S R value of all examined rivers. The Hei river is also noteworthy as having the highest average a(355). The observation of marked terrigenous feature in the Zoige rivers could be attributed to the abundant terrigenous CDOM (DOM) input from the Zoige peatland, where large amounts of peat and humus are stored. Meanwhile, the flat terrain characteristics along the Zoige rivers allow a higher DOC leaching degree of the soil and more DOC release from the bank or sediment into the water system [33]. Furthermore, owing to the high elevation (over 3 km) and low temperature in the Zoige region, the decrease of microbial activity contributes to the accumulation of high-molecular-weight and aromatic DOM in the rivers [17,33]. Seasonal variation of CDOM was significant (p < 0.01) in the Zoige rivers, with higher values observed in June and lower values observed in November (Figure 3c,d). This finding supports the former inference that Zoige rivers accept more terrestrial derived DOC in June. Unlike the Zoige Rivers, the alpine-gorge rivers showed no significant difference in the seasonal variations of a(355) (Figure 3a,b), indicating that the increased precipitation and runoff in June may dilute the terrestrial input of DOC (CDOM).
CDOM/DOC in the rivers was calculated as a proxy for the relative contribution of CDOM to DOC, and covered the range from 0 to1.737 L·mg −1 ·m −1 among all sampling sites (Figure 3e-h). The Bai River exhibited the highest mean CDOM/DOC value of 0.61 ± 0.43 L·mg −1 m −1 , followed by the Hei River (0.43 ± 0.18 L·mg −1 ·m −1 ), the Zagunao River (0.19 ± 0.28 L·mg −1 ·m −1 ), and the Min River (0.08 ± 0.08 L·mg −1 ·m −1 ), which was indicative of the great contribution of non-chromophoric organic matter to DOC in the alpine-gorge rivers and the important contribution of CDOM to DOC in the Zoige Rivers. For the alpine-gorge rivers, the spatial and temporal variabilities of CDOM/DOC were similar to that of CDOM absorption, indicating that the fluctuation of non-chromophoric DOM concentrations along sites may be quite stable.
Water 2020, 12, x FOR PEER REVIEW 7 of 16 no significant difference in the seasonal variations of a(355) (Figure 3a,b), indicating that the increased precipitation and runoff in June may dilute the terrestrial input of DOC (CDOM). CDOM/DOC in the rivers was calculated as a proxy for the relative contribution of CDOM to DOC, and covered the range from 0 to1.737 L·mg −1 ·m −1 among all sampling sites (Figure 3e-h). The Bai River exhibited the highest mean CDOM/DOC value of 0.61 ± 0.43 L·mg −1 m −1 , followed by the Hei River (0.43 ± 0.18 L·mg −1 ·m −1 ), the Zagunao River (0.19 ± 0.28 L·mg −1 ·m −1 ), and the Min River (0.08 ± 0.08 L·mg −1 ·m −1 ), which was indicative of the great contribution of non-chromophoric organic matter to DOC in the alpine-gorge rivers and the important contribution of CDOM to DOC in the Zoige Rivers. For the alpine-gorge rivers, the spatial and temporal variabilities of CDOM/DOC were similar to that of CDOM absorption, indicating that the fluctuation of non-chromophoric DOM concentrations along sites may be quite stable.

SUVA 254 and S R
It has been shown that the SUVA 254 are positively correlated with the percentage aromaticity of DOM as an effective tool to determine DOC characteristics [14]. The SUVA 254 values of the examined sites ranged from 0.17 L·mg C −1 ·m −1 in the Min River to 7.19 L·mg C −1 ·m −1 in the Bai River (Figure 4a-d). The Bai River had the highest average SUVA 254 level (2.53 ± 1.81 L·mg C −1 ·m −1 ), followed by the Hei (1.77 ± 0.78 L·mg C −1 ·m −1 ), Zagunao (0.81 ± 0.53 L·mg C −1 ·m −1 ), and the Min rivers (0.45 ± 0.31 L·mg C −1 ·m −1 ). The Zagunao, Min, and Hei rivers all showed that SUVA 254 values in June were higher than that in November (p < 0.05), indicating a larger input of terrestrial DOM during the wet season.
Water 2020, 12, x FOR PEER REVIEW 8 of 16

SUVA254 and SR
It has been shown that the SUVA254 are positively correlated with the percentage aromaticity of DOM as an effective tool to determine DOC characteristics [14]. The SUVA254 values of the examined sites ranged from 0.17 L·mg C −1 ·m −1 in the Min River to 7.19 L·mg C −1 ·m −1 in the Bai River (Figure 4ad). The Bai River had the highest average SUVA254 level (2.53 ± 1.81 L·mg C −1 ·m −1 ), followed by the Hei (1.77 ± 0.78 L·mg C −1 ·m −1 ), Zagunao (0.81 ± 0.53 L·mg C −1 ·m −1 ), and the Min rivers (0.45 ± 0.31 L·mg C −1 ·m −1 ). The Zagunao, Min, and Hei rivers all showed that SUVA254 values in June were higher than that in November (p < 0.05), indicating a larger input of terrestrial DOM during the wet season.  The spectral ratio S R (calculated as the ratio of S 275-295 to S 350-400 ) is suggested to be related to the source and quality of DOM, with higher values indicative of a decrease in molecular weight and aromaticity [12,13]. S R values tested in our sites ranged from 0 to 11.53, with the majority of the values ranging from 2 to 5 (Figure 4e-h). The order of the average S R level in the rivers was the Bai River (3.25 ± 0.55) ≤ the Hei River (3.36 ± 0.58) < the Zagunao River (4.26 ± 2.82) < the Min River (4.28 ± 1.92). No significant difference was observed in the seasonal variation of S R in the Zagunao river, the Bai River, and the Hei river (p > 0.05), although the alpine-gorge rivers had relatively higher S R values in November, which suggests that S R and SUVA 254 have different sensitivities towards DOM aromaticity.
Meanwhile, it is of interest that although the Bai River and the Hei River had similar average S R and average a(355) levels, slightly higher SUVA 254 levels were detected in the Bai River. Combined with the presence of a higher average a(355)/DOC level in the Bai River than that in the Hei River (Figure 3g,h), we can infer that the Hei River has a larger proportion of non-chromophoric DOM which is of low-aromaticity. In addition, since S 275-295 and S 350-400 are corelated with different molecular formulae of DOM [52], their ratio (i.e., S R ) could be affected by the molecular formulae indicated by any of these spectral slopes. Wagner et al. [52] also found that S 275-295 provided similar information to other indices, including SUVA 254 , while the molecular formulae negatively corelated with S 350-400 were more oxidized and more aromatic than molecular formulae negatively corelated with S 275-295 . The fact that S R provides different information from SUVA 254 on some occasions may be induced by the complex aromatic structures of DOM.

Spatial Fluctuations Along Rivers
Concluding from the research conducted by Liu et al. [33] and our results, the wide fluctuation ranges of CDOM and CDOM properties along the Bai River can be attributed to the unequal input intensity of terrigenous sources, while the autochthonous impacts on the Zagunao River and the Min River are more significant. However, exhibiting extremely high DOC concentration in June but normal a(355) and low SUVA 254 values (Figures 2c, 3c and 4c), B4 is also influenced by autochthonous or anthropogenic sources. Furthermore, the fluctuations of a(355) and S R among the sampling sites of the Hei River were relatively stable, which may be due to the similar topographic, climatic, and hydrological conditions of these sites.

Correlations Between DOC and CDOM
Positive and significant relationships between the CDOM absorption coefficients and DOC concentrations have been observed in many inland water systems like rivers, lakes, and reservoirs in recent years, indicating the possibility of estimating DOC from in-situ or remote monitoring of CDOM absorption [4,13,17,40,53]. However, this does not apply to different river systems under various spatial and temporal backgrounds, and decoupled CDOM-DOC relationships were more common in our study ( Figure 5), noting that significant correlation was only observed in November for the Zagunao river (r = −0.929, p = 0.003). It was reported that river systems with CDOM parameters indicative of autochthonous or anthropogenic, or photochemically degraded allochthonous sources had weaker CDOM-DOC relationships than other typical rivers systems in America [13]. In addition, autochthonous (i.e., microbial/algae) produced DOM is easily removed from DOC pools through microbial activities in the water systems [18,54], and CDOM or high-aromaticity DOM are suggested to be more sensitive to photo-irradiation than bulk DOC in waters [55][56][57], dissociating CDOM from DOC to some extent. It was also reported that highly polluted tributaries derived from anthropogenic sources could weaken the CDOM-DOC relationships in the main stream of the Songhua River in Northeast China [48]. Thus, considering that the fluctuation of non-chromophoric DOM is quite stable along the alpine-gorge rivers (will be discussed later), the disassociated CDOM-DOC relationships may be attributed to the impact of autochthonous sources on DOC pools in these rivers. However, as for the Zoige rivers with distinct terrigenous feature, the presence of disassociated CDOM-DOC relationships may have resulted from the intensive photochemical degradation of terrestrial DOM (CDOM), which is supported by research on bog water, arctic lakes and rivers, and even estuarine rivers, suggesting that the role of photodegradation on terrestrial DOM is higher than previously thought and an increased fraction of low-aromaticity and low-molecular-weight organic matter is produced during the process [56][57][58].
Water 2020, 12, x FOR PEER REVIEW 10 of 16 terrestrial DOM (CDOM), which is supported by research on bog water, arctic lakes and rivers, and even estuarine rivers, suggesting that the role of photodegradation on terrestrial DOM is higher than previously thought and an increased fraction of low-aromaticity and low-molecular-weight organic matter is produced during the process [56][57][58].

Correlations Among DOC and Light Absorption Characteristics
Although positive DOC-SUVA254 correlation was evidenced in some cases, there was overall negative correlation between DOC and SUVA254, which was significant in November in the Zagunao River (p < 0.01), in June in the Bai River and the Hei River (p < 0.05, p < 0.01, respectively) ( Table 1). The significantly negative correlation between SUVA254 and DOC indicates that the increase of DOC in the river is mainly induced by low-aromaticity components, for example, microbial or algae derived organic matter (in the alpine-gorge rivers) and some photochemically produced organic matter (in the Zoige rivers).
Significantly positive correlation between CDOM and SUVA254 existed in many cases (especially in the Bai river), suggesting that the increase of aromaticity in DOC was closely related to the input of CDOM. However, no CDOM-SUVA254 correlation was found in the Hei River, suggesting that river CDOM concentration was not the only critical factor of DOC aromaticity. The SUVA254-CDOM/DOC relationship tended to be significant and strong on many occasions and was somehow dependent on the SUVA254-CDOM relationship. However, the independence of the SUVA254-CDOM/DOC relationship from the SUVA254-CDOM relationship observed in June in the Hei River (r = 0.993, p < 0.01 and r = −0.071, p > 0.05, respectively) indicates that the aromaticity percentage of DOC is not only determined by CDOM but can also be diluted by non-chromophoric DOM. Regardless of the few occasions exhibiting not significant CDOM-CDOM/DOC relationships (e.g., the Hei River), significantly positive CDOM-CDOM/DOC relationships were common across the four rivers, and were independent of spatiotemporal differences. The simultaneous presence of these two results (the dependence of SUVA254-CDOM/DOC on SUVA254-CDOM relationship and significantly positive CDOM-CDOM/DOC relationship) demonstrates that the fluctuation of non-chromophoric DOM is quite stable along the rivers. Combining the discussed relationships of the parameters in the Hei River, we inferred that the Hei River may be profoundly affected by anthropogenic sources, as the fluctuation of non-chromophoric DOM in the Hei River is not as stable as that in the alpinegorge rivers and the Bai River.

Correlations Among DOC and Light Absorption Characteristics
Although positive DOC-SUVA 254 correlation was evidenced in some cases, there was overall negative correlation between DOC and SUVA 254 , which was significant in November in the Zagunao River (p < 0.01), in June in the Bai River and the Hei River (p < 0.05, p < 0.01, respectively) ( Table 1). The significantly negative correlation between SUVA 254 and DOC indicates that the increase of DOC in the river is mainly induced by low-aromaticity components, for example, microbial or algae derived organic matter (in the alpine-gorge rivers) and some photochemically produced organic matter (in the Zoige rivers).
Significantly positive correlation between CDOM and SUVA 254 existed in many cases (especially in the Bai river), suggesting that the increase of aromaticity in DOC was closely related to the input of CDOM. However, no CDOM-SUVA 254 correlation was found in the Hei River, suggesting that river CDOM concentration was not the only critical factor of DOC aromaticity. The SUVA 254 -CDOM/DOC relationship tended to be significant and strong on many occasions and was somehow dependent on the SUVA 254 -CDOM relationship. However, the independence of the SUVA 254 -CDOM/DOC relationship from the SUVA 254 -CDOM relationship observed in June in the Hei River (r = 0.993, p < 0.01 and r = −0.071, p > 0.05, respectively) indicates that the aromaticity percentage of DOC is not only determined by CDOM but can also be diluted by non-chromophoric DOM. Regardless of the few occasions exhibiting not significant CDOM-CDOM/DOC relationships (e.g., the Hei River), significantly positive CDOM-CDOM/DOC relationships were common across the four rivers, and were independent of spatiotemporal differences. The simultaneous presence of these two results (the dependence of SUVA 254 -CDOM/DOC on SUVA 254 -CDOM relationship and significantly positive CDOM-CDOM/DOC relationship) demonstrates that the fluctuation of non-chromophoric DOM is quite stable along the rivers. Combining the discussed relationships of the parameters in the Hei River, we inferred that the Hei River may be profoundly affected by anthropogenic sources, as the fluctuation of non-chromophoric DOM in the Hei River is not as stable as that in the alpinegorge rivers and the Bai River.

DOC Classification
All of the samples were divided into two main classes (A, B), and each class was subdivided into two subclasses (expressed as Type 1, Type 2, Type 3, and Type 4) ( Table 2). Samples collected from the Zoige rivers (except for B1) belonged to class A, exhibiting significant terrigenous feature in a(355) and HIX, while samples collected from the alpine-gorge rivers were classified as class B, with a(355) and HIX indicative of autochthonous feature. In class A, Type 1 waters (including most samples from the Hei River) exhibited higher DOC concentrations, lower SUVA 254 values, and higher TDN concentrations in comparison with Type 2 waters, indicating that Type 1 waters may be intensively influenced by anthropogenic sources, as increased TDN may result from wastewater discharge and fertilizer use in the river basins [43]. In class B, Type 3 waters had lower a(355), higher S R , and lower HIX values than Type 4 waters, which demonstrates the autochthonous feature in Type 3 waters was more significant than the average level in the alpine-gorge rivers.

Summary and Conclusions
DOC concentration, CDOM absorption coefficient a(355), SUVA 254 , the spectral slope ratio S R , and a(355)/DOC were examined in four rivers (the Zagunao, Min, Bai, and Hei rivers) of the northeastern Qinghai-Tibetan Plateau in June and November of 2018. Distinct spatiotemporal patterns were evidenced in the distribution of these parameters. The following results were obtained: (1) the average DOC level (15.04 mg·L −1 ) of the rivers was relatively high in the study periods, and the browning level was particularly high in June in the Zoige rivers, indicating that a large amount of terrestrial DOC was transferred to the surface waters in the Zoige region; (2) the Zoige rivers stand out as having higher average DOC concentration, a(355), SUVA 254 , a(355)/CDOM, and lower S R value than the alpine-gorge rivers, which is indicative of significant terrigenous characteristics. The temporal distribution patterns of DOC and CDOM were also different between the Zoige rivers and the alpine-gorge rivers; (3) differing from other studies reporting a significant linear relationship between the CDOM absorption coefficient and DOC concentration, disassociated CDOM-DOC relationships were common in our study. Therefore, deriving the DOC flux from CDOM parameters may not be feasible in these rivers; (4) combining the results of DOC classification and correlations of DOC and CDOM properties, we conclude that the Hei River may have experienced a profound influence by anthropogenic sources.
Significant warming and wetting trends were observed on the Qinghai-Tibetan Plateau over the past few decades [59], and if these trends continue, the terrestrial input of DOM into the plateau's freshwaters will be greatly increased, which may ultimately affect the carbon cycle in the aquatic systems at both the regional and global scales. Since the fate of DOM and CO 2 emissions in waters is related to photochemical and microbial degradation processes, and many studies have pointed out that the increase of terrigenous DOM will enhance the CO 2 emissions in various water systems [17,30,57,60], long-term monitoring of the DOC concentration and its optical characteristics with higher time resolution is needed to track the carbon dynamics in the aquatic systems in the Zoige alpine wetlands. Our findings also indicate the intensive anthropogenic influence on aquatic DOM in this region (especially in the Hei River). In order to further reveal the impacts of global change on the DOC fate and carbon cycle in these aquatic ecosystems, more attention should be given to the sources, dynamics, characteristics, and photochemical/microbial degradation mechanisms of aquatic DOM.