jh：青藏高原腹地纳木错流域化学风化研究取得新进展-凯发k8官网手机客户端

　　岩石化学风化从大气中吸收co2并将其封存于风化产物中，构成碳循环的重要环节，也是“碳中和”的重要途径。受“新生代青藏高原隆升导致全球降温”假说的影响，青藏高原是化学风化研究的热点地区。当前，关于喜马拉雅山脉南坡和青藏高原东部地区的化学风化研究较为充分，但青藏高原内陆地区的化学风化状况仍不清楚，特别是冰川的存在对流域化学风化具有何种影响还存在争论。这阻碍了人们对于青藏高原对全球碳循环影响的理解。

　　我所环境变化与多圈层过程团队邬光剑研究员等选取高原腹地的纳木错流域为研究对象，系统采集了纳木错流域的大气降水、雪冰、河水和河流泥沙等样品，综合运用正演模型和端元混合模型量化了该区域的化学风化过程，探讨了纳木错流域化学风化的季节变化特征、冰川作用对化学风化的影响以及纳木错各子流域化学风化的时空差异及控制因素。

　　图1 纳木错各子流域的位置、岩性与采样点分布

　　选取你亚曲（编号22)代表纳木错东、西两岸地势平缓、流域面积较大的河流，选取曲嘎切（编号9)代表纳木错南岸地形陡峻、流域面积较小的河流（图1）。研究结果显示，你亚曲的河水溶质来自大气降水、硅酸盐岩风化、碳酸盐岩溶解和硫化物氧化的年均值分别为5%、12%、65%和18%，季节差异小；上述来源对曲嘎切河水溶质贡献的年均值分别为10%、35%、34%和20%，季节差异显著：硅酸盐岩风化的贡献在季风期低于非季风期，而碳酸盐岩溶解呈相反变化趋势（图2）。你亚曲与曲嘎切的碳酸盐岩风化速率（cwr）与硅酸盐岩风化速率（swr）均显示了十分显著的季节变化，季风期的风化速率远高于非季风期（4~8倍）。径流是这两个典型子流域化学风化速率季节变化的控制因素。

　　图2 纳木错你亚曲(a)和曲嘎切(b)子流域化学风化过程的季节变化

　　选择地理条件与曲嘎切相似但无冰川分布的曲嘎琼流域（编号8)为研究点，与具备山地冰川发育的曲嘎切流域进行对照研究，探讨冰川作用对化学风化的影响。结果显示，无冰川分布的曲嘎琼流域的cwr和swr分别为6.9 t/km2/y和3.8 t/km2/y，而有冰川分布的曲嘎切流域cwr和swr分别为9.7 t/km2/y和6.7 t/km2/y，曲嘎切的化学风化速率是曲嘎琼的1.4~1.8倍，说明该区域的冰川作用促进了化学风化，冰川消融和磨蚀带来的更大径流量和更强物理侵蚀速率是导致冰川区化学风化强度更大的主要原因。

　　纳木错地区季风期的化学风化速率远高于非季风期，非季风期因化学风化产生的溶质通量对全年的贡献可忽略不计。通过季风期对纳木错24个子流域进行定期采样，较为精准地估算了整个纳木错流域的化学风化速率。研究表明纳木错各子流域的cwr存在很强的空间异质性，最高值与最低值分别为55 t/km2/y和9.2 t/km2/y。cwr与岩性系数呈显著正相关，与平均海拔呈显著负相关，说明岩性对纳木错流域cwr的影响较大。swr的空间分异不如cwr显著，最高值与最低值分别为7.9 t/km2/y和2.3 t/km2/y，径流是影响swr最为重要的因素。纳木错南岸8个冰川子流域的化学风化速率是岩性相似非冰川流域的1.3倍，再次证实了该地区的冰川作用促进了化学风化。地形是导致子流域化学风化强度呈现空间差异的更深层次原因（图3)。从流域尺度上看，纳木错流域cwr与swr平均值分别为24.3 t/km2/y和4.5 t/km2/y，相比于青藏高原的边缘地区，高原腹地的化学风速率处于较低水平，尤其是硅酸盐岩风化。

　　图3 纳木错各个小流域化学风化速率与环境因子的关系

　　图中数字代表相关系数，实线代表显著相关，蓝色代表负相关，红色代表正相关

　　上述研究成果近期陆续发表于国际地学期刊《journal of hydrology》，题目分别为“seasonal variation of chemical weathering and its controlling factors in two alpine catchments, nam co basin, central tibetan plateau”、“glaciation enhanced chemical weathering in a cold glacial catchment, western nyaingêntanglha mountains, central tibetan plateau”和“small-catchment perspective on chemical weathering and its controlling factors in the nam co basin, central tibetan plateau”，第一作者为我所博士后于正良，通讯作者为邬光剑研究员。该研究获得国家自然科学基金项目(41725001)、中国科学院a类战略性先导科技专项(xda20060201)和第二次青藏高原综合科学考察研究(2019qzkk0201)等的联合资助。