全球冰川徑流中硅的研究進展

劉 莎， 李向應， 楊船洋， 韓添丁， 井哲帆， 朱永華

（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年來，冰蓋的冰儲量持續減少，冰川普遍退縮，凍土層厚度逐漸減?。?］。在冰川加速消融背景下，夾帶著泥沙的融水徑流向下游輸送了大量的生物活性元素（如Si、C、Fe、N 和P），當這些元素進入水庫、湖泊和海洋后［7］會影響生態系統和生物地球化學循環［11-21］?；钚許i元素可被硅藻等硅質海洋生物吸收利用并形成細胞結構，細胞死亡后將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 μ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%、入?？偭康?.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 的變化范圍最大，加拿大西部/美國冰川最?。粊喼藓透窳晏m冰蓋流域的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］?；谕凰胤逐s模型，基巖風化速率與δD30Si 負相關且δD30Si 會隨冰川覆蓋度而變化，說明現代和冰期-間冰期尺度上的氣候變化會影響陸源的δ30Si［13］。此外，冰川徑流的δ30Si小于非冰川徑流，非冰川徑流的δ30Si 值與基巖接近［15，45，74，85-88］；冰期的δ30Si 小于間冰期，說明冰川變化會影響δ30Si 的時空分布［42，44］?，F代海洋的δ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平均值最大、阿拉斯加冰川最?。▓D3）。冰島冰川的δD30Si 變化范圍最大、斯瓦爾巴冰川最小，斯堪的納維亞冰川δA30Si 的變化范圍最大、格陵蘭冰川最?。ǜ奖?）。δ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 的生物地球化學循環乃至氣候變化具有重要意義?？臻g分析結果表明，冰川徑流中的DSi 濃度大于冰山，這與冰山中泥沙含量較少有關；溢出冰川徑流中ASi 的濃度大于DSi，說明格陵蘭冰蓋主要以ASi的形式輸出Si通量。盡管目前已取得一些研究資料和成果，但仍存在較多不足，需加強以下研究：

（1）冰川Si 的野外監測急需加強。目前對冰川Si的研究剛剛起步，監測資料稀疏/不連續且監測點的空間分布不均，導致對其時空變化格局及其通量評估等科學問題的研究仍是空白。因此急需在典型冰川流域一個完整消融期內獲取各形態Si 的高頻次、多類型樣品，綜合分析冰川Si 濃度的時間變化過程（日、月和季節）和空間分布特征，有助于深入認識冰川消融與Si釋放過程的關聯機制。

（2）冰川Si 的遷移轉化過程研究。冰川輸出的Si進入河道以后可能會受吸附、解吸附、沉淀等物理化學作用影響，導致進入下游生態系統的Si 數量急劇減?。淮送猓拥乐蠸i 通過河口進入湖泊/海洋生態系統的過程中，水質的改變可能也會導致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）