馬思遠(yuǎn) 劉曉燕 董國濤 高云飛

摘要：1999年以來黃河河龍區(qū)間西北片入黃沙量減少約94%，分析植被變化及其驅(qū)動力對科學(xué)認(rèn)識該區(qū)來沙銳減原因具有重要意義。利用遙感調(diào)查數(shù)據(jù)，分析了該區(qū)1978-2017年林草植被變化程度和變化過程，并對這一時段植被變化的主要驅(qū)動力進(jìn)行了研究。結(jié)果表明：過去40a該區(qū)植被一直處于不斷改善的過程中，其中1998-2017年的貢獻(xiàn)接近80%；2016年研究區(qū)林草地植被蓋度平均達(dá)59.7%、區(qū)域林草覆蓋率為48.8%，均約為 1978年的3倍；2000年以前，人類對林草植被的破壞是該區(qū)林草覆蓋狀況惡化的決定性因素；近十幾年來，人類對植被的破壞很小，降雨是影響林草植被狀況的關(guān)健因素。

關(guān)鍵詞：植被；變化；驅(qū)動力；河龍區(qū)間；黃土高原

中圖分類號：S714.7；TV882.1 文獻(xiàn)標(biāo)志碼：A doi：10.3969/j.issn.1000-1379.2018.05.022

1 研究區(qū)概況

黃河河口鎮(zhèn)至龍門區(qū)間西北片（簡稱河龍間西北片）屬黃土高原水土流失最嚴(yán)重的地區(qū)，包括皇甫川、孤山川、窟野河和禿尾河等支流，涉及內(nèi)蒙古和陜西兩省（區(qū)）的準(zhǔn)格爾、伊金霍洛、東勝、府谷、神木、佳縣、榆陽等7個縣（旗、區(qū)）[1]，總面積19 444km2。受人類活動干擾較少的1954-1979年研究區(qū)4條支流把口斷面輸沙量為2.47億t/a，1999年以來人黃沙量急劇減少，較1954-1979年減少94%，引起人們高度關(guān)注。

研究區(qū)多年平均降水量390mm，其中4-9月降雨量344mm，地表土壤有黃綿土、風(fēng)積沙、砒砂巖和沙礫石等多種類型，林草植被類型主要為草本植物和灌木。土地利用類型主要為植被蓋度不同的灌草地，耕地面積占總面積的比例只有20%，主要分布在沿黃的黃土丘陵區(qū)和河川，其中梯田面積占總面積的比例不足1%[2]。

眾所周知，土壤、地形和植被是決定區(qū)域水土流失程度的關(guān)鍵下墊面因素，其中植被是最活躍、最容易受人類干擾的因素。有關(guān)研究表明，20世紀(jì)80年代以來黃土高原地區(qū)氣候變化趨勢明顯[3]；植被生長與氣候變化具有密切的相關(guān)性，特別是在干旱半干旱地區(qū)，水資源是影響生態(tài)恢復(fù)和生態(tài)建設(shè)的根本因素[4]；張翀等認(rèn)為，降水量增加促使黃土高原地區(qū)植被改善、干旱程度降低[5]；沈斌等指出，生長季內(nèi)植被NDVI與同期氣溫和降水量均成高度正相關(guān)，時滯尺度為1個月，植被NDVI對月平均氣溫及降水響應(yīng)最為強(qiáng)烈[6]；植被覆蓋的時空變化與人類活動有密切的關(guān)系，它深刻反映了人類活動的過程[7-10]。不過，現(xiàn)有文獻(xiàn)很少聚焦在河龍間西北片，且研究時段大多截至2010年前后，對河龍間西北片植被變化的原因也鮮見深入分析。鑒于研究區(qū)耕地不多、土地利用類型主要為灌草地的特點，本文對該區(qū)植被變化及其驅(qū)動力進(jìn)行分析，為科學(xué)認(rèn)識人黃泥沙銳減原因提供參考。

2 數(shù)據(jù)采集與處理

（1）植被數(shù)據(jù)。植被數(shù)據(jù)全部來源于衛(wèi)星遙感影像，利用空間分辨率為30～56m的衛(wèi)星遙感影像，結(jié)合地面樣方調(diào)查，提取相關(guān)支流1978年、1998年、2016年的土地利用數(shù)據(jù)和林草地植被蓋度數(shù)據(jù)[2]。為闡明流域林草植被覆蓋狀況，引入“林草覆蓋率”概念，即林草葉莖正投影面積占易侵蝕區(qū)面積的比例，其中“易侵蝕區(qū)面積”指流域內(nèi)剔除地面坡度小于30的河川地、平原區(qū)、建設(shè)用地、石山區(qū)后的山丘區(qū)土地面積，主要由不同植被蓋度的林草地（荒草地）、坡耕地和梯田組成，該區(qū)域顯然是入河泥沙的主要來源地。

利用衛(wèi)星遙感影像，得到了研究區(qū)1981-2017年的植被歸一化指數(shù)（NDVI）。NDVI是遙感影像的近紅外波段和紅外波段反射率的比值，它與植被長勢、生物量、蓋度和葉面積指數(shù)等有極強(qiáng)相關(guān)性，可以在相當(dāng)程度上反映植被覆蓋狀況，且有逐年數(shù)據(jù)、能部分消除輻射條件變化對反演參數(shù)的影響[6]。研究采用的各年NDVI值均為年最大值，數(shù)據(jù)源為第三代GIMMS NDVI數(shù)據(jù)集（NDV13g，v1.0）。

（2）氣象數(shù)據(jù)。降雨數(shù)據(jù)主要源自黃河水文年鑒，個別氣象數(shù)據(jù)源自國家氣象科學(xué)數(shù)據(jù)共享服務(wù)平臺，氣溫數(shù)據(jù)全部源自國家氣象科學(xué)數(shù)據(jù)共享服務(wù)平臺。

（3）社會經(jīng)濟(jì)數(shù)據(jù)。利用1981-2016年《陜西區(qū)域統(tǒng)計年鑒》《榆林統(tǒng)計年鑒》《中國縣（市）社會經(jīng)濟(jì)統(tǒng)計年鑒》《中國分縣農(nóng)村經(jīng)濟(jì)統(tǒng)計概要》和《神木縣國民經(jīng)濟(jì)和社會發(fā)展統(tǒng)計公報》等，結(jié)合實地問卷調(diào)查，得到研究區(qū)國內(nèi)生產(chǎn)總值（GDP）、總?cè)丝?、農(nóng)業(yè)人口和耕地面積等數(shù)據(jù)。

3 植被的時空變化分析

基于遙感影像提取的研究區(qū)3個典型年份的林草地植被蓋度和林草覆蓋率見表1。由表1可知，無論是林草地植被蓋度還是流域的林草覆蓋率，2016年的數(shù)值都達(dá)到1978年的3倍左右。1978年，研究區(qū)林草地植被蓋度平均為20.0%，其中蓋度≤30%的低蓋度林草地面積占比高達(dá)90%；而在2016年，林草地植被蓋度平均為59.7%，其中蓋度≤30%的低蓋度林草地面積占比只有1%、蓋度≥50%的高蓋度林草地面積占比約為78%。

圖1為研究區(qū)1978年、2016年林草覆蓋率空間格局，可以看出，1978年研究區(qū)林草覆蓋率大都不足20%，其中窟野河中部和禿尾河上中游的風(fēng)沙區(qū)甚至不足10%，而2016年大部分地區(qū)的林草覆蓋率為40%～60%，只有降雨和土壤條件均較差的西北部、耕地面積占比大的禿尾河下游仍不足40%。事實上，禿尾河下游2016年林草地的植被蓋度已超過60%，有些區(qū)域林草覆蓋率不足40%的原因主要是耕地面積較大、林草地面積的占比較小。

由表1可知，近40 a研究區(qū)林草覆蓋狀況的巨大改善主要發(fā)生在1998年以后，前20a（1978-1998年）的改善幅度僅占1/4，這與1981-2017年研究區(qū)NDVI的變化過程（見圖2）大體一致。由圖2可以看出，研究區(qū)植被一直處于不斷改善過程中，但2005年以前改善幅度不大，2006年以后改善幅度較大。土地利用情況分析表明[2]，研究區(qū)耕地不多，林草地面積占比20世紀(jì)70年代末為70%、2016年約為82%。因此，圖2揭示的1981-2017年NDVI變化情況可以大體反映研究區(qū)林草植被覆蓋變化情況。

研究區(qū)1978年以來的植被覆蓋變化過程可大體劃分成3個階段：一是1978-1985年，二是1986-2005年，三是2006-2017年。顯然，1986-2005年是植被的緩慢改善期，2006年以來是快速改善期。

4 近年來植被改善的原因分析

降雨和氣溫是影響植被生長的關(guān)鍵氣候因素，在植被生長發(fā)育的關(guān)鍵期，降雨越多、氣溫越高植被生長越茂盛（相應(yīng)的NDVI值越大）。在氣溫低、降雨少的河龍間西北片，4-9月降雨對植被生長發(fā)育非常重要[4]。由圖3可以看出：研究區(qū)4-9月降雨1981-1996年略偏枯，1997-2011年偏枯，2012-2016年明顯偏豐；4-9月平均氣溫與降雨變化有所不同，20世紀(jì)后半期研究區(qū)氣溫一直處于緩慢上升狀態(tài)，其中1998-2016年氣溫較70年代以前偏高2.1℃；2012年以來的降雨和氣溫條件均有利于研究區(qū)植被生長。有關(guān)研究結(jié)果表明，水熱條件有利于植被生長是近年來黃土高原植被改善的主要原因[11]。值得注意的是，降雨明顯偏枯的2006-2011年恰恰是研究區(qū)植被覆蓋狀況快速改善的時期（見圖2），這說明2006年以來除降雨因素外，還有其他因素促使植被覆蓋狀況大幅改善。

研究區(qū)1981-2017年不同時段NDVI與4-9月降雨量的關(guān)系見圖4。由圖4可以看出，1981年以來NDVI與4-9月降雨量關(guān)系線的斜率越來越大，其中：NDVI與1981-1995年降雨量幾乎無關(guān)，表明該時段降雨對植被的促生作用基本上被人類對植被的破壞所抵消，擴(kuò)耕、過度放牧和砍伐等人類活動是該時段植被覆蓋狀況惡化的決定性因素；2006年以來NDVI與4-9月降雨量之間成顯著正相關(guān)，相關(guān)系數(shù)高達(dá)0.77，表明該時期人類對植被的破壞活動很少，降雨條件是影響植被覆蓋狀況的關(guān)鍵因素。

以神木縣為典型，分析NDVI和人均GDP的變化過程，見圖5。2004年以前，神木縣人均GDP雖然一直在增加但增幅不大，2004-2011年快速增加，2011年以來基本穩(wěn)定，NDVI值與人均GDP相呼應(yīng)，在2006年以后增幅較大。這樣的良好呼應(yīng)關(guān)系，揭示了以煤炭產(chǎn)業(yè)為主的經(jīng)濟(jì)拉動促使越來越多的農(nóng)民進(jìn)城或進(jìn)礦務(wù)工，進(jìn)而使人類對植被的干擾破壞越來越少、植被自然修復(fù)的環(huán)境越來越好。

為證明上述推論的合理性，2017年上半年對研究區(qū)陜西省境內(nèi)的神木、府谷和佳縣等地進(jìn)行了社會調(diào)查，結(jié)果見表2。由表2可以看出，目前三縣農(nóng)民有80%左右外出務(wù)工，留守在農(nóng)村的基本上是60歲以上的老人。當(dāng)?shù)卮迕裾f，這樣的外出務(wù)工狀況已經(jīng)持續(xù)10多a。實際務(wù)農(nóng)人口減少，使很多耕地撂荒，牛羊數(shù)量大幅減少，很多村甚至不再養(yǎng)殖。

表2反映的情況在研究區(qū)內(nèi)蒙古境內(nèi)更為突出，2003年以來在農(nóng)牧民進(jìn)城務(wù)工的同時，內(nèi)蒙古相關(guān)旗（區(qū)）還實施了生態(tài)移民，因而現(xiàn)在不少村莊幾乎無人居住。

綜上分析認(rèn)為，研究區(qū)近年來植被的大幅改善，主要原因是大量農(nóng)牧民外出務(wù)工，減少了對植被的干擾和破壞，為植被自然修復(fù)創(chuàng)造了良好環(huán)境。此外，4-9月降雨偏豐也為近年植被恢復(fù)提供了有利條件。

5 結(jié)論

（1）20世紀(jì)70年代末以來，河龍間西北片林草植被不斷改善，到2016年林草地植被蓋度已由20.0%提高到59.7%，林草覆蓋率由14.0%提高到48.8%。2005年以前，該區(qū)植被改善程度不大，林草植被大幅改善主要發(fā)生在2006年以來。

（2）20世紀(jì)90年代及其以前，降雨對林草植被的促生作用基本上被人類破壞所抵消，人類的擴(kuò)耕、過度放牧和砍伐等活動是該區(qū)植被覆蓋狀況惡化的決定性因素。近十幾年來，80%的農(nóng)村人口進(jìn)城、進(jìn)礦務(wù)工，大幅減少了對土地的干擾和破壞，進(jìn)而使降雨成為改善植被覆蓋狀況的關(guān)鍵因素。

參考文獻(xiàn)：

[1]黃河上中游管理局.黃河流域水土保持圖集[M].北京：地震出版社，2012；92-93.

[2]劉曉燕.黃河近年水沙銳減成因[M].北京：科學(xué)出版社，2016：99-113，361-362.

[3]劉憲鋒，任志遠(yuǎn).西北地區(qū)植被覆蓋變化及其與氣候因子的關(guān)系[J].中國農(nóng)業(yè)科學(xué)，2012，45（10）：1954-1963.

[4]吳麗麗，任志遠(yuǎn)，張翀.陜北地區(qū)植被指數(shù)對水熱條件變化的響應(yīng)及其時滯分析[J].中國農(nóng)業(yè)氣象，2014，35（1）：103-108.

[5]張翀，任志遠(yuǎn)，李小燕.黃土高原植被對氣溫和降水的響應(yīng)[J].中國農(nóng)業(yè)科學(xué)，2012，45（20）：4205-4215.

[6]沈斌，房世波，余衛(wèi)國，等.NDVI與氣候因子關(guān)系在不同時間尺度上的結(jié)果差異[J].遙感學(xué)報，2016，20（3）；481-490.

[7]易浪，任志遠(yuǎn)，張種，等.黃土高原植被覆蓋變化與氣候和人類活動的關(guān)系[J].資源科學(xué)，2014，36（1）：166-174.

[8]信忠保，許炯心，鄭偉.氣候變化和人類活動對黃土高原植被覆蓋變化的影響[J].中國科學(xué)，2007，37（11）：1504-1514.

[9]徐小任，徐勇.黃土高原地區(qū)人類活動強(qiáng)度時空變化分析[J].地理研究，2017，36（4）：661-672.

[10]ZHAO J，YANG Y，ZI-IAO Q，et al.Effects of EcologicalRestoration Projects on Changes in Land Cover：A caseStudy on the Loess Plateau in China[J].ScientificReports，2017，7：44496.

[11]張翀，任志遠(yuǎn)，韋振鋒.近12年來黃土高原植被覆蓋對年內(nèi)水熱條件的響應(yīng)[J].資源科學(xué)，2013，35（10）：2017-2023.