夏季強降雨期間千島湖有機碳的時空分布特征及影響因素

鄭 清，黎云祥，朱廣偉，俞 潔，朱夢圓，許 海，王裕成，劉明亮

1. 西華師范大學環境科學與工程學院，四川 南充 637002

2. 中國科學院南京地理與湖泊研究所，湖泊與環境國家重點實驗室，江蘇 南京 210008

3. 浙江省生態環境監測中心，浙江 杭州 310012

4. 杭州市淳安生態環境監測站，浙江 杭州 311700

5. 杭州市生態環境科學研究院，浙江 杭州 310005

湖泊是內陸水體的重要組成部分，是流域內有機碳等物質的最終匯集場所，也是大氣CO2的重要交換場所、全球碳循環的重要組成部分[1-4]. 水體有機碳是湖庫總碳的重要組成部分，主要有溶解性有機碳(DOC)和顆粒有機碳(POC)兩種基本賦存形式[5-7].其中DOC是主要存在形式[8]，能通過礦化作用釋放CO2和CH4[3]，也可為異養微生物提供碳源和能源[9]，也是天然水體中消毒副產物的重要前驅物[10]，有機碳與生物的生命過程、初級生產力等密切相關[11-14]. 鑒于此，研究湖庫有機碳的時空分布特征及影響因素不僅對量化湖庫碳循環具有重要意義，也可為深入了解湖泊生態系統動態過程提供科學依據[15-16].

湖泊有機碳的來源主要有外來和內生兩種方式[17].其中，外源有機碳主要通過降雨形成的地表徑流進入水體，而內源有機碳主要源于浮游植物、大型水草的生長合成[18]. 有研究表明，湖泊水體中有機碳的時空變化很大程度上受其來源的影響[19]，國內外大量學者研究了湖庫有機碳的時空分布及影響因素. 例如，Park等[20]研究了韓國兩個大型水庫有機碳的來源和時空特征，結果表明，氣象和水文決定了Paldang湖有機碳的狀態，而原位生成和轉化是影響Chungju湖有機碳的決定因素；Chow等[16]對翡翠水庫有機碳的長期趨勢和動態研究結果也表明，氣候和水文因素是河流DOC動態的主要驅動力；范志偉等[21]研究表明，三峽庫區DOC主要受降水和溫度的影響；王秀君等[17]研究表明，博斯騰湖DOC濃度變化主要受外部因素影響，而夏季POC濃度受浮游生物的影響較大；楊平等[22]研究表明，文武砂水庫浮游植物和外源輸入分別是引起庫區有機碳濃度季節性和空間性變化的重要因素；盧曉漩等[23]研究表明，桂林五里峽水庫DOC在豐水期低于枯水期且表層大于底層，這是受Chla和DIC濃度影響的結果. 但由于大多數研究區域面積大、水體情況復雜，且采樣頻次低、點位少，可能難以準確反映有機碳濃度的時空分布等情況.

在全球變暖，極端氣候事件發生的頻率和強度顯著提高的背景下[24]，湖庫有機碳濃度的時空格局可能會發生深遠改變[25]. 有研究表明，夏季風降雨是控制韓國大型深水水庫-昭陽湖有機碳分布的重要因素[26]，DOC濃度主要由溫度和降雨量決定[16]. 此外，由于極端降雨和高溫，貧營養水庫近年來也發生了藍藻水華[27]. 因此，對貧營養水庫而言，研究極端降水對其時空格局的影響更為關鍵，因為這類水體有機碳一般以外源輸入為主. 外來碳會顯著影響湖泊碳收支，導致外來負荷遠超過本地負荷[26]，因此研究極端天氣對水庫有機碳通量變化和收支平衡的影響非常有必要. 國內外大量研究人員估算了水體有機碳通量[28-30]，但也存在采樣頻次低、點位少等情況，并且對于有機碳賦存量研究主要集中在沉積物碳庫方面[31]，而對水相碳儲量的研究較少.

千島湖是我國長三角地區最大的水資源儲存庫，是杭州市的飲用水源地，集發電、防洪、旅游等多功能的典型清澈深水山谷型水庫[32]，水質保護要求高.但截至目前，有關千島湖水體有機碳的研究較少[33-34].1996年特大洪水導致千島湖出現極為突出的“污染峰”[35]，對水庫水質造成了極大的影響，而2020年6-7月千島湖流域又出現了建庫以來(60年)最大的強降雨事件，這更加印證了研究極端降水對千島湖以及類似水體有機碳影響研究的必要性. 該研究擬利用2020年5-8月在千島湖布設的100個監測點位有機碳監測數據，結合同期水質理化和氣象水文數據，分析了千島湖有機碳時空分布特征及影響因子，并著重分析強降雨對大型水庫有機碳濃度、通量和儲量的影響，以期為類似地理環境背景下水庫的營養鹽管控和碳效應評估提供數據支撐.

1 材料與方法

1.1 研究區域概況

千島湖，又名新安江水庫，位于浙江省淳安縣(29°22'N～29°50'N、118°36'E～119°14'E)，是我國長三角地區最大的戰略水源地及杭州市飲用水水源[36]，具有發電、防洪、旅游等多種功能，流域植被覆蓋率較高[37]. 千島湖地處亞熱帶季風氣候區，氣候溫暖，雨量充沛，年均氣溫約17.3 ℃，年均降水量約為1 733 mm[36]，從新安江入口(街口斷面)到大壩縱長150 km，水面最寬處約50 km, 平均水深31.13 m，最大水深100 m，在水位達到108 m時，水面面積為573.33 km2，庫容為178.6×108m3. 千島湖的主要入湖河流有25條，包括新安江、東源港、武強溪等，其中最大入流為新安江，約占入庫徑流的2/3. 湖區徑流集中在5-9月，占全年徑流量的77%[38]. 千島湖雖然全湖平均水質良好，但仍面臨水體富營養化風險，入湖河道、局部庫灣出現過藍藻水華[39]，部分庫灣春季硅藻生物量較高，優良水質的長期維持面臨挑戰.

1.2 樣品采集與分析

在千島湖全湖布設100個點位進行采樣(見圖1).采樣時間分別為2020年5月27-28日、6月24-25日、7月21-22日和8月24-25日. 在100個點位均分別使用塞氏盤和手持式測深儀(SM-5A, 美國Speedtech公司)現場測量透明度(SD)和水深(WD)；使用重力采水器分別采集表層(0.5 m)、次表層(葉綠素最大層)、中層(溫躍層位置)和底層(距離底泥2 m)各1.5 L水樣(該研究數據均使用每個點位各層濃度的加權平均濃度值，垂向分布處除外)，低溫避光保存，帶回實驗室用于后續分析. 5-8月樣品數分別為326、341、353和344個.

圖 1 2020年5-8月千島湖采樣點的分布Fig.1 Distribution of sampling sites in Qiandaohu Reservoir from May to August in 2020

水樣帶回室內，12 h內使用Whatman GF/F玻璃纖維膜(直徑45 mm、孔徑0.7 μm)過濾(每個樣品均通過未稱重和已烘干稱重的濾膜過濾). 原水、濾膜及濾后水立刻冷藏，及時進行測定和分析. 分析指標包括TOC、DOC、POC、Chla、CODMn、TN、TP、SS、OSS和ISS濃度.

TOC、DOC濃度采用高溫燃燒氧化法、使用multi N/C 2 100分析儀(德國耶拿公司)測定，即在水樣中添加2 mol/L鹽酸酸化(控制pH＜2)，將無機碳酸鹽分解成CO2以去除無機碳，注入高溫燃燒管(最高溫度可達950 ℃)中，分別對未過濾和過濾的樣品進行TOC和DOC濃度測定，二者之差為POC濃度.

水質指標測定參考文獻[40]，GF/F濾膜上過濾而得的藻類細胞中Chla濃度使用熱乙醇提取、分光光度法測定(波長665 nm和750 nm)；CODMn濃度采用高錳酸鉀氧化滴定法測定；TN、TP濃度采取堿性過硫酸鉀消解，分別采用紫外分光光度法(波長210 nm)和鉬銻抗顯色分光光度法(波長700 nm)測定；GF/F濾膜濾得物的SS濃度測定利用105 ℃烘干重量法測得(4 h)，烘干后濾膜再用馬弗爐450 ℃灼燒4 h，稱重，損失部分為OSS濃度，剩余灰分為ISS濃度.

1.3 數據分析

1.3.1湖區劃分

根據水庫的水文形態特征及不同區域流域土地利用情況，將千島湖劃分為6個區域(見圖1)：安徽區(AH，n=5)、西北區(NW，n=17)、中心湖區(C，n=12)、東南區(SE，n=29)、東北區(NE，n=12)和西南區(SW，n=25)，其中AH、SW、NE三個區為河流區，NW為過渡區. 由于流經安徽區的新安江支流占新安江水庫的60%[38]，依次沿AH、NW、C和SE方向變化是新安江主流向.

1.3.2降水類型劃分

為探明降水強度對千島湖有機碳的影響，根據氣象部門關于降水強度的等級劃分標準，將日累計降水量25～50 mm定義為大雨、50～100 mm定義為暴雨、100～250 mm定義為大暴雨.

1.3.3有機碳通量和有機碳儲量的計算

為探討千島湖有機碳收支變化情況，根據式(1)計算有機碳通量(FTOC、FDOC、FPOC)：式中：F為有機碳通量，t；k為常數；ci為對應點位的有機碳濃度，mg/L；Di為25條主要入庫河流和出庫流的流量，m3/s.

為計算千島湖全庫水體有機碳儲量，利用ArcGIS 10.3軟件以全湖100個點位為中心，將新安江水庫劃分為100個泰森多邊形，計算出每個區域面積，用式(2)計算有機碳儲量(RTOC、RDOC、RPOC).

式中：R為每個月千島湖全庫的水體有機碳儲量，t；Si為點位i的泰森多邊形面積，m2；Hi為點位i的水深，m.

1.3.4相關分析和空間插值

使用SPSS 23.0軟件進行有機碳和各參數之間的相關性分析(Spearman)以及非參數檢驗(Kruskal-Wallis)差異性；使用ArcGIS 10.3軟件進行Kriging插值，繪制有機碳濃度空間分布圖；利用Origin 2018、SPSS 23.0及Excel 2010軟件進行制圖和統計分析.

2 結果與討論

2.1 水文氣象及水質參數的動態特征

2.1.1水文氣象條件的時間變化

5月采樣時間距離第一次大雨(日降雨量43 mm)間隔22 d；6月采樣時間距離第一次暴雨(84 mm)、大暴雨(123.5 mm)分別間隔22和5 d；7月采樣時間距離第二次大暴雨(121.9 mm)間隔14 d，且采樣前十余天只有少量降雨；8月幾乎無降雨〔見圖2(a)〕.

入庫流量對降雨事件響應迅速，幾乎同步升高或滯后于降雨事件1～2 d. 受極端強降雨影響，千島湖7月8日平均出庫流量達7 090 m3/s〔見圖2(b)〕. 如圖2(c)所示，千島湖2020年5-8月的平均氣溫分別為23.0、25.9、26.9和30.1 ℃，總降雨量分別為240.4、528.8、509.7和25.8 mm，6月、7月的總降雨量約占歷史平均全年降雨量的2/3.

2.1.2水質參數的時間變化

整體上6月后各水質參數(TN和SD除外)濃度都明顯(P＜0.01)升高，7月Chla、SS、OSS、ISS和TP濃度均處于研究期間的最高水平. 其中，5-7月SS和ISS濃度不斷升高(P＜0.05)，8月顯著降低(P＜0.01)；SD則相反，5-7月逐漸降低而8月升高. 此外，與7月相比，8月各水質參數(CODMn和SD除外)濃度都明顯(P＜0.01)降低(見表1).

2.1.3水質參數的空間變化

5-8月TN、TP、SS和ISS濃度具有相似的空間變化趨勢，即安徽區到東南區逐漸下降(見圖3). 7月上述參數濃度下降較快且濃度明顯高于其余3個月，安徽區顯著較高，分別是東南區的2、3、9和17倍. SD相反，安徽區到東南區具有上升趨勢，河流區(東北區、西南區)開始降低. CODMn濃度在安徽區到東南區有一定的下降趨勢，OSS和Chla濃度處于波動變化并且變化趨勢相似.

2.2 有機碳濃度的時空變化

2.2.1有機碳濃度的時間分布特征

圖 2 2020年5—8月千島湖日降雨量、出入庫流量總降雨量和平均氣溫的變化情況Fig.2 Variation of daily rainfall, flow (i.e., inflow and outflow), total rainfall and mean temperature in Qiandaohu Reservoir from May to August in 2020

表 1 2020年5—8月千島湖各水質參數的變化Table 1 Variation of water quality parameters in Qiandaohu Reservoir from May to August in 2020

研究期間TOC、DOC濃度最大值出現在7月，且顯著高于其余3個月(P＜0.01)，7月后顯著降低(P＜0.01)，而POC濃度有逐漸升高趨勢(見圖4)；5-7月3種有機碳雖呈逐漸升高趨勢，但5月、6月濃度相差不大(P＞0.05). 研究期間TOC濃度范圍為1.19～3.94 mg/L，平均值為(2.06±0.48) mg/L；DOC濃度范圍為0.91～3.25 mg/L，平均值為(1.73±0.39) mg/L；POC濃度范圍為0～1.35 mg/L，平均值為(0.33±0.21) mg/L.

千島湖夏季DOC濃度平均值(1.78 mg/L)低于密云水庫(2.61 mg/L)[9]、水布埡水庫(2.22 mg/L)[7]、三峽庫區(3.86 mg/L)[21]、Tingalpa淺水水庫(暴雨導致雨季月均值為14.3 mg/L)[41]、東湖(7.31 mg/L)[2]等湖泊水庫；POC濃度(0.33 mg/L)也遠小于博斯騰湖(0.71 mg/L)[17]、巢湖(最低平均濃度1.38 mg/L)[42]等湖泊，這可能與千島湖良好的水質狀況和高森林覆蓋率以及高水土保持等條件有關[37].

2.2.2有機碳濃度的空間分布特征

垂向分布上，5-8月TOC、DOC、POC濃度的垂直分布都具有表層＞中層＞底層的規律(見圖5).

圖 3 2020年5—8月千島湖各水質參數的空間變化特征Fig.3 Spatial variation characteristics of water quality parameters in Qiandaohu Reservoir from May to August in 2020

圖 4 2020年5—8月千島湖有機碳濃度的月變化特征Fig.4 Monthly variation characteristics of organic carbon concentration in Qiandaohu Reservoir from May to August in 2020

水平分布上，5-7月有機碳濃度高值沿新安江入庫水流向逐漸降低(見圖6)，整體具有河流區到湖泊區逐漸降低趨勢，并且5-7月有機碳濃度高值在全庫的分布范圍逐漸擴大. 其中，5月DOC濃度高值沿新安江主線路輸移到三潭島，6月擴散到大壩，7月擴散到全庫. 此外，5-7月POC濃度高值的輸移距離始終小于DOC濃度. 8月DOC濃度高值主要集中在東南區、東北區和西南入庫區，POC濃度高值主要集中在中心區、東北區和東南區.

2.3 有機碳通量和碳儲量

圖 5 2020年5—8月千島湖有機碳濃度的垂向分布特征Fig.5 Vertical distribution characteristics of organic carbon concentration in Qiandaohu Reservoir from May to August in 2020

新安江是千島湖全庫25條主要入庫河流有機碳碳通量的重要貢獻河流，研究期間3種有機碳通量平均值都約占全庫總入庫碳通量的69%，并且與總入庫碳通量的變化趨勢幾乎一致，呈5-7月逐漸升高而8月降低的特征〔見圖7(a)〕. 具體來說，5-8月總入庫FTOC分別為2 863、9 796、11 122和268 t，5-7月FTOC分別是8月的11、36、41倍，FDOC分別為2 252、8 064、9 300和217 t，FPOC分別為610、1 732、1 823和52 t.

全庫水體有機碳儲量如圖7(b)所示，研究期間3種有機碳儲量平均值分別為44 611、38 452和6 159 t，RDOC占RTOC的86%，結合DOC濃度與POC濃度的比值(約為5)，表明DOC是千島湖水體有機碳的主要組成部分. 5-8月RTOC分別為37 300、41 724、59 399和40 023 t，RDOC分別為33 089、36 428、52 874和31 419 t，RPOC分 別為4 212、5 296、6 525和8 604 t.5-8月總入庫FTOC分別占RTOC的7.7%、23.5%、18.7%和0.7%，6月和7月入庫通量均約占當月全庫儲量的1/5，分別是8月的35和28倍.

2.4 千島湖有機碳濃度時空分布的影響因素

Spearman相關分析表明，千島湖DOC、POC濃度均與Chla、CODMn、SS、OSS和TP濃度呈極顯著正相關(P＜0.01，見表2和圖8)，說明這些水質參數是影響千島湖有機碳濃度和分布的關鍵因素. 在千島湖中，Chla和OSS濃度主要表征浮游植物的生物量，因此Chla、OSS濃度均與DOC、POC濃度呈極顯著線性正相關(P＜0.01)，表明浮游植物內源生產是千島湖有機碳的關鍵來源[43]. 千島湖的研究結果符合傳統認知，即浮游植物的內源釋放是湖庫有機碳的主要來源之一[9]. 鑒于此，千島湖有機碳濃度的時空分布受浮游植物光合作用強度的影響. 因此，可以直接或間接影響浮游植物生物量的環境條件(如TP)，也會成為影響千島湖有機碳濃度時空分布的顯著因子之一.

另一方面，湖庫水體有機碳也來自陸源有機物輸入，即降雨對流域的沖刷作用也是影響湖庫有機碳的關鍵因素[7]. 對千島湖而言，DOC、POC濃度與ISS濃度均呈顯著正相關(P＜0.05) (見表2)，而已有研究[38]發現，千島湖ISS濃度主要與泥沙含量有關，這主要受到河流入庫等外源輸入影響. 綜合說明千島湖有機碳濃度與外源性有機物輸入具有顯著關聯性. 另外，DOC、POC濃度均與SD呈顯著線性負相關(P＜0.01)(見圖8). 實質上，SD反應了藻類濁度(即Chla濃度)和非藻類濁度(如ISS濃度)的綜合效應[44]，因此有機碳濃度與SD的反相關，反映出千島湖有機碳濃度的變化受內源藻類生長和外源輸入的綜合影響.

降雨是外源性有機碳輸入的重要驅動力，尤其是強降雨過程帶來的流域沖刷作用[21,26]. 筆者研究發現，強降雨過程將大量陸源有機物沖入千島湖，對有機碳濃度的時空分布產生較大影響. 降水對不同區域的影響不同，對入庫河流區的沖刷影響較大，對湖泊庫區影響較小. 具體來說，5-7月千島湖有機碳濃度的高值都出現在新安江入庫河流附近水域，而遠離新安江入庫口區域有機碳濃度相對較低(如大壩、中心湖區等)(見圖6). 相比較而言，POC濃度輸移距離始終小于DOC濃度(見圖6)，主要是由于POC通過坡面流與泥沙耦合運移，而DOC通過坡面流和地下流與水耦合運移，導致二者有不同的流路[45]. 此外，5-7月有機碳濃度的高值隨湖流輸移分別到三潭島、大壩和全庫(見圖6)，說明月間高值的輸移距離不同，這可能是由于采樣與降雨事件的時間間隔及降雨強度存在差異導致的[46].

圖 6 2020年5—8月千島湖有機碳濃度的空間分布特征Fig.6 Spatial distribution characteristics of organic carbon concentration in Qiandaohu Reservoir from May to August in 2020

同時，該研究還發現，距大暴雨事件較近的6月，DOC、POC濃度與Chla濃度的相關系數均小于5月(見表2)，這可能與強降水導致水體SD降低，以及水體穩定性較差等不利于浮游植物生長有關[17]. 相對地，在基本無降雨的8月，遠離河口的中心區與東南區有機碳濃度相對較高，這可能與無降雨影響水動力條件(如相對較低的流速)，以及高溫利于藻類生長產生內源釋放有關[22]. 另外，淳安縣居民區生活污水的排放也可能是促進上述區域藻類生長及有機碳濃度進一步增加的重要原因[34]. 綜上，千島湖有機碳濃度的時空分布主要與藻類生長過程和流域外源輸入過程有關. 其中，浮游植物生長過程受營養鹽、透明度、溫度等多種因素的綜合影響，而流域外源輸入主要與降水有關，尤其是強降雨過程.

2.5 降雨對千島湖有機碳濃度和通量的影響

圖 7 2020年5—8月千島湖有機碳通量和碳儲量隨月份的變化Fig.7 Monthly variation of organic carbon flux and carbon storage in Qiandaohu Reservoir from May to August in 2020

表 2 2020年5—8月有機碳濃度與環境因子的相關性分析Table 2 Correlation analysis between organic carbon concentration and environmental factors from May to August in 2020

圖 8 2020年5—8月有機碳濃度與水質參數的相關性分析(n=400)Fig.8 Correlation analysis between organic carbon concentration and water quality parameters from May to August in 2020(n=400)

調查期間，5月采樣前發生了4次大雨、多次小雨，6月臨近采樣前發生了暴雨和大暴雨，7月采樣前十余天發生過暴雨和大暴雨，8月幾乎無降雨. 筆者研究表明，有機碳濃度對不同降雨強度的響應明顯不同. 具體來說，千島湖非雨季時DOC的背景濃度為1.24 mg/L(筆者所在課題組未發表的數據)；相比較而言，5月DOC濃度增幅為28%(1.59 mg/L)，這與千島湖以往研究中30次暴雨事件的DOC濃度平均值相差較小(1.55 mg/L)[33]，說明持續性大雨將大量外源有機質沖進水體[21,26]. 與5月相比，6月千島湖DOC濃度增幅 (1.3%，P＞0.05)較小. 究其原因可能是：①千島湖流域的水土保持率較高[37]；②大量雨水的稀釋作用[47]以及水動力條件不利于內源生產；③水力停留時間短，外來有機碳難以在湖泊中滯留[20]. 值得一提的是，同樣發生了強降水的7月，DOC濃度卻顯著升高(40%，P＜0.01)，這可能是由于強降雨帶來的地表有機碳和營養鹽擴散到全庫(見圖7)[46]，且采樣前十余天只有少量降雨，水動力條件和相對高溫有利于浮游植物生長，從而引起碳滯留和大量內源有機碳的產生(DOC、POC濃度與Chla濃度的相關性均大于6月，見表2).此外，該現象可能也與降雨強度和降雨頻次對流域水土保持的嚴重破壞存在某個臨界點有關. 相對地，幾乎無降雨的8月DOC濃度顯著低于7月(P＜0.01)，這可能是由于降雨后更高的溫度往往會提高微生物活性，從而導致有機碳的分解大于累積[17].

分析降雨對有機碳通量的影響可能更具有科學和現實意義[2,48]. 強降雨期間，來自流域的外來碳會顯著影響湖泊碳收支[26]. 強降雨及其伴隨的高流量是千島湖有機碳通量升高的關鍵控制因子[48]，并且高入庫碳通量會對全庫水體碳儲量產生強烈沖擊. 雖然6月有機碳濃度與5月相比僅略微增加，但6月的總入庫TOC通量是5月的3.4倍且達到顯著水平(P＜0.01)，這是由于6月更大的入湖流量所致[18]. 從數值比例來看，有暴雨及大暴雨的6-7月入庫TOC通量分別是8月的36和41倍，都約占當月全庫水體碳儲量的1/5，所占比例分別是未降雨月份(8月)的35和28倍；由于強降雨的影響，月間出庫碳通量差異也較大，7月3種有機碳的總入庫通量與出庫通量差值均為負值，說明開閘泄洪使水庫階段性表現出有機碳碳源功能(以往研究中千島湖常作為DOC碳匯)[33]. 相比較而言，7月新安江入庫和出庫DOC通量分別是千島湖以往研究中年入庫和年出庫通量的39%和1.43倍[33]，這與國外一些研究中風暴事件導致森林流域、溪流等DOC輸出占全年的24%～50%[49-50]相當.但是，千島湖有機碳對強降雨的響應程度遠超英國泥炭地集水區[48]、青藏高原泥炭地排水溝[45]、西江[18]等地的有機碳對洪水過程的響應.

氣候變化背景下，極端降雨事件可能會更加頻發[51].同樣需要關注的是，極端降雨帶來的碳、氮、磷是同步輸入的(見表1)，這會為藻類生長提供大量的生源要素，加之降雨過后的高溫以及平靜的水動力條件(一般夏季)，十分有利于藻類異常增殖，這可能導致貧營養水庫發生水華[27]，從而引起藻毒素、異味物質等多種經濟生態負面效應. 此外，強降水會帶來更多的外源有機碳，這會增加水處理成本和消毒副產物的形成風險[10]. 綜上，筆者認為，人們需要額外關注極端氣候事件對水庫碳循環和水生態的影響，同時需要構建相關應急防控體系. 強降雨事件對千島湖水體有機碳濃度、分布、通量和儲量均影響較大，探究其變化規律及影響因素對于全球碳循環的認識與研究具有重要意義.

3 結論

a) 2020年5-8月千島湖全湖水體TOC、DOC和POC濃度的平均值分別為2.06、1.73和0.33 mg/L，隨著強降雨開始，5-7月TOC、DOC濃度呈逐漸上升趨勢，而雨量急劇下降的8月，其濃度均顯著下降；水平分布上，5-7月有機碳濃度高值具有河流區到湖泊區逐漸降低的趨勢.

b) 新安江入庫碳通量約占全庫25條主要河流總入庫碳通量的69%；強降雨及其伴隨的高流量是有機碳通量升高的關鍵控制因子，并且對全庫水體碳儲量產生強烈沖擊. 降雨期間5-7月總入庫FTOC分別是8月的11、36和41倍；6月、7月入庫碳通量都約占當月全庫水體碳儲量的1/5，所占比例分別8月的35和28倍；RDOC約占RTOC的86%，DOC是千島湖水體有機碳的主要儲存形式.

c) 浮游植物內源生產及外源輸入均是千島湖有機碳的重要來源，有機碳濃度的時空分布受水文和氣象因素的綜合影響，其中強降雨過程是關鍵的驅動因素.