王永鵬, 楊鵬年, 周 龍, 徐生武, 鄧曉雅, 馮思陽(yáng)
(1.新疆農(nóng)業(yè)大學(xué) 水利與土木工程學(xué)院, 新疆 烏魯木齊 830052; 2.新疆水利工程安全與水災(zāi)害防治重點(diǎn)實(shí)驗(yàn)室,新疆 烏魯木齊 830052; 3.塔里木河流域干流管理局, 新疆 庫(kù)爾勒 841000; 4.中國(guó)水利水電科學(xué)研究院, 北京 100044)
塔里木河(以下簡(jiǎn)稱“塔河”)下游屬極度干旱區(qū),蒸散是該區(qū)域地下水的主要排泄方式,準(zhǔn)確計(jì)算蒸散量是把握塔河下游生態(tài)需水量與植被修復(fù)面積的關(guān)鍵。解析歷次輸水后植被耗水時(shí)空分布及其變化特征,對(duì)認(rèn)清下游的生態(tài)水文過(guò)程與開(kāi)展生態(tài)調(diào)度具有重要的作用。
蒸散是植被耗水的主要方式,近年來(lái),有大量研究利用蒸散估算植被耗水量[1]。傳統(tǒng)的蒸散量計(jì)算方法局限于點(diǎn)值,在向流域尺度擴(kuò)展時(shí)其精度明顯下降;遙感技術(shù)在估算區(qū)域蒸散發(fā)領(lǐng)域經(jīng)過(guò)多年的探索研究取得了長(zhǎng)足的發(fā)展[2]。Lu等[3]基于植被修復(fù)背景下對(duì)中國(guó)干旱半干旱地區(qū)植被耗水進(jìn)行了研究,發(fā)現(xiàn)干旱地區(qū)植被耗水量遠(yuǎn)超過(guò)區(qū)域降水量,這是導(dǎo)致地下水位大范圍持續(xù)下降的原因。李寶富等[4]基于地表能量平衡原理,通過(guò)建立SEBAL模型結(jié)合土地利用數(shù)據(jù),確定了塔河干流區(qū)蒸散量,得出植被耗水與土地利用類(lèi)型、植被覆蓋度等密切相關(guān)的重要結(jié)論。李霞等[5]通過(guò)Landsat系列數(shù)據(jù),提取歸一化植被指數(shù)(NDVI)、土壤調(diào)整植被指數(shù)(SAVI)和簡(jiǎn)單比值植被指數(shù)(SR),結(jié)合地面實(shí)測(cè)土壤含水量和葉面積指數(shù),建立了科爾沁地區(qū)植被系數(shù)(Kc)模型,驗(yàn)證了草甸和沙丘試驗(yàn)區(qū)Kc模型的可行性。目前,塔河下游植被耗水量的研究以水量平衡法[6-7]、潛水蒸發(fā)法[8]和面積定額法等[9]為主,其均以有限監(jiān)測(cè)地下水埋深數(shù)據(jù)、單一植被需水定額或潛水蒸發(fā)系數(shù)為基礎(chǔ),計(jì)算結(jié)果空間分辨率不高,且未能揭示植被耗水過(guò)程和機(jī)理,尚不足以作為確定塔河下游生態(tài)需水量的依據(jù)。結(jié)合干旱地區(qū)天然植被覆蓋度、水分脅迫等特征以及涉及耗水過(guò)程的植被耗水量研究在塔河下游有待探討。
遙感數(shù)據(jù)具有高效、量大、經(jīng)濟(jì)等特點(diǎn),且隨著地理信息系統(tǒng)等技術(shù)的發(fā)展,實(shí)現(xiàn)植被耗水周期性監(jiān)測(cè),客觀描述水文循環(huán)過(guò)程成為可能[10]。Penman-Monteith公式法(以下簡(jiǎn)稱P-M法)應(yīng)用能量平衡方程結(jié)合空氣動(dòng)力學(xué)原理,充分考慮了各氣象參數(shù),物理意義明確,是國(guó)際糧農(nóng)組織(FAO)推薦的定義和計(jì)算參照蒸散發(fā)的唯一標(biāo)準(zhǔn)方法,被廣泛應(yīng)用于植被蒸散領(lǐng)域[11];其在裸土、植被全覆蓋[12]及草甸、沙丘等[5]不同下墊面情況均有較好的應(yīng)用效果,且適用于任何時(shí)間尺度的蒸散量研究,能合理地揭示植被耗水過(guò)程和機(jī)理。
基于此,本文以Landsat系列遙感影像、氣象數(shù)據(jù)及土壤含水率等數(shù)據(jù)為基礎(chǔ),采用改進(jìn)的P-M法,結(jié)合GIS技術(shù),將塔河下游天然植被覆蓋度與P-M法中參考作物的差異以植被系數(shù)Kc在柵格尺度上進(jìn)行修正,將土壤含水率的差異以土壤水分限制系數(shù)Ks在不同土地利用類(lèi)型上進(jìn)行修正,在30 m柵格尺度上計(jì)算塔河下游天然植被耗水量。分析其時(shí)空分布及變化趨勢(shì),以探討生態(tài)輸水效應(yīng),為今后塔河下游制定生態(tài)輸水策略及修復(fù)目標(biāo)提供依據(jù)。
研究區(qū)地處塔河下游大西海子水庫(kù)—阿拉干斷面的其文闊爾河和老塔河(以下簡(jiǎn)稱“兩河”)之間狹長(zhǎng)的沖積平原區(qū)域(見(jiàn)圖1),地處于新疆維吾爾自治區(qū)南部尉犁縣和若羌縣境內(nèi),位于東經(jīng)87°00′—88°40′,北緯40°00′—40°50′之間。年降水量?jī)H有17.4~42.0 mm,多年平均蒸發(fā)量高達(dá)2 500~3 000 mm,屬于典型的大陸性干旱氣候[13]。塔河下游水源補(bǔ)給方式主要通過(guò)大西海子水庫(kù)下泄生態(tài)水,淺層土壤含水率受潛水埋深影響大,天然植被生長(zhǎng)受水分脅迫嚴(yán)重。天然植被主要由依靠潛水生長(zhǎng)的非地帶性植被構(gòu)成,其物種單一,覆蓋度低,以廊道式分布在兩河沿線,主要有胡楊、檉柳、黑果枸杞、鈴鐺刺、蘆葦、花花柴、脹果甘草等,構(gòu)成喬木、灌木和草本植物群落[8]。2000—2020年,塔河下游已實(shí)施了21次生態(tài)輸水,河岸帶地下水位普遍抬升了3~5 m,垂直主河道地下水顯著抬升范圍超過(guò)1 000 m[14]。以其文庫(kù)勒湖為代表的多個(gè)自然漫溢區(qū)蘆葦?shù)炔荼局参镩L(zhǎng)勢(shì)良好,且出現(xiàn)胡楊及檉柳等幼苗萌發(fā)現(xiàn)象,塔河下游退化的生態(tài)格局出現(xiàn)了逆轉(zhuǎn)。
圖1 塔里木河下游土地利用類(lèi)型及監(jiān)測(cè)井位置
2.1.1 氣象數(shù)據(jù) 氣象資料從國(guó)家氣象數(shù)據(jù)中心下載(http:∥data.cma.cn)。使用距研究區(qū)最近的鐵干里克氣象站日尺度數(shù)據(jù),包括逐日降雨量、平均氣溫、平均風(fēng)速、日照時(shí)數(shù)和平均相對(duì)濕度等。對(duì)于部分未監(jiān)測(cè)到的數(shù)據(jù),采用鄰近時(shí)段數(shù)據(jù)內(nèi)插補(bǔ)齊。
2.1.2 遙感數(shù)據(jù) 遙感數(shù)據(jù)來(lái)源于美國(guó)地質(zhì)調(diào)查局(https:∥earthexplorer.usgs.gov/)。篩選2000,2005,2010,2015,2020年7—11月云量覆蓋度低的Landsat系列遙感影像,利用ENVI5.3等軟件對(duì)獲取的遙感影像進(jìn)行影像裁剪、輻射定標(biāo)、大氣校正等預(yù)處理。
2.1.3 生態(tài)輸水量和地下水位數(shù)據(jù) 輸水量資料和地下水位數(shù)據(jù)均由塔里木河流域管理局提供,包括2000—2020年生態(tài)輸水期間,大西海子水庫(kù)及以下英蘇、喀爾達(dá)依—博孜庫(kù)勒、阿拉干3個(gè)典型水文斷面的地表徑流量和地下水位數(shù)據(jù)。
2.2.1 植被耗水量估算方法 本文采用改進(jìn)Penman-Monteith法[15],計(jì)算塔河下游非完全覆蓋、受水分脅迫下的植被耗水量。計(jì)算公式如下:
ETc=KsKcET0
(1)
式中:ETc為植被耗水強(qiáng)度(mm/d);Ks為土壤水分限制系數(shù);Kc為植被系數(shù); ET0為參考蒸散強(qiáng)度(mm/d)。
考慮到植被耗水強(qiáng)度的差異,基于各類(lèi)土地利用面積和典型植被的耗水強(qiáng)度,通過(guò)累加的方式確定植被耗水總量。計(jì)算公式如下:
(2)
式中:VWC為植被耗水總量(m3);An為各類(lèi)土地利用面積(m2); ETcn為第n天的典型植被蒸散強(qiáng)度(mm/d);n為典型植被生長(zhǎng)期日序數(shù)。
2.2.2參考蒸散發(fā)強(qiáng)度(ET0) 參考作物蒸散發(fā)強(qiáng)度ET0根據(jù) FAO Penman-Monteith公式[15]計(jì)算得到,通過(guò)站點(diǎn)氣象數(shù)據(jù)計(jì)算日尺度ET0,通過(guò)累加的方式確定植被各生長(zhǎng)期ET0。
2.2.3 土地利用類(lèi)型劃分及植被系數(shù)(Kc) 參考中國(guó)植被圖編輯委員會(huì)2001年編纂完成的中國(guó)植被區(qū)劃圖及白元[9]、劉新華等[8]對(duì)塔河下游的研究,基于遙感影像,采用ENVI軟件中監(jiān)督分類(lèi)的方法確定塔河下游土地利用類(lèi)型。選取各類(lèi)土地利用的典型植被,根據(jù)FAO-56推薦的方法[15],結(jié)合中國(guó)科學(xué)院發(fā)布的《中國(guó)動(dòng)植物物候觀測(cè)年報(bào)》第8號(hào)(1982年)[16]中相關(guān)區(qū)域的植被物候信息,將植物生長(zhǎng)周期劃分為生長(zhǎng)初期、發(fā)育期、生長(zhǎng)中期和生長(zhǎng)后期4個(gè)階段(表1)。
表1 典型植被生長(zhǎng)階段劃分
植被系數(shù)Kc是表征植被耗水能力的重要參數(shù),Kc曲線反應(yīng)了植被全生育周期的耗水過(guò)程。生長(zhǎng)初期植被系數(shù)(Kcnin)主要由土壤自身持水能力、濕潤(rùn)時(shí)間間隔、大氣蒸發(fā)能力等因素決定;生長(zhǎng)中期、末期植被系數(shù)(Kcmid,Kcend)主要由植被覆蓋度、有效覆蓋度等決定;發(fā)育期和生長(zhǎng)晚期分別通過(guò)Kcini與Kcmid及Kcmid與Kcend線性內(nèi)插獲取[15]。生長(zhǎng)中期植被覆蓋度通過(guò)植被長(zhǎng)勢(shì)旺盛的7—9月的Landsat系列影像解譯獲取,生長(zhǎng)末期通過(guò)11月的影像獲取。
2.2.4 土壤水分限制系數(shù)(Ks) 土壤水分限制系數(shù)Ks是表征植被生長(zhǎng)水分盈缺的關(guān)鍵指標(biāo),以植被根系層土壤含水量等確定。Ks參考FAO-56推薦的非完全灌溉下、受水分脅迫的天然植被耗水量確定方法[15]。其中所需的植被根系深度參考張麗等[17]的研究確定,田間持水率θFC采用Saxton[18]的方法確定,凋萎系數(shù)θWP參考喬照華[19]的方法,土壤黏粒含量等參數(shù)通過(guò)土壤顆粒分析結(jié)果確定,不同根深的土壤含水率θ通過(guò)現(xiàn)場(chǎng)烘干試驗(yàn)測(cè)定。
土壤顆粒分析結(jié)果、含水率及各類(lèi)土地利用的Ks計(jì)算結(jié)果分別見(jiàn)表2—4。
表2 土樣顆粒分析結(jié)果相關(guān)參數(shù)
表3 土壤體積含水率(θ)隨埋深分布
表4 不同土地利用類(lèi)型的Ks值
研究區(qū)土地利用類(lèi)型分布及面積變化趨勢(shì)(圖2)發(fā)現(xiàn),間隔5 a的天然植被面積呈明顯增加趨勢(shì),河岸帶荒漠面積大量向天然植被轉(zhuǎn)變,但兩河區(qū)域土地利用類(lèi)型整體仍以荒漠為主。2000—2020年平均荒漠面積為1 406.03 km2,占兩河區(qū)域面積比例為66.11%;林草地主要分布在主河道沿線及漫溢區(qū),林地面積為124.03,比例為5.83%;草地面積為380.16 km2,比例為17.87%。水體受既有的間歇性輸水方式影響,以大西海子水庫(kù)庫(kù)區(qū)為主,平均水體面積為79.10 km2,比例為3.72%。而2020年水庫(kù)庫(kù)區(qū)解譯結(jié)果為草地,其原因是2020年對(duì)大西海子水庫(kù)部分區(qū)域進(jìn)行了維護(hù)修繕工作,庫(kù)區(qū)無(wú)水,水草生長(zhǎng)旺盛;耕地則集中分布于第二師34團(tuán),面積為137.53 km2,比例6.47%。21 a來(lái),荒漠與天然植被面積呈顯著變化(圖3),其中,荒漠面積多年以48.14 km2/5 a(占多年平均荒漠面積的3.42%)的速率遞減,至2020年,其面積比例降至61.65%;天然植被面積以62.35 km2/5 a(占多年平均天然植被面積的12.37%)的速率增長(zhǎng),比例上升至31.39%。
圖2 2000—2020年塔里木河下游土地利用類(lèi)型分布
圖3 2000—2020年塔里木河下游
不同時(shí)期各類(lèi)土地利用的植被覆蓋度值(表5)顯示,2000—2020年林地平均覆蓋度為0.324,草地0.175,荒漠僅有0.037,天然植被自然生長(zhǎng)下覆蓋度整體偏低,長(zhǎng)勢(shì)較差。21 a來(lái)塔河下游天然植被覆蓋度整體呈上升趨勢(shì),其中,林草地增長(zhǎng)趨勢(shì)最為顯著,增長(zhǎng)速率分別為0.042/5 a,0.027/5 a;而荒漠區(qū)域增長(zhǎng)相對(duì)較少,增長(zhǎng)速率僅為0.002/5 a。
表5 塔里木河下游植被覆蓋度Fc值
根據(jù)塔河下游植被覆蓋度、高度和氣象條件等實(shí)際情況,對(duì)Penman-Monteith法中標(biāo)準(zhǔn)條件加以修正,確定了各生長(zhǎng)階段的植被系數(shù)Kc。Kc曲線(圖4)顯示,天然植被生長(zhǎng)中期時(shí)段長(zhǎng)、耗水強(qiáng)度高,其耗水主要集中于中期。2000—2020年,塔河下游耕地Kc多年平均為1.08,天然植被則受限于其覆蓋度,整體較小,多年平均林地為0.85,草地為0.44,荒漠植被為0.21。不同土地利用類(lèi)型耗水強(qiáng)度(表6)差異明顯,生長(zhǎng)中期多年平均耗水強(qiáng)度耕地為3.88 mm/d,林地為3.43 mm/d,草地為1.25 mm/d,荒漠植被為0.45 mm/d;整體呈現(xiàn)耕地耗水強(qiáng)度最高,其次為林地,且在多年份保持一致。這與正常森林與農(nóng)作物耗水強(qiáng)度的差異相悖,表明塔河下游天然植被自然生長(zhǎng)下長(zhǎng)期受限于水分脅迫。就年際間植被同一生長(zhǎng)階段的耗水強(qiáng)度變化而言(圖5),生長(zhǎng)初期保持平穩(wěn)狀態(tài),發(fā)育期、生長(zhǎng)中期及晚期耗水強(qiáng)度波動(dòng)較大,整體呈上升趨勢(shì)。
圖4 塔里木河下游植被系數(shù)Kc曲線
圖5 2000—2020年塔里木河下游天然植被
塔河下游植被耗水量空間分布(圖6)顯示,植被耗水量空間分布差異較大,2015以前植被耗水量高值區(qū)主要集中在連片耕地,2015年及以后主河道沿線及局部漫溢區(qū)出現(xiàn)高值;表明經(jīng)過(guò)21 a的生態(tài)輸水,河岸帶天然植被修復(fù)狀況良好。植被耗水量變化程度顯示:植被耗水增長(zhǎng)量空間差異較懸殊,植被耗水量明顯增長(zhǎng)區(qū)域主要集中在主河道沿線及自然漫溢區(qū),其文庫(kù)勒湖、博孜庫(kù)勒、庫(kù)木吐格3處自然漫溢區(qū)植被耗水量上升最為顯著,有向濕地轉(zhuǎn)變傾向,荒漠區(qū)耗水量變幅則較小。天然植被耗水量(表7—8)顯示:多年平均天然植被耗水量為1.90×108m3,其中林地6.48×107m3,草地4.63×107m3,荒漠植被為7.90×107m3。多年來(lái)天然植被耗水量以2.44×107m3/5 a(占多年平均天然植被耗水量的12.82%)的速率增長(zhǎng),其耗水量增大是面積和耗水強(qiáng)度雙增長(zhǎng)的結(jié)果;荒漠植被耗水量2000—2005年上升,2005年以后下降,其變化趨勢(shì)是荒漠面積減小和耗水強(qiáng)度增加共同作用的結(jié)果。
圖6 2000—2020年塔里木河下游植被耗水量空間分布及變化程度
表7 2000—2020年塔里木河下游各類(lèi)土地年耗水量 108 m3
研究區(qū)除大西海子水庫(kù)下泄生態(tài)水外無(wú)其他補(bǔ)給來(lái)源,大西海子水庫(kù)下泄斷面與阿拉干斷面之間的水量差為河道消耗量,其主要轉(zhuǎn)化項(xiàng)為地下水補(bǔ)給量、植被蒸散量及水面蒸發(fā)量。
由于河道水面面積較小且為間歇性過(guò)水,故河道消耗量主要以前二項(xiàng)為主。以5 a為步長(zhǎng),將2000—2020年分為4個(gè)輸水時(shí)段,天然植被耗水量與河道總消耗量對(duì)比(始末年植被耗水量平均值×5)(表8)顯示,不同輸水時(shí)段植被耗水量占河道總消耗量比例相差較大,其原因在于枯水年下泄量小,地下水被植被大量消耗,表現(xiàn)為植被耗水量占河道總消耗量的比例偏高;豐水年生態(tài)輸水優(yōu)先填補(bǔ)枯水年導(dǎo)致的地下水虧空,導(dǎo)致結(jié)果偏低,這也體現(xiàn)了地下水的調(diào)蓄功能。采取水量累加的方式可避免植被耗水對(duì)生態(tài)輸水的滯后效應(yīng),至2020年,植被耗水累積量為3.99×109m3,河道消耗累積量為5.85×109m3,植被耗水占河道消耗的68.23%,已為河道消耗的主要轉(zhuǎn)化項(xiàng)。
表8 2000—2020年塔里木河下游天然植被耗水量與河道消耗量對(duì)比
本文通過(guò)GIS技術(shù)將改進(jìn)的P-M法應(yīng)用于30 m空間柵格尺度,實(shí)現(xiàn)植被耗水量逐像元計(jì)算,確定出的植被耗水量轉(zhuǎn)化比例為68.23%。與鄧銘江等[20]采用水均衡原理確定的塔河下游植被耗水轉(zhuǎn)換比例51.55%;劉新華等[9]基于2005年中巴資源衛(wèi)星數(shù)據(jù),采用潛水蒸發(fā)法和面積定額計(jì)算的轉(zhuǎn)化比例44.14%有一定差異。其客觀原因表現(xiàn)在大西海子水庫(kù)—阿拉干段天然植被密度高于阿拉干—臺(tái)特瑪湖段;且隨著生態(tài)輸水次數(shù)的增加,下游植被得到改善,植被耗水量增大,使得后期輸水植被耗水轉(zhuǎn)化比例上升。本文較前述二者的研究在數(shù)據(jù)源的選取與計(jì)算方法上均有了較大的改進(jìn),所得柵格尺度的植被耗水量應(yīng)更符合塔河下游植被耗水現(xiàn)狀,結(jié)果具有較高的可信度。
隨著輸水次數(shù)的增多,塔河下游天然植被耗水量的增加愈發(fā)顯著,尤以植被覆蓋度增加的貢獻(xiàn)為主。目前,天然植被耗水量多年平均為1.90×108m3/a,仍小于兩河之間平均河道消耗量2.79×108m3/a,林草地生態(tài)用水仍有較大的提升空間。下游草地多以1 a生為主,其萌新對(duì)外界環(huán)境變化較敏感,對(duì)生態(tài)輸水響應(yīng)快,至2020年面積占比已上升至24.67%,已得到較好的恢復(fù);林地則是以胡楊為代表的非地帶性多年生喬木,在干旱區(qū)生態(tài)系統(tǒng)中的地位不可代替,其擴(kuò)展則受制于現(xiàn)有胡楊林的樹(shù)勢(shì)及地下水位,其面積增長(zhǎng)速率明顯緩于草地,至2020年比例僅為6.71%,是中遠(yuǎn)期的修復(fù)目標(biāo),也是塔河下游生態(tài)修復(fù)的關(guān)鍵。
既有的現(xiàn)狀輸水方式,已使河岸帶地下水位呈頂托態(tài)勢(shì),地下水轉(zhuǎn)化比例下降,造成河岸帶無(wú)效潛水蒸發(fā)增大。地勢(shì)低洼的自然漫溢區(qū),長(zhǎng)期滯水,導(dǎo)致生態(tài)貢獻(xiàn)度較低的蘆葦?shù)戎参锟焖僭黾樱哟罅藷o(wú)效耗散;未過(guò)水區(qū)域地下水則只能通過(guò)側(cè)滲補(bǔ)給,補(bǔ)給速率受限于水位梯度,兩河區(qū)域地下水位抬升不均衡。后期輸水可基于地形,通過(guò)劃分漫溢區(qū),采用輪滲輸水方式避免地下水頂托,可減小無(wú)效蒸發(fā)以合理調(diào)配生態(tài)用水。
本研究中也發(fā)現(xiàn),塔河下游距離河道2 km以外區(qū)域天然植被覆蓋度極低,胡楊和檉柳以單株?duì)顟B(tài)分布,受柵格尺度效應(yīng)影響大[21];若能采用5 m等更高空間分辨率的遙感數(shù)據(jù),則可彌補(bǔ)30 m尺度研究中局部信息丟失的不足。
(1) 本文將改進(jìn)的P-M法系統(tǒng)地應(yīng)用于塔河下游,在30 m空間柵格上進(jìn)行了天然植被耗水量的計(jì)算,確定出下游各類(lèi)土地利用上植被系數(shù)Kc林地0.85,草地0.44,荒漠植被0.21,土壤水分限制系數(shù)分林地0.65,草地0.46,荒漠植被0.35,是塔河下游植被修復(fù)研究與生態(tài)輸水效應(yīng)評(píng)價(jià)的一項(xiàng)基礎(chǔ)成果。
(2) 經(jīng)過(guò)21次生態(tài)輸水后,塔河下游兩河區(qū)域天然植被面積由2000年的376.16 km2(占兩河區(qū)域17.69%)上升至2020年的504.19 km2(31.39%);林草地覆蓋度分別由0.20,0.11提高至0.37,0.22,天然植被耗水量由1.25×108m3/a增長(zhǎng)至2.22×108m3/a;土地利用格局變化趨勢(shì)和天然植被改善狀況對(duì)生態(tài)輸水響應(yīng)良好,生態(tài)輸水效應(yīng)顯著。
(3) 2000—2020年兩河區(qū)域多年平均天然植被耗水量為1.90×108m3/a,其中林地耗水量為6.50×107m3/a,草地為4.60×107m3/a,荒漠為7.90×107m3/a,荒漠區(qū)無(wú)效耗散占比高;且受限于既有現(xiàn)狀輸水方式固化的制約,植被耗水量主要集中于主河道沿線及自然漫溢區(qū),加劇了水面無(wú)效蒸發(fā),生態(tài)輸水策略有必要根據(jù)耗水方式的變化而進(jìn)行相應(yīng)的改進(jìn)。