新疆科克蘇濕地植被生物量時空分布特征及與水文連通的關(guān)系

湯子同,李興麗,劉華兵,李謙維,高俊琴

北京林業(yè)大學(xué)生態(tài)與自然保護(hù)學(xué)院,北京 100083

濕地生態(tài)系統(tǒng)擁有豐富的物種和較高的生產(chǎn)力,發(fā)揮著巨大的生態(tài)環(huán)境效益。植被生物量作為衡量生態(tài)系統(tǒng)植物生產(chǎn)力和健康狀況的重要指標(biāo),對濕地生態(tài)系統(tǒng)碳儲存功能和物質(zhì)循環(huán)具有關(guān)鍵影響[1]。研究表明,濕地植被生物量積累受水文過程調(diào)控,尤其是水文連通度、水位及其變化等[2]。其中,水文連通度是指在水文循環(huán)各要素內(nèi)部和各要素之間,物質(zhì)、能量及生物以水為媒介進(jìn)行遷移和傳遞的順暢程度[3]。研究表明,水文連通通過影響水分的供給、養(yǎng)分的供應(yīng)及植物物種多樣性等,直接或間接地影響濕地植物群落的生長狀況,從而影響植被生物量的大小[4]。此外,水文連通還通過調(diào)整植被群落結(jié)構(gòu)與分布格局,進(jìn)一步影響植被生物量及碳庫分布格局[5—6]。因此,維持適宜的水文連通有助于促進(jìn)濕地內(nèi)外物質(zhì)和能量的交換,推動養(yǎng)分循環(huán),從而促進(jìn)植物生長,增強(qiáng)濕地生態(tài)系統(tǒng)的固碳能力[7—8]。然而,水文連通程度過高可能不利于某些植物的生長,從而影響植被群落生物量的積累,降低生態(tài)系統(tǒng)碳輸入及儲存功能。因此,研究水文連通和植被生物量之間的定量關(guān)系對于濕地生態(tài)系統(tǒng)的保護(hù)和管理具有重要的理論和實(shí)踐意義。以往研究更多從定性角度描述水文連通對植物、水生動物或微生物的分布或遷移的影響[9],定量研究植被生物量的時空變化,而很少揭示水文連通與植被生物量之間的定量關(guān)系,且不同區(qū)域和生態(tài)系統(tǒng)的適宜水文連通度可能存在差異,因此需開展?jié)竦厮倪B通與植被生物量的關(guān)系研究,為濕地保護(hù)管理以及水文連通調(diào)控提供科學(xué)支撐。

新疆科克蘇濕地位于額爾齊斯河與其支流克蘭河交匯的三角地帶,包括河流、沼澤(森林沼澤和草本沼澤)、灘地等多種濕地類型[10]。受阿爾泰山冰雪融水補(bǔ)給、克蘭河水壩蓄水[11]以及降水的影響,科克蘇濕地水體斑塊及水文連通變化較強(qiáng),具有明顯的季節(jié)性特征,對植物生物量積累和植被碳庫具有重要影響。目前有關(guān)科克蘇濕地的研究主要集中在水文情勢變化[12]以及植物物種或群落分布上[13—14],而水文連通與植被生物量的關(guān)系還很不明確。因此,為了闡明科克蘇濕地植被地上生物量的時空分布特征及其與水文連通的關(guān)系,本研究基于landsat-8影像和野外實(shí)測數(shù)據(jù)反演科克蘇濕地生長季不同月份的植被地上生物量,并通過提取和計算濕地斑塊的水文連通度,建立水文連通與植被地上生物量及植被碳庫之間的量化關(guān)系,為科克蘇濕地的植被保護(hù)管理以及水文連通調(diào)控提供科學(xué)支撐。

1 研究區(qū)概況

科克蘇濕地國家級自然保護(hù)區(qū)位于新疆阿勒泰市西南部(47°28′31″— 47°40′9″N, 87°9′12″—87°34′59″E),克蘭河與額爾齊斯河交匯處[13](圖1)。保護(hù)區(qū)總面積約307 km2,海拔高度476—796 m,年平均溫度4.9 ℃,年平均降水量和蒸發(fā)量分別為112.6 mm和2000 mm,屬于溫帶大陸性干旱半干旱氣候[13]?？瓶颂K濕地國家級自然保護(hù)區(qū)是新疆北部最大的沼澤濕地,植物資源豐富,白柳、吐倫柳、布爾津柳、灰毛柳等僅分布于此,額爾齊斯河木蓼、雪白睡蓮、阿勒泰菱角等為該區(qū)特有物種[10]?？瓶颂K濕地作為額爾齊斯河國際生物廊道的重要區(qū)段,為水生野生動物和河谷植物提供了良好的棲息和繁衍環(huán)境[14]。

圖1 科克蘇國家級自然保護(hù)區(qū)范圍及實(shí)驗(yàn)樣點(diǎn)分布圖Fig.1 Map of Kekesu national nature reserve and sample point distribution

2 數(shù)據(jù)與方法

2.1 影像數(shù)據(jù)來源與處理

本研究使用的遙感影像為Landsat-8 OLI數(shù)據(jù),數(shù)據(jù)來源于美國USGS地球資源觀測與科學(xué)中心(http://earthexplorer.usgs.gov/)?？瓶颂K濕地具有明顯的季節(jié)性特征,一般在6月上旬上游克蘭河托克孜塔勞段攔河閘開閘放水,對科克蘇濕地進(jìn)行生態(tài)補(bǔ)水;8月科克蘇濕地積水逐漸消退,此時牧民會進(jìn)行植被刈割[15],為牲畜提供草料,解決冬春飼草不足的問題。因此,6—8月植被地上生物量與濕地水體斑塊變化較為明顯。本研究選取2022年6月、7月和8月影像進(jìn)行水體斑塊提取、水文連通計算以及植被地上生物量反演。使用ENVI5.3中的Radiometric Calibration工具及FLAASH工具對影像進(jìn)行輻射定標(biāo)和大氣校正,根據(jù)研究區(qū)范圍對影像進(jìn)行剪裁。

2.2 濕地水體斑塊提取

采用WIW法以及WI2015水體指數(shù)相結(jié)合對科克蘇濕地的水體斑塊(包含被地表植被覆蓋的沼澤)進(jìn)行提取。WIW法[16]公式如下:

Landsat8:WIW=NIR≤0.1735且SWIR2≤0.1035;

WI2015水體指數(shù)[17]:

WI2015=1.7204+171GREEN+3RED-70NIR-45SWIR1-71SWIR2;

計算完成后在IDL中采用OATU大津法識別WI2015水體指數(shù)閾值,利用閾值對濕地水體進(jìn)行提取。根據(jù)研究區(qū)已有Landsat全色影像進(jìn)行人工目視校正,校正后的水體在ArcGIS 10.8中生成100個隨機(jī)點(diǎn),利用新疆科克蘇濕地已有資料,采用二分類混淆矩陣對已提取斑塊精度進(jìn)行評價。提取結(jié)果的總體解譯精度在90%以上,Kappa系數(shù)為78.99%,提取精度較高。

2.3 濕地水文連通度計算

根據(jù)景觀連通度原理,計算可能連通性指數(shù)(probability of connectivity index,PC)及斑塊重要性指數(shù)dPC。PC可以評估整體的水文連通可能程度,但它缺乏評估單個水體斑塊重要性的能力。因此,為了識別單個斑塊的水文連通度,用dPCi來表示單個斑塊i與其他斑塊的連通度[18—19]。

2.3.1閾值選取

連通度指數(shù)的計算需要確定濕地水體斑塊間連通的距離閾值[20]。距離閾值是保證斑塊連通的最小值,當(dāng)斑塊間的距離大于閾值時,則認(rèn)為斑塊間不連通;當(dāng)斑塊間的距離小于或等于閾值時,則認(rèn)為斑塊間是連通的[21]。距離閾值是通過整體連通性指數(shù)(Integral Index of Connectivity, IIC)和可能連通性指數(shù)(Probability of Connectivity Index,PC)的變化計算確定, IIC和PC計算公式如下[18,22]:

為減小隨機(jī)性,合理確定穩(wěn)定的距離閾值,本文選取三幅不同年份相同月份的影像進(jìn)行IIC和PC的計算。當(dāng)兩個連通性指數(shù)(IIC和PC)不隨距離產(chǎn)生較大變化時(即斜率接近于0),可以確定為穩(wěn)定的距離閾值。結(jié)果顯示,研究區(qū)域不同年份的整體連通性指數(shù)IIC和可能連通性指數(shù)PC呈現(xiàn)一致的變化趨勢,IIC和PC隨距離閾值的增加而增加(圖2)。距離閾值為50—400 m時,IIC和PC迅速增長,此區(qū)間內(nèi)斑塊的連通性不穩(wěn)定,易受閾值影響。距離閾值為400—600 m時,IIC和PC增長但趨近于平穩(wěn),區(qū)域內(nèi)斑塊連通性較穩(wěn)定,受閾值變化影響較小。距離閾值大于600 m時,IIC和PC保持平穩(wěn),整個生境內(nèi)的斑塊被認(rèn)為是連通的[21]。綜上,將計算水文連通度的距離閾值定為600 m。

圖2 不同距離閾值下IIC和PC值變化Fig.2 Changes of IIC and PC values at different distanceIIC:整體連通性指數(shù) Integral index of connectivity;PC:可能連通性指數(shù) Probability of connectivity index

2.3.2水文連通度指數(shù)計算

利用Conefor Sensinode 2.6軟件計算水文連通度指數(shù)。使用ArcGIS 10.8軟件中的 Conefor 插件,生成計算連通度指數(shù)所需的節(jié)點(diǎn)和連接文件,將節(jié)點(diǎn)和連接文件在Conefor Sensinode 2.6中打開,輸入閾值距離以及研究區(qū)面積,計算水文連通指數(shù)。斑塊重要性指數(shù)dPC可表示某個斑塊對區(qū)域整體水文連通度的貢獻(xiàn)程度,計算公式如下[18]:

dPCi=100×(PC-PC′)/PC

式中,PC為研究區(qū)原本連通性指數(shù),PC′為去除單個斑塊i后的連通性指數(shù),dPCi用來指示單個斑塊i對整個區(qū)域連通性的貢獻(xiàn)程度,同時,dPCi可在一定程度上反映斑塊i與其他斑塊的連通性[18—19],以dPC值計算的水文連通度可反映濕地斑塊之間連通的概率[18—19]。

2.4 遙感影像反演植被生物量及植被碳庫計算

2.4.1 植被生物量測定

于2022年7月進(jìn)行了濕地植被生物量樣品采集和測定,保證了樣品采集時間與遙感影像反演時間相一致。實(shí)地樣點(diǎn)的選擇遵循植被類型分布的均勻性和代表性以及交通可達(dá)性,共設(shè)置30個樣點(diǎn),每個樣點(diǎn)設(shè)置3個重復(fù)樣方,樣方大小為1 m×1 m。樣點(diǎn)分布如圖1所示。采集時,將樣方內(nèi)植物地上部分沿地面貼根部剪下,共采集90個植物地上樣品,并隨機(jī)選擇5個樣方進(jìn)行地下生物量采集,用于計算地下/地上的比值。植物樣品帶回實(shí)驗(yàn)室,放入 65℃ 的恒溫烘箱中烘干72 h至恒重后稱重,所得重量即為植被地上生物量。

2.4.2遙感影像反演植被生物量數(shù)據(jù)驗(yàn)證

選取常用的8種植被指數(shù)NDVI、RVI、SAVI、MSAVI、DVI、EVI、GNDVI進(jìn)行計算,用實(shí)測生物量與8種植被指數(shù)進(jìn)行相關(guān)性分析,其中NDVI與SAVI兩種指數(shù)與生物量存在顯著的正相關(guān)關(guān)系(表1)。

表1 植被指數(shù)計算公式及生物量相關(guān)性分析Table 1 Calculation formulas of vegetation indices and correlation analysis of biomass

NDVI與SAVI的指數(shù)計算結(jié)果如圖3所示。采用同緯度地帶已有的SAVI指數(shù)生物量回歸模型與和NDVI指數(shù)生物量回歸模型[30—31](表2)對科克蘇濕地地上生物量進(jìn)行反演,并根據(jù)實(shí)測樣本對反演結(jié)果進(jìn)行驗(yàn)證,利用均方根誤差(RMSE)和決定系數(shù)(R2)判定植被生物量反演模型的精度,計算公式如下[32]:

表2 兩種植被指數(shù)反演模型Table 2 The models of the vegetation index

圖3 NDVI與SAVI植被指數(shù)Fig.3 NDVI and SAVI vegetation indexNDVI:歸一化植被指數(shù) Normalizaed difference vegetation index;SAVI:土壤調(diào)整植被指數(shù) Soil-adjusted vegetation index

式中,n為樣本數(shù)量,yi和yi′分別為實(shí)測植被生物量和遙感反演的植被生物量,y′為實(shí)測植被生物量的平均值。

擬合結(jié)果如圖4所示,SAVI指數(shù)構(gòu)建的模型和NDVI指數(shù)模型精度都較高(R2>0.91),但SAVI模型擬合線與1∶1線差距較小,因此選用SAVI指數(shù)構(gòu)建的生物量模型對整個研究區(qū)的植被地上生物量進(jìn)行反演。

圖4 兩種模型模擬地上生物量與實(shí)測地上生物量擬合效果Fig.4 Relationships between measured and estimated biomass by two models

2.4.3植被總生物量及碳庫估算

根據(jù)樣點(diǎn)實(shí)地采集測定的植物平均根冠比估算植被總生物量[33]:

Biomassi=AGB+BGB

式中,AGB為植被地上生物量,BGB為植被地下生物量,由AGB除以植物平均根冠比獲得。

濕地植被碳庫采用國際上常用的植被總生物量乘以碳轉(zhuǎn)化系數(shù)估算[34]:

Ci=Biomassi×vi

式中Ci為濕地植被碳庫,Biomassi為植被總生物量,vi為植被碳轉(zhuǎn)化系數(shù),取值0.41[35]。

3 結(jié)果與分析

3.1 科克蘇濕地水體斑塊分布

科克蘇濕地斑塊分布有明顯的季節(jié)性特征(圖5)。6月的濕地水體斑塊面積較大,占保護(hù)區(qū)面積的63.12%,除了北部和東南部外,大多數(shù)區(qū)域處于淹水狀態(tài);7月初水體大幅退去,小面積水體斑塊零散分布在整個保護(hù)區(qū)中;7月末水體斑塊持續(xù)減少,北部的闊克蘇村、蒙古灣及保護(hù)區(qū)的核心區(qū)依舊分布著水體斑塊,其它區(qū)域大多無水體斑塊分布;8月水體斑塊面積僅占研究區(qū)面積的6.27%,只有主河道和北部的闊克蘇村和蒙古灣依舊處于淹水狀態(tài)。

圖5 科克蘇濕地不同月份濕地水體斑塊分布Fig.5 Distribution of surface water bodies in different months in Kekesu wetland

3.2 科克蘇濕地植被生物量時空分布

科克蘇濕地植被生物量具有明顯的時空變化特征(圖6)。6月植被生物量較高的區(qū)域分布在東北部的闊克蘇農(nóng)場,地上生物量平均值為226 g/m2,總生物量平均值為809 g/m2。7月初植被生物量較高區(qū)域分布在西部的阿克鐵熱克村及東部的巴勒喀木斯村,地上生物量平均值為915 g/m2,總生物量平均值為3275 g/m2。7月末植被生物量較高區(qū)域主要分布在西部及南部的阿勒特拜,地上生物量平均值為836 g/m2,總生物量平均值為2991 g/m2。8月植被生物量較高區(qū)域分布在東北部,地上生物量平均值為598 g/m2,總生物量平均值為2142 g/m2。

圖6 科克蘇濕地不同月份植被地上生物量空間分布Fig.6 Spatial distribution of aboveground biomass in different months in Kekesu wetland

7月地上生物量達(dá)到生長季最高值,該時段科克蘇濕地的植被總生物量為1.09×109kg,最大總生物量為4832 g/m2。

3.3 科克蘇濕地植被地上生物量及植被碳庫與水文連通的關(guān)系

科克蘇濕地植被地上生物量與水文連通度呈現(xiàn)非線性關(guān)系(圖7)。水文連通度0—0.6范圍內(nèi),植被地上生物量隨水文連通度增加而增加;水文連通度0.6—1范圍內(nèi),植被地上生物量隨水文連通度增加而降低,植被地上生物量在水文連通度0.6左右達(dá)到最大。

圖7 科克蘇濕地水文連通與植被地上生物量的關(guān)系Fig.7 The relationships between hydrological connectivity and aboveground biomass in Kekesu wetland

科克蘇濕地植被碳庫較高的區(qū)域分布在額爾齊斯河道及克蘭河河道兩側(cè)的濕地區(qū)域以及保護(hù)區(qū)東部部分區(qū)域,北部及東南部部分區(qū)域植被碳庫較低(圖8)。保護(hù)區(qū)內(nèi)植被總碳庫為4.5×1011kg C?？瓶颂K濕地植被碳庫與水文連通呈現(xiàn)非線性關(guān)系,與地上生物量和水文連通關(guān)系類似,在水文連通度0.6左右,植被碳庫達(dá)到最大,在1197—1406 g C/m2之間。

圖8 科克蘇濕地植被碳庫分布及水文連通與植被碳庫的關(guān)系Fig.8 Distribution of plant carbon pool and relationship between hydrological connectivity and plant carbon pool in Kekesu wetland

4 討論

4.1 科克蘇濕地植被地上生物量

科克蘇濕地植被生物量隨季節(jié)變化呈現(xiàn)先增加后減少的趨勢。5—6月是科克蘇濕地降水量和徑流量較大的2個月[36],6月初上游攔河閘放水,科克蘇濕地進(jìn)入淹水狀態(tài),此時大部分植物處于生長初期,其生長會受到淹水限制。7月氣溫升高,濕地水分迅速蒸發(fā),導(dǎo)致濕地水體斑塊面積顯著減少,受水澇影響的植物快速生長[37],植被地上生物量逐漸增加,達(dá)到0—1350 g/m2,這與神祥金等[38]對中國草本沼澤反演的植被地上生物量密度范圍一致。濕地水分的持續(xù)蒸散發(fā)使得8月份只有主河道依舊處于淹水狀態(tài),大部分地區(qū)不再積水,此時科克蘇濕地會進(jìn)行刈割,為牲畜提供草料,因此植被地上生物量下降。

科克蘇濕地植被生物量呈現(xiàn)聚集分布的空間分布模式。額爾齊斯河及克蘭河河道兩側(cè)的濕地和東部的闊克蘇農(nóng)場為高生物量區(qū)域;北部阿熱勒齊及闊克蘇村和東南部薩爾胡松鄉(xiāng)為低生物量區(qū)域。這一空間分布格局的形成,很大程度上是濕地植被對生境條件長期適應(yīng)尤其是對水文條件適應(yīng)的結(jié)果。高生物量區(qū)域的水體斑塊分布較多,水分較為充足,多為水生及濕生植被,如蘆葦(Phragmitesaustralis)、兩棲蓼(Persicariaamphibia)、水蓼(Persicariahydropiper)、水蔥(Schoenoplectustabernaemontani)、荸薺(Eleocharisdulcis)、花藺(Butomusumbellatus)等,蓋度大,生物量高。低生物量區(qū)域水體斑塊較少,水分相對匱乏,植被類型多為堿蓬(Suaedaglauca)、鹽生車前(Plantagosalsa)、小獐毛(Aeluropuspungens)、蒼耳(Xanthiumstrumarium)、披堿草(Elymusdahuricus)、車軸草(Galiumodoratum)等,生物量相對較低。

4.2 水文連通對科克蘇濕地植被地上生物量及植被碳庫的影響

水文連通與植被生物量的關(guān)系表明,濕地水文連通度對植被生物量及植被碳儲存功能具有顯著影響。我們發(fā)現(xiàn),科克蘇濕地植被地上生物量和植被碳庫隨著水文連通度呈現(xiàn)先增加后減小的趨勢,在水文連通度0.6左右時,植物生物量和植被碳庫達(dá)到最高。研究表明,水文連通可以通過影響植被群落結(jié)構(gòu)與分布格局,改變植被地上生物量[4,39],且適宜的水文連通可在時間和空間尺度上有效促進(jìn)碳儲存[40],其能加強(qiáng)植物的光合作用,促進(jìn)光合產(chǎn)物向地上或地下部分分配,從而碳庫增加。同時,適宜的水文連通能加強(qiáng)物質(zhì)和能量循環(huán),有利于植物生長。水文連通度過高時,水動力的加強(qiáng)不利于植物幼苗萌發(fā)以及生長[41];同時水位超過大多數(shù)植物的耐受范圍,植物受缺氧脅迫,光合作用能力下降;過多的水分還會導(dǎo)致土壤中有害細(xì)菌和真菌過度繁殖,引起根腐病,嚴(yán)重的情況下會導(dǎo)致植物死亡[42]。科克蘇濕地水文連通度0.6左右時植被地上生物量及植被碳庫最高,這與Wang 等[41]的研究中幼苗移植策略的水文連通對植被影響的結(jié)果一致;Zhang[43]等的研究中水文連通與植被的關(guān)系也呈現(xiàn)拋物線關(guān)系,但水文連通度在0.4左右時,植被覆蓋度達(dá)到峰值,這可能是因?yàn)檠芯繀^(qū)地理位置的不同以及植被類型的不同而導(dǎo)致的差異。因此,在科克蘇濕地水文調(diào)控時,可考慮維持適宜的水文連通度范圍,維持適宜的水位或土壤水分條件,從而促進(jìn)植被生物量積累和碳儲存功能提升,維持濕地生態(tài)系統(tǒng)的穩(wěn)定性。

5 結(jié)論

(1)科克蘇濕地斑塊分布有明顯的時空特征。6月的濕地水體斑塊面積最大,占研究區(qū)面積的63.12%,大多數(shù)區(qū)域處于淹水狀態(tài);7月濕地水體斑塊面積持續(xù)減小,呈零散分布在整個保護(hù)區(qū)中;8月水體斑塊面積僅占研究區(qū)面積的6.27%,只有主河道和北部的闊克蘇村和蒙古灣依舊處于淹水狀態(tài)。

(2)科克蘇濕地植被生物量呈現(xiàn)聚集分布的空間分布模式。額爾齊斯河及克蘭河河道兩側(cè)的濕地和東部的闊克蘇農(nóng)場為高生物量區(qū)域,北部阿熱勒齊及闊克蘇村和東南部薩爾胡松鄉(xiāng)為低生物量區(qū)域。7月生物量達(dá)到最高值,植被總生物量為1.09×109kg。8月受到刈割的影響,生物量降低。

(3)濕地水文連通度對植被生物量及植被碳儲存功能具有顯著影響?？瓶颂K濕地水文連通與植被地上生物量及植被碳庫之間呈現(xiàn)拋物線關(guān)系,水文連通度0.6左右時植被地上生物量和植被碳庫最大,有利于植被碳固定,植被總碳庫達(dá)到4.5×1011kg C。本研究揭示了科克蘇濕地植被生物量的時空分布特征,建立了科克蘇濕地水文連通度與植被地上生物量及植被碳庫的量化關(guān)系,明確了適宜的水文連通度對植被生物量積累存在促進(jìn)作用,可為濕地水文連通調(diào)控和植被碳儲存功能提升提供有效參考。