張 萌,杜福光,張瀟文

(唐山師范學(xué)院 資源管理系,河北 唐山 063000)

0 引言

青藏高原是亞洲內(nèi)陸高原,有亞洲水塔之稱,是中國(guó)最大、世界海拔最高的高原,南起喜馬拉雅山脈南緣,北至昆侖山、阿爾金山脈及祁連山北緣,西部為帕米爾高原與喀喇昆侖山脈,東及東北部與秦嶺山脈西段及黃土高原相接,草原面積占總面積的60%,占全國(guó)總面積的32%,主要分布于青藏高原中部與東部地區(qū)[1-4],而高原的降水量變化對(duì)該地的草原生態(tài)具有重要影響。草原放牧業(yè)是青藏高原牧民的重要經(jīng)濟(jì)來(lái)源之一,也是青藏高原生態(tài)系統(tǒng)中的重要一環(huán),做好青藏高原草原及草甸的生態(tài)保護(hù)及合理利用具有十分重要的意義。很多學(xué)者對(duì)青藏高原降水來(lái)源及時(shí)空變化特征進(jìn)行了研究。黃鵬等[5]探討了青藏高原降水的時(shí)空變化特征,提出高原全年降水量為西北—東南向的均勻帶狀分布。李虎等[6]分析了主要的水來(lái)源,得出南部邊界輸送的水汽對(duì)青藏高原有巨大的貢獻(xiàn)。但目前對(duì)青藏高原植被分布與降水時(shí)空變化關(guān)系的研究不多,為了研究二者的對(duì)應(yīng)關(guān)系,利用青藏高原生長(zhǎng)季降水分布數(shù)據(jù)及NDVI數(shù)據(jù),采用皮爾遜相關(guān)系數(shù),計(jì)算降水量與NDVI的相關(guān)性,得出降水對(duì)青藏高原不同地區(qū)植被的影響程度,有助于揭示與反映青藏高原植被-水汽的影響與互動(dòng)關(guān)系,進(jìn)一步促進(jìn)青藏高原的生態(tài)保護(hù)。

1 資料方法

1.1 數(shù)據(jù)來(lái)源及預(yù)處理

NDVI數(shù)據(jù)來(lái)源及處理:數(shù)據(jù)來(lái)源于國(guó)家青藏高原科學(xué)數(shù)據(jù)中心整編的青藏高原NDVI數(shù)據(jù)集(國(guó)家青藏高原科學(xué)數(shù)據(jù)中心 (tpdc.ac.cn)),源數(shù)據(jù)來(lái)自MODIS產(chǎn)品,經(jīng)過(guò)數(shù)據(jù)格式轉(zhuǎn)換、投影、重采樣等預(yù)處理流程,現(xiàn)有格式為T(mén)IFF格式,投影為Krasovsky_1940_Albers投影,數(shù)據(jù)空間分辨率1000 m。選取2016—2020年的NDVI數(shù)據(jù),利用ArcGIS統(tǒng)一空間坐標(biāo)系。

降水?dāng)?shù)據(jù)來(lái)源與預(yù)處理:數(shù)據(jù)來(lái)源于國(guó)家科技基礎(chǔ)條件平臺(tái)—國(guó)家地球系統(tǒng)科學(xué)數(shù)據(jù)中心(http://www.geodata.cn)。國(guó)家地球系統(tǒng)科學(xué)數(shù)據(jù)中心1901—2021年中國(guó)1公里分辨率(0.008333°)逐月降水量數(shù)據(jù)集。該產(chǎn)品具有長(zhǎng)時(shí)間序列、較高分辨率、高精度的特征,可為我國(guó)歷史降水格局變化分析、陸地生態(tài)系統(tǒng)模型氣候因子驅(qū)動(dòng)及其他區(qū)域尺度與降水相關(guān)的生態(tài)、環(huán)境、資源等變量的格局分析提供數(shù)據(jù)支持。選用該數(shù)據(jù)集中2016—2020年覆蓋全國(guó)范圍的降水?dāng)?shù)據(jù),使用Arc GIS軟件對(duì)數(shù)據(jù)進(jìn)行格式轉(zhuǎn)換、數(shù)據(jù)裁剪、投影轉(zhuǎn)換、重采樣等處理。

1.2 相關(guān)性分析

采用皮爾遜(Pearson)相關(guān)系數(shù),分析降水與草地歸一化指數(shù)的相關(guān)性,計(jì)算公式如下:

(1)

式中,r為相關(guān)系數(shù);x、y為各相關(guān)變量的平均值。相關(guān)系數(shù)r的取值范圍為[-1,1],其值越偏近于1說(shuō)明變量之間的相關(guān)性越高。

2 研究?jī)?nèi)容

平均降水量分布。將2016—2020年全國(guó)月均降水量數(shù)據(jù)導(dǎo)入Arc GIS軟件中,利用裁剪工具,按照青藏高原的矢量圖形進(jìn)行裁剪,利用局部像元統(tǒng)計(jì)分析計(jì)算出2016—2020年的平均降水量,渲染成圖。

生長(zhǎng)季降水量分布。將2016—2020年全國(guó)月均降水量數(shù)據(jù)導(dǎo)入Arc GIS軟件中,根據(jù)青藏高原矢量圖進(jìn)行裁剪,利用像元統(tǒng)計(jì)工具分別計(jì)算2016—2020年每年5—9月的累計(jì)降水,渲染成圖。

生長(zhǎng)季降水與NDVI相關(guān)性計(jì)算。將NDVI數(shù)據(jù)與計(jì)算好的生長(zhǎng)季降水?dāng)?shù)據(jù)導(dǎo)入Arc GIS軟件中,統(tǒng)一坐標(biāo)系并改為投影坐標(biāo),重采樣數(shù)據(jù)并捕捉柵格,所有數(shù)據(jù)整理完成后根據(jù)公式計(jì)算相關(guān)系數(shù),利用柵格計(jì)算器求平均,再根據(jù)公式分別求出分子、分母,利用柵格計(jì)算器計(jì)算分子與分母,求出每個(gè)柵格點(diǎn)的相關(guān)系數(shù),渲染成圖。

3 結(jié)果與分析

3.1 平均降水量空間分布

由圖1可知,2016—2020年平均降水量最大值為2405 mm,分布于青藏高原南部地區(qū);最小值為2 mm,分布于青藏高原西部地區(qū),空間上呈現(xiàn)由東南向西北逐漸遞減。

圖1 2016—2020年平均降水量空間分布Fig.1 Spatial distribution of average precipitation from 2016 to 2020

3.2 生長(zhǎng)季降水的空間分布

由圖2可知,2016—2020年青藏高原生長(zhǎng)季(5—9月)降水整體相對(duì)穩(wěn)定,呈現(xiàn)東南多西北少的情況。2016年5—9月,累積降水量最高值為22 361 mm,最低值為10 mm。2017年5—9月,青藏高原累計(jì)降水量最高值為24 848 mm,最低值為37 mm。2018年5—9月,青藏高原累計(jì)降水量最高值為24 196 mm,最低值為47 mm。2019年5—9月,青藏高原累計(jì)降水量最高值為24 440 mm,最低值為1 mm。2020年5—9月,青藏高原累計(jì)降水量最大值為28 283 mm,最低值為42 mm。其中2020年相比其他年份,中部地區(qū)的降水有所減少。2016—2019年西部部分地區(qū)降水與周?chē)忻黠@差異,而在2020年該差異明顯減少。

圖2 2016—2020年生長(zhǎng)季累計(jì)降水空間分布Fig.2 Spatial distribution of cumulative precipitation in growing season from 2016 to 2020

3.3 NDVI的空間變化

由圖3可知,NDVI在空間上呈現(xiàn)西北低東南高的情況,南部地區(qū)5年內(nèi)變化較為明顯,其中2016年、2018年、2020年南部地區(qū)NDVI較低,而2017年、2019年南部地區(qū)的NDVI較高。

圖3 2016—2020年NDVI空間分布Fig.3 Spatial distribution of NDVI from 2016 to 2020

3.4 相關(guān)性計(jì)算結(jié)果及分析

通過(guò)計(jì)算可知,甘肅南部地區(qū)、青海東北部地區(qū)、新疆南部及西部地區(qū)與西藏西部地區(qū)的降水與NDVI的相關(guān)性較為明顯,即降水對(duì)該地植被影響作用較強(qiáng)。西藏中部地區(qū)相關(guān)性性中等偏高,降水與NDVI的相關(guān)性適中。青海南部、西藏東部、四川西北部的相關(guān)系數(shù)較低,降水與NDVI的相關(guān)性較弱,即降水對(duì)該地區(qū)植被生長(zhǎng)的影響較小,更多受到溫度濕度等因素的影響較大,但并非降水對(duì)其毫無(wú)影響。生態(tài)保護(hù)要著重從降水角度對(duì)相關(guān)系數(shù)較高的地區(qū)進(jìn)行保護(hù)與治理,如在旱季適時(shí)對(duì)降水與NDVI相關(guān)性較強(qiáng)的地區(qū)進(jìn)行灌溉,在相關(guān)性較弱的地區(qū)要多關(guān)注其他氣候因素。

4 植被生態(tài)保護(hù)措施

通過(guò)對(duì)青藏高原降水與NDVI的相關(guān)性分析發(fā)現(xiàn),甘肅南部地區(qū)、青海東北部地區(qū)、新疆南部地區(qū)、新疆西部地區(qū)、西藏西部地區(qū)的相關(guān)系數(shù)較高,降水與NDVI的相關(guān)性較為明顯,需更加關(guān)注這些地區(qū)因降水不充分導(dǎo)致的生態(tài)環(huán)境惡化問(wèn)題。可采用以下方法保護(hù)植被與生態(tài):適當(dāng)灌溉青藏高原極寒地區(qū),不宜開(kāi)展大規(guī)模灌溉,積極開(kāi)展牧區(qū)草原節(jié)水灌溉,探索草原節(jié)水灌溉新模式。采用機(jī)械動(dòng)力淋灑技術(shù),利用蓄積山上的雪融水或泉水引入防滲襯砌渠,在常壓下采用自流噴灌,無(wú)壓力時(shí)利用拖拉機(jī)為機(jī)械動(dòng)力,加裝水泵、貯水箱、水管及淋灑器等設(shè)備給水加壓,牽引灌溉設(shè)備,移動(dòng)淋灑灌水。利用步行式噴灌技術(shù)在有地下水源與電力能源的地方打井開(kāi)源,水泵提水,管道輸水至步行式噴灌設(shè)備,動(dòng)力移動(dòng)噴灌。采用集雨滴灌技術(shù),將雨水集蓄與滴灌相結(jié)合,在牧民居住相對(duì)集中的地方把溫室暖棚、牲畜防寒棚、人居屋舍作為集雨場(chǎng)收集雨水,采用滴灌或滲灌等微型節(jié)水灌溉方法對(duì)網(wǎng)圍欄內(nèi)的牧草進(jìn)行補(bǔ)充灌溉。這些節(jié)水灌溉措施可解決水源難以保障、能源動(dòng)力不足等問(wèn)題,噴灑均勻,移動(dòng)方便,滿足開(kāi)源節(jié)流與節(jié)水灌溉的需要。