陳曉霞，馮俊人△，廖禮彬，石福孫，王健順，高武超，孫飛達(dá)*

(1.四川農(nóng)業(yè)大學(xué) 草業(yè)科學(xué)系，四川 成都 611130；2.中國科學(xué)院成都生物研究所，四川 成都 610041；3.蘭州大學(xué)草地農(nóng)業(yè)科技學(xué)院，甘肅 蘭州 730020)

川西北牧區(qū)作為青藏高原東南緣和橫斷山脈的一部分，是我國重要的高寒牧區(qū)，也是四川省重要的畜牧生產(chǎn)基地，在生態(tài)保護(hù)和畜牧業(yè)發(fā)展方面意義突出[1-3]。不同畜種的地理分布格局，反映了生態(tài)環(huán)境對家畜的影響及不同畜種對環(huán)境的適應(yīng)程度，揭示了特定地域的家畜生態(tài)特征[4]。牦牛和藏綿羊作為在川西北地區(qū)生活的優(yōu)勢畜種，是草業(yè)系統(tǒng)物流與能流中轉(zhuǎn)站動(dòng)物生產(chǎn)層的重要組成部分，其與草地互作是放牧生態(tài)系統(tǒng)物質(zhì)循環(huán)的關(guān)鍵驅(qū)動(dòng)力之一[5-8]。畜群結(jié)構(gòu)和規(guī)模直接影響到草地生態(tài)系統(tǒng)的第二性生產(chǎn)力，間接地影響著經(jīng)濟(jì)價(jià)值的生產(chǎn)[9]。合理的畜群結(jié)構(gòu)對草原畜牧業(yè)發(fā)展具有重要意義，有利于畜群再生產(chǎn)和擴(kuò)大再生產(chǎn)、正常周轉(zhuǎn)的同時(shí)取得最大的經(jīng)濟(jì)效益；畜群結(jié)構(gòu)不合理，會導(dǎo)致畜牧業(yè)生產(chǎn)出欄率低、周轉(zhuǎn)慢和效能低[10-12]。因此，了解川西北典型牧區(qū)牦牛和藏綿羊的數(shù)量變化及畜群結(jié)構(gòu)特征，明晰川西北畜群存在空間分布差異的原因，對促進(jìn)青藏高原高寒草地的可持續(xù)利用和畜牧業(yè)健康發(fā)展具有重要意義。

家畜作為重要的土地利用者，其與自然、經(jīng)濟(jì)、政治、社會環(huán)境之間存在錯(cuò)綜復(fù)雜的聯(lián)系[13]。多年來，許多學(xué)者在家畜地理空間分布與氣候因素方面進(jìn)行了大量的研究工作[14-16]，但利用空間異質(zhì)性探測工具，并結(jié)合人文和自然因素的分析相對較少。鑒于此，本研究利用地理探測器模型結(jié)合降雨、溫度、高程、坡度、土地覆蓋類型、草地面積、植被類型、植被覆蓋度等自然因素，并考慮人口數(shù)量、人口密度、出欄管理等人為因素的影響，對川西北放牧家畜存在空間分異的原因進(jìn)行探究。地理探測器模型( http: //www.geodetector.org /)最初被用于探尋地理空間分區(qū)因素對疾病風(fēng)險(xiǎn)的影響，之后在自然科學(xué)和環(huán)境科學(xué)研究中逐漸普及，是探測空間分異性，以及揭示分異驅(qū)動(dòng)力的一組統(tǒng)計(jì)學(xué)方法，與其他空間異質(zhì)性探測工具相比地理探測器具有更高的解釋效率[17-20]。其核心思想是：如果某個(gè)自變量對某個(gè)因變量有重要影響，那么自變量和因變量的空間分布應(yīng)該具有相似性[21]。利用該模型既可以探測數(shù)值型數(shù)據(jù)，也可以探測定性數(shù)據(jù)的一大優(yōu)勢，對川西北放牧家畜空間分布的影響因素進(jìn)行定量分析，對川西北家畜空間尺度下分布差異歸因的探究具有重要意義。

1 材料與方法

1.1 研究區(qū)域概況

阿壩藏族羌族自治州位于青藏高原東段，地理位置處于東經(jīng)100°30′～104°27′，北緯 30°35′～34°19′之間；幅員84 200 km2，現(xiàn)轄馬爾康1個(gè)市，阿壩、若爾蓋、紅原、壤塘、松潘 5個(gè)牧區(qū)縣，金川、小金、汶川、理縣、茂縣、九寨溝、黑水7個(gè)半農(nóng)半牧區(qū)縣，共有54個(gè)鎮(zhèn)、165個(gè)鄉(xiāng)；總體地形是西北高，東南低，自西北向東南傾斜，西北部分海拔在3 500～4 000 m之間，屬于丘狀高原，而南部為高山峽谷，海拔780～6 250 m；全州屬于高原季風(fēng)氣候，光照充沛，晝夜溫差大，無霜期短，冬春季節(jié)空氣干燥，多陣性大風(fēng)，旱、霜、雪、低溫、大雪各類災(zāi)害性天氣頻繁；本文主要研究區(qū)域?yàn)榘沃莸娜魻柹w、阿壩、紅原、壤塘、松潘5個(gè)典型牧業(yè)縣，該區(qū)域位于阿壩州西北部的丘狀高原地帶，屬大陸高原性氣候，四季氣溫?zé)o明顯差別，冬季嚴(yán)寒漫長，夏季涼寒濕潤，年平均氣溫0.8～4.3 ℃，年均降雨量600～800 mm之間[22]。若爾蓋、阿壩、紅原、壤塘、松潘5縣作為阿壩州主要的畜牧生產(chǎn)基地，面積44 621 km2，占全州幅員面積的54%，天然草地面積32 000 km2，占全州天然草地面積的71%。據(jù)統(tǒng)計(jì)2016年末若爾蓋、阿壩、紅原、壤塘、松潘五縣的戶籍總?cè)丝跀?shù)為327 700，占全州人口的35.6%；5個(gè)牧業(yè)縣牦牛和綿羊混合存欄量為1 953 300 混合頭，占整個(gè)阿壩州地區(qū)牛羊混合存欄量的70.64 %[23]。

1.2 影響因子的確定

從自然和人為兩個(gè)因素探究川西北總放牧家畜、牦牛、綿羊空間分布的原因，并分別從兩類影響因素下選取代理變量。自然因素對于家畜空間格局的作用機(jī)理體現(xiàn)家畜適應(yīng)性，基于已有家畜地理空間格局研究的理論和實(shí)證分析，選取年均降雨量(X1)、年均溫(X2)、高程(X3)、坡度(X4)、草地面積(X5)土地覆蓋類型(X6)、植被類型(X7)、植被覆蓋度(X8)等8個(gè)因子作為自然因素代理變量。人口和家畜的出欄管理等人為因素能間接反映經(jīng)濟(jì)情況，與家畜分布密切相關(guān)，所以選用人口數(shù)量(X9)、人口密度(X10)、出欄率(X11)等作為人為因素的代理變量(見圖1-3)。

圖2 地理探測器中自然影響因素的空間分布Fig. 2 Distribution of natural factors in geographical detector

圖3 地理探測器中人為影響因素的空間分布Fig. 3 Distribution of human factors in geographical detector

1.3 數(shù)據(jù)來源及處理

本研究分別以川西北總的放牧家畜數(shù)量及不同畜種牦牛和綿羊的數(shù)量為因變量，以8個(gè)自然因素代理變量和3個(gè)人為因素代理變量為自變量數(shù)據(jù)，使用王勁峰等[21]提出的地理探測器模型，分析川西北放牧家畜數(shù)量和畜種結(jié)構(gòu)存在空間分布差異的原因。

1.3.1 家畜空間分布 作為因變量的家畜數(shù)據(jù)來源于阿壩州農(nóng)業(yè)畜牧局提供的2017年家畜鄉(xiāng)鎮(zhèn)單位數(shù)據(jù)，利用屬性數(shù)據(jù)空間化技術(shù)實(shí)現(xiàn)家畜分布格局圖的構(gòu)建。

1.3.2 自然因素 自然因素代理變量年降雨量、年均溫?cái)?shù)據(jù)使用2006-2016年中國地面累積月值數(shù)據(jù)集( 2006-2016年) 由中國氣象科學(xué)數(shù)據(jù)共享服務(wù)平臺 ( http: /data.cma.cn /site /index.html) 提供，利用Kring和CoKring方法進(jìn)行插值處理；高程數(shù)據(jù)來源于地理空間數(shù)據(jù)云(http://www.gscloud.cn/)提供的30 M分辨率的數(shù)字高程數(shù)據(jù)，坡度信息是利用Arc GIS10.5中的3D分析模塊從高程數(shù)據(jù)中提取的獲得；鄉(xiāng)鎮(zhèn)尺度草地面積數(shù)據(jù)來源于阿壩州農(nóng)業(yè)畜牧局2017年家畜統(tǒng)計(jì)年鑒；土地覆蓋類型數(shù)據(jù)使用的是美國國家航空航天局(https://search.earthdata.nasa.gov/)提供的MCD12Q1產(chǎn)品數(shù)據(jù)；植被類型數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心提供的1∶1 000 000植被類型空間分布數(shù)據(jù)；植被覆蓋度利用遙感數(shù)據(jù)MOD13Q1合成，該產(chǎn)品是美國航空航天局提供的陸地3級標(biāo)準(zhǔn)數(shù)據(jù)產(chǎn)品，空間分辨率為 250 m×250 m，時(shí)間分辨率為 16 d，本文選擇2016年23景MODISNDVI數(shù)據(jù)，利用最大合成法 ( Maximum Value Composite，MVC) 合成年最大 NDVI，然后用像元二分法進(jìn)行植被覆蓋度的估算，估算公式如下：

(1)

式中：NDVI為像元實(shí)際值； NDVIsoil為地表裸土覆蓋所對應(yīng)的 NDVI值；NDVIveg為地表全部植被所覆蓋對應(yīng)的NDVI值，在運(yùn)用像元二分模型估算植被覆蓋度時(shí)，通常參照經(jīng)驗(yàn)在給定置信區(qū)間內(nèi)確定NDVIsoil和NDVIveg的值[24]。本研究采用的置信區(qū)間取值范圍是3%～95%，選取累計(jì)頻率為3%的 NDVI值作為NDVIsoil，累計(jì)頻率為95%的NDVI 值NDVIsoil，然后進(jìn)行植被覆蓋度的估算。

1.3.3 人為因素 人為因素代理變量人口密度、人口數(shù)量、出欄率為鄉(xiāng)鎮(zhèn)尺度數(shù)據(jù)，均來源于2017年中國縣域統(tǒng)計(jì)年鑒(縣市卷)和阿壩州農(nóng)業(yè)畜牧局2017年家畜統(tǒng)計(jì)年鑒，利用ArcGIS10.5實(shí)現(xiàn)數(shù)據(jù)空間化。

1.4 研究方法

地理探測器模型要求自變數(shù)據(jù)為類型變量，所以除了土地覆蓋類型和植被類型外，其他代理變量均使用自然斷點(diǎn)法進(jìn)行重分類，利用ArcGIS 10.5中的漁網(wǎng)工具將研究區(qū)域劃分為3 km×3 km的網(wǎng)格，研究區(qū)域采樣網(wǎng)格點(diǎn)4876個(gè)，用于將因變量Y和自變量X屬性值相匹配，然后將X和Y值導(dǎo)入地理探測器模型進(jìn)行運(yùn)算。運(yùn)算原理公式表達(dá)如下：

(2)

2 結(jié)果與分析

2.1 川西總放牧家畜空間分布的影響因子分析

使用地理探測器模型對放牧家畜空間分布的影響因素進(jìn)行分析，得出各因子對家畜空間分布影響程度排序?yàn)椋喝丝跀?shù)量X9>年均溫X2>草地面積X5>出欄率X11-1>坡度X4>年降雨量X1>植被類型X7>人口密度X10>高程X3>土地覆蓋類型X6>植被覆蓋度X8，其中人口數(shù)量的q值最高為0.48，年均溫和草地面積次之分別為0.42和0.35，植被覆蓋度和土地覆蓋類型的q值最小分別為0.027和0.069。當(dāng)兩種影響因子同時(shí)作用于因變量Y時(shí)，交互作用對Y的解釋作用增強(qiáng)，其中草地面積∩年均溫、人口數(shù)量∩家畜出欄率、草地面積∩人口數(shù)量q值最大分別為0.809、0.706、0.759。

2.2 川西牦?？臻g分布的影響因子分析

各影響因子對牦牛的空間分布影響程度排序?yàn)椋喝丝跀?shù)量X9>草地面積X5>年均溫X2>植被類型X7>人口密度X10>牦牛出欄率X11-2>坡度X4>土地覆蓋類型X6>高程X3>年降雨量X1>植被覆蓋度X8。人口數(shù)量和草地面積是牦?？臻g分布最主要的影響因素，其q值分別為0.45和0.42，年均溫和植被類型次之分別為0.30和0.20，年降雨量和植被覆蓋度q值最低為0.07和0.04，對牦牛的空間分布的影響最??；從兩因素的交互作用來看，自變量的代理變量之間兩兩交互對牦?？臻g分布的影響大于單一變量的作用，人口數(shù)量和草地面積分別與其他因子互作，q值均高于0.45，牦牛出欄率和人口數(shù)量、草地面積之間的交互對牦?？臻g分布的影響最大，q值分別為0.61和0.64；人口密度∩草地面積以及年均溫∩草地面積對牦?？臻g分布的影響次之，q值分別為0.586和0.576，植被覆蓋度和土地覆蓋類型除了與人口數(shù)量和草地面積交互作用較強(qiáng)以外，與其他因子的交互作用相對較弱q值小于0.35，影響力度最小。

圖4 各影響因素對總的放牧家畜空間分布的解釋力程度Fig. 4 Explanatory power of various influencing factors on the spatial distribution of total grazing livestock

圖5 各影響因素交互作用對總的放牧家畜分布的影響Fig. 5 Interactive impact of various influencing factors on the spatial distribution of total grazing livestock

2.3 川西綿羊空間分布的影響因子分析

從q值來看，降雨和溫度這兩種氣象因子對綿羊在川西北空間分布的影響程度最大，其q值分別為0.30和0.29，其次是人文因素的代理變量人口數(shù)量和綿羊出欄率，q值分別為0.22和0.17，植被覆蓋度和坡度對綿羊的空間分布影響最小，q值為0.007和0.004；從交互作用來看，草地面積和年均溫、人口數(shù)量的交互作用對綿羊空間分布的影響最大，q值為0.72和0.61，年降雨量分別與人口數(shù)量、年均溫交互作對綿羊空間分布的影響次之，q值分別為0.59和0.51，說明草地面積、年降雨量、年均溫等自然因素變量對川西北綿羊的空間分布影響較大于人口數(shù)量因素。

圖6 各影響因素對牦?？臻g分布的解釋力程度Fig. 6 Explanatory power of various influencing factors on the spatial distribution of yak

圖7 各影響因素交互作用對牦牛分布的影響Fig. 7 Interactive impact of various influencing factors on the spatial distribution of yak

圖8 各影響因素對綿羊空間分布的解釋力程度Fig. 8 Explanatory power of various influencing factors on the spatial distribution of sheep

圖9 各影響因素交互作用對綿羊分布的影響Fig. 9 Interactive impact of various influencing factors on the spatial distribution of sheep

3 討 論

家畜數(shù)量隨著人口的增加和家庭規(guī)模的擴(kuò)大而呈增加趨勢[25-26]，本研究也驗(yàn)證了這一觀點(diǎn)，對比其他影響因子，發(fā)現(xiàn)人口數(shù)量對川西北放牧家畜空間分布的影響起著主導(dǎo)作用。川西北地區(qū)屬于青藏高原高寒脆弱區(qū)，氣候是影響該地區(qū)畜牧業(yè)生產(chǎn)的重要因素，溫度是限制草地生產(chǎn)力的主要生態(tài)因子之一[27]，本研究表明氣象因子溫度對川西北典型牧業(yè)縣放牧家畜的空間分布影響突出，且遠(yuǎn)大于降雨對家畜分布的影響，說明溫度和人口數(shù)量一樣也是限制該地區(qū)放牧家畜對生存環(huán)境適應(yīng)性選擇的主要因子。草原占地球陸地面積的一半以上,其植物生產(chǎn)占全球總量的30 %～35 %，是家畜生產(chǎn)的主要基地，養(yǎng)活著世界上約1/3 的人口[28]。在川西北牧區(qū)家畜生產(chǎn)是牧民主要的收入來源，劉媛媛[29]根據(jù)牧業(yè)部門提供的數(shù)據(jù)發(fā)現(xiàn)在紅原和若爾蓋地區(qū)理論載畜數(shù)量完全不能滿足牧民基本的生產(chǎn)生活所需，隨著人口的增加，牧戶一味追求數(shù)量型發(fā)展，導(dǎo)致草畜矛盾突出。因此，筆者認(rèn)為草地面積對家畜空間分布的解釋程度低于人口數(shù)量，也反映了隨著社會經(jīng)濟(jì)的發(fā)展草地資源不能滿足該地區(qū)人均基本的需要，存在超載的事實(shí)。

牦牛作為川西北的主要放牧家畜，空間分布跟總放牧家畜的分布具有相似特征，從較大空間尺度范圍來看，川西北地區(qū)人口分布對牦牛及總的放牧家畜空間分布的解釋程度大于其他因素，這與黃敬蜂等[30]認(rèn)為的牦牛常年放牧，其分布主要取決于海拔高度、氣候及天然草場的植被類型，而受人為干擾較少的觀點(diǎn)有所不同。與各影響因素對總放牧家畜及牦?？臻g分布的解釋程度存在差異，綿羊在川西北地區(qū)主要分布于降雨量相對較少的若爾蓋縣地區(qū)，其空間分布的異質(zhì)性主要?dú)w因于降雨和溫度，其次是人口數(shù)量和出欄率。部分研究表示溫度和降雨與羊的幼畜成活率[31]、腐蹄病的發(fā)病率具有一定的關(guān)系[32-33]。本研究認(rèn)為降雨量過多、溫度過高，增加羊只患病率，提高養(yǎng)殖成本和管理難度，是制約川西北藏綿羊分布范圍的重要因素。

從各影響因素對家畜空間分布的交互作用來看，各影響因素對牦牛、綿羊、總的家畜空間格局的影響表現(xiàn)為協(xié)同增強(qiáng)，這與胡江等[34]提出的家畜空間分布受多種因子交互影響的觀點(diǎn)相類似。

本文以2016年川西北典型牧業(yè)縣鄉(xiāng)鎮(zhèn)尺度的放牧家畜總存欄數(shù)、牦牛存欄數(shù)、綿羊存欄數(shù)、人口數(shù)量、人口密度、出欄率、植被覆蓋度和2006-2016年十年間的年降雨量和年均溫及該地區(qū)的高程、坡度、植被類型、土地覆蓋類型、草地面積等數(shù)據(jù)為基礎(chǔ)，有效地分析了各因子對川西北家畜空間分布的影響。但影響家畜空間分布的因素較多，以鄉(xiāng)鎮(zhèn)為單元的數(shù)據(jù)獲取比較困難，所以家畜分布和少數(shù)影響因素的分析僅立足于2016年的部分?jǐn)?shù)據(jù)，對于其他影響因素的補(bǔ)充和更長時(shí)間尺度的分析，有待進(jìn)一步探討。

5 結(jié) 論

(1) 各影響因素對家畜空間分布的影響具有差異，牦牛作為川西北的主要放牧家畜，跟總放牧家畜的分布具有相似特征。單因子方面，人口數(shù)量是影響川西北總的放牧家畜和牦?？臻g分布的主要人為因子，年均溫和草地面積是主要的自然因子；

(2) 總體來看，自然因素對川西北地區(qū)綿羊空間分布的影響大于人為因素，其中自然因素代理變量降雨和溫度的解釋程度最大；

(3) 家畜存在空間分布差異是多種因素互相作用的結(jié)果，各影響因子之間的交互作用對川西北放牧家畜的影響力表現(xiàn)為協(xié)同增強(qiáng)。