青藏高原氣候變化對(duì)農(nóng)牧民開墾的影響

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

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

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

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

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

青藏高原是生態(tài)脆弱區(qū),農(nóng)牧交錯(cuò)帶草地的墾殖會(huì)導(dǎo)致草地退化和沙漠化,并且不可逆轉(zhuǎn)。研究氣候變化背景下青藏高原農(nóng)牧民是否會(huì)采取開墾耕地的適應(yīng)策略,不僅僅涉及當(dāng)?shù)剞r(nóng)牧民的生計(jì)和氣候變化適應(yīng)策略,也和區(qū)域的可持續(xù)發(fā)展息息相關(guān)。對(duì)于青藏高原氣候變化背景下的墾殖,目前雖然有少量區(qū)域尺度的模型模擬和遙感調(diào)查,但一直缺少對(duì)農(nóng)牧民是否開墾以及開墾原因的實(shí)證研究,牧民如何平衡氣候變化帶來的開墾機(jī)遇、風(fēng)險(xiǎn)以及人口壓力,仍然缺少認(rèn)識(shí)。本文擬以青藏高原的“一江兩河”地區(qū)、河湟谷地以及川西地區(qū)的壤塘縣為研究區(qū),調(diào)查氣候變化和勞動(dòng)力析出背景下農(nóng)牧民的耕地開墾情況,并定量研究氣候變化對(duì)農(nóng)牧民開墾耕地的影響。論文由以下幾個(gè)部分組成,第二部分介紹了研究區(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,占全國(guó)總面積的26.8%,平均海拔4383 m。該地區(qū)輻射強(qiáng),日照時(shí)間長(zhǎng),氣溫低且日較差大,積溫少；干濕分明,多夜雨；冬季干冷漫長(zhǎng),大風(fēng)多；夏季溫涼多雨,冰雹多。青藏高原產(chǎn)糧區(qū)主要包括西藏自治區(qū)、青海省和川西地區(qū),因此本文選擇這3個(gè)地區(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í)數(shù)為2800—3300 h,年降水量在251.7—580.0 mm之間,蒸發(fā)量大,約為降水量的4—5倍。該地區(qū)下轄拉薩、山南、日喀則3個(gè)地市的18個(gè)縣(市、區(qū)),農(nóng)牧業(yè)發(fā)達(dá),農(nóng)作物以青稞、小麥、油菜為主,牲畜包括牦牛、犏牛、奶牛、馬和綿羊。據(jù)2017西藏統(tǒng)計(jì)年鑒可知,2016年,該地區(qū)耕地面積達(dá)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í)數(shù)在2600—3000 h之間,年平均溫度在5—9℃之間,年平均降水量在252—535 mm之間。該地區(qū)下轄西寧市、海東市、海北藏族自治州、黃南藏族自治州和海南藏族自治州的17個(gè)區(qū)縣,農(nóng)牧業(yè)較發(fā)達(dá),農(nóng)作物以小麥、青稞、玉米、油菜和枸杞為主,牲畜主要包括綿羊、山羊和豬。據(jù)2017青海統(tǒng)計(jì)年鑒可知,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)作物實(shí)際播種面積18.42 km2,其中糧食作物播面積為15.79 km2,總產(chǎn)量3114 t；經(jīng)濟(jì)作物播面積1.59 km2,總產(chǎn)量264 t；農(nóng)村經(jīng)濟(jì)總收入26891.51萬元,農(nóng)牧民人均純收入6413元。

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

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

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

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

2.1 模型的選擇

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

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

(1)

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

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

2.2 變量的設(shè)定

2.2.1解釋變量的選取

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

雖然農(nóng)牧民對(duì)氣候的感知與實(shí)測(cè)數(shù)據(jù)不一定一致[32- 33],但諸多研究表明農(nóng)牧民氣候感知是影響其響應(yīng)行為的關(guān)鍵因素[34]。當(dāng)農(nóng)牧民感知?dú)夂蜃兓鶐淼娘L(fēng)險(xiǎn)和適應(yīng)能力較低時(shí),農(nóng)牧民參與適應(yīng)行為的可能性降低；相反,當(dāng)農(nóng)牧民感知?dú)夂蜃兓瘡?qiáng)烈時(shí),農(nóng)牧民便會(huì)采取適應(yīng)策略[35]。本文將農(nóng)牧民氣候感知分為三類：季節(jié)性感知、氣象感知、風(fēng)險(xiǎn)感知,并從這三個(gè)角度來選取了農(nóng)牧民感知?dú)夂蛑笜?biāo)。

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

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

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

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

2.2.2解釋變量間共線性檢驗(yàn)

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

3 結(jié)果分析

3.1 農(nóng)牧民感知?dú)夂蜃兓?/h3>

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

表2 研究區(qū)樣本農(nóng)牧民感知?dú)夂蜃兓闆r/戶

圖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)牧民實(shí)際開墾情況

研究區(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)牧交錯(cuò)帶,為了滿足其人口不斷地增長(zhǎng)與牲畜數(shù)量增加而造成的糧食需求,農(nóng)牧民會(huì)開墾較大的面積。河湟谷地雖然開墾戶數(shù)最少且開墾面積最小,但戶均開墾面積最大,為0.19 hm2。

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

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

樣本農(nóng)牧民大多在2005—2016年發(fā)生開墾行為,占總開墾農(nóng)戶的60.83%。壤塘縣在2010—2016年間開墾戶數(shù)最多,高達(dá)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àn)槿丝趬毫Χ_墾(主要是為了自食和養(yǎng)殖牲畜)。同時(shí),農(nóng)牧民開墾地主要種植燕麥和青稞兩種作物(表6),這也足以說明開墾的原因是人口壓力。在壤塘縣,有54戶開墾農(nóng)牧民種植了飼料性作物,而在一江兩河地區(qū),開墾農(nóng)牧民主要種植糧食作物。

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

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

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

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

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

3.5 計(jì)量結(jié)果分析

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

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

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

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

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

戶主受教育水平與農(nóng)牧民開墾行為呈顯著負(fù)相關(guān)關(guān)系。這主要是因?yàn)槭芙逃皆礁叩霓r(nóng)牧民有更大的機(jī)會(huì)去從事非農(nóng)工作,從而會(huì)有更高的收入。

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

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

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

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

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

4 討論

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

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

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

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

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

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

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

5 結(jié)論

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