近30 a雅魯藏布江流域高寒濕地動態變化及其對氣候變化的響應

劉　冬，王　濤，沈渭壽，林乃峰，鄒長新①

(1.環境保護部南京環境科學研究所，江蘇 南京　210042；2.南京信息工程大學地理與遙感學院，江蘇 南京　210044)

?

近30 a雅魯藏布江流域高寒濕地動態變化及其對氣候變化的響應

劉冬1，王濤2，沈渭壽1，林乃峰1，鄒長新1①

(1.環境保護部南京環境科學研究所，江蘇 南京210042；2.南京信息工程大學地理與遙感學院，江蘇 南京210044)

摘要：利用1980、1990、2000和2010年4期Landsat、環境一號衛星遙感影像提取雅魯藏布江(以下簡稱“雅江”)流域濕地面積動態變化信息,并結合實地調查和近30 a氣象、水文、積雪數據研究濕地動態變化對氣候變化的響應特征。結果表明，2010年雅江流域高寒濕地面積占流域總面積的3.73%,其中以沼澤濕地面積最大(59.46%),其次為河流(30.19%)、湖泊(10.13%)和人工濕地(0.23%);近30 a間雅江流域湖泊濕地和河流濕地面積總體上呈增加趨勢,而沼澤濕地面積呈降低趨勢;近30 a來雅江流域年均溫、年均最高溫和年均最低溫均呈顯著增加趨勢,降水量波動較大,呈輕微增加趨勢但不顯著,平均相對濕度呈不顯著下降趨勢。年均氣溫和年均最高氣溫這2個因素與雅江流域湖泊面積變化相關性顯著,河流濕地面積、雅江5大支流干流河長的增長分別與年均溫、年均最高溫、年均最低溫相關性顯著,沼澤濕地面積與該區域氣溫以及平均相對濕度相關性顯著,開發大壩、水庫等國家建設戰略是雅江流域人工濕地面積增加的主要驅動力。

關鍵詞：高寒濕地；雅魯藏布江；氣候；響應

氣候變化能夠顯著影響濕地的各種生態過程和濕地生產系統的生產力,是控制濕地動態變化的動因[1]。作為全球海拔最高的一個獨特地域單元,青藏高原是我國高原濕地的最主要分布區,區域內的濕地不僅是中國長江、黃河和瀾滄江等大江大河的發源地,還是眾多珍稀物種的棲息地和候鳥遷徙地,具有生態蓄水、水源涵養和氣候調節等重要的生態功能[2-3]。目前,國內針對高寒濕地的研究還相對較少,對青藏高原高寒濕地區域尺度的研究主要集中于若爾蓋高原濕地[4-5]、拉薩拉魯濕地及拉薩河流域濕地[6]、黃河與長江等三江源源區濕地[7]以及藏北典型湖泊[8]。在已有的相關研究中,對高寒濕地定義的理解也不盡相同,對濕地類型的劃分也因研究尺度及研究區域的不同而出現差異。關于西藏高原上大面積分布的濕地生態系統對氣候變化的響應機制問題,目前尚沒有野外實驗和遙感等直接數據能夠給出確切答案。

考慮到氣候變化對西藏高原生態系統的影響和西藏高原生態環境的脆弱性、戰略性,選擇雅魯藏布江(以下簡稱“雅江”)濕地為研究對象,通過大量野外實地調查,以1980、1990、2000和2010年4期遙感影像為數據源,采用人機目視解譯方法提取雅江流域濕地的類型數據,對其空間分布及近30 a的演變趨勢進行分析。雅江是西藏最大的河流,是西藏的主要淡水來源和水汽通道[9],雅江上中游人口稀少,研究湖泊、河流和沼澤等水系的面積變化趨勢能夠很好地反映氣候變化對高寒濕地生態系統的影響,深入認識青藏高原的濕地動態變化特征以及對全球氣候變化的響應關系。同時雅江流域下游還跨越溫帶和亞熱帶濕潤氣候區,下游河谷地區是西藏自治區最重要的經濟社會發展區域[10],分析雅江全流域各種類型濕地的分布及動態變化對于西藏當地居民的水源供給和生存環境十分必要。

1研究區概況

雅江是一條重要的國際河流,發源于西藏高原西南部、喜馬拉雅山脈北麓海拔5 200 m的杰馬央宗冰川,流域平均海拔超過4 000 m,是世界上海拔最高的大河[11]。雅江是雅江流域橫跨西藏的拉薩、日喀則、山南地區的大部分和阿里、那曲、昌都地區的一小部分,涉及西藏41個縣(市)。雅江流域地勢西高東低,海拔垂直差異大。根據水文特點及河谷地貌特征,以昂仁與拉孜縣界、桑日與加查縣界為分界點,將雅江流域劃分為上游區、中游區和下游區3個生態分區,其中雅江上游地區平均海拔4 600 m以上,以高山谷地地貌為主,屬高原亞寒帶半干旱氣候區[10];雅江中游地區在地貌上屬于藏南谷地,主要為溫暖半干旱氣候[12]; 雅江下游地區屬高原溫帶濕潤、半濕潤氣候區和山地亞熱帶濕潤氣候區,氣候溫暖濕潤,是生物物種分化變異、生物區系最豐富的地方[13]。

2材料與方法

2.1數據源

高寒濕地的解譯選用了4期遙感影像數據:1980年的美國陸地衛星Landsat MSS影像21景,1990年美國陸地衛星Landsat TM影像21景,2000年美國陸地衛星Landsat ETM+影像21景,2010年環境一號衛星CCD影像6景,影像采集時間為2011年6月。1980年的MSS數據分辨率為78 m,1990、2000和2010年的TM/ETM+/環境一號衛星數據分辨率為30 m。

雪深數據來自1979—2010年被動微波遙感SSM/I逐日雪深資料,下載于寒區旱區科學數據中心(http:∥westdc.westgis.ac.cn/),空間分辨率為25 km。

氣象數據是從國家氣象部門獲取的1980—2010年雅江流域地區8個氣象站(江孜、日喀則、澤當、拉薩、林芝、波密、當雄和嘉黎)逐日氣象實際觀測資料,包括平均氣溫、最高氣溫、最低氣溫、降水量、平均相對濕度和年均蒸發量,其中年均蒸發量數據截至2005年。

以奴各沙、羊村和奴下作為代表雅江流域上、中、下游地區的水文站,這3個水文站數據資料較詳實,基本能夠反映全流域的徑流狀況。采用這3個水文站1980—2010年實測逐月天然徑流量系列資料探討雅江各流域段年徑流量變化的基本特征。數據由西藏自治區水文局提供。

研究采用的其他數據包括西藏30 m分辨率的DEM高程數據、西藏2000年1∶10萬土地利用類型圖和1∶100萬青藏高原植被類型圖等。

2.2研究方法

2.2.1濕地分類

參考《關于特別是作為水禽棲息地的國際重要濕地公約(Ramsar公約)》，將雅魯藏布江流域濕地分為河流濕地、湖泊濕地、沼澤濕地和人工濕地4類(表1)。

2.2.2濕地信息提取

遙感圖像預處理主要包括幾何校正、鑲嵌裁剪、最佳波段組合選擇和圖像增強。利用EARDAS軟件對影像進行預處理,采用最小二乘法分別對1980、1990、2000和2010年遙感影像進行幾何精校正,像元重采樣采用最鄰近法或者雙線性差值法,從而消除不同遙感影像空間分辨率不同所造成的影響,使遙感影像幾何精校正誤差不超過1～2個像元。

表1雅魯藏布江流域高寒濕地分類

Table 1Types of the alpine wetlands in the Yarlung Zangbo River Valley

分類系統一級二級三級類型特征解譯標志天然濕地河流濕地河流永久性或暫時性的河流及其支流,以及河流泛濫淹沒的兩岸地勢平坦地區,包括河灘及泛濫河谷采用432波段合成,河道呈黑色或藍褐色,紋理清晰,呈狹長彎曲線狀分布,河流兩旁的河灘濕地沿河道兩側呈不規則形狀分布,有較清晰邊界湖泊濕地湖泊面積大于1km2,常年或季節性積水的湖泊采用432波段合成,湖泊呈深藍色,顏色細膩,表面光滑,水陸界線清楚,色調均勻沼澤濕地苔草沼澤濕地草丘和水坑相間分布,地表長期處于積水狀態,植物生長密集,平均蓋度在90%左右,以青藏苔草為優勢種TM/ETM影像采用753波段合成,MSS/環境衛星影像采用431波段合成,呈綠色或灰綠色不規則斑塊狀圖案,多可見密密麻麻的小面積水體點綴分布嵩草沼澤濕地草丘和水坑相間分布,地表長期處于積水狀態,植物生長密集,平均蓋度在90%左右,以藏北嵩草為優勢種TM/ETM影像采用753波段合成,MSS/環境衛星影像采用431波段合成,呈暗灰綠色不規則斑塊狀圖案,紋理較粗糙人工濕地水庫水庫能攔水、蓄水,人類開發建造并供人類利用的水域采用432波段合成,棱角分明、形狀規則的幾何形邊界的較小面積水域稻田稻田人類開發種植水稻的區域采用753波段合成,呈亮綠色,小且成片分布在河流或河灘周圍的農田區域

研究組通過大量野外實地調查,根據雅江流域各類型濕地景觀的光譜和幾何特征,借助 1∶10萬西藏土地利用類型圖、1∶100萬青藏高原植被類型圖和30 m分辨率的青藏高原DEM高程圖,建立雅江流域濕地分類解譯標志(表1),在EARDAS軟件的支持下對遙感影像進行人機交互目視解譯。在濕地信息提取中,河流寬度不小于2個像素(60～120 m),湖泊面積不小于3×3個像素(0.09～0.36 km2)。在解譯1980、1990和2010年濕地數據時以2000年的解譯結果作為參考,在解譯時進行對比及檢查。

2.2.3濕地空間動態格局變化分析

采用變化動態度作為確定濕地不同類型不同時期變化幅度的分析指標,其計算公式為

(1)

式(1)中,D為濕地變化動態度,%;Ua為起始年濕地面積,km2;Ub為終結年濕地面積,km2;t為間隔時間,a。

2.2.4氣候分析

采用氣候傾向率法,以氣候要素為因變量,時間為自變量,建立2個變量之間的一元回歸(或N階曲線模擬)的趨勢方程,計算模擬序列與原序列的相關系數及傾向變化率等,并比較各參量的變化特征,對不同流域段氣溫、降水量、濕度和蒸發量等的變化趨勢進行對比分析;計算各年份氣溫、降水量、濕度和蒸發量等氣象因子值與研究時段各因子平均值的差值,作為代表各因子的距平,研究雅江流域氣候的歷史演變規律和年代際變化特點。

2.2.5雪深數據提取及分析

利用ArcGIS 9.3軟件將被動微波SSM/I反演的雪深數據統一轉換為Albers投影,通過柵格運算獲取月尺度平均雪深數據,然后基于雅江流域區劃對月尺度雪深數據進行裁切,提取雅江流域的雪深信息。采用平均值法統計雅江流域內所有像元的逐年雪深數據,得到雅江流域地區雪深年際序列。采用一元線性回歸分析其變化趨勢,建立雪深(y)與時間序列(年份x)的一元線性回歸方程。

2.2.6灰色關聯度分析

以灰色系統2要素歷史數據序列之間的關聯度來表征高寒濕地面積動態與氣候變化之間的密切程度。關聯度的計算公式為

(2)

(3)

式(2)～(3)中,Rom為子序列m與母序列o的關聯度;N為數據序列長度,即數據個數;Rom(t)為子序列m與母序列o在t時刻的關聯系數;Δmin和Δmax分別為各時刻2個序列絕對差中的最小值和最大值;Δom(t)為t時刻2個序列的絕對差;ρ為分辨系數,ρ∈(0,1),一般取0.1～0.5。

采用灰度關聯分析時,對資料做如下處理:首先根據前文4個時期濕地面積變化動態度求逐年濕地面積;其次,對所選取的氣候因子和濕地面積原始數據進行標準化處理;該研究中分辨系數取值0.5,對標準化數據進行關聯度分析,得到濕地面積變化與各驅動因子的灰色關聯度。

3結果與分析

3.1雅江濕地分類及動態特征

由表2可知,2010年雅江流域濕地總面積為9 033.53 km2,占流域總面積的3.73%。其中,沼澤濕地面積最大,占濕地總面積的59.46%;其次是河流濕地,占濕地總面積的30.19%;湖泊濕地面積為915.46 km2,占流域內濕地總面積的10.13%;流域內人工濕地分布較少,僅占濕地總面積的0.23%。

從1980—2010年4個研究時期雅江流域不同濕地類型的面積動態特征(表2)可知,1980、1990、2000和2010年濕地總面積分別為8 331.18、8 504.04、8 404.28和9 033.53 km2,濕地面積總體呈上升趨勢。雅江流域是青藏高原典型高寒沼澤濕地的主要分布區,但過去30 a間沼澤濕地面積一直呈降低趨勢,總共減少81.74 km2。

2010年雅魯藏布江流域濕地空間分布見圖1。雅江流域上游地區濕地總面積為6 195.39 km2,占全流域濕地總面積的68.58%。其中河流、湖泊和沼澤這3類濕地占該區域濕地總面積的90%以上,是上游分布最廣的濕地類型。中游地區濕地面積占全流域濕地總面積的17.0%,河流和沼澤濕地分別占60.53%和37.35%。下游地區是流域內濕地面積最少的區域,其中沼澤和河流濕地分別占45.76%和38.53%;水稻田人工濕地僅在下游地區的墨脫和察隅縣境內有零星分布,面積為14.67 km2。

表21980—2010年雅魯藏布江流域濕地面積變化

Table 2Dynamic change in the area of wetlands in the Yarlung Zangbo River Valley during 1980-2010

年份濕地面積/km2動態變化度/%河流湖泊沼澤人工濕地河流湖泊沼澤人工濕地19802188.97681.855452.677.6919902361.90709.325422.5010.310.790.40-0.063.4120002334.79660.775393.4515.26-0.11-0.68-0.054.8120102726.81915.465370.9420.331.683.85-0.043.32

圖1　2010年雅魯藏布江流域濕地空間分布

1980—2010年雅江流域各時段湖泊數量變化見表3。湖泊濕地在近30 a總體呈不斷增加趨勢,總增幅為34.36%。1980—1990年湖泊濕地面積從681.85增加到709.32 km2,變化動態度為0.40%,面積大于1 km2的湖泊數從167增長到172個;1990—2000年湖泊濕地面積出現下降,變化動態度為-0.68%,湖泊數從172降低到162個,減少的12個湖泊都是面積小于10 km2的小型湖泊;2000—2010年湖泊濕地面積重新出現大幅度的增長趨勢,變化動態度為3.85%,湖泊數從162增加到199個(表3)。

河流濕地面積在近30 a間總體呈增長趨勢,總面積增加537.84 km2,增幅為24.57 %。但1990—2000年河流面積出現階段性小幅萎縮,面積減少27.11 km2。2000年以后,河流面積又恢復增長趨勢,2000—2010年增加392.02 km2。由圖2可見,雅江流域內5個最大支流的河長在近30 a間也出現不斷增加趨勢,拉薩河干流1980年河長為333.17 km,到2010年增長29.69%;尼洋河、年楚河、帕隆藏布和多雄藏布的干流河長在近30 a間也分別增長90.58、50.39、75.93和95.81 km。

總體而言,雅江流域人工濕地面積呈現穩步增加趨勢,面積在近30 a增加12.64 km2,增幅達164.37%,且該變化主要發生在1990—2000年,這個時段的變化動態度高達4.81%。

表31980—2010年雅魯藏布江流域湖泊數量變化

Table 3Change in the number of lakes in the Yarlung Zangbo River Valley during 1980-2010

湖泊面積/km2湖泊數量1980年1990年2000年2010年<18687821041～1075777287>10～304666>30～502222總計167172162199

圖2　1980—2010年雅魯藏布江流域

3.2雅江流域氣候變化特征

由圖3可見,1980—2010年雅江流域氣溫呈顯著上升趨勢,近30 a雅江流域年均溫氣候傾向率為0.50 ℃·(10 a)-1,其中雅江中游地區增溫最為顯著,氣候傾向率為0.50 ℃·(10 a)-1,這與已有報道的青藏高原從20世紀60年代以來呈逐漸升溫態勢的研究結論一致[14]。從年均溫的距平來看,1980—1994年期間各年份年均溫以負距平為主,即普遍低于整個流域多年年均溫的平均值,而在1995年后氣溫明顯升高,以正距平為主。年均溫的年際變化特征為:1980年代為相對低溫階段,1990年代中后期進入相對高溫階段,2000—2005年升溫不明顯,2005年之后為近30 a來最溫暖時期。氣溫總體上逐漸上升,上、中、下游地區年均溫分別大約以0.45、0.62和0.41 ℃·(10 a)-1的速率上升,明顯高于全國和全球的增溫速率,尤其是海拔3 000 m以上的上游及中游地區升溫最快,高于青藏高原近40 a以來的增溫速率〔0.44 ℃·(10 a)-1〕[15],可見在整個青藏高原地區雅江流域的氣候變暖趨勢也非常明顯。

近30 a雅江流域年均最高溫從14.33升高到16.24 ℃,增幅為13.29%,年均最低溫從0.28升高到1.96 ℃,也增長近7倍。與年均溫的距平趨勢相似,年均最高溫和最低溫自1995年之后普遍以正距平為主,高于全流域年均溫的平均值,年均最高溫的極值出現在2009年的下游地區,達17.74 ℃。

近30 a雅江流域年均降水量波動較大,呈現增加趨勢但不顯著,全流域年均降水量的氣候傾向率為9.70 mm·(10 a)-1,這與已有文獻中近30 a青藏高原降水大體上以增加趨勢為主,但總體顯著性水平不高的研究結論相一致[14,16]。近30 a雅江下游地區平均年降水量為794.34 mm,明顯高于上游(354.92 mm)和中游地區(436.50 mm)。全流域多年平均降水量呈現從偏多到偏少的周期性波動特征,這一特征在5 a滑動平均序列圖上反映更清楚:從1980 年代中期開始,曲線一直呈上升趨勢,1990年代初期有所降低,從1999年開始逐漸上升,之后又呈現降低趨勢,負距平開始增多。

全流域年均相對濕度波動較大,總體上隨時間變化呈現下降趨勢,但降幅較小,未通過0.05水平的顯著性檢驗,氣候傾向率為-0.43%·(10 a)-1。而下游地區年均相對濕度顯著高于上、中游地區(P<0.01)。1980—1990年年均相對濕度的距平值多為負值,低于多年平均值;1990—2000年年均相對濕度較高,距平值也大多高于多年平均值;進入2000年代后前5 a年均相對濕度呈現高值,后5 a情況則恰好相反,距平值為負值。

1980—2005年雅江流域年均蒸發量呈顯著降低趨勢(P<0.01),平均氣候傾向率為-235.31 mm·(10 a)-1,符合青藏高原總體年均蒸發量的變化趨勢,但低于藏北高原的年均蒸發量氣候傾向率〔-617 mm·(10 a)-1〕[17-18]。下游地區年均蒸發量多年平均值(1 455 mm)明顯低于上游(2 138 mm)和中游地區(2 177 mm)。1980—1995年年均蒸發量距平值幾乎為正值,均高于多年平均值;而1995—2000年年均蒸發量呈波動變化,正距平逐漸減少,負距平開始增多;2000年之后年均蒸發量距平值均低于多年平均值,表明年均蒸發量自2000年起明顯降低。

圖3　1980—2010年雅魯藏布江流域主要氣象因子變化

3.3雅江流域河流徑流量和雪深的變化特征

3個測站天然徑流的5 a滑動趨勢線見圖4。

圖4　1980—2010年雅魯藏布江干流

由圖4可見，奴各沙、羊村和奴下3 個水文站的歷年徑流變化具有相同的變化周期,即整個雅江干流區在1980年代初期徑流持續下降,1983年達最低值,1984年開始波動上升,1990年代進入相對豐水期,2000年達最高值,隨后呈波動性降低趨勢。從總體上看,徑流年際變化存在一定的周期性波動,2個豐水時間段為1980年代后期至1990年代初期和2000年以后。2000年以后上、中、下游各站徑流量總體呈快速增長趨勢。1998—2000年3個測站年徑流量均劇增,這可能是受1998年大洪水影響所致。

通過分析雅江流域的年平均雪深變化(圖5～6)發現,1980年代年平均雪深的氣候傾向率為-0.10 cm·(10 a)-1,1980—1984年雅江流域年平均雪深呈顯著下降趨勢(P<0.05),1985—1988年年平均雪深呈波動性下降趨勢,距平普遍為負值。1989—1998年年平均雪深呈顯著上升趨勢(P<0.05),距平值在1991年之后也升高為正值,10 a間年平均雪深的氣候傾向率為0.18 cm·(10 a)-1,與整個青藏高原1980—1995年年平均雪深呈顯著線性上升趨勢〔約為0.06 cm·(10 a)-1〕的研究結論一致[19-20],但高于青藏高原的上升幅度。1998年雅江流域平均雪深達到峰值0.25 cm,與降水峰值年份相同。1998—2007年積雪深度呈波動下降趨勢,10 a間年平均雪深的氣候傾向率為-0.05 cm·(10 a)-1。

圖5　1980—2007年雅魯藏布江流域年平均雪深變化

圖6　1980—2007年雅魯藏布江流域年平均雪深距平序列

4雅江流域高寒濕地對氣候變化的響應

4.1河流濕地

雅江流域上、中、下游高寒濕地面積與自然氣候因子的灰色關聯度見表4。由表4可知,年均溫和年均最高溫與河流面積的關聯性均達顯著水平(P<0.05)。近30 a雅江流域氣溫總體呈升高趨勢,而河流濕地面積總體上也呈增加趨勢。由于雅江流域高寒的氣候特征,現代冰川多發育于此,冰雪融水也成為河流、湖泊等濕地水體的重要補給水源。氣候升高會引起冰川萎縮,冰川融水產生的冰川徑流直接匯入其下段的河流、湖泊等濕地水體,對這些水體的徑流變化產生影響[20-21]。

1980年代雅江流域年平均雪深呈顯著下降趨勢,而同期冰川也呈現退縮狀態[22],這說明冰川末端融化引起冰雪融水增多,雪水流入附近的河流后使得河流徑流量增加。前文研究結果顯示,從1984年開始河流徑流量呈上升趨勢,1980—1990年河流面積增加172.93 km2。 而進入2000年后,2000—2010年雅江流域多年平均氣溫為7.69 ℃,年均溫距平幾乎普遍呈正值(圖3),為近30 a來最溫暖時期。冰川退縮進一步加劇,年平均雪深也再次呈現緩慢降低趨勢,使得融雪量增加,冰雪融水變多。以融水補給為主的尼洋河、易貢藏布等支流以及干流中游河段徑流量增加,上、中、下游各站徑流量在近10 a間快速增長,與1980年代相比分別增加28.4%、34.9%和21.7%(圖4)。因此,冰雪融水和降水量的增加使得河流面積在2000—2010年出現快速增長,面積動態變化度為1.68%(表2),而同期年均溫的氣候傾向率也高達1.20 ℃·(10 a)-1,說明氣溫升高最迅速的10 a也是河流面積增加最為劇烈的時期。

表4雅魯藏布江流域上、中、下游高寒濕地面積與自然氣候因子的灰色關聯度

Table 4Gray relational analysis of the areas of alpine wetlands in the upper, middle and lower reaches of the Yarlung Zangbo River Valley with natural meterological factors

流域濕地類型關聯度年均溫年均最高溫年均最低溫年均相對濕度年降水量年蒸發量上游河流0.8684*0.9346*0.61540.8262*0.68690.8018湖泊0.8723*0.8835*0.67570.82630.72230.8040*沼澤0.8281*0.9103*0.61660.8233*0.66900.8246人工濕地0.8119*0.8771*0.65080.7931*0.6639*0.8114*中游河流0.8652*0.8977*0.56490.81920.69640.7417湖泊0.8284*0.8990*0.56180.8512*0.71730.7993沼澤0.8310*0.8904*0.60390.8666*0.75330.8418*人工濕地0.7957*0.7428*0.62110.71600.68440.6239下游河流0.7042*0.6854*0.7476*0.68530.63520.6225湖泊0.7693*0.8292*0.62780.72150.65350.7079沼澤0.8023*0.8606*0.65860.74060.57330.6824人工濕地0.67080.66010.70420.65700.71030.6611

*表示P<0.05。

4.2湖泊濕地

由表4可知,年均溫和年均最高溫這2個因素與雅江上、中、下游地區湖泊面積均呈顯著關聯關系(P<0.05),說明氣溫是與湖泊面積動態變化的最主要的關聯因子。雅江流域內較大的湖泊一般均發育在河谷的開闊地段,而小型湖泊則分布在各支流的源頭附近,是冰川在長期進退消融過程中對地表局部挖蝕、堆積的結果,因此冰川作用在很大程度上影響著湖泊水情[7,20]。近30 a雅江流域年平均雪深呈現先降低后升高又降低的波動變化趨勢,由此引起的融雪雪水量也呈現先增加后降低又增加的變化趨勢,這與湖泊面積的動態變化規律一致(表2),而發生冰川退縮的區域附近也常常伴隨出現小型湖泊(面積<10 km2)的擴張。雅江上游的湖泊面積變化除與氣溫顯著關聯之外,還與年均蒸發量顯著關聯(P<0.05)。這是因為對于上游地區封閉的內陸湖泊[7]來說,其水源的主要輸出途徑是蒸發,而降水量的增加和蒸發量的減少成為湖泊面積在近30 a總體呈持續增長的主要原因。另一方面,通過轉移矩陣計算,筆者發現1980年以來有91.78%～97.86%的沼澤濕地轉變為小型湖泊,這可能是因為冰川退縮的加快引起了以冰川融水補給為主的湖泊水位上漲,湖泊面積增大,導致周邊沼澤濕地被淹沒,使其陸續轉變成小型湖泊。

4.3沼澤濕地

對沼澤濕地而言,年均溫、年均最高氣溫和年均相對濕度這3個因子與其關聯度大多達顯著水平(P<0.05)。由于氣溫升高導致冰雪融水增加,沼澤濕地轉變為湖泊濕地;此外,不斷升高的氣溫使得與沼澤濕地息息相關的多年凍土層不斷消融[23],從而加速了沼澤濕地的退化。另一方面,高原濕地消耗水分的主要途徑是蒸散發,而與蒸散發直接相關的一個重要因子是相對濕度[24]。近30 a間雅江流域上、中、下游年平均相對濕度有不同程度下降(圖3),說明該流域氣候呈干熱化發展趨勢。相對濕度的降低意味著植物生長季內空氣濕度飽和差增大,這必然會促使生長季內蒸散發量增加,從而為沼澤濕地水分的不斷喪失創造條件;同時較低的相對濕度也促使植被的蒸騰作用加劇,導致土壤墑情惡化[25]。上述原因都對植被生長構成危害,這也可能是高原沼澤草甸植被發生退化、沼澤萎縮的原因之一。

4.4人工濕地

從整個雅江流域看,人工濕地變化主要受人類活動干擾,與特定氣象因子的關聯性大多不顯著。但上中游的人工濕地面積變化與氣溫顯著關聯(表4),這可能是由于隨著氣溫升高上中游的稻田人工濕地種植面積擴大。近30 a間雅江流域人工濕地面積增幅達164%,且人工濕地大幅增加主要發生在1990—2000年。1991年開始,國家對雅江中游及其支流拉薩河、年楚河實施“一江兩河”開發,重點扶持種植業、農業和公路、大壩、水庫等的建設[13]。該項目的實施是1990年代以來雅江流域人工濕地面積增加的主要驅動力。

5結論

通過野外調查和遙感目視解譯相結合的方法分析了1980—2010年雅江流域濕地的變化特征,同時利用同期氣候、徑流和雪深數據討論了其與雅江濕地面積變化的關聯關系?？傮w上來看,沼澤、河流和湖泊濕地構成了雅江流域濕地的主體,占雅江濕地總面積的99.77%,而水稻田和水庫等人工濕地分布較少,僅占濕地總面積的0.23%。近30 a來雅江流域年均溫、年均最高溫和年均最低溫總體呈顯著增加趨勢,年均蒸發量呈顯著下降趨勢,降水量波動較大,年均相對濕度呈不顯著下降趨勢。雅江流域高寒濕地的動態變化與氣候變化密切相關,尤其是年均溫、年均最高溫和年均相對濕度。近30 a來河流濕地面積、雅江5大支流的干流河長增加主要是受年均溫和年均最高溫升高影響所致;湖泊面積總體呈增加趨勢,主要與氣溫升高有關;沼澤濕地面積持續降低,這與區域氣溫升高以及年均相對濕度降低顯著關聯;水庫、水稻田等人工濕地面積增加則主要是受人類活動影響。

由于該研究首次嘗試摸清雅江全流域的濕地動態變化,對于較大范圍內衛星遙感影像的使用存在一定的局限性,各類衛星資料的獲取時間不在同一個季相,使得獲取的濕地變化情況可能存在一定誤差。另外,由于缺乏雅江流域更多站點的水文觀測數據和湖泊周圍冰雪融水的變化情況資料,對于冰川變化引起的湖泊、河流變化以及它們之間的相互作用等問題尚需進一步的深入研究。

參考文獻：

[1]SHI Y,GAO X J,ZHANG D F,etal.Climate Change Over the Yarlung Zangbo-Brahmaputra River Basin in the 21st Century as Simulated by a High Resolution Regional Climate Model[J].Quaternary International,2011,244(2):159-168.

[2]孫鴻烈,鄭度.青藏高原形成演化和發展[M].廣州:廣東科學技術出版社,1998:155-169.

[3]關志華,陳傳友,區裕雄,等.西藏河流與湖泊[M].北京:科學出版社,1984:176-182.

[4]白軍紅,歐陽華,崔保山.近40年來若爾蓋高原高寒濕地景觀格局變化[J].生態學報,2008,28(5):2245-2252.

[5]BAI J H,LU Q Q.WANG J J,etal.Landscape Pattern Evolution Processes of Alpine Wetlands and Their Driving Factors in the Zoige Plateau of China[J].Journal of Mountain Science,2013,10(1):54-67.

[6]ZHANG Y L,WANG C L,BAI W Q.Distribution of Alpine Wetland in Lhasa River Basin,China[J].Journal of Geographical Sciences,2010,20(3):375-388.

[7]SARKAR A,GRAG R D,SHARMA N Y.RS-GIS Based Assessment of River Dynamics of Brahmaputra River in India[J].Journal of Water Resource and Protection,2012,4(2):63-72.

[8]林乃峰,沈渭壽,張慧,等.近35 a西藏那曲地區湖泊動態遙感與氣候因素關聯度分析[J].生態與農村環境學報,2012,38(3):231-237.

[9]GUO Y P,WANG C H.Trends in Precipitation Recycling Over the Qinghai-Xizang Plateau in Last Decades[J].Journal of Hydrology,2014,517:826-835.

[10]SHEN W S,LI H D,SUN M,etal.Dynamics of Aeolian Sandy Land in the Yarlung Zangbo River Basin of Tibet,China From 1975 to 2008[J].Global and Planetary Change,2012,86/87:37-44.

[11]湯奇成,熊怡.中國河流水文[M].北京:科學出版社,1998:15-19.

[12]舒儉民,王家驥,高吉喜,等.雅魯藏布江中部地區生物多樣性保護規劃[J].生態與農村環境學報,1997,13(2):1-5.

[13]LI H D,SHEN W S,ZOU C X,etal.Spatio-Temporal Variability of Soil Moisture and Its Effect on Vegetation in a Desertified Aeolian Riparian Ecotone on the Tibetan Plateau,China[J].Journal of Hydrology,2013,479:215-225.

[14]吳紹洪,尹云鶴,鄭度,等.青藏高原近30年氣候變化趨勢[J].地理學報,2005,60(1):3-11.

[15]LI X J,XU L Y.An Evaluation of Ecological Losses From Hydropower Development in Tibet[J].Ecological Engineering,2015,76:178-185.

[16]CHEN H.Assessment of Hydrological Alterations From 1961 to 2000 in the Yarlung Zangbo River,Tibet[J].Ecohydrology & Hydrobiology,2012,12(2):93-103.

[17]LIU X M,ZHENG H X,ZHANG M H,etal.Identification of Dominant Climate Factor for Pan Evaporation Trend in the Tibetan Plateau[J].Journal of Geographical Sciences,2011,21(4):594-608.

[18]杜軍,邊多,鮑建華.藏北高原蒸發皿蒸發量及其影響因素的變化特征[J].水科學進展,2008,19(6):786-791.

[19]馬麗娟,秦大河.1957—2009年中國臺站觀測的關鍵積雪參數時空變化特征[J].冰川凍土,2012,34(1):1-11.

[20]LIN Z J,NIU F J,LIU H,etal.Hydrothermal Processes of Alpine Tundra Lakes,Beiluhe Basin,Qinghai-Tibet Plateau[J].Cold Regions Science and Technology,2011,65(3):446-455.

[21]姚檀棟,李治國,楊威,等.雅魯藏布江流域冰川分布和物質平衡特征及其對湖泊的影響[J].科學通報,2010,55(20):2072-2078.

[22]YAO T D,THOMPSON L,YANG W,etal.Different Glacier Status With Atmospheric Circulations in Tibetan Plateau and Surroundings[J].Nature Climate Change,2012,2:663-667.

[23]YANG Z P,OUYANG H,XU X L,etal.Effects of Permafrost Degradation on Ecosystems[J].Acta Ecologica Sinica,2010,30(1):33-39.

[24]牛沂芳,李才興,習曉環,等.衛星遙感檢測高原湖泊水面變化及與氣候變化分析[J].干旱區地理,2008,31(2):284-290.

[25]WANG G X,WANG Y B,LI Y S,etal.Influences of Alpine Ecosystem Responses to Climatic Change on Soil Properties on the Qinghai-Tibet Plateau,China[J].Catena,2007,70(3):506-514.

(責任編輯： 許素)

Dynamic of the Alpine Wetlands and Its Response to Climate Change in the Yarlung Zangbo River Valley in Recent 30 years.

LIUDong1,WANGTao2,SHENWei-shou1,LINNai-feng1,ZOUChang-xin1

(1.Nanjing Institute of Environmental Sciences, Ministry of Environmental Protection, Nanjing 210042, China；2.School of Geography and Remote Sensing, Nanjing Universitos of Information Science & Technology, Nanjing 210044, China)

Abstract:The knowledge about dynamic of the alpine wetlands in the Yarlung Zangbo River Valley and its response to climate change would help understand in depth variation of the alpine ecosystem of the Qinghai Tibet Plateau and its response to global climate change. Information dynamics of the wetlands, in area, in the Yarlung Zangbo River Valley was extracted from the remote sensing images and data transmitted from the Landsat and Environmental Satellite #1 in 1980, 1990, 2000 and 2010, and collated with the data obtained in field surveys and the recent 30 years of metereological, hydrological and snow-accumulation data to explore dynamics of wetlands and its response to climate change. Results show that in 2010 alpine wetlands accounted for 3.73% of the whole valley in area, with swampy wetlands being the largest in area and accounting for 59.46% of the total area of the wetlands, followed by riverine wetlands (30.19%), lacustrine wetlands (10.13%) and man-made wetlands (0.23%). In the recent 30 years riverine and lacustrine wetlands have been on a rising trend in area, while swampy wetlands on a reverse trend. The Yarlung Zangbo River Valley has been on a rising trend in annual mean temperature, mean annual maximum temperature, and mean annual minimum temperature, fluctuating greatly or on a slightly, but not apparently, rising trend in precipitation, and a declining trend, though not apparent, in mean relatively humidity. The changes in annual mean temperature and mean annual maximum temperature were significantly related to the changes in the area of the lakes in the valley; the increases in the area of the riverine wetlands and in the length of the main streams of the five major tributaries of the Yarlung Zangbo River were to the changes in annual mean temperature, mean annual maximum temperature, and mean annual minimum temperature; and the area of the swampy wetlands was to the air temperature and annual mean relative humidity. The country′s strategy of construction, like building dams and reservoirs is the major driving force of the increase in the area of man-made wetlands in the river valley.

Key words:alpine wetland;Yarlung Zangbo River;climate；response

作者簡介：劉冬(1984—),女,河南洛陽人,助理研究員,博士,主要從事氣候變化與生態安全方面的研究。E-mail: liudong@nies.org

DOI：10.11934/j.issn.1673-4831.2016.02.012

中圖分類號：X22；X37

文獻標志碼：A

文章編號：1673-4831(2016)02-0243-09

通信作者①E-mail： zcx@nies.org

基金項目：國家“十二五”科技支撐計劃(2012BAC19B06-1);國家自然科學基金青年基金(41501596)；江蘇省普通高校研究生科研創新計劃(KYLX15_0869)

收稿日期：2015-03-26