全球冰川径流中硅的研究进展

刘 莎， 李向应， 杨船洋， 韩添丁， 井哲帆， 朱永华

（1.河海大学水文水资源学院，江苏南京210098； 2.陕西省地表系统与环境承载力重点实验室，陕西 西安 710127； 3.西北大学城市与环境学院，陕西 西安 710127； 4.中国科学院西北生态环境资源研究院冰冻圈科学国家重点实验室，甘肃 兰州 730000）

0 引言

硅（Si）是自然界极为常见的元素，在地壳中的丰度仅次于氧［1］。Si 极少溶于水，大多与氧结合形成SiO2或其他硅酸盐矿物［2］。硅酸盐矿物是最常见的造岩矿物，在化学风化过程中会释放出Si。Si 元素是硅质生物体必需的营养元素［3-6］，在海洋生态系统平衡中扮演重要角色。Si 的迁移转化是维持生态系统平衡的关键因素之一，其生物地球化学循环过程与碳循环和气候变化关系密切。Si 循环的主要过程为硅酸盐矿物的化学风化和硅酸盐沉积物的固结成岩，中间过程包括Si 的生物地球化学传输和转化以及与碳等其他元素的相互作用［7］。大坝和水库修建以及耕地灌溉等人类活动通过增加水体的滞留时间也会影响Si循环［8］。

IPCC 指出，大气温室气体浓度已达到过去80万年以来的最高水平，全球地表及浅层海水的温度迅速升高，海平面也加速上升［9］。气候变暖对冰冻圈的影响最为显著［10］。近20年来，冰盖的冰储量持续减少，冰川普遍退缩，冻土层厚度逐渐减薄［9］。在冰川加速消融背景下，夹带着泥沙的融水径流向下游输送了大量的生物活性元素（如Si、C、Fe、N 和P），当这些元素进入水库、湖泊和海洋后［7］会影响生态系统和生物地球化学循环［11-21］。活性Si元素可被硅藻等硅质海洋生物吸收利用并形成细胞结构，细胞死亡后将Si 沉降到海底参与Si 的生物地球化学循环过程（图1）；硅藻的呼吸作用会固定大气CO2，因此活性Si 的增加会促进硅藻繁殖、影响海洋初级生产力和全球碳循环，进而对气候变化产生反馈作用［22-27］。近期研究表明，冰川很可能是海洋中生物活性Si 的重要来源之一［19］。与其他来源相比，冰川Si 的活性高，很可能是海洋初级生产力的重要限制因子［28-29］。随着气候持续变暖，从冰川中释放的Si大大增加，进而可能显著影响下游生态系统和碳循环，最终影响气候变化。然而，已有研究主要关注可被生物直接利用的可溶性Si（DSi，Dissolved Si）、忽视了悬浮泥沙附着的颗粒态Si（ASi，Amorphous Si）的输出潜力和贡献，因此大大低估了冰川Si 的释放率及其对Si循环的影响，导致针对Si循环收支平衡的重建研究具有较大误差［5，30-31］。此外，Si稳定同位素比率（δ30Si）是一种重要工具，可用于追踪Si 的来源、迁移转化和沉积等循环过程，可为现代及历史时期的生物地球化学循环过程研究提供可靠资料［7］。

图1 冰川消融与生态系统和碳循环的关系Fig.1 Relationship between glacier melt and aquatic ecosystems and carbon cycle

涉及冰川Si 的研究始于20 世纪90 年代。初期关注流域硅酸盐的风化机制及其对溶质释放的影响，同时指出硅酸盐风化可能影响万年尺度上的碳循环［22，32-40］。本世纪以来开始关注冰川消融对Si浓度/通量及δ30Si 的影响。随着研究的深入，逐渐认识到冰川在Si 循环中的角色可能比以前认为的更加重要。随着δ30Si 测定技术的发展，Si 循环研究也相继展开，目前已获得一些岩芯、基岩及径流的δ30Si 资料，同时指出冰川可能影响δ30Si 的时空分布并提出了可能的影响机理［7，13，20，41-44］。然而，专门针对冰川Si 的监测研究仍较少，主要集中在格陵兰冰盖和一些冰川流域［19，40，44-46］，已有资料的时间和空间连续性非常不足。本文综述了冰川径流中各种形态Si 浓度/通量及δ30Si 的数量级，对比分析了各形态Si 浓度/通量及δ30Si 的空间变化特征，并指出了存在的问题与亟待加强的研究方向。

1 研究方法回顾

1.1 实验分析方法

冰川径流中的Si（简称冰川Si）由DSi 和ASi 组成。DSi是过滤后滤液中包含的Si，ASi是滤膜附着的泥沙中包含的Si。研究方法一般包括样品采集、实验分析和数据处理。

（1）样品采集。采样前先将样品瓶和过滤装置清洗干净；采样点尽可能靠近冰川，将采样瓶置于水面以下并使瓶口朝着来水方向采样，采样前先用母液冲洗采样瓶，采样时间尽可能覆盖消融期；采样后用0.45 μm 孔径滤膜过滤样品，起初的滤液用于冲洗过滤装置和样品瓶，随后用作DSi样品、滤膜附着的泥沙为ASi 样品［19］；在野外将样品避光冷藏/冷冻保存［44，47］。

（2）实验分析。通常应用流动注射分析仪、原子吸收分光光度计、电感耦合等离子体质谱仪和直读光谱仪等测定Si浓度［19，25，34-35，39，43，46，48-51］。一般用标准钼酸比色法［44，47］和分光光度计法［23，52］测定DSi 浓度。在测定ASi之前需要对样品预处理和碱性提取。一般步骤为：先将样品风干，后称取定量样品并装入HDPE瓶中，随后加入Na2CO3溶液，接着在水浴中加热，然后在不同时段移取定量溶液并分别置于微量离心管中，最后冷藏保存；实验分析时，先移取定量冷藏溶液至离心管中，再加入盐酸溶液稀释并中和多余的碱溶液［19，24，44］。Si 的稳定同位素有28Si、29Si和30Si［53］，其稳定同位素比率常用δ30Si 表示，计算公式为δ30Si样品=［（30Si/28Si）样品/（30Si/28Si）标准-1］×1000，这里的标准为国际标准物质NBS-28［44，54-55］。常用高分辨率多接收电感耦合等离子体质谱仪（如Nu Plasma HR-MC-ICP-MS）的干燥等离子模式测定DSi（δD30Si）和ASi（δA30Si）的Si 同位素［13，56-59］。在测定δA30Si之前要对样品进行预处理：先称取定量样品并置于层流罩中风干，然后在风干样品中加入NaOH溶液并在水浴中持续加热，随后向加热的溶液中加入硝酸，最后将混合溶液稀释-过滤-离心分离［44］。

（3）数据处理。冰川径流中DSi 和ASi 平均浓度与年通量、δ30Si 平均值的计算方法一般为算数平均或流量加权平均，冰川流域DSi 和ASi 的年产量则根据年通量与流域面积获取。

1.2 资料获取方法

文献查阅主要有两种查阅方式：（1）找到目前最新发表、相关性最高的文章，再从这些文章的参考文献中查找可能相关的文献。（2）主要从Research Gate、中国知网、Science Direct、Springer 等数据库搜索关键词“冰川”、“冰川径流”、“冰川融水”、“Si”、“Si 同位素”。通过两种方式共查阅了约200篇文献，其中有冰川径流中DSi 浓度资料的文献有33篇，然而有ASi和Si同位素资料的文献仅各3篇，且均为近五年新发表的文献。

从查阅的文献中收集并整理了全球冰川径流中DSi 和ASi 的浓度、通量与产量数据，方法如下：（1）文献中已给出的浓度、通量与产量数据不进行其他处理（单位换算除外）。（2）对于单个冰川，径流中浓度资料的样品获取地点应尽可能靠近冰川末端，舍去距离冰川末端太远的资料；若通量和产量知其一且已知流域面积，则根据“通量=产量×面积”这一关系得到未知的通量或产量数据；若已知浓度与径流量，则根据“通量=浓度×径流量”得到通量数据。（3）区域/全球平均浓度/产量及其标准偏差由区域/全球范围内已有资料冰川径流中Si 浓度/产量计算均值和标准偏差得到。

2 研究动态

2.1 硅浓度和通量研究

目前，涉及DSi 的研究区主要集中在格陵兰冰盖、冰岛、斯瓦尔巴、斯堪的纳维亚、阿拉斯加、欧洲中部、加拿大北极、加拿大西部/美国、南极/亚南极和亚洲流域，涉及ASi 的研究集中在格陵兰冰盖流域且资料较少（图2）。

图2 全球冰川Si监测点的空间分布以及径流中DSi（a）和ASi（b）浓度的空间变化（DSi资料来源见附表1，ASi资料引自Hawkings等［19，23，66］）Fig.2 Spatial distribution of monitoring points about Si，and mean concentrations of dissolved Si（a）and amorphous Si（b）from glacial runoff globally（datas of dissolved Si are shown in Supplement table 1，amorphous Si are from Hawkings et al［19，23，66］）

冰山的DSi 和ASi 资料都仅分布在格陵兰冰盖周围。数据显示冰山中DSi 和ASi 的平均浓度分别为8.0 μmol·L-1（＜0.02～20 μmol·L-1）和49.2 μmol·L-1（47.9～50.5 μmol·L-1）（附表1～2），可见ASi 浓度显著大于DSi 浓度，可能是冰山中DSi 的溶解度小造成的。冰山中的DSi 浓度很低，如2014 年Godthåbs 峡湾冰山的浓度低于0.02 μmol·L-1［60］，其他一些冰山的DSi 浓度甚至可能低于仪器检测限［19，61］。此外，富含泥沙的冰山中DSi 浓度会显著大于洁净冰山［19］。由格陵兰冰盖周围冰山的DSi和ASi 平均浓度以及年径流量（612 km-3·a-1）［62］，得到冰山的DSi 和ASi 通量分别为0.05×105mmol·a-1和0.03×106mmol·a-1［19］。根据冰山周周叶绿素浓度的分布和南大洋冰山底部硅藻群落的分布，已有研究发现冰山也是海洋中ASi的重要来源之一［19，63］。

对DSi 的研究广泛分布在格陵兰冰盖，对ASi的研究较少且集中分布在格陵兰冰盖，南极冰盖的研究极少（图2）。格陵兰冰盖径流中的DSi 平均浓度为27.7 μmol·L-1，显著大于冰山（8.0 μmol·L-1）；其中格陵兰Russell 冰川流域的DSi 浓度最小（5.5 μmol·L-1），Kuannersuit冰川流域的DSi浓度最大（114 μmol·L-1）（附表1~2）。DSi 浓度具有季节变化：在冰川消融初期（5 月初），径流主要（＞49%）源于“慢速流”，这时DSi 浓度较高；到5 月末，随着“慢速流”比例（12%～36%）减小，DSi 浓度相应减小且变化范围增大；到7月份，由于径流主要来自冰面且“慢速流”比例（＜18%）继续减小，此时的DSi浓度更低［35］。Aciego 等［25，64-65］发现，与北大西洋表层水相比，冰川径流中较高的DSi 浓度可能与气候变暖导致的冰川径流增加有关。DSi 浓度具有空间变化，不同的采样位置浓度会有显著差异：由于冰下水-岩作用强烈，冰下的DSi浓度比冰川河水中的浓度要大［24，32，36］。比如，南极Whillans 冰下湖中的DSi浓度为130～210 μmol·L-1［39］，大于全球大多数冰川径流（附表1）。此外，Meire 等［46］发现冰川融水进入峡湾会增大表层水的DSi 浓度，富氮的深层水上涌会改变Si/N 的比例，从而会提高下游硅藻等浮游植物的初级生产力。格陵兰冰盖径流中的ASi平均浓度为145.9 μmol·L-1（9.87～392 μmol·L-1），溢出冰川的ASi 浓度显著大于冰山（49.2 μmol·L-1），这种差异与单位质量泥沙中ASi的含量密切相关［45］。对比溢出冰川中的DSi 和ASi 浓度，发现ASi 浓度显著大于DSi。然而，Kiattuut Sermiat 冰川冰下径流中 的DSi 浓 度（22.2 μmol·L-1）大于ASi 的浓 度（9.87 μmol·L-1）［55］，这可能与该冰川流域冰下化学风化作用较弱有关。

格陵兰冰盖流域的DSi 平均通量变化范围为0.09～660 mmol·a-1，显著大于冰山；区域内最大、最小DSi 通量相差4 个数量级（附表1）。不同来源融水对DSi 通量的贡献不同，如Godthåbs 峡湾冰面融水、冰下融水和冰山的DSi 通量分别为825 mmol·a-1、130 mmol·a-1和93 mmol·a-1，说明DSi 主要来自冰面融水且冰川融水（冰面融水+冰下融水）的DSi通量远大于冰山［46］。Hawkings等［19，46，62］基于2012年Leverett 冰川的中值浓度（21.2 μmol·L-1）和模拟径流量（437 km3·a-1），估算得到格陵兰和南极冰盖的DSi 通量相当，均为0.01×106mmol·a-1，可与北极一些大河的DSi 通量相比较。比如，格陵兰及南极冰盖的DSi 通量与西伯利亚的Mackenzie 河相当（0.02×106mmol·a-1）［67］。随着气候持续变暖，冰川径流及其输送的DSi 通量增大，预计21 世纪末全球冰川释放的DSi 通量增加20%～160%［44］。据估算，格陵兰和南极冰盖的ASi 通量分别为0.16×106mmol·a-1和0.01×106mmol·a-1［19］。发现格陵兰冰盖的ASi 通量是DSi 通量的十几倍，说明格陵兰冰盖径流主要以ASi的形式输出Si元素。冰盖的Si通量（DSi+ASi，约0.2×106mmol·a-1）与大气沉降（0.5×106mmol·a-1）、地下水（0.6×106mmol·a-1）和深海热液（0.6×106mmol·a-1）的通量相当，约占陆源输入量的3%、入海总量的1.8%［19，31］。模型（ICE-6G_C）模拟结果显示，在较大时间尺度上，间冰期的Si 输出量可能会更大；格陵兰冰盖流域在末次间冰期输入到海洋的Si 通量可达5.5×106mmol·a-1，这与古河流的输出量（5.5×106～5.8×106mmol·a-1）相当［19］。在气候持续变暖背景下，入海Si 通量增大会影响全球Si循环，因此在建立全球Si收支平衡时应考 虑 冰 川 的 重 要 性［19，30-31，68］。格 陵 兰 冰 盖 流 域 的DSi 平均产量的变化范围为195～9 180 kg·km-2·a-1（附表1），Leverett 冰川流域的ASi 产量为3.6×104kg·km-2·a-1［44］。DSi 产量可能主要与流域气候和基岩性质相关。

在冰川流域（冰岛、斯瓦尔巴、斯堪的纳维亚、阿拉斯加、欧洲中部、加拿大北极、加拿大西部/美国、亚洲），径流中DSi浓度的相关研究较粗浅，大多只是列出数据，且暂未发现ASi 的相关研究报道。收集的资料显示，这些流域冰川径流中的DSi 浓度较大，平均浓度为26.5 μmol·L-1（3.03～94.8 μmol·L-1），略小于冰盖（27.7 μmol·L-1）、但显著大于冰山（8.0 μmol·L-1）；其中，斯瓦尔巴Ebbabreen 冰川的DSi 浓度最小（3.03 μmol·L-1），冰岛Solheimajokull冰川最大（94.8 μmol·L-1）。冰川的DSi平均通量变化范围为0.03～72.6 mmol·a-1，最大、最小通量相差3 个数量级；产量变化范围为57～4 025 kg·km-2·a-1（附表1）。

值得注意的是，DSi 和ASi 并不相互孤立，两者可以互相转化。冰川消融会促进基岩风化，但风化释放的ASi 结构松散且含水量高［2，38，68-72］。虽然ASi的生物活性不高，但其溶解会释放DSi［31，73-74］；在低温环境下，弱碱性水体中ASi 的溶解度远大于淡水且在河口的溶解速率更快［19，69，75］。冰川源ASi 易溶解的原因为：冰川侵蚀会破坏矿物结构、降低ASi的化学稳定性［69-70，76］；实验提取的ASi就是泥沙中化学性质活跃的Si［19］；碱金属和碱土金属会促进ASi 溶解［77］。此外，DSi 的生物利用过程也会促进ASi 溶解，底栖生物也会分解泥沙并将Si 等营养物质运输回透光带，这可能会影响更长时间尺度上的Si循环［31，73-74］。

综上所述，就DSi平均浓度的空间分布而言，在全球尺度上，DSi 浓度的大小次序为冰盖＞冰川＞冰山，DSi 的产量为冰川＞冰盖（附表2）。冰盖的DSi浓度稍大于冰川，可能与冰盖下部玄武岩的广泛分布有关［46］；冰川的DSi浓度也较大，说明冰川区基岩风化对DSi 的贡献大；冰山的DSi 浓度最小，这与冰山的含沙量较低有关［19］。由于冰川分布区域广，不同流域间的气候和基岩组成以及消融强度等的差异显著，因而冰川的DSi 变化范围较大；对于冰盖，底部的侵蚀作用会导致输出的DSi浓度增加。冰川的DSi 产量大于冰盖，说明冰川的侵蚀风化作用整体上强于冰盖。在区域尺度上，冰岛冰川的DSi 平均浓度最大、加拿大北极和加拿大西部/美国冰川最小，相差1 个数量级；冰岛较高浓度的DSi 可能与玄武岩的广泛分布且风化效率较高有关［68］。格陵兰冰盖DSi 的变化范围最大，加拿大西部/美国冰川最小；亚洲和格陵兰冰盖流域的DSi 变化范围明显大于其他流域，很可能与这两个流域的数据密集有关（附表2）。亚洲冰川的DSi平均产量最大、斯瓦尔巴冰川最小，二者相差2个数量级（附表2）。亚洲较高的DSi 产量可能与所在流域的降水量大、气温高且化学风化强度高有关［78］，斯瓦尔巴较低的DSi 产量与气温和化学风化速率均较低有关。在流域尺度上，格陵兰Kuannersuit冰川的DSi平均浓度最大、斯瓦尔巴Ebbabreen 冰川最小，二者相差2 个数量级；格陵兰Godthåbs 峡湾的DSi 产量最大、斯瓦尔巴Scott Turnerbreen 冰川最小，亦相差2 个数量级（附表1）。亚洲南部冰川的DSi 浓度较大［图2（a）］，这可能与降水较多且侵蚀作用较强有关［78］。格陵兰西部冰川的DSi 平均浓度大于东部冰川，这可能与大气环流和东/西海岸的洋流影响有关［79］。其他流域之间的DSi 浓度没有明显的空间变化规律。就ASi 平均浓度的分布而言，现有资料仅限于格陵兰冰盖的个别流域，不具有空间分析价值，故不作讨论。

冰川Si 的浓度的空间变化很大，这与采样方法、径流量和基岩类型等因素有关：

（1）采样方法。采样点的选择及其距冰川的远近会导致Si 浓度的差异，这主要与融水的稀释作用、河道内Si 的迁移转化和其他源的贡献有关。采样季节及频率会影响Si浓度，原因为Si浓度的季节变化比较明显。采样后的过滤操作也会影响Si 浓度，这主要与过滤与采样的间隔时间、滤膜孔径及材质有关。推测不同材质的滤膜也会引起Si 浓度的差异，目前尚未见报道。

（2）径流量。流量对Si 浓度的影响可归因于融水侵蚀与稀释的综合作用。ASi 的变化比较复杂，具体影响机制难以确定。虽然流量增大和含沙量增加会使ASi 浓度增大，但ASi 的溶解和融水稀释又会减小ASi 浓度。在大流域内，因为排水系统复杂、排水通道相互隔离、冰下水体滞留时间长，所以洪峰期的流量减小且洪峰出现的时间延后，这也会影响DSi和ASi的浓度变化［35，45］。

（3）基岩特征。基岩是影响Si 浓度的重要因子，原因为不同基岩的物理化学性质和风化效率有较大差异。比如，花岗岩和片麻岩流域的DSi 浓度低，玄武岩流域的DSi浓度高［50］。

2.2 硅同位素研究

δ30Si 可以识别陆地风化过程，尤其是河流系统的风化［7，13，15，54］。在研究历史时期海洋中的Si 浓度、硅藻对Si的利用水平以及Si来源的变化时，可以利用海洋沉积物中生物活性Si 的δ30Si 来反演Si 的浓度。冰川径流中δ30Si 的变化过程反映了生物地球化学循环过程的分异，可用于评估硅酸盐岩的风化强度、解析风化产物的再溶解过程、识别排水系统的特征［45］。当前针对冰川径流中δ30Si 的研究主要集中在格陵兰冰盖、冰岛、斯瓦尔巴、斯堪的纳维亚和阿拉斯加流域（图3）。

图3 全球冰川径流中δD30Si和δA30Si的空间变化。a指阿拉斯加，b指格陵兰冰盖，c指冰岛，d指斯瓦尔巴，e指斯堪的纳维亚（数据来源见附表3）Fig.3 Spatial distribution and mean values of δD30Si（a）and δA30Si（b）from glacial runoff globally.a refers to Alaska，b refers to Greenland Ice Sheet，c refers to Iceland，d refers to Svalbard，e refers to Scandinavia（datas are shown in Supplement table 3）

格陵兰冰盖径流中δD30Si 的平均值为0.036‰（-0.36‰～0.41‰），δA30Si 的 平 均 值 为-0.467‰（-0.67‰～0.467‰）（附表2～附表3）。Leverett冰川径流中δD30Si 和δA30Si 的流量加权平均值分别为-0.25‰和-0.22‰，相比其他时期，洪峰期的δD30Si减小，而δA30Si 无明显的季节变化［44］。值得注意的是，冰盖径流中的δD30Si 小于基岩的δ30Si［0.00‰（±0.07‰）］，表明冰下风化作用是其变化的主要驱动因子［40］。

冰岛、斯瓦尔巴、斯堪的纳维亚和阿拉斯加流域冰川径流中δD30Si一般大于δA30Si，与格陵兰冰盖一致，这说明基岩风化和溶解-沉降过程引起了同位素分馏［44，80-82］。尽管地理位置和风化壳组成差异显著，但风化条件造成的δ30Si 分馏最为重要［48，83-84］。基于同位素分馏模型，基岩风化速率与δD30Si 负相关且δD30Si 会随冰川覆盖度而变化，说明现代和冰期-间冰期尺度上的气候变化会影响陆源的δ30Si［13］。此外，冰川径流的δ30Si小于非冰川径流，非冰川径流的δ30Si 值与基岩接近［15，45，74，85-88］；冰期的δ30Si 小于间冰期，说明冰川变化会影响δ30Si 的时空分布［42，44］。现代海洋的δD30Si 呈增大趋势，可能与末次冰期以来冰川径流增加有关［7，55］。

在全球尺度上，冰川和冰盖的δD30Si 平均值分别 为0.225‰ 和0.036‰、δA30Si 平 均 值 分 别 为-0.4‰和-0.673‰（附表4），δD30Si 的变化范围分别为-0.58‰～0.78‰和-0.36‰～0.41‰、δA30Si 的变化范围分别为-0.86‰～-0.05‰ 和-0.67‰～0.467‰。在区域尺度上，阿拉斯加冰川的δD30Si平均值最大、格陵兰冰川最小，冰岛冰川的δA30Si平均值最大、阿拉斯加冰川最小（图3）。冰岛冰川的δD30Si 变化范围最大、斯瓦尔巴冰川最小，斯堪的纳维亚冰川δA30Si 的变化范围最大、格陵兰冰川最小（附表4）。δD30Si 变化范围的差异可能与数据量的多少有关。在流域尺度上，冰岛Soheimajokull 冰川的δD30Si 最 大、Langjokull 冰 川 最 小，二 者 相 差1.36‰；冰岛Eyjabakkajokull 冰川的δA30Si 最大、阿拉斯加Lemon 冰川最小，二者相差0.43‰（附表3，图3）。D30Si 的 变 化 范 围 大 于δA30Si，这 可 能 与δD30Si的分馏过程先于δA30Si有关，具体机理还不清楚。引起δD30Si 和δA30Si 差异的原因尚不清楚。可见，冰川δ30Si 的空间变化大，这可能与排水系统、水文路径、水体滞留时间、流域规模（冰川覆盖度）及采样方法等因素有关，还有待进一步研究［40，44-45，84，89］。

3 总结和展望

冰川Si 是全球Si循环的重要组成部分，对于重建现代和历史时期的气候环境记录、评估冰川影响Si 的生物地球化学循环乃至气候变化具有重要意义。空间分析结果表明，冰川径流中的DSi 浓度大于冰山，这与冰山中泥沙含量较少有关；溢出冰川径流中ASi 的浓度大于DSi，说明格陵兰冰盖主要以ASi的形式输出Si通量。尽管目前已取得一些研究资料和成果，但仍存在较多不足，需加强以下研究：

（1）冰川Si 的野外监测急需加强。目前对冰川Si的研究刚刚起步，监测资料稀疏/不连续且监测点的空间分布不均，导致对其时空变化格局及其通量评估等科学问题的研究仍是空白。因此急需在典型冰川流域一个完整消融期内获取各形态Si 的高频次、多类型样品，综合分析冰川Si 浓度的时间变化过程（日、月和季节）和空间分布特征，有助于深入认识冰川消融与Si释放过程的关联机制。

（2）冰川Si 的迁移转化过程研究。冰川输出的Si进入河道以后可能会受吸附、解吸附、沉淀等物理化学作用影响，导致进入下游生态系统的Si 数量急剧减小；此外，河道中Si 通过河口进入湖泊/海洋生态系统的过程中，水质的改变可能也会导致Si被“清除”。因此，应抓住冰川Si 发生明显变化的关键节点，刻画Si在河道内的迁移转化过程和运移规律。

（3）冰川Si 的通量变化规律研究。冰川输出的Si 会影响下游生态系统的生物群落结构与功能，对冰川Si 通量变化规律的研究可以更合理地评估进入下游生态系统的Si 数量，并探究不同量级的冰川Si 通量水平对生态系统的影响程度，为冰川源Si 的重要性评价提供借鉴意义。

（4）冰川Si 的来源解析研究。冰川径流中的Si除了来自冰川消融，还可能源于降水、地下水、土壤水等，此外，河道下游数公里乃至更远的河段/湖泊/海洋中是否有冰川Si 的输入尚不清楚，为此需要借助同位素手段来识别Si 的具体来源和不同源的贡献大小，有助于更好地理解冰川消融影响Si 循环的方式、评估冰川在全球Si循环中的重要性。

附表1 全球冰川径流中DSi和ASi的平均浓度（标准偏差）和通量及其变化范围Attached table 1 Mean concentrations and fluxes of DSi and ASi as well as their ranges for glacial runoff globally

续附表1

附表2 全球不同冰川区DSi的平均浓度（标准偏差）及平均产量（数据来源见附表1）Attached table 2 Mean concentrations and yields of DSi from glacial regions globally（data sources are shown in Supplement table 1）

附表3 全球冰川径流中δD30Si和δA30Si平均值（标准偏差）及变化范围Attached table 3 Mean δD30Si，δA30Si values and standard deviations in glacial runoff globally

附表4 全球不同冰川区δD30Si和δA30Si的平均值（标准偏差）及变化范围（数据来源见附表2）Attached table 4 Mean δD30Si，δA30Si values and standard deviations in glacial region globally（data sources are shown in Supplement table 2）