楊 柳,閻建忠,王 盼,王 宏

西南大學(xué)資源環(huán)境學(xué)院, 重慶 400715

青藏高原是全球氣候變化研究熱點區(qū)域[1- 2],是全球氣候變化的啟動區(qū)和敏感區(qū)[3- 4]。該地區(qū)在1980年后發(fā)生了暖突變[5],持續(xù)變暖趨勢明顯,1971—2011年期間達到了0.039℃/a,增溫幅度顯著高于我國其他地區(qū),且增溫以冬季最為明顯,主要增溫區(qū)位于高海拔地區(qū)[6- 9]。該地區(qū)降水量的變化呈波動式增加,增長率為0.67 mm/a,且1986年以后降水增加更為顯著[10- 11]。同時,降水量空間分布差異十分顯著且季節(jié)分配不均[12]。隨著氣溫、降水的變化,該地區(qū)在20世紀(jì)80年代中后期氣候條件由暖干型向暖濕型突變[13- 14]。此外,該地區(qū)極端高溫事件(熱晝、熱夜和持續(xù)暖期)頻次呈增加態(tài)勢,極端低溫事件(冷晝、冷夜、持續(xù)冷期、霜日和冰凍日)頻次顯著降低,且下降率高于增長率[2,15- 16]。

以變暖為主的氣候變化使當(dāng)?shù)責(zé)崃抠Y源得以改善,從而使青藏高原適宜耕種的土地面積增加。張戈麗等運用種植適宜性模型,發(fā)現(xiàn)西藏適宜耕種面積在1970—2000年間增加了870 km2,適宜耕種土地向原有耕地上限擴張,擴張區(qū)域主要分布在農(nóng)牧交錯帶上[17]。利用遙感影像和GIS工具,一些研究發(fā)現(xiàn)近30多年來青藏高原耕地面積擴張明顯,且農(nóng)牧交錯帶是耕地增加的主要區(qū)域。楊春艷等[18]利用4期遙感影像數(shù)據(jù),采用土地利用動態(tài)度、土地利用區(qū)位指數(shù)和土地區(qū)域差異模型,對1980—2010年西藏耕地面積時空變化進行研究,結(jié)果表明,耕地面積總體呈增加態(tài)勢,相對增加較快的是昌都地區(qū),而昌都地區(qū)是西藏的主要牧區(qū)之一。擺萬奇等[19]運用同樣的方法研究了1976—2011年西藏拉薩河流域耕地的時空變化,結(jié)果也表明耕地增加主要來自草地,且耕地擴張不斷向更高海拔和更大坡度的范圍推進。同樣地,Shen等[20]研究了西藏南部地區(qū)的作物氣候生產(chǎn)力和耕地分布的時空變化,結(jié)果表明氣候變化使該地區(qū)物候作物氣候生產(chǎn)力提高,耕地面積增加且向高海拔地區(qū)和更北的地方擴張；Wu等[21]通過對天?？h的研究也可以發(fā)現(xiàn)該地區(qū)耕地以每年2.34%的速率遞增,且增加的耕地主要來源于草地。因此,青藏高原的農(nóng)牧民有極大的可能性采取開墾策略以適應(yīng)氣候變化。

然而,農(nóng)牧民也有極大的可能性不選擇開墾,主要的原因在于：第一,氣候變化帶來的負面影響會加劇青藏高原農(nóng)業(yè)的脆弱性,導(dǎo)致作物單產(chǎn)下降[22]。主要體現(xiàn)在：極端高溫事件頻發(fā),導(dǎo)致鼠蟲草害等自然災(zāi)害和干旱發(fā)生頻率增加；高溫不僅增加了作物病蟲害的范圍和強度[23- 24],還使得作物灌漿期縮短,農(nóng)作物抗寒能力減弱[25]；氣溫升高加快了土壤有機質(zhì)的分解,促使土壤肥力下降[26]。第二,隨著城市化進程的加快,青藏高原同時面臨勞動力析出和勞動力價格上漲的問題,農(nóng)牧民更愿意追求生計多樣化,從而導(dǎo)致耕地撂荒現(xiàn)象明顯[27- 28]。第三,在青藏高原農(nóng)區(qū)耕地大多處于飽和狀態(tài),無地可墾[29]。第四,開墾耕地需要較高的費用(如追加灌溉、平整土地等投資)。農(nóng)牧民的行為是風(fēng)險規(guī)避型[30],為降低風(fēng)險,他們往往選擇集約化利用耕地,而不是開墾。

青藏高原是生態(tài)脆弱區(qū),農(nóng)牧交錯帶草地的墾殖會導(dǎo)致草地退化和沙漠化,并且不可逆轉(zhuǎn)。研究氣候變化背景下青藏高原農(nóng)牧民是否會采取開墾耕地的適應(yīng)策略,不僅僅涉及當(dāng)?shù)剞r(nóng)牧民的生計和氣候變化適應(yīng)策略,也和區(qū)域的可持續(xù)發(fā)展息息相關(guān)。對于青藏高原氣候變化背景下的墾殖,目前雖然有少量區(qū)域尺度的模型模擬和遙感調(diào)查,但一直缺少對農(nóng)牧民是否開墾以及開墾原因的實證研究,牧民如何平衡氣候變化帶來的開墾機遇、風(fēng)險以及人口壓力,仍然缺少認識。本文擬以青藏高原的“一江兩河”地區(qū)、河湟谷地以及川西地區(qū)的壤塘縣為研究區(qū),調(diào)查氣候變化和勞動力析出背景下農(nóng)牧民的耕地開墾情況,并定量研究氣候變化對農(nóng)牧民開墾耕地的影響。論文由以下幾個部分組成,第二部分介紹了研究區(qū)的情況和調(diào)查方法；第三部分是模型選擇與變量設(shè)定；第四部分是研究結(jié)果；最后是討論和結(jié)論。

1 研究區(qū)概況與數(shù)據(jù)來源

1.1 研究區(qū)概況

青藏高原(73°19′—104°47′E,26°00′—39°47′N)幅員面積約250×104km2,占全國總面積的26.8%,平均海拔4383 m。該地區(qū)輻射強,日照時間長,氣溫低且日較差大,積溫少；干濕分明,多夜雨；冬季干冷漫長,大風(fēng)多；夏季溫涼多雨,冰雹多。青藏高原產(chǎn)糧區(qū)主要包括西藏自治區(qū)、青海省和川西地區(qū),因此本文選擇這3個地區(qū)的典型區(qū)域作為研究區(qū)(圖1)。

圖1 研究區(qū)圖及調(diào)研路線Fig.1 Location of the study area & map of inspection line

“一江兩河”地區(qū)(87°00′—92°35′E,28°20′—31°20′N)是指西藏自治區(qū)雅魯藏布江中游及其支流拉薩河和年楚河中下游的河谷地區(qū),位于藏南湖盆谷地,幅員面積逾6.57×104km2(占西藏自治區(qū)總面積的5.47%),海拔多在2700—4200 m之間。氣候?qū)儆诟咴瓬貛Ъ撅L(fēng)半干旱氣候,年均氣溫在4.7—8.3℃之間,全年日照時數(shù)為2800—3300 h,年降水量在251.7—580.0 mm之間,蒸發(fā)量大,約為降水量的4—5倍。該地區(qū)下轄拉薩、山南、日喀則3個地市的18個縣(市、區(qū)),農(nóng)牧業(yè)發(fā)達,農(nóng)作物以青稞、小麥、油菜為主,牲畜包括牦牛、犏牛、奶牛、馬和綿羊。據(jù)2017西藏統(tǒng)計年鑒可知,2016年,該地區(qū)耕地面積達16.17×104hm2,占西藏總耕地面積的67.92%；農(nóng)作物播種面積為16.47×104hm2,占西藏農(nóng)作物播種總面積的63.85%；農(nóng)業(yè)產(chǎn)值為81.92×104萬元,占西藏農(nóng)業(yè)總產(chǎn)值的50.43%；糧食產(chǎn)量為93.78×104t,約占西藏糧食產(chǎn)量的72.07%。

河湟谷地(100°—103°E,35°—38°N)地處青藏高原東北緣、青海省的東北部,位于黃河及其支流湟水河河谷地區(qū)域,幅員面積3.5×104km2(占青海省總面積的4.85%),海拔在1689—5218 m之間。該區(qū)域?qū)儆诘湫偷母珊?、半干旱大陸性氣?年日照時數(shù)在2600—3000 h之間,年平均溫度在5—9℃之間,年平均降水量在252—535 mm之間。該地區(qū)下轄西寧市、海東市、海北藏族自治州、黃南藏族自治州和海南藏族自治州的17個區(qū)縣,農(nóng)牧業(yè)較發(fā)達,農(nóng)作物以小麥、青稞、玉米、油菜和枸杞為主,牲畜主要包括綿羊、山羊和豬。據(jù)2017青海統(tǒng)計年鑒可知,2016年,該地區(qū)耕地面積約為41.82×104hm2,約占青海省總耕地面積的71.07%；糧食產(chǎn)量約為82.74×104t,約占青海省糧食總產(chǎn)量的79.98%；農(nóng)業(yè)總產(chǎn)值約為116.71億元,約占青海省農(nóng)業(yè)總產(chǎn)值的51.94%。

壤塘縣(100°30′57″—101°29′16″E,31°28′14″—32°41′38″N)是川西北牧區(qū)的重要組成部分,地處青藏高原東南邊緣、川西北高原的西北部,位于大渡河上游。該地區(qū)幅員面積6606 km2,平均海拔3280 m以上。該地區(qū)氣候?qū)儆诟吆箨懶约撅L(fēng)氣候,年均氣溫在2.2—7.5℃之間,年日照數(shù)2607.4—3036 h,年蒸發(fā)量1132.4 mm,年降水量在666.9—790 mm之間,雨熱同期,干濕季節(jié)分明；無霜期短。該地區(qū)耕地18.82 km2,占全縣幅員面積的0.28%；天然草原面積4189 km2,占全縣幅員面積的63.41%。該地區(qū)養(yǎng)殖的牲畜主要包括牦牛、藏羊、馬、黃牛和藏豬,農(nóng)作物主要以青稞、小麥、油菜為主。據(jù)2016年阿壩州年鑒可知,2015年,農(nóng)作物實際播種面積18.42 km2,其中糧食作物播面積為15.79 km2,總產(chǎn)量3114 t；經(jīng)濟作物播面積1.59 km2,總產(chǎn)量264 t；農(nóng)村經(jīng)濟總收入26891.51萬元,農(nóng)牧民人均純收入6413元。

1.2 數(shù)據(jù)來源

本文所用數(shù)據(jù)來自于三輪調(diào)查。課題組分別于2015年7月—8月、2016年7月—8月和2017年7月—8月對“一江兩河”地區(qū)、壤塘縣和河湟谷地進行入戶調(diào)查,采用分層隨機抽樣的方法。課題組首先在每個調(diào)研區(qū)進行為期7天的預(yù)調(diào)研,結(jié)合農(nóng)牧民回答情況,修改完善問卷,然后開展正式調(diào)研。首先,與當(dāng)?shù)剞r(nóng)牧局進行座談,結(jié)合當(dāng)?shù)厣鐣?jīng)濟、交通及農(nóng)業(yè)狀況確定樣本縣和樣本鄉(xiāng)。其次,通過與負責(zé)農(nóng)牧業(yè)的鄉(xiāng)長商談,在每個鄉(xiāng)的農(nóng)牧交錯帶上隨機抽取了3—5個半農(nóng)半牧村。入戶問卷調(diào)查分為兩步。第一步,與村長等村干部及相關(guān)知情人士進行半結(jié)構(gòu)式訪談,充分掌握該村的主要情況,主要包括氣候變化情況、自然災(zāi)害、耕地變化情況、作物種植情況、牲畜養(yǎng)殖情況等。第二步,采用參與式農(nóng)村調(diào)查評估法中的半結(jié)構(gòu)式訪談進行農(nóng)牧民問卷調(diào)查。在一江兩河地區(qū)和壤塘縣,由于調(diào)研區(qū)域幅員面積較大,農(nóng)牧戶主要居住在帳篷中,居住較分散,訪談難度較大,因此,以隨機抽樣的方式在每個村調(diào)查10戶左右；在河湟谷地,由于基礎(chǔ)建設(shè)較為完善,農(nóng)牧民分布較為集中,因此在每個村調(diào)查20戶左右。為了避免調(diào)查過程中的失誤,在預(yù)調(diào)研前對每個調(diào)查員進行了培訓(xùn)。如,在問及農(nóng)牧民感知氣候情況時,調(diào)查員給予多項選項(升高、降低、不變)供農(nóng)牧民選擇；由于青藏高原在20世紀(jì)80年代后期發(fā)生暖突變,因此在問及農(nóng)牧民開墾情況的時間點均在1980年以后。為確保問卷質(zhì)量,調(diào)查對象多為戶主,且年齡多處于40—50之間。每戶問卷調(diào)查時間控制在1—2 h之內(nèi)(主要為了避免農(nóng)牧民反感,時間過長農(nóng)牧民不愿配合),同時聘請3名藏族學(xué)生作為語言翻譯。

問卷調(diào)查內(nèi)容主要包括：人力資本(家庭人口、受教育水平、年齡等)、自然資本(經(jīng)營的耕地塊數(shù)等)、物質(zhì)資本(牲畜數(shù)量等)、金融資本(政府補貼、采集蟲草收入、非農(nóng)收入等)、社會資本(住宅距集市距離、可否獲得借貸機會等)、農(nóng)牧民感知氣候情況(主要是近10—30年的情況)和農(nóng)牧民開墾情況等。野外調(diào)查在問及農(nóng)牧民開墾原因時除去了與氣候變化相關(guān)的選項(如氣候變化會使適宜耕種面積擴大等具有導(dǎo)向型的選項),這是因為農(nóng)牧民的受教育程度普遍相對較低,對氣候變化產(chǎn)生的相應(yīng)結(jié)果認識較為薄弱。經(jīng)清理,共獲得861份有效問卷,考慮到計量分析的必要,且河湟谷地同村農(nóng)戶相似性較大,故每村隨機選取10戶作為研究對象。最終,共獲得605份問卷,其中一江兩河地區(qū)188份,壤塘縣169份,河湟谷地248份。

2 模型選擇與變量設(shè)定

2.1 模型的選擇

本文旨在基于可持續(xù)生計框架研究氣候變化背景下農(nóng)牧民的開墾決策。由于因變量為離散的二值選擇型變量,即開墾或不開墾,故采用Probit模型進行計量分析?；诖?本文將計量模型設(shè)定如下：

P(Yi=1)=Pi=Φ{α+β(Xi)+γ(Zi)+ε}

(1)

式中,Y是耕地利用方式/農(nóng)牧民適應(yīng)策略指標(biāo)；X為農(nóng)牧民感知氣候變化指標(biāo)；Z為控制變量矩陣,即各項生計資本指標(biāo),包括自然資本、人力資本、金融資本、物質(zhì)資本和社會資本,ε是具有正態(tài)分布特征的隨機擾動項。

考慮到估計結(jié)果的穩(wěn)健性,本文應(yīng)用Robust進行校正；由于青藏高原具有復(fù)雜多樣的地形條件,且經(jīng)濟條件也有明顯的差異,具體表現(xiàn)在地區(qū)之間的差異,故以地區(qū)(即“一江兩河”地區(qū)、壤塘縣和河湟谷地)為單位來解決空間自相關(guān)問題,并用Cluster來進行校正。運用軟件STATA 12.0進行統(tǒng)計、計量分析。

2.2 變量的設(shè)定

2.2.1解釋變量的選取

氣候變化通過對生計資本的直接影響,進而影響農(nóng)牧民所采取的適應(yīng)策略[31],故本文以五大類生計資本(人力資本、金融資本、社會資本、自然資本和物質(zhì)資本)為控制變量進行研究。

雖然農(nóng)牧民對氣候的感知與實測數(shù)據(jù)不一定一致[32- 33],但諸多研究表明農(nóng)牧民氣候感知是影響其響應(yīng)行為的關(guān)鍵因素[34]。當(dāng)農(nóng)牧民感知氣候變化所帶來的風(fēng)險和適應(yīng)能力較低時,農(nóng)牧民參與適應(yīng)行為的可能性降低；相反,當(dāng)農(nóng)牧民感知氣候變化強烈時,農(nóng)牧民便會采取適應(yīng)策略[35]。本文將農(nóng)牧民氣候感知分為三類：季節(jié)性感知、氣象感知、風(fēng)險感知,并從這三個角度來選取了農(nóng)牧民感知氣候指標(biāo)。

除了氣候變化的影響,在日益增長的人口壓力影響下,青藏高原的農(nóng)牧民不得不開墾以滿足對飼草的需求[29,36- 37]。氣候變暖導(dǎo)致灌木叢入侵高山草甸,影響了牧草的生長[38]。且在放牧規(guī)模不斷擴大,國家禁牧、休牧、輪牧政策的影響下,致使牧草供給十分短缺[39]。同時,陶娟平等[40]認為人口數(shù)量對土地開墾起著促進作用,人口增加會擴大土地開墾范圍。故選取家庭總?cè)丝诤蜕髷?shù)量作為人口壓力的指標(biāo)。

青藏高原越來越多的農(nóng)民接受和從事非農(nóng)業(yè)活動[41],這就導(dǎo)致在農(nóng)忙季節(jié)的勞動力短缺。而且,我國農(nóng)村轉(zhuǎn)移的勞動力以受教育程度較高的青壯年為主[42- 43],這就導(dǎo)致參與農(nóng)作的勞動力呈現(xiàn)老齡化和女性化的趨勢[44- 45],勞動力質(zhì)量下降[46],開墾的可能性就會降低。故選取了人均打工月收入為農(nóng)戶家庭參與非農(nóng)工作的指標(biāo)。

采集蟲草作為青藏高原特有的行業(yè),需要考慮采集蟲草與農(nóng)作之間也存在勞動力競爭的關(guān)系,因此,選取了指標(biāo)采集蟲草收入指標(biāo)。具體地,各研究變量的描述性統(tǒng)計見表1。

表1 研究變量描述性統(tǒng)計(n=605)

2.2.2解釋變量間共線性檢驗

為防止解釋變量間存在多重共線性問題,本文利用Pearson相關(guān)系數(shù)、容忍度(Tolerance)和方差膨脹因子(VIF)來進行檢驗。根據(jù)Pearson相關(guān)性分析結(jié)果顯示,居住地距集鎮(zhèn)距離和政府農(nóng)業(yè)生產(chǎn)性補貼之間的相關(guān)系數(shù)絕對值最高,為0.3505<0.8；政府農(nóng)業(yè)生產(chǎn)性補貼的容忍度最低,為0.7105>0.1,而其方差膨脹因子最高,為1.41<10。綜合以上三類檢驗指標(biāo),說明解釋變量之間不存在多重共線性問題,不會對模型的分析造成影響。

3 結(jié)果分析

3.1 農(nóng)牧民感知氣候變化

除了降水情況和極端低溫天氣的感知,研究區(qū)樣本農(nóng)牧民感知情況與氣象站數(shù)據(jù)基本相符(表2、圖2、圖3)。農(nóng)牧民對極端高溫天氣感知最為準(zhǔn)確,干旱增多,而對于極端低溫天氣感知不明顯；對季節(jié)性變化感知較為明顯,夏季時間延長,冬季時間縮短；對氣溫變化的感知能力強于對降水變化的感知,94.31%的農(nóng)牧民感知氣溫上升,僅有21.60%的農(nóng)牧民感知降水增多,這一結(jié)果與趙雪雁和侯向陽等的研究結(jié)果一致[47- 48]。絕大多數(shù)農(nóng)牧民感知氣候條件仍是“暖干型”氣候?？赡艿脑蛟谟谝韵聝蓚€方面：一方面,雖然氣象數(shù)據(jù)顯示每年降雨量增加,但當(dāng)?shù)亟涤觐l次降低、降雨強度增加且多為夜雨,因此農(nóng)牧民感知對降水變化感知不明顯；另一方面,農(nóng)戶的氣候感知受其記憶時段和參考基準(zhǔn)影響,通過實際的降雨量數(shù)據(jù)我們可以發(fā)現(xiàn)每年的平均降雨量不穩(wěn)定,呈波動式增加態(tài)勢,且波動幅度較大。

表2 研究區(qū)樣本農(nóng)牧民感知氣候變化情況/戶

圖2 1970—2013年青藏高原溫度變化圖Fig.2 Temperature change of the Tibetan Plateau during 1970—2013

圖3 1970—2013年青藏高原降水量變化圖Fig.3 Precipitation change of the Tibetan Plateau during 1970—2013

3.2 農(nóng)牧民實際開墾情況

研究區(qū)樣本農(nóng)牧民開墾行為不明顯,大多開墾面積較小。調(diào)研區(qū)域?qū)嶋H開墾的農(nóng)牧民共有120戶(占總樣本數(shù)的19.80%),開墾總面積為19.46 hm2。其中,開墾面積在0.03—0.13 hm2之間的農(nóng)牧民有79戶,開墾面積大于0.13 hm2的有32戶,開墾面積小于0.03 hm2的有9戶。

壤塘縣樣本農(nóng)牧民開墾行為最明顯,河湟谷地戶均開墾面積最大(表3)。壤塘縣的開墾戶數(shù)最多且開墾面積最大,戶均開墾面積為0.18 hm2。壤塘縣地處農(nóng)牧交錯帶,為了滿足其人口不斷地增長與牲畜數(shù)量增加而造成的糧食需求,農(nóng)牧民會開墾較大的面積。河湟谷地雖然開墾戶數(shù)最少且開墾面積最小,但戶均開墾面積最大,為0.19 hm2。

表3 研究區(qū)樣本農(nóng)牧民實際開墾情況

3.3 農(nóng)牧民開墾與人口壓力

樣本農(nóng)牧民大多在2005—2016年發(fā)生開墾行為,占總開墾農(nóng)戶的60.83%。壤塘縣在2010—2016年間開墾戶數(shù)最多,高達44戶,這可能與壤塘縣2010—2011年間人口急劇增多有關(guān)(表4、圖4)。同樣地,在2010年之后河湟谷地的人口開始下降,也就幾乎沒有開墾行為了。一江兩河地區(qū)樣本農(nóng)牧民主要集中在1995—2009年開墾,在這幾年該地區(qū)的人口也處于直線上升狀態(tài)。

表4 研究區(qū)樣本農(nóng)牧民開墾年份情況

圖4 1995—2016年青藏高原農(nóng)村人口變化圖Fig.4 The change of rural population of the Tibetan Plateau in 1995—2016

研究區(qū)樣本農(nóng)牧民開墾是人口壓力(主要表現(xiàn)在家庭人口和牲畜數(shù)量的增加)導(dǎo)致的結(jié)果。問卷調(diào)查表明有88戶(占總開墾農(nóng)戶的73.33%,表5)是因為人口壓力而開墾(主要是為了自食和養(yǎng)殖牲畜)。同時,農(nóng)牧民開墾地主要種植燕麥和青稞兩種作物(表6),這也足以說明開墾的原因是人口壓力。在壤塘縣,有54戶開墾農(nóng)牧民種植了飼料性作物,而在一江兩河地區(qū),開墾農(nóng)牧民主要種植糧食作物。

表5 研究區(qū)樣本農(nóng)牧民因人口壓力而開墾的情況

表6 研究區(qū)樣本農(nóng)牧民開墾地種植作物情況/戶

3.4 農(nóng)牧民打工及其年收入情況

由于勞動力在農(nóng)業(yè)和非農(nóng)活動之間存在競爭關(guān)系,開墾農(nóng)戶的打工比例明顯低于未開墾農(nóng)戶(表7)。從區(qū)域的差異來看,壤塘縣的打工比例和打工年收入都是最低。這主要是因為壤塘縣的農(nóng)牧民的文化程度普遍偏低,多為文盲或者小學(xué)學(xué)歷(占總?cè)丝跀?shù)的92.90%),且大多數(shù)人不會說普通話,只會說藏語,這就導(dǎo)致了當(dāng)?shù)氐霓r(nóng)牧民很難非農(nóng)就業(yè)。即使能夠非農(nóng)就業(yè),其工資也處于較低水平的狀態(tài)。此外,他們?nèi)匀槐в兄灰酗埑杂幸麓?就不愿意到其他地方去打工,或冒風(fēng)險進入其他非農(nóng)行業(yè)。該區(qū)域幅員面積較大,每個農(nóng)牧民都擁有自己的草場,農(nóng)牧民即使想從事非農(nóng)行業(yè),但距集鎮(zhèn)距離較遠,交通不便利；且由于非農(nóng)勞動力市場信息不對稱,農(nóng)牧民無法獲得非農(nóng)就業(yè)信息,從而導(dǎo)致農(nóng)牧民非農(nóng)工作的可獲得性低。相反,河湟谷地的打工比例在50%左右,打工年收入也是最高的。這主要是因為,河湟谷地第二三產(chǎn)業(yè)較其他兩個地區(qū)要發(fā)達得多,農(nóng)牧民的受教育水平較高、普通話好。

表7 研究區(qū)樣本農(nóng)牧民打工及其年收入情況

3.5 計量結(jié)果分析

本文使用STATA 12.0統(tǒng)計軟件對影響農(nóng)牧民開墾的因素進行回歸分析,模型1僅納入了農(nóng)牧民感知氣候變化情況,模型2僅納入了家庭人均打工年收入情況,模型3僅納入了人口壓力相應(yīng)指標(biāo)(即家庭總?cè)丝诩吧髷?shù)量),計量結(jié)果如表8所示。

農(nóng)牧民感知冬季持續(xù)時間的整體變化趨勢與農(nóng)牧民開墾行為呈顯著正相關(guān)關(guān)系。當(dāng)農(nóng)牧民感知冬季持續(xù)時間越長時,他們會存儲更多的糧食和飼料以滿足越冬的需求,因此農(nóng)牧民會傾向于開墾土地。

農(nóng)牧民感知降水的整體變化趨勢與農(nóng)牧民開墾行為呈顯著負相關(guān)關(guān)系?？赡艿慕忉屖怯捎谇嗖馗咴^缺乏灌溉設(shè)施,當(dāng)農(nóng)牧民感知降水減少,農(nóng)作物的產(chǎn)量就會降低,農(nóng)牧民就不得不開墾以保障糧食產(chǎn)量。

人均打工年收入與農(nóng)牧民開墾行為呈顯著負相關(guān)關(guān)系?？赡艿慕忉尀椋旱谝?有非農(nóng)業(yè)收入的農(nóng)牧民不需要從開墾更多的耕地來獲得更多的收入。第二,開墾需要較多的勞動力,由于家庭勞動力從事非農(nóng)工作,從而限制了他們對農(nóng)業(yè)(尤其是開墾)的投入。第三,由于青藏高原不存在農(nóng)業(yè)勞動力市場,農(nóng)牧民往往在農(nóng)忙季節(jié)進行勞動力交換,因而不可能雇傭農(nóng)業(yè)工人。

家庭總?cè)丝谂c農(nóng)牧民開墾行為呈顯著正相關(guān)關(guān)系。家庭總?cè)丝跀?shù)量越多,對糧食的需求也就越大,因此農(nóng)牧民傾向于開墾土地以獲得足夠多的糧食。

戶主受教育水平與農(nóng)牧民開墾行為呈顯著負相關(guān)關(guān)系。這主要是因為受教育水平越高的農(nóng)牧民有更大的機會去從事非農(nóng)工作,從而會有更高的收入。

撫養(yǎng)比與農(nóng)牧民開墾行為呈顯著正相關(guān)關(guān)系。當(dāng)家庭中需要撫養(yǎng)的人員比重越大,為了滿足非勞動力人員的口糧需求,就增加了開墾耕地的概率。

是否借貸與農(nóng)牧民開墾行為呈顯著負相關(guān)關(guān)系。農(nóng)牧民借錢大多是為了治病。當(dāng)農(nóng)牧民發(fā)生借貸時,就表明他們的健康狀況較差,這就造成了農(nóng)牧民家庭的勞動力短缺。開墾耕地需要大量的勞動力來實施并管理,勞動力的短缺就使得他們無法開墾更多的耕地。

居住地距集鎮(zhèn)距離與農(nóng)牧民開墾行為呈顯著負相關(guān)關(guān)系。主要原因在于距離集鎮(zhèn)較遠的地方其海拔相對較高,適宜開墾的耕地較少,因此農(nóng)牧民開墾的可能性較小。

表8 影響農(nóng)牧民開墾的因素估計結(jié)果

***、**、*分別表示估計系數(shù)在1%、5%、10%的統(tǒng)計水平上顯著；()中的數(shù)值代表Robust標(biāo)準(zhǔn)誤差

4 討論

以變暖為主的氣候變化使青藏高原適宜耕種的土地面積增加,已有的研究采用氣候模型和遙感調(diào)查,發(fā)現(xiàn)農(nóng)牧交錯帶是耕地增加的主要區(qū)域[17- 21]。本文利用青藏高原三個主要產(chǎn)糧區(qū)的農(nóng)戶調(diào)查數(shù)據(jù),證實了氣候變化驅(qū)動了耕地的開墾。受氣候變化影響,在有地可墾且非農(nóng)化程度不高的地區(qū),農(nóng)牧民開墾才被作為一種主要的適應(yīng)策略[29]。在尼日爾,受頻繁干旱和強風(fēng)的影響,土壤肥力低下,但由于長期高溫且有可開墾的土地,農(nóng)戶往往擴大種植面積來彌補作物產(chǎn)量下降[49]。氣候變暖,為我國北方地區(qū)創(chuàng)造了開墾的有利條件,Dong等[50]發(fā)現(xiàn)氣候變化影響著耕地的變化,尤其是在內(nèi)蒙古東部和東北地區(qū)。在內(nèi)蒙古,由于草場退化,農(nóng)牧交錯帶上的農(nóng)牧民會開墾天然草地[27,39]。

在研究區(qū),人口壓力也是開墾耕地的重要原因。氣候變暖使當(dāng)?shù)氐哪敛萃嘶F(xiàn)象明顯,加之受當(dāng)?shù)卣琳叩挠绊?農(nóng)牧民擁有的草場的面積較以往少得多。較大的家庭規(guī)模和牲畜數(shù)量增加了當(dāng)?shù)剞r(nóng)牧民人口壓力,農(nóng)牧民不得不采取生計多樣化或開墾措施。然而,當(dāng)?shù)剞r(nóng)牧民受教育程度較低,且語言不通,往往阻礙了其生計多樣化。同時,雖然調(diào)研數(shù)據(jù)表明農(nóng)牧民借貸現(xiàn)象較為普遍(占總樣本農(nóng)牧民的65.50%),但多是為了維持日常生活。這與Yan等[41]研究相似,在高原地區(qū),阻礙農(nóng)牧民從事非農(nóng)活動的因素在于以下四點：沒有技能,文盲,缺乏資金,缺乏社會網(wǎng)絡(luò)。在其他地區(qū),也有類似發(fā)現(xiàn)。在撒哈拉地區(qū),Vliet等[29]發(fā)現(xiàn)與氣候變化、社會經(jīng)濟因素以及土地產(chǎn)權(quán)相比,人口的增長是耕地擴張的主要原因。在印度,在人口增長與糧食生產(chǎn)矛盾日益加劇的情況下,農(nóng)戶不斷地開墾土地以滿足需要[36]。在內(nèi)蒙古,由于草地質(zhì)量下降且放牧規(guī)模擴大,農(nóng)牧民也會采取開墾舉措[39]。

農(nóng)牧民非農(nóng)工資越高,其開墾的可能性就越低。這與Clay等[51]、Holden等[52]的研究一致,勞動力析出導(dǎo)致非農(nóng)勞動力與農(nóng)業(yè)勞動力之間存在競爭,且較高的非農(nóng)收入降低了農(nóng)牧民從事農(nóng)牧業(yè)的積極性,從而降低了農(nóng)牧民的開墾意愿。Goldstein等的研究表明[53]青藏高原的農(nóng)牧民已經(jīng)越來越能夠接受高收入的非農(nóng)工作,而且這種趨勢將會繼續(xù)。

在氣候變化背景下,一江兩河地區(qū)將開墾作為適應(yīng)策略的較少。這與諸多研究相似,在氣候變化背景下,以農(nóng)業(yè)為主導(dǎo)的地區(qū)由于存在諸多的非農(nóng)就業(yè)機會或耕地處于飽和狀態(tài),農(nóng)戶更多地選擇了耕地集約化和生計多樣化適應(yīng)策略[54- 56]。發(fā)現(xiàn)壤塘縣是開墾戶數(shù)多、開墾面積最大的區(qū)域,而河湟谷地開墾最少。不難理解,壤塘縣處于農(nóng)牧交錯帶上,農(nóng)牧民將畜牧業(yè)作為主要的生計活動,因此,為了滿足牲畜養(yǎng)殖的需要,他們會開墾更多的土地。而河湟谷地的第二、第三產(chǎn)業(yè)較發(fā)達,農(nóng)牧民愿意從事更高收入的非農(nóng)活動,且受國家禁墾管制,他們開墾最少。

在其他國家和地區(qū),氣候變化下農(nóng)牧民更多地采取耕地集約化策略或從事生計多樣化。如在尼泊爾朱姆拉區(qū),氣候向暖干型變化,加之受資產(chǎn)基礎(chǔ)和社會服務(wù)較差、適宜種植土地面積銳減等因素的限制,大多數(shù)農(nóng)戶選擇種植多樣化的作物(如種植抗旱作物和蔬菜),且傾向于從事高收入的非農(nóng)活動,此外,中上等經(jīng)濟條件的農(nóng)戶還會撂荒原有耕地并購買有灌溉設(shè)施的土地[55]。在伊朗,在過去50年中經(jīng)歷了近27次干旱,為維持并提高農(nóng)業(yè)生產(chǎn)率,農(nóng)戶不得不增加對耕地以及水資源的投資,具體表現(xiàn)在減少種植面積,改善灌溉設(shè)施、挖掘或深化水井、平整土地等[56]。在俄羅斯,農(nóng)戶大量種植抗霜凍和抗病性作物、使用天氣預(yù)報以及增加對耕地的投資(如增加溫室保護措施,引進新技術(shù)以保證灌溉,購置農(nóng)用設(shè)備,增加有機肥和防蟲劑的使用等)等[57]。在加拿大,降水量下降、蒸發(fā)蒸騰量大幅度增加、極端事件(熱浪和干旱等)更加頻繁,大多數(shù)農(nóng)戶采取的措施是使用省工型工具、改善灌溉條件、購買農(nóng)業(yè)保險等[54,58]。其中,俄羅斯和加拿大是氣候變化下受益的國家(主要表現(xiàn)在適宜耕種面積的增加)[59],與其他國家相比,這些地區(qū)的農(nóng)戶都更傾向于選擇耕地集約化適應(yīng)策略。

青藏高原是生態(tài)脆弱區(qū)[60],土地開墾對該地區(qū)造成的影響是不可逆的,因此,我們必須考慮土地開墾的生態(tài)效應(yīng)。諸多研究表明,土地開墾是沙漠化的一個重要原因[61- 62]。Wang等[63]和Zhou[64]通過對青藏高原東北邊緣地區(qū)的研究發(fā)現(xiàn),人口壓力促進了耕地的擴張,但是耕地的擴大和過度放牧導(dǎo)致草地退化現(xiàn)象嚴(yán)重,進而加劇了當(dāng)?shù)厣鷳B(tài)系統(tǒng)的脆弱性。同時,也有研究表明國家實施的禁牧輪牧休牧政策使土地沙化情況得到了改善[65]。因此,有必要控制農(nóng)牧民放牧規(guī)模以及開墾行為。從調(diào)研數(shù)據(jù)可以發(fā)現(xiàn)僅有120戶(占總樣本數(shù)的19.80%)選擇開墾,且大多農(nóng)牧民開墾面積較小,多在0.03—0.13 hm2之間。然而,今后人口壓力是否會增大,是否會促進耕地開墾,其環(huán)境效應(yīng)如何,這些問題仍需開展持續(xù)的研究。

在野外調(diào)查中,與農(nóng)牧局、鄉(xiāng)鎮(zhèn)、村組以及農(nóng)戶討論了耕地開墾的監(jiān)管問題及解決措施。結(jié)合多個尺度的意見以及本文的結(jié)果,為了防止耕地開墾,提出如下建議：(1)需要為當(dāng)?shù)剞r(nóng)牧民提供非農(nóng)就業(yè)機會,加快生計多樣化的步伐；(2)提高灌溉水平,推行集約化農(nóng)業(yè)；(3)政府和村集體對于集體所有的草地要加強監(jiān)管,防止開墾；(4)對于貧困牧民,結(jié)合精準(zhǔn)扶貧措施,提高其生活補助標(biāo)準(zhǔn),降低其對耕地和草地資源的依賴。

5 結(jié)論

對于青藏高原氣候變化背景下的墾殖,目前雖然有少量區(qū)域尺度的模型模擬和遙感調(diào)查,但一直缺少對農(nóng)牧民是否開墾以及開墾原因的實證研究。本文通過對青藏高原的3個典型產(chǎn)糧區(qū)的調(diào)查數(shù)據(jù)分析,證實了農(nóng)牧民開墾是氣候變化和人口壓力共同作用的結(jié)果。研究還發(fā)現(xiàn)該地區(qū)農(nóng)牧民非農(nóng)收入與農(nóng)牧民的開墾行為呈顯著負相關(guān)關(guān)系。隨著城鎮(zhèn)化不斷地加快,農(nóng)牧民對非農(nóng)工作的認可度越來越高,加之政府的管制,使得開墾現(xiàn)象并不普遍。由于農(nóng)牧民的開墾對青藏高原的生態(tài)環(huán)境極有可能造成不可逆的負面影響,應(yīng)該提供更多的非農(nóng)就業(yè)機會,促進農(nóng)業(yè)集約化,加強監(jiān)管,以降低開墾的可能性。