摘要 基于2000～2010年土地利用/覆被變化、植被覆蓋度、坡耕地分布等數(shù)據(jù)，分析陜甘寧退耕還林/草工程區(qū)耕地轉(zhuǎn)換為林地或草地狀況、不同坡耕地退耕情況及退耕區(qū)植被覆蓋度變化等，從而揭示該區(qū)域退耕還林/草工程的生態(tài)成效。結(jié)果表明：①2000～2010年，陜甘寧退耕還林/草工程區(qū)耕地面積減少了9 501.12 km2，減少比例占工程區(qū)總面積的1.41%，其中2 168.80 km2和6 753.48 km2分別轉(zhuǎn)變?yōu)榱值睾筒莸?，分別占減少耕地面積的22.83%和71.08%，以退耕還草為主說明工程實(shí)施過程中考慮了干旱/半干旱區(qū)水熱條件限制。②2000～2010年，研究區(qū)地表植被覆蓋度增加了8.37個百分點(diǎn)，呈現(xiàn)較明顯的增長趨勢，其中耕地轉(zhuǎn)林地區(qū)域增加了14.63個百分點(diǎn)，耕地轉(zhuǎn)草地區(qū)域增加了18.06個百分點(diǎn)，達(dá)到甚至超過了工程區(qū)林草覆蓋率增加的規(guī)劃目標(biāo)。③2000～2010年，陜甘寧退耕還林/草工程區(qū)內(nèi)大于15°坡耕地面積減少了1 800.02 km2，其中15°～25°坡耕地轉(zhuǎn)換為林地和草地面積為408.68 km2和10 48.04 km2，大于25°坡耕地轉(zhuǎn)為林地和草地面積分別為149.83 km2和152.36 km2。

關(guān)鍵詞 陜甘寧地區(qū)；退耕還林/草；土地利用變化；植被覆蓋度；生態(tài)成效

中圖分類號 S725.7 文獻(xiàn)標(biāo)識碼 A 文章編號 0517-6611（2014）17-05524-05

Abstract Based on data of land use/coverage change， vegetation coverage， slope land distribution during 2000-2010， the conversion of cropland to forest or grassland condition， the vegetation coverage change were analyzed， so as to reveal ecological effectiveness of GrainforGreen project in ShanxiGansuNingxia region. The results are shown as follows.① From 2000 to 2010， farmland in the ShanxiGansuNingxia Region decreased by 9 501.12 km2， with the reducing proportion of 1.41% accounted for the total area of the study area. Converted from farmland to forest and grassland were 2 168.80 km2 and 6 753.48 km2， which accounted for 22.83% and 71.08% of reduced farmland， respectively.②From 2000 to 2010， the vegetation coverage was improved progressively， the average coverage of the whole study area increased by 8.37 percentage points. Among them， the coverage of regions with farmland to forest increased by 14.63 percentage points， and the coverage of regions with farmland to grassland increased by 18.06 percentage points.③In resent ten years， the area of farmland with slope greater than 15°decreased by 1 800.02 km2. In which， farmland with slops between 15°-25° slope， 408.68 km2 and 1 048.04 km2 were converted to forest and grassland respectively. The farmland with slops more than 25°， with 149.83 km2 and 152.36 km2 transferred to forest and grassland respectively.

Key words ShanxiGansuNingxia Region； Grain for Green Program； Land use change； Vegetation coverage； Ecological effectiveness

黃土高原是我國水土流失最嚴(yán)重的區(qū)域，超過60%的土地存在土壤質(zhì)量退化問題，年平均土壤流失量為2 000～2 500 t/km2，土質(zhì)山坡道路路面土壤侵蝕模數(shù)高達(dá)3.49萬～9.02萬t2/km2 [1-2]。地形破碎、土壤疏松、植被稀疏以及不合理的土地利用結(jié)構(gòu)是黃土高原地區(qū)水土流失嚴(yán)重和生態(tài)環(huán)境退化的主要原因[3]。隨著生態(tài)環(huán)境的日益惡化，土地資源的合理利用受到越來越多人們的重視[4]。從1999年開始，我國開始實(shí)施退耕還林還草工程，以此來減少水土流失和土壤侵蝕等生態(tài)環(huán)境問題[5-6]。陜甘寧地區(qū)的主體位于黃土高原中部丘陵溝壑區(qū)[7]，是我國退耕還林/草工程的典型區(qū)[8]。準(zhǔn)確掌握退耕還林/草對區(qū)域生態(tài)環(huán)境的影響，對退耕還林還/草工程的可持續(xù)發(fā)展具有極其重要的意義。

近年來，國內(nèi)外學(xué)者在退耕還林/草工程的生態(tài)、經(jīng)濟(jì)、社會效應(yīng)等方面開展了大量研究，得出土地利用、植被覆蓋狀況發(fā)生了較大變化，水土流失有減輕趨勢的結(jié)論[9-10]。許多研究利用遙感和GIS技術(shù)分析退耕還林/草區(qū)域的歸一化植被指數(shù)（NDVI）、植被覆蓋度、土地利用變化等指標(biāo)，結(jié)果表明退耕前后的地表植被覆蓋度呈現(xiàn)逐步上升的趨勢[11-15]。曹世雄[6]研究表明陜北地區(qū)地表植被覆蓋度由1998年的34.6%增至2005年的49.2%；崔曉臨[16]得出2000～2009年陜西省植被覆蓋增加區(qū)域面積占陜西省總面積的27.8%；宋富強(qiáng)[17]研究表明陜北地區(qū)7°～15°、15°～25°及25°～35°坡度植被明顯改善面積分別占總改善面積的39.91%、25.81%、2.28%；李雙雙[7]研究表明2000～2009 年陜甘寧地區(qū)植被覆蓋呈現(xiàn)明顯增加趨勢。

然而，退耕還林/草工程規(guī)劃的“陡坡耕地基本退耕還林，嚴(yán)重沙化耕地基本得到治理”是否實(shí)現(xiàn)，工程區(qū)林草覆蓋率是否增加了4.5%，工程治理地區(qū)的生態(tài)狀況是否得到較大改善等問題并沒有明確答案。因此，該研究選擇退耕還林/草工程典型區(qū)——陜甘寧3省作為研究區(qū)，通過土地利用變化、植被覆蓋度等數(shù)據(jù)，分析研究區(qū)耕地轉(zhuǎn)換為林地或草地狀況，不同坡度退耕情況以及退耕區(qū)植被覆蓋度變化等，以期評價陜甘寧地區(qū)退耕還林工程的生態(tài)效應(yīng)。

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

1.1 研究區(qū)概況

陜甘寧地區(qū)的主體位于黃土高原中部丘陵溝壑區(qū)，是黃河中上游生態(tài)環(huán)境重點(diǎn)建設(shè)地區(qū)，橫跨半濕潤、半干旱和干旱3個氣候帶，南北差異較大。研究區(qū)年均氣溫3.7～13.8 ℃，年平均降水量在173.9～814.5 mm之間[18]，自東南向西北減少，全年降水主要集中在夏季，汛期降水量占全年的55%～74%[19]，且大部分區(qū)域年降水量在400 mm以下，但潛在蒸散量在1 000 mm以上，土壤侵蝕模數(shù)為10 000～20 000 t/km2，局部地區(qū)達(dá)30 000 t/km2[20]，是我國水土流失最嚴(yán)重的地區(qū)。區(qū)域內(nèi)自然環(huán)境復(fù)雜，生態(tài)條件多樣，植物資源種類繁多，植被分布具有明顯的地帶性，從南到北依次分布北亞熱帶常綠闊葉林、暖溫帶落葉闊葉林、森林草原和溫帶草原。在研究區(qū)域內(nèi)植被有良好的攔截泥沙作用，當(dāng)植被度達(dá)到90%時可減少徑流90% 以上[21-23]。

在大空間尺度上，具體實(shí)施退耕還林的工程區(qū)邊界難以界定。根據(jù)已有研究的經(jīng)驗(yàn)[6，24-25]以及“全國生態(tài)環(huán)境十年變化遙感調(diào)查項(xiàng)目”需求，將納入退耕還林工程的縣作為工程區(qū)，因此該研究的工程區(qū)邊界為陜甘寧3省實(shí)施該工程的210個縣的縣界，包括陜西省104個縣、甘肅省 87個縣、寧夏自治區(qū)19個縣。

1.2 土地利用/覆被變化數(shù)據(jù)

利用“全國生態(tài)環(huán)境十年變化遙感調(diào)查項(xiàng)目”提供的2000、2005、2010年3期30 m空間分辨率的土地利用/覆蓋數(shù)據(jù)，將其重新分類為林地、草地、濕地、耕地、人工表面、其他等6類一級生態(tài)系統(tǒng)類型數(shù)據(jù)。該土地利用/覆蓋數(shù)據(jù)通過對陸地衛(wèi)星TM/ETM+影像進(jìn)行輻射傳輸校正和大氣校正，利用1 ∶5萬地形圖進(jìn)行控制點(diǎn)校正，利用各省提供的地面點(diǎn)位信息建立標(biāo)志庫，采用自動分類輔以人工干預(yù)進(jìn)行遙感解譯得到分類精度為86%。采用式（1）的轉(zhuǎn)移矩陣數(shù)學(xué)模型對研究區(qū)生態(tài)系統(tǒng)類型間的相互轉(zhuǎn)移進(jìn)行分析。

1.3 植被覆蓋度

植被覆蓋度的差異體現(xiàn)了植被動態(tài)變化過程的時間特征，也是區(qū)域生態(tài)系統(tǒng)環(huán)境變化的重要標(biāo)志[26]。植被覆蓋度年均值能夠反映植被的年際變化，根據(jù)像元二分模型理論，可以認(rèn)為1個像元的NDVI值是由綠色植被部分貢獻(xiàn)的信息與無植被覆蓋部分貢獻(xiàn)的信息組合而成。該文采用16 d最大值合成的2000～2010年逐年MODIS/NDVI時間序列數(shù)據(jù)，并通過下式估算植被覆蓋度：

Fc= NDVI-NDVIsoil NDVIveg-NDVIsoil （2）

式中，F(xiàn)c為1個像元內(nèi)的植被覆蓋度；NDVI為該像元上的歸一化植被指數(shù)；NDVIveg是純植被像元的NDVI值；NDVIsoil是完全無植被覆蓋像元的NDVI值。參考常規(guī)值并結(jié)合研究區(qū)特點(diǎn)，利用土地利用/覆蓋圖和NDVI圖進(jìn)行疊加分析，查看每個類型的NDVI值，最終將NDVIsoil表示全裸土的NDVI值取0.1，NDVIveg表示全植被的NDVI值取0.9。

1.4 坡耕地分布數(shù)據(jù)

收集研究區(qū)90 m分辨率的SRTMDEM數(shù)據(jù)，通過對90 mDEM進(jìn)行重采樣，得30 m分辨率的DEM。通過30 mDEM生成坡度圖，并對坡度重新分類。根據(jù)退耕還林政策規(guī)定，大于25°的陡坡耕地必須無條件地退耕，15°～25°的陡坡耕地有選擇地退耕。景可與鄭粉莉等[22]認(rèn)為黃土高原區(qū)大于15°陡坡耕地應(yīng)該全部退耕。因此，該文將坡度劃分為15°～25°和大于25°。2類坡度數(shù)據(jù)分別與2000、2010年的耕地分布數(shù)據(jù)進(jìn)行疊加分析，得到不同坡度范圍內(nèi)的耕地面積以及不同坡度范圍內(nèi)2000～2010年耕地變化狀況。

2 結(jié)果與分析

2.1 2000～2010年陜甘寧地區(qū)土地利用變化

由圖1可知，①2000～2010年，陜甘寧地區(qū)林地面積持續(xù)增加，凈增加面積為2 171.28 km2，增加比例占陜甘寧退耕還林/草工程區(qū)總面積674 669.44 km2的0.32%。其中2000～2005年凈增加了1 940.01 km2，2005～2010年凈增加了231.27 km2。②2000～2010年，草地面積持續(xù)增加，凈增加面積為6 483.37 km2，增加比例為0.96%。其中2000～2005年凈增加了5 189.32 km2，

2005～2010年凈增加了294.05 km2。③2000～2010年，濕地面積凈增加了204.71 km2，增加比例為0.03%。其中2000～2005年凈增加了144.15 km2，2005～2010年凈增加了60.56 km2。④2000～2010年，耕地面積持續(xù)減少，凈減少面積為9 501.11 km2，減少比例為1.41%。其中2000～2005年減少了7 508.7 km2，2005～2010年減少了1 992.41 km2。⑤2000～2010年人工表面凈增加了1 305.67 km2，增加比例為0.19%。其中2000～2005年凈增加了611.41 km2，2005～2010年凈增加了694.26 km2。⑥2000～2010年，其他面積減少了663.92 km2，減少比例為0.1%。其中2000～2005年凈減少了376.19 km2，2005～2010年凈減少了287.73 km2（表1）。

從圖3可以看出，陜甘寧地區(qū)10年平均植被覆蓋度分布不均勻，跨度比較大，從南往北植被覆蓋度逐漸降低。陜西省平均植被覆蓋度最高，甘肅省西北部平均植被覆蓋度相對較低，最低至10%。利用最小二乘法分析陜甘寧地區(qū)2000～2010年植被覆蓋度的變化趨勢，可以看出工程區(qū)10年間植被覆蓋總體呈現(xiàn)好轉(zhuǎn)態(tài)勢，特別是陜北地區(qū)的寶塔區(qū)、延長縣、延川縣、子長縣、宜川縣、安塞縣、吳堡縣和清澗縣，植被增加趨勢高達(dá)5%/a，部分區(qū)域植被覆蓋度呈現(xiàn)下降趨勢，如甘肅省的玉門市、金川區(qū)、永昌區(qū)、肅南縣，陜西省的未央?yún)^(qū)、雁塔區(qū)、高陵縣，以及寧夏自治區(qū)的金鳳區(qū)、永寧縣和賀蘭縣，減少的趨勢高達(dá)6%/a。

2.3 陜甘寧退耕還林/草區(qū)生態(tài)成效

2.3.1 陜甘寧退耕還林/草工程區(qū)退耕還林/草面積變化。

2000～2005年，陜甘寧退耕還林/草工程區(qū)內(nèi)退耕面積7 508.70 km2。其中，退耕還林的面積為1 924.94 km2，主要分布在陜西省、甘肅省的東部和南部以及寧夏南部，退耕還草面積5 424.39 km2，主要分布在陜西省的北部，甘肅省的東部和南部。2005～2010年，退耕面積1 992.41 k m2。其中，退耕還林的面積為256.53 km2，主要分布在陜西省、甘肅省的東部和南部，退耕還草面積1 389.33 km2，主要分布在寧夏自治區(qū)的南部。整個10年內(nèi)，退耕面積9 501.11 km2。其中，退耕還林的面積為2 168.8 km2，主要分布在陜西省、甘肅省的東部和南部。退耕還草面積6 753.48 km2，主要分布在陜西省的北部，甘肅省的東部和南部，寧夏自治區(qū)的西部和南部（圖4）。

2.3.2 退耕還林/草區(qū)域植被覆蓋度變化。

如圖5所示，2000～2010年，陜甘寧退耕還林/草工程區(qū)內(nèi)耕地轉(zhuǎn)換為林地區(qū)域的植被覆蓋度基本上呈現(xiàn)出逐年增長的趨勢，年平均植被覆蓋度由2000年的65.55%增至2010年的80.18%，高于整個區(qū)域的增速。其中，2000～2005年，退耕還林區(qū)域的平均植被覆蓋度為46.13%，2005～2010年平均植被覆蓋度為49.94%，說明退耕還林/草工程使植被覆蓋度得以快速增加。從各省情況來看，退耕還林區(qū)域

前、后5年平均植被覆蓋度增速最快的是陜西省，增加了8.2個百分點(diǎn)；其次是甘肅省，增加了6.99個百分點(diǎn)；最后是寧夏自治區(qū)，增加了1.16個百分點(diǎn)。

陜甘寧退耕還林/草工程區(qū)內(nèi)耕地轉(zhuǎn)換為草地區(qū)域的植被覆蓋度由2000年的42.93%增至2010年的61.00%。其中，2000～2005年，退耕還草區(qū)域的平均植被覆蓋度為49.84%，2005～2010年平均植被覆蓋度為

58.09%。從各省情況來看，退耕還草區(qū)域

前、后5年平均植被覆蓋度增速最快的是陜西省，增加了11.46個百分點(diǎn)；其次是甘肅省，增加了5.79個百分點(diǎn)；最后是寧夏自治區(qū)，增加了3.03個百分點(diǎn)。

2.3.3 不同坡度退耕還林/草狀況分析。

2000年，陜甘寧退耕還林/草工程區(qū)大于15°坡耕地面積17 873.52 km2，其中，坡度在15°～25°之間的坡耕地面積14 369.78 km2；大于25°的坡耕地面積3 503.74 km2。2010年，陜甘寧退耕還林/草工程區(qū)大于15°坡耕地面積減少至16 073.5 km2，其中，15°～25°之間的坡耕地面積12 879.55 km2，大于25°的坡耕地面積3 193.95 km2。

2000～2010年陜甘寧退耕還林/草工程區(qū)內(nèi)大于15°坡耕地減少了1 800.02 km2，坡耕地轉(zhuǎn)換為林地面積為558.5 km2，坡耕地轉(zhuǎn)換為草地面積1 200.4 km2。其中，15°～25°耕地減少面積為1 490.23 km2，轉(zhuǎn)換為林地和草地面積分別為408.68 km2、1 048.04 km2；大于25°耕地減少了309.79 km2，轉(zhuǎn)換為林地和草地面積為149.83 km2和152.36 km2（圖6）。

3 結(jié)論與討論

該文利用土地利用、植被覆蓋度等數(shù)據(jù)，分析陜甘寧退耕還林/草工程區(qū)土地利用變化、耕地轉(zhuǎn)換為林地或草地狀況，不同坡度退耕情況以及退耕區(qū)植被覆蓋度變化等，得出主要結(jié)論如下。

（1）2000～2010年，林地凈增加面積為2 171.28 km2，增加比例占工程區(qū)總面積的0.32%；草地凈增加面積為6 483.37 km2，增加比例為0.96%；濕地面積凈增加了204.71 km2，增加比例為0.03%；耕地面積持續(xù)減少，凈減少面積為9 501.11 km2，減少比例為1.41%；人工表面凈增加了1 305.67 km2，增加比例為0.19%；其他面積減少了663.92 km2，減少比例為0.1%。

（2）2000～2010陜甘寧退耕還林/草工程區(qū)植被覆蓋度呈明顯的增長趨勢，增長速率為每年0.837個百分點(diǎn)。從各省情況來看，前、后5年相比增速最快的是陜西。

（3）2000～2010年，陜甘寧退耕還林/草工程區(qū)耕地轉(zhuǎn)為林地和耕地轉(zhuǎn)為草地的面積分別為2 168.8、6 753.48 km2。從分省情況來看，前后5年相比，平均植被覆蓋度增速最快的均是陜西省；其次是甘肅省，最后是寧夏自治區(qū)。

（4）退耕還林/草工程規(guī)定25°以上坡耕地全部退耕，而在黃土區(qū)可以考慮將15°以上的坡耕地實(shí)施退耕[22]。2000～2010年，陜甘寧退耕還林/草工程區(qū)內(nèi)大于15°坡耕地減少了1 800.02 km2，15°～25°之間的坡耕地轉(zhuǎn)換為林地和草地面積為面積分別為408.68 km2和1 048.04 km2；大于25°的耕地轉(zhuǎn)換為林地和草地面積為149.83 km2和152.36 km2。

4 結(jié)語

退耕還林/草工程使得植被覆蓋度得以逐年提高，生態(tài)環(huán)境持續(xù)改善。但對于陜甘寧退耕還林/草工程區(qū)內(nèi)降水量在400 mm以下的地方，是適宜退耕還草、還林還是還灌木，與水分條件有直接關(guān)系。當(dāng)降雨量在400 mm以下，樹木難以成活，但水資源豐富的地區(qū)樹木也可成活[8]。干旱半干旱生態(tài)脆弱區(qū)大規(guī)模的人工造林可誘發(fā)嚴(yán)重的水資源短缺，導(dǎo)致人工林地植被覆蓋度減小以及對植物種類的多樣性存在負(fù)面影響[6，12]。生態(tài)環(huán)境的改善不僅受氣候條件的驅(qū)動，同時也受人為活動的驅(qū)動。不僅要考慮自然環(huán)境因素，還要考慮當(dāng)?shù)厝藗兊囊庠负托袨?，如果?dāng)?shù)厝嗣癫⒉恢С猪?xiàng)目的進(jìn)

行，那么該項(xiàng)目將難以長久存在。農(nóng)戶樂意接受退耕還林項(xiàng)目提供的糧食補(bǔ)償，但如果項(xiàng)目結(jié)束停止補(bǔ)助會再次墾荒種糧[27]。因此退耕還林/草項(xiàng)目要獲得持續(xù)的成功，當(dāng)?shù)剞r(nóng)戶的參與支持至關(guān)重要。

參考文獻(xiàn)

[1] 馬巖，陳利頂，虎陳霞.黃土高原地區(qū)退耕還林工程的農(nóng)戶響應(yīng)與影響因素——以甘肅定西大牛流域?yàn)槔齕J].地理科學(xué)，2008，28（1）：34-39.

[2] 曹世雄，陳莉，高旺盛.黃土丘陵區(qū)路面種草對水土保持的影響及成本效益分析[J].農(nóng)業(yè)工程學(xué)報(bào)，2006，22（8）：73-76.

[3] FU B J，CHEN L D.Agricultural landscape spatial pattern analysis in the semi-arid area of the Loess Plateau，china[J].China Journal of Arid Environment，2000，44（3）：291-303.

[4] 楊存建，劉紀(jì)遠(yuǎn)，張?jiān)鱿?，?遙感和GIS支持下的中國退耕還林還草決策分析[J].遙感學(xué)報(bào)，2002，6（3）：205-210.

[5] UCHIDA E，XU J T，ROZELLE S.Grain for green：cost-effectiveness and sustainability of China’s conservation set-aside program[J].Land Economics，2005，81（18）：247-264.

[6] 曹世雄，陳軍，陳莉，等.退耕還林項(xiàng)目對陜北地區(qū)自然與社會的影響[J].中國農(nóng)業(yè)科學(xué)，2007，40（5）：972-979.

[7] 李雙雙，延軍平，萬佳.近10 年陜甘寧黃土高原區(qū)植被覆蓋時空變化特征[J].地理學(xué)報(bào)，2012，67（7）：960-970.

[8] 劉紀(jì)遠(yuǎn)，張?jiān)鱿?，徐新良，?21 世紀(jì)初中國土地利用變化的空間格局與驅(qū)動力分析[J].地理學(xué)報(bào)，2009，64（12）：1411-1420.

[9] 汪邦穩(wěn)，楊勤科，劉志紅，等.延河流域退耕前后土壤侵蝕強(qiáng)度的變化[J].中國水土保持科學(xué)，2007，5（4）：27-33.

[10] 謝紅霞，李銳，楊勤科，等.退耕還林（草）和降雨變化對延河流域土壤侵蝕的影響[J].中國農(nóng)業(yè)科學(xué)，2009，42（2）：569-576.

[11] 王朗，傅伯杰，呂一河，等.生態(tài)恢復(fù)背景下陜北地區(qū)植被覆蓋的時空變化[J].應(yīng)用生態(tài)學(xué)報(bào)，2010，21（8）：2109- 2116.

[12] 杜靈通，李國旗.基于SPOTVGT 的寧夏鹽池縣近8 年生態(tài)環(huán)境動態(tài)監(jiān)測[J].北京林業(yè)大學(xué)，2008，30（5）：46-51.

[13] 王情，劉雪華，呂寶磊.基于SPOT-VGT數(shù)據(jù)的流域植被覆蓋動態(tài)變化及空間格局特征——以淮河流域?yàn)槔齕J].地理科學(xué)進(jìn)展，2012，32（2）：270-277.

[14] 楊亞娟，劉政鴻.陜北地區(qū)退耕還林十年植被恢復(fù)變化分析——以子長縣為例[J].中國水土保持，2013（2）：54-56.

[15] 孫智輝，雷延鵬，卓靜，等.延安北部丘陵溝壑區(qū)退耕還林（草）成效的遙感監(jiān)測[J].生態(tài)學(xué)報(bào)，2010，30（23）：6555-6562.

[16] 崔曉臨，趙娟，白紅英，等.基于MODIS NDVI的陜西省植被覆蓋變化分析[J].西北農(nóng)林科技大學(xué)學(xué)報(bào)：自然科學(xué)版，2013，41（5）：87-93.

[17] 宋富強(qiáng)，邢開雄，劉陽，等.基于MODIS/NDVI的陜北地區(qū)植被動態(tài)監(jiān)測與評價[J].生態(tài)學(xué)報(bào)，2011，31（2）：354-363.

[18] 馬玉平.陜甘寧地區(qū)退耕還林還草的氣候和農(nóng)業(yè)效應(yīng)模擬研究[D].北京：中國氣象科學(xué)研究院，2011.

[19] LI S S，YAN J P，LIU X Y，et al.Response of vegetation restoration to climate change and human activities in Shaanxi-Gansu-Ningxia Region[J].J Geogr Sci，2013，23（1）：98-112.

[20] 李世東.干熱干旱河谷區(qū)和黃土丘陵溝壑區(qū)退耕還林還草模式初步研究[J].北京林業(yè)大學(xué)學(xué)報(bào)，2002，24（3）：35-38.

[21] SUN W Y，SHAO Q Q，LIU J Y.Soil erosion and its response to the changes of precipitation and vegetation cover on the Loess Plateau[J].J Geogr Sci，2013，23（6）：1091-1106.

[22] 景可，鄭粉莉.黃土高原植被建設(shè)的經(jīng)驗(yàn)教訓(xùn)與前景分析[J].水土保持研究，2004，11（4）：25-27.

[23] 許炯心.降水-植被耦合關(guān)系及其對黃土高原侵蝕的影響[J].地理學(xué)報(bào)，2006，61（1）：58-65.

[24] LIU J G，LI S X，OUYANG Z Y，et al.Ecological and socioeconomic effects of China’spolicies for ecosystem services[J].PNAS，2008，105（28）：9477-9482.[25] CAO S X，SUN G，ZHANG Z Q，et al.Greening China naturally[J].AMBIO，2011，40：828-831.

[26] 陳學(xué)兄，張小軍，陳永貴，等.陜西省1998～2008年植被覆蓋度的時空變化研究[J].武漢大學(xué)學(xué)報(bào)：信息科學(xué)版，2013，38（6）：674-678.

[27] 曹世雄，陳莉，余新曉.陜北農(nóng)民對退耕還林的意愿評價[J].應(yīng)用生態(tài)學(xué)報(bào)，2009，20（2）：426-434.