徐敏云,高立杰,李運起*
(1.河北農(nóng)業(yè)大學動物科技學院,河北 保定071000;2.北京大學城市與環(huán)境學院,北京100871)
家畜放牧是天然草地最主要的利用方式,合理放牧對于草地利用的可持續(xù)性以及家畜生產(chǎn)和草地生態(tài)系統(tǒng)服務至關重要[1-2]。草地放牧管理的目標是維持家畜牧草需求量和草地飼料供應量的平衡,獲取最大的經(jīng)濟效益[3]??茖W地確定草地載畜量和放牧率,實現(xiàn)草畜平衡,有利于保護草地生態(tài)環(huán)境和促進草地畜牧業(yè)持續(xù)發(fā)展[4-5]。目前,確定載畜量的方法在家畜采食量、牧草產(chǎn)量、放牧區(qū)域選擇、放牧利用率、水源距離調(diào)整以及補飼等逐漸趨向標準化[5]。但相關參數(shù)計算和載畜量估算研究中存在著爭議和分歧[6-9]。本文依據(jù)國內(nèi)外研究進展,綜述草地載畜量校正參數(shù)、計算方法與調(diào)整策略,以期為草地放牧管理提供參考依據(jù)。
草地載畜量是在不危害草地的情況下,草地能夠負荷的多年最大平均放牧率[10]。草地載畜量由多個因素決定,主要取決于草地初級生產(chǎn)力[11]。草地總生物量是確定載畜量的基礎,以可食牧草產(chǎn)量比用總產(chǎn)量更能準確計算草地載畜量,但對地面測量或遙感數(shù)據(jù)獲取而言都存在實際困難[12]。此外,地面測量對于大面積草地而言并不現(xiàn)實[13]。因此,應合理使用模擬模型,包括植被-氣候模型、過程模型、遙感估算模型等,獲取草地產(chǎn)量[13-14]。過程模型參數(shù)需要長期的數(shù)據(jù);如果校準方法不當,遙感估算的精度難以滿足估算要求[14-15]。植被-氣候模型由于使用方便,在草地產(chǎn)量估測上應用廣泛[16]。有研究[14]在內(nèi)蒙古阿拉善草地,采用 Miami模型、Schuur模型、Chikugo模型、Beijing模型、Synthetic模型、Classification Indices-based模型等6個植被-氣候模型估測草地產(chǎn)量。草地初級生產(chǎn)力主要受控于氣候因素,尤其是降水及其季節(jié)分布[17]。全球尺度上,有基于蒸散-初級生產(chǎn)力[18],及年降水、年均溫度與凈初級生產(chǎn)力(NPP)[19]的關系模型,Oesterheld等[20]還建立了溫帶草地載畜量與草地初級生產(chǎn)力關系模型。
降水是影響草地產(chǎn)量和植被構成的最主要因素,尤其對干旱、半干旱生態(tài)系統(tǒng),由于其年內(nèi)和年際降水變率大,干旱抑制草地植被生長,降水是影響草地植被生產(chǎn)力的主要因素[21],從而影響載畜量。但僅利用降水指標估計地上生物量進而估算載畜量過于簡單,由于地形、地貌差異,土壤養(yǎng)分的有效性劇烈變化,影響地上生物量,也是區(qū)域范圍內(nèi)決定載畜量的重要因素[22]。土壤因素對草地植被生產(chǎn)力的影響等同于降水,甚至超過降水的效應[23]。在降水量小于250mm的地區(qū),有效水分是限制產(chǎn)量的主要因素;但當降水量大于250mm,土壤氮含量則成為主導產(chǎn)量的關鍵[23]。即使降水量差別不大,土壤類型的差異也會造成產(chǎn)量和載畜量的巨大差異[23]。土壤有效水取決于土壤水分滲漏、土壤持水量以及徑流和蒸發(fā)的強度和速率,這些都受地形坡度、植被覆蓋、土壤類型等因素影響,從而影響草地載畜量[23]。Noy-Meir[24]提出IST(inverse soil texture)假說,干旱、半干旱地區(qū),粗質(zhì)土的初級生產(chǎn)力要高于細質(zhì)土,進而影響草地載畜量。Fritz和Duncan[22]證實,土壤質(zhì)量是降水之外,影響非洲熱帶稀樹草地載畜量的重要因素。在干旱、半干旱地區(qū),物種豐度也是影響牧草產(chǎn)量從而決定載畜量的因素,低產(chǎn)生產(chǎn)系統(tǒng),隨著物種豐度提高,載畜量也相應提高[22]。
1.2.1 草地利用率定義 草地畜牧業(yè)生產(chǎn)中,放牧管理存在3種可能性,過牧、輕牧和適度放牧。過牧常會導致牧草及土壤肥力的下降,優(yōu)良牧草的減少。輕牧常引起放牧不均勻,形成草地斑塊,還會造成牧草浪費和畜牧業(yè)生產(chǎn)效益下降。因此,草地放牧的關鍵是在畜牧業(yè)生產(chǎn)發(fā)展和草地生態(tài)安全之間取得平衡[12]。相關綜述[25-27]總結并實例介紹草地對家畜放牧的反應,及家畜對草地的影響。放牧草地系統(tǒng)中,保持適量的地上部分生物量是牧草開花、結實不受影響,及草地放牧后恢復生長、維持草地植被成分和生產(chǎn)能力的關鍵。草本植物具有一定的耐牧性,但也有限度,一旦放牧超過牧草忍耐的臨界點,草地植被失去活力,生產(chǎn)能力下降,最終導致死亡[28]。因此,采用合理的草地利用率、留茬高度對保持草地生態(tài)系統(tǒng)穩(wěn)定性非常重要[29-30]。草地利用率取決于草地條件、放牧系統(tǒng)類型、放牧利用季節(jié)及干旱程度[31]。利用率值的大小對草地載畜量的影響很大[30]。有研究表明,45%的草地利用率比30%的草地利用率,估算的草地載畜量增加1倍[30]。
利用率是被家畜、野生動物和昆蟲利用數(shù)量占當年草地產(chǎn)量的比率,草地產(chǎn)量可以是單一物種,也可以是多個物種,或指全部植被[32-33]。Parker和Glendening[29]把利用率定義為:不會導致優(yōu)良牧草密度和活力下降,避免土壤侵蝕和地表徑流的放牧強度。美國林業(yè)協(xié)會最初用長期監(jiān)測的草地利用率來表示放牧率[34],用以預警和指導調(diào)整草地管理措施,避免草地退化[35]。從更廣的定義講,利用率是在一個放牧時期內(nèi),家畜采食草地飼草的比率[12],既可以作為放牧策略用以確定安全載畜量,也可以作為放牧技術用于根據(jù)環(huán)境條件變化和牧草產(chǎn)量調(diào)整放牧率[12]。作為放牧策略,利用率通常是固定的,不會隨著草地產(chǎn)量的變化而變化[12]。作為放牧技術,草地放牧家畜的數(shù)量隨草地產(chǎn)量的動態(tài)而變化,但在公共草地上,不適合根據(jù)草地等級調(diào)整利用率[36]。
FAO定義了草地安全利用率,即不導致草地退化的情況下,草地放牧采食的最大比率[37]。美國草地管理協(xié)會1999年發(fā)布了一個聲明,承認草地利用率和留茬技術,對草地管理的意義和價值,并將草地利用率定義為:實現(xiàn)草地管理目標,維持或提高草地生產(chǎn)能力,草地牧草當年生長量的利用程度[32]。由于地形、水源地、放牧季節(jié)和其他管理因素的影響,草地利用率差異很大[31]。氣候變化引起草地產(chǎn)量波動,多年生草地降水波動范圍大致在平均降水的-30%~25%之間,草地產(chǎn)量也會在平均產(chǎn)量±30%范圍內(nèi)波動,草地利用率應進行適時調(diào)整[31]。當降水量偏離年均降水50%或更多,就無法給出可靠的利用率[31]。
1.2.2 確定草地利用率的方法 草地利用率的確定方法,主要有經(jīng)驗估計和實際測量兩類[34]。草地安全利用率必須確保草地處于較好的狀態(tài)[38]。Holechek[31]建立了確定草地利用率的程序。實際測量法主要根據(jù)放牧試驗及刈割試驗[39]。放牧試驗法通常耗費財力、歷時較長,無法推廣。經(jīng)驗回顧法成本較低[38],但經(jīng)驗回顧法確定的草地安全利用率一般低于草地實際利用率[12]。此外,還有基于牧草高度和重量函數(shù)關系的植被高度-重量法[40],但由于不同年份的生長差異,可能會有10%~25%的誤差[40]。另一個計算利用率的方法是加權平均法,根據(jù)不同地塊的利用率和面積,得到加權系數(shù),根據(jù)各地塊的相對利用率計算整塊草地的利用率[35,41-42]。大多數(shù)確定利用率的方法都是適用于多年生草地,并不適用于灌叢草地、一年生草地[36]。景觀外貌法(landscape appearance method)[33]采用投影技術,廣泛應用于其他草地[36]。一年生草地的利用率[12,38],草地生長量可以用Aussie-GRASS模型估計[43],采用實際測量法進行糾正。
關鍵種(key-plant)和關鍵場(key-area)在確定利用率時十分有用[36]。關鍵種通常是草地植被中適口性好、蓋度大、產(chǎn)量高的物種,一般草地上存在1~3個關鍵種[27],關鍵種并不是單純指適口性,適口性最好的物種也叫嗜食種(ice cream plants),往往豐度較低,并不能反應整個草地的利用情況[33]。確定關鍵種的利用率,比較其他物種和關鍵種的適口性差異,能夠確定整個草地的利用率。由于地形差異、飼草適口性差別,草地利用不均勻[31,44],關鍵場不應位于飲水點、道路、河床等附近,這些區(qū)域雖然放牧利用率高,但面積較小,并不能反映整個草地的利用情況[33]。根據(jù)草地狀況,分別選擇代表性的地塊作為關鍵場[5]。
1.2.3 草地利用率經(jīng)驗法則 美國林業(yè)服務機構,首次提出維持草地生產(chǎn)的草地安全利用率為15%~20%(體積百分率)[34]。草地利用率的經(jīng)驗法則是采食一半、保留一半(take half-leave half),這也是草地管理多年遵守的標準。采食的50%中,只有25%被家畜采食,另外25%被踐踏、棄食、昆蟲和其他動物采食等,或者由于分解而消失。經(jīng)驗法則其實只有1/4的草地牧草被家畜采食[5]。采食一半、保留一半的經(jīng)驗法則對濕潤草地和一年生草地是適用的[31]。對于干旱半干旱草地,草地利用率一般為35%[5]。大多數(shù)研究表明,美國西部干旱、半干旱草地,35%的利用率比較適宜,并不會導致草地灌叢化[44]。Galt等[5]認為,對大部分美國西部草地而言,要避免長期的牧草匱乏和草地退化,25%的草地利用率是適宜的[5]。高于25%的草地利用率,一旦干旱發(fā)生時,由于牧民不愿意縮減飼養(yǎng)規(guī)模,必然會導致草地退化和經(jīng)濟損失[5]。在非均衡放牧管理的草地,草地利用率通常會超過80%[45]。根據(jù)對草地生產(chǎn)力、家畜生產(chǎn)及草地畜牧業(yè)經(jīng)濟收益的現(xiàn)有文獻,當不確定草地植被狀況及草地產(chǎn)量時,設置初始放牧率,灌叢草地利用率為30%、干旱草地40%、濕潤草地50%、一年生草地55%[31]。研究表明,由于家畜踐踏、野生動物采食及牧草風化等,實際利用率通常會超過設定利用率,因此設置利用率時至少把初始利用率降低5%[5]。
1.2.4 草地利用率應用上的爭議 由于草地利用率的概念、構成及計算草地利用率的時間跨度和方法并不明確,其在草地管理上的應用也有很大爭議[12],計算草地載畜量時,廣泛采用草地年產(chǎn)量、季節(jié)產(chǎn)量和現(xiàn)存量,家畜利用的部分是僅指家畜采食的部分,還是應包括由于踐踏、脫落等部分,也不明確。利用率的計算沒有考慮放牧后牧草的再生,如果當前放牧利用在生長季初期,放牧后再生對草地利用率就很重要,利用率的計算明顯存在錯誤[34]。草地利用率計算周期一般指的是12個月,但其開始和結束的時間也差別很大[12]。不是基于全年生長量計算的利用率,沒有考慮牧草的再生,應稱為相對利用率或季節(jié)利用率[33]。另一個不確定性,放牧試驗采用的牧場面積較小,家畜分布較為均勻,而實際放牧由于草地異質(zhì)性,家畜分布不均勻[12]。因此,真正的問題不是如何測量和計算,而是如何解釋這些數(shù)據(jù)[34]。Scarnecchia[46]認為,利用率的概念是錯誤的。Mace[47]表示根據(jù)當年未采食牧草及下年草地生長量確定草地利用率是無法理解的。
由于體型大小、類型和生產(chǎn)水平等差異,家畜的采食量差異很大,不同放牧家畜的采食量很難量化,尤其是當存在補飼的情況下[48-49]。在載畜量估算中,為了計算方便,通常是把不同種類、不同體型大小、不同年齡的家畜統(tǒng)一轉換為“標準家畜”(standard animal),以方便載畜量的估算。一個標準家畜的采食量稱為一個家畜單位(animal unit,AU)。家畜單位也稱為家畜等價系數(shù)或家畜指數(shù)(animal index)[50]。基于家畜放牧壓力和飼料需求,把某一種類的畜群和標準家畜相比,根據(jù)其與標準家畜單位當量(animal unit equivalent)比值,確定其家畜單位。家畜單位當量用于近似確定家畜對草地的影響程度[50]。家畜單位轉換系數(shù)是特定種類家畜飼草需求量與標準家畜飼草需求量的比值。和家畜單位相關的概念,還有家畜-天 (animal unit day,AUD)、家畜-月 (animal unit month,AUM)、家畜-年(animal unit year,AUY),分別指一個家畜單位放牧1天、1月和1年的牧草采食量[25]。
1.3.1 家畜單位法(牛單位法) 家畜單位的定義不盡相同,最近的研究較為系統(tǒng)的評述了家畜單位的概念[11,51-53]。Scarnecchia和Kothmann[54]綜述了家畜單位這一術語的發(fā)展變化歷程,實踐上應用最多的標準家畜是牛。
早在1907年,美國林業(yè)管理部門就使用“牛-天”表示家畜的日食量,是現(xiàn)在使用的“家畜-天”的最早形式[51]。Sampson[25,55]使用牛單位表示載牧量,5只羊等同于1個牛單位。Pickford[56]把家畜單位定義為成年牛一年的采食量,相當于5只羊或1匹馬的年采食量,等同于家畜-年(AUY)[32]。Stoddart和Smith[57]對家畜單位的定義和Sampson[55]的牛單位一致。12年后,Stoddart和Smith[26]將家畜單位定義為454kg的成年母牛及其犢牛或其他活重達到454kg的家畜,為一個家畜單位。Voisin[58]定義家畜單位是活重為500kg的家畜。Vallentine[59]提出家畜單位是454kg的干母牛的維持或妊娠需求。Robert和Richard(2001)[60]認為一個家畜單位是454kg的成年母牛及其犢牛(未明確犢牛的月齡),日食干草11.8kg或日食量為其體重的2.6%。Scarnecchia和Kothmann[54]定義家畜單位為日食量11.8kg干物質(zhì)/d)。目前,最為廣為接受的家畜單位概念,是活重為454kg的成年母牛及其哺乳犢牛(不超過6個月齡),日食干草11.8kg[32],或日食量占其體重的2.6%[61]。
Holechek等[44]給出了不同種類家畜轉換為家畜單位標準和方法。對牛而言,由于體重差異,每100磅體重增減對應0.1個家畜單位變化[62]。體重小于900磅的牛,其折合家畜單位的計算公式為(體重+100)/1000,體重大于1100磅的牛,其折合家畜單位的計算公式為(體重-100)/1000;超過3個月齡的犢牛直至400磅左右斷奶,由于牧草比母乳更能滿足其營養(yǎng)需求,折合0.3個家畜單位。
1.3.2 羊單位法 目前,采用羊作為標準家畜的國家主要有澳大利亞、新西蘭和中國。澳大利亞的羊單位(dry sheep equivalent,DSE)是體重為45kg的2歲美利奴綿羊的維持飼草需求量[63]。新西蘭目前的載畜量家畜轉換體系也采用羊單位(ewe equivalent,EE),以體重為54.5kg的母羊(帶有一個羔羊)作為基本家畜單位[64]。新西蘭家畜單位最初是Fawcett和Paton[65]提出的家畜單位(stock unit,SU),1965年,Coop[64]以SU 為基礎,基于家畜消化有機物質(zhì)(DOM)需求,提出了標準羊單位(standard ewe),標準羊單位DOM 年需求量為370kg[66],根據(jù)牧草62%的平均消化率,干物質(zhì)需求量為595kg[64]。中國采用綿羊單位[67],1個羊單位是活重為40kg的綿羊及其哺乳羔羊,每天的飼草(青草)需求量為5.0~7.5kg。
1.3.3 熱帶家畜單位法 熱帶草地常用熱帶家畜單位(tropical livestock unit,TLU),指的是活重為250kg的家畜的飼草需求量。熱帶牛單位(tropical cattle unit,TCU)使用的較少,指的是活重為175kg的牛的飼草需求量。還有小反芻家畜單位,指的是活重為25kg的家畜的飼草需求量[68]。
1.3.4 家畜平均體重法 家畜平均體重法(average animal weight,AAW),不考慮家畜的種類、品種,采用統(tǒng)一轉換系數(shù)0.02667,與家畜體重相乘,計算家畜的飼草需求量。這個系數(shù)是根據(jù)母牛及其哺乳犢牛的代謝需求計算得來,一頭母牛及其哺乳犢牛1d的飼草需求量約為其體重的2.667%。Holechek等[44]認為家畜日均采食干物質(zhì)量應約為其體重的2%。依據(jù)體重法估算放牧率時,一般綿羊采食量為其體重的3.0%~3.5%,山羊采食量為其體重的4.0%~4.5%[60]。馬和驢雖然是單胃動物,但其盲腸發(fā)達,單位體重的飼草需求量高于反芻家畜,飼草采食量一般占其體重的3%[31]。雖然,不同研究給出不同的數(shù)字,但大多數(shù)體重法基本上以體重的2.6%計算家畜的飼草需求量。
1.3.5 家畜單位應用上的缺陷 家畜單位具有一定的模糊性,作為評價放牧壓力或家畜采食量的單位,和載畜量這個評價目標多元性的變量并不相容[11]。家畜單位當量僅用以表達一個家畜單位的飼料需求,并不是家畜采食量的單位[51,69],也不是轉化率[52]。家畜采食量是家畜-草地關系及環(huán)境因素的函數(shù),家畜單位僅和影響家畜飼料需求的因素有關,這些因素包括家畜的代謝體重、妊娠與否、泌乳階段,并不包括牧草或環(huán)境特征[51]。家畜單位作為衡量家畜營養(yǎng)需求的單位,不應涉及植物-動物互作關系[51-52,69],這樣才便于嚴格界定放牧水平和載畜量。家畜單位不能既表示家畜的飼料需求,又表示采食量[52]。
1.4.1 家畜代謝體重(metabolic body size,MBS)與家畜折算系數(shù) 家畜單位通常以家畜活重為基礎,在大范圍的草地上,只考慮活重或許是充分的,但基于活重的家畜單位沒有對家畜生長階段、體型等做合理的區(qū)分,也沒有考慮生理階段,如哺乳和犢牛年齡對采食量的影響。即使是同一種類的家畜,由于其品種、生長階段也差異很大。根據(jù)管理經(jīng)驗,家畜泌乳期的蛋白和能量需求要比非泌乳時期高33%和50%,因此計算家畜的飼草需求量,應依據(jù)家畜的代謝體重[70]。但Hutton[71]報道妊娠對家畜的攝食量影響很小。Voisin[58]對家畜單位的定義就是基于代謝體重為基礎。代謝體重糾正了家畜體型對活重的影響,考慮了家畜活重及表面積,可以降低由于體形差異導致計算飼料需求量時的誤差[72]。反芻家畜的代謝體重,是其活重的0.75次方,表達式為MBS=W0.75(式中,MBS是家畜代謝體重,W是家畜活重,kg),根據(jù)代謝體重的家畜單位(AU)表達為AU=(W/450)0.75(W是家畜活重,kg)[73]。
1.4.2 采食重疊及家畜單位 Cook[74]和Smith[75]先后提出混群放牧(mixed livestock grazing)概念并應用于草地管理。草地混群放牧能夠有效提高草地畜產(chǎn)品產(chǎn)量,有利于草地資源的管理和保護,但增大了草地管理的難度和草地維護成本[76]。Holechek等[77]提出了混群放牧提高草地穩(wěn)定性的機理,即家畜混群放牧使得草地植被能夠被均衡利用?;烊悍拍劣捎诩倚蟮倪x擇性采食,能夠充分利用草地的禾本科草類、雜草類和灌木植被,提高了牧草利用效率。與此同時,由于一些種類的飼草并不為家畜喜食,有的牧草對家畜的適口性要好于其他飼草,因此混群放牧存在資源分異和采食重疊現(xiàn)象[50-51,78]。在資源有限的情況下,資源利用重疊常導致競爭,造成資源利用分異。資源分異還能夠解釋為什么采食重疊度高的物種能夠共存[79]。由于物種的協(xié)同進化,采食重疊現(xiàn)象往往只存在于資源充足的生態(tài)系統(tǒng)[80]。同時,采食重疊也和資源缺乏相聯(lián)系[81]。例如,在半干旱森林草地,山羊、綿羊和?;烊悍拍粒捎诓墒持丿B,山羊和牛的資源競爭比較大[78]?;诩倚篌w重或代謝體重的家畜單位轉換系數(shù)沒有考慮采食重疊,計算載畜量存在誤差[28]。因此,混群放牧,計算草地載畜量時,家畜單位的計算應賦予采食重疊一定的權重,采用加權計算法消除采食重疊對計算載畜量帶來的誤差[51,82]。
關于采食重疊,有多種計算指數(shù)[83],Sorenson’s相似指數(shù)[84]是第一個定量計算相似度的指數(shù),也是計算畜群采食重疊的重要方法。Kulcyznski提出的相似度指數(shù)用以計算牧草在家畜采食中的重疊度[85]。Dunbar[86]提出了評估采食重疊程度的簡單指數(shù)。Taylor[87]基于食管瘺技術提出了計算家畜選擇性采食的喜食指數(shù)。Roux[88]制定了計算采食重疊相似度的方法。Schoener指數(shù)[89]用于計算大型偶蹄動物的采食重疊度。關于采食和生態(tài)位重疊,還有生態(tài)位重疊指數(shù)公式[90]、2種及3種家畜間的采食重疊度指數(shù)公式[91]。
考慮了采食重疊,家畜單位調(diào)整公式為[84,92]:
式中,X是家畜單位,DO是采食重疊度(%),M是家畜活重。
根據(jù)采食重疊度家畜單位后,由于引入了植物因子,可能會混淆家畜單位和家畜單位當量的界限,不利于載畜量的計算和分析[51]。但是Hobbs和Carpenter[50]堅持認為,飼草的可獲得性、家畜生理和解剖學特征是家畜的食性及采食重疊的基礎[93]。
關于載畜量計算的文獻較少[94],生態(tài)學家主要基于降水和草地產(chǎn)量[95]、草地利用率、可食牧草比例[47]以及氣候、土壤、植被、家畜和人類的影響[96-97]來估算草地載畜量。最基本的計算方法是,先計算生長季末的草地產(chǎn)量,利用草地利用率等校正系數(shù)校正,除以家畜的年平均需求量,得到草地理論載畜量[47]。Vallentine[41]列舉了初始載畜量法、歷史載畜量法、草地產(chǎn)量估計法、植被利用率法、牧草產(chǎn)量比較、能量需求法、牧草密度法等7種確定載畜量/放牧率的方法。
目前,草地載畜量的確定主要歸為牧草產(chǎn)量法、家畜牧草采食量法等2類方法。牧草產(chǎn)量法,采用放牧或刈割試驗,或者利用降水或土壤濕度估測草地產(chǎn)量,根據(jù)草地可持續(xù)利用率得到校正載畜量[7]。在降水量低于700 mm的草地,水分是草地產(chǎn)量的決定因子,載畜量一般利用降水量結合水分年際變化計算[7]。家畜牧草需要量法,程序和牧草產(chǎn)量法基本一致,但家畜牧草需要量法要考慮家畜的營養(yǎng)需求,因此除了考慮牧草產(chǎn)量外,還應該把牧草質(zhì)量和營養(yǎng)物質(zhì)產(chǎn)量包括在內(nèi)[7]。有學者[98]提出,除了考慮牧草產(chǎn)量、營養(yǎng)物質(zhì)含量外,還應考慮家畜對牧草的轉化效率。草地放牧是一個動態(tài)的過程,因此,草地歷史放牧資料和牧民長期的放牧經(jīng)驗在確定載畜量的地位無可替代[99],不同時期實際放牧家畜數(shù)量、草地利用率、草地趨勢動態(tài)及降水資料對計算載畜量比較重要[5]。Johnston等[38]基于科學基準和當?shù)氐姆拍翆嵺`,提出了評估草地載畜量的生態(tài)框架。
載畜量的確定需要測定草地產(chǎn)量,草地產(chǎn)量測定一般在未放牧草地上進行,與不放牧草地相比,放牧刺激植物葉片的生長,而家畜更喜采食葉片,造成得出的數(shù)據(jù)與實際放牧草地相比存在一些誤差。草地產(chǎn)量一般在生長季末生物量最大的時候進行,而在實際放牧草地上,放牧后的草地牧草存在繼續(xù)生長、枯死、凋落,或被昆蟲和其他家畜采食,這是生長季末測產(chǎn)法無法獲得的。因此,生長季多次測產(chǎn)比生長季末測產(chǎn)得到的數(shù)據(jù)更準確。此外,準確計算草地載畜量,需要多年的測產(chǎn)數(shù)據(jù)。由于在牧草產(chǎn)量、家畜采食量及牧草利用率等測定上,存在著誤差和主觀性[47],使得計算出來的載畜量無法做到準確[30],有時是實際放牧率的3~8倍[47]。
根據(jù)以上參數(shù)估算的草地載畜量還要根據(jù)草地可利用程度(例如地形復雜、生態(tài)敏感區(qū)、自然保護區(qū)等天然屏障使家畜無法或不能采食)及管理需要(如調(diào)整家畜放牧布局)進行調(diào)整,最常用的校正因子有牧草損失率、草地利用率及放牧有效性。調(diào)整后的載畜量(carrying capacity)叫載牧量(grazing capacity)[31]。
牧草的適口性差異很大,由于牧草老化,以及家畜放牧造成的畜蹄踐踏和糞便污染會降低牧草采食率。放牧有效性因素主要有地形坡度、水源距離、植被蓋度等,Holechek[31]首次提出了根據(jù)水源距離、地形坡度調(diào)整載畜量的方法。就坡度而言,坡度0~10%不需要調(diào)整,坡度11%~30%減少載畜量30%,坡度31%~60%減少載畜量60%,坡度大于60%家畜無法采食。和牛相比,山羊和綿羊,由于體型小、靈活、攀爬能力高,能夠更好利用崎嶇地形上的草地[31]。坡度>45°的草地對綿羊而言,利用難度也較高,但坡度<45°的草地計算綿羊的載畜率時不需要調(diào)整[100]。牛的取水范圍大約為3km,而綿羊和山羊由于不是每天需要飲水,可以利用距離水源地3km以外的草地[100-101]。
放牧制度也影響草地載畜量,與連續(xù)放牧相比,輪牧由于使草地具有休養(yǎng)生息的機會,使草地載畜量提高20%,而短期放牧、控制放牧提高載畜量30%~50%。此外,“自上而下”還是“自下而上”的社會經(jīng)濟政策,也是決定載畜量的重要因素[102]。
草地放牧管理的目標是維持家畜牧草需求量和草地飼料供應量的平衡,獲取最大的經(jīng)濟效益[3]。草地基況理論長期以來是放牧管理的主要依據(jù)[103],強調(diào)草食家畜和草地資源之間的生物反饋,主張采取保守放牧策略保證草地生態(tài)系統(tǒng)的穩(wěn)定性。在保守放牧策略之下,草地利用采取輕度至中度的放牧強度[99],草地利用率一般固定在2/3的生態(tài)載畜量水平上[104-105]。但對于降雨量低、降水變率高的牧區(qū),由于氣候的劇烈波動及不可預測性,并不適合采取固定草地利用率。而合理的放牧率是草地生態(tài)系統(tǒng)主要動力[39]。因此,草地利用率應根據(jù)氣候條件、草地類型、草地基況等進行調(diào)整,而不應簡單遵循經(jīng)驗法則。
20世紀50年代,家畜單位的定義只表示家畜活重。直到1974年,美國草地管理學會的家畜單位定義才既包含家畜體重,也包括家畜日采食量[51]。此外,實際生產(chǎn)管理中,很少有畜牧業(yè)從業(yè)者定時監(jiān)測家畜重量進行飼養(yǎng)管理,大多以家畜體況評分(body condition score,BCS)作為管理的依據(jù)[106]。僅僅依靠家畜體重,而不考慮家畜體況來估測牧草需求,常常會帶來一些錯誤[60]。只有家畜體況評分達到5分,依據(jù)家畜體重計算采食量才是可靠的[60]。家畜單位除了考慮代謝體重、采食重疊等因素外,還應考慮家畜體況,并在不同轉換方法之間尋求統(tǒng)一轉換尺度。
草地載畜量計算中,草地產(chǎn)量法受草地類型、植物種類、牧草成熟度等因素的影響,不能反映特定時段和全年草畜平衡的真實狀況和動態(tài)變化[107],估算結果通常與家畜營養(yǎng)需求法估算結果并不一致。因此,應采用2種方法相結合估算草地載畜量,根據(jù)估算數(shù)較低者確定草地載畜量有利于草地可持續(xù)利用和保護。為了調(diào)控草原各生產(chǎn)要素以獲得最優(yōu)的產(chǎn)品和為社會提供最優(yōu)質(zhì)服務[108],實現(xiàn)草地的可持續(xù)利用,確定草地載畜量除了采用數(shù)量化方法外,還應結合放牧經(jīng)驗、草地屬性、歷史氣象資料、草地產(chǎn)量等資料[38],確定特定地區(qū)的草地載畜量[38]。
[1]Bestelmeyer B T,Estell R E,Havstad K M.Big questions emerging from a century of rangeland science and management[J].Rangeland Ecology & Management,2012,65(6):543-544.
[2]Hardin G.The tragedy of the commons[J].Science,1968,162:1243-1248.
[3]Abel N O J.Comments on PDN paper 29aland degradation,stocking rates and conservation policies in the communal rangelands of Botswana and Zimbabwe[A].In:Abel N O J,Blaikie P M.Pastoral Development Network Paper(29d)[M].London:Overseas Development Institute,1990.
[4]李青豐,齊智鑫.草畜平衡管理系列研究(4)—草畜平衡計算軟件開發(fā)[J].草業(yè)科學,2011,28(11):2042-2045.
[5]Galt D,Molinar F,Navarro J,etal.Grazing capacity and stocking rate[J].Rangelands,2000,22(6):7-11.
[6]汪詩平.天然草原持續(xù)利用理論和實踐的困惑[J].草地學報,2006,14(2):188-192.
[7]Leeuw P N D,Tothill J C.The concept of rangeland carrying capacity in sub-saharan Africa:myth or reality[A].In:Behnke R H,Scoones I,Kerven C.Range Ecology at Disequilibrium:New Models of Natural Variability and Pastoral Adaptation in African Savannas[M].London,UK:Overseas Development Institute,1993.
[8]賈幼陵.關于草畜平衡的幾個理論和實踐問題[J].草地學報,2005,13(4):265-268.
[9]董世魁,江源,黃曉霞.草地放牧適宜度理論及牧場管理策略[J].資源科學,2002,24(6):35-41.
[10]Holechek J L,Hess K J.Grazing lands:prices,value,and the future[J].Rangelands,1996,18(3):102-105.
[11]Scarnecchia D L.Concepts of carrying capacity and substitution ratios:a systems viewpoint[J].Journal of Range Management,1990,43(6):553-555.
[12]Hunt L P.Safe pasture utilisation rates as a grazing management tool in extensively grazed tropical savannas of northern Australia[J].The Rangeland Journal,2008,30(3):305-315.
[13]Scurlock J M O,Johnson K,Olson R J.Estimating net primary productivity from grassland biomass dynamics measurements[J].Global Change Biology,2002,8(8):736-753.
[14]Lin H L,Zhang Y J.Evaluation of six methods to predict grassland net primary productivity along an altitudinal gradient in the Alxa rangeland,western Inner Mongolia,China[J].Grassland Science,2013,59(2):100-110.
[15]Lin H L,Zhao J,Liang T G,etal.A classification indices-based model for net primary production(NPP)and potential productivity of vegetation in China[J].International Journal of Biomathematics,2012,5(3):1-12.
[16]Grosso S D,Parton W,Stohlgren T,etal.Global potential net primary production predicted from vegetation class,precipitation,and temperature[J].Ecology,2008,89(8):2117-2126.
[17]Houerou H N L,Bingham R L,Skerbek W.Relationship between the variability of primary production and the variability of annual precipitation in world arid lands[J].Journal of Arid Environments,1988,15(1):1-18.
[18]Rosenzweig M L.Net primary productivity of terrestrial communities:Prediction from climatological data[J].The American Naturalist,1968,102:67-74.
[19]Lieth H.Modeling the primary productivity of the world[A].In:Lieth H,Whittaker R H.Primary Productivity of The Biosphere[M].New York,USA:Springer-Verlag,1975.
[20]Oesterheld M,Sala O E,McNaughton S J.Effect of animal husbandry on herbivore-carrying capacity at regional scal[J].Nature,1992,356(19):234-236.
[21]Behnke R H,Scoones I.Rethinking range ecology:implications for rangeland management in Africa[A].In:Behnke R H,Scoones I,Kerven C.Range Ecology at Disequilibrium,New Models of Natural Variability and Pastoral Adaptation in African Savannas[M].London,UK:Overseas Development Institute and International Institute for Environment and Development,1993.
[22]Fritz H,Duncan P.On the carrying capacity for large ungulates of African savanna ecosystems[J].Proceedings:Biological Sciences,1994,256:77-82.
[23]Hocking D,Mattick A.Dynamic carrying capacity analysis as tool for conceptualizing and planning range management improvements,with a case study from India[A].In:Overseas Development Institute.Pastoral Development Network.Papers,34c[M].London:Overseas Development Institute,1993.
[24]Noy-Meir I.Desert ecosystems:environment and producers[J].Annual Review of Ecology and Systematics,1973,4(1):25-51.
[25]Sampson A W.Range Management:Principles and Practices[M].New York:John Wiley and Sons,Inc.,1952.
[26]Stoddart L A,Smith A D.Range Management(2nd edition)[M].New York:McGraw-Hill Book Co.,1955.
[27]Hedrick D W.Proper utilization—A problem in evaluating the physiological response of plants to grazing use:A Review[J].Journal of Range Management,1958,11(1):34-43.
[28]Heady H F,Child R D.Rangeland Ecology and Management[M].Colorado:WestviewPress,Inc.,1994.
[29]Parker K W,Glendening G E.General Guide to Satisfactory Utilization of The Principal Southwestern Range Grasses[M].Tucson,US:Res.Note 104,Department of Agriculture,F(xiàn)orest Service,Southwestern Forest and Range Experiment Station,1942.
[30]Bartels G B,Perrier G K,Norton B E.The applicability of the carrying capacity concept to Africa[A].In:Leeuw P N,Tothill J C.The Concept of Rangeland Carrying Capacity in Sub-saharan Africa:Myth or Reality[M].London:Pastoral Development Network 29b,Overseas Development Institute,1990.
[31]Holechek J L.An approach for setting the stocking rate[J].Rangelands,1988,10(1):10-14.
[32]Society for Range Management.A Glossary of Terms Used in Range Management:A Definition of Terms Commonly Used in Range Management[M].Denver,Colo:Society for Range Management,1989.
[33]Interagency Technical Reference.Utilization Studies and Residual Measurements[R].Cooperative Extension Service,U.S.D.A.Forest Service,U.S.D.A.Natural Resource Conservation Service,U.S.D.I.Bureau of Land Management,1999.
[34]Heady H F.Methods of determining utilization of range forage[J].Journal of Range Management,1949,2(2):53-63.
[35]Valentine K A.Distance from water as a factor in grazing capacity of rangeland[J].Journal of Forestry,1947,45(10):749-754.
[36]Smith L,Ruyle G,Maynard J,etal.Principles of Obtaining and Interpreting Utilization Data on Rangelands(AZ1375)[M].Tucson,Arizona:Arizona Cooperative Extension,College of Agriculture and Life Sciences,The University of Arizona,2007.
[37]FAO.Guidelines:Land Evaluation for Extensive Grazing[S].Rome:Soil Resources,Management and Conservation Service,Soils Bulletins Series,No.58,1988.
[38]Johnston P W,McKeon G M,Day K A.Objective‘safe’grazing capacities for south-west Queensland Australia:development of a model for individual properties[J].The Rangeland Journal,1996,18(2):244-258.
[39]Holechek J L,Pieper R D,Herbel C H.Range Management:Principles and Practices(6th Edition)[M].New Jersey:Prentice Hall,2010.
[40]Cook C W,Taylor K,Harris L E.The effect of range condition and intensity of grazing upon daily intake and nutritive value of the diet on desert ranges[J].Journal of Range Management,1962,15(1):1-6.
[41]Vallentine J F.Grazing capacity inventory[A].In:Vallentine J F.Grazing Management[M].New York:Academic Press,1990.
[42]陳全功.關鍵場與季節(jié)放牧及草地畜牧業(yè)的可持續(xù)發(fā)展[J].草業(yè)學報,2005,14(4):29-34.
[43]Carter J O,Hall W B,Brook K D,etal.Aussie GRASS:Australian grassland and rangeland assessment by spatial simulation[A].In:Hammer G,Nicholls N,Mitchell C.Applications of Seasonal Climate Forecasting in Agricultural and Natural Ecosystems-The Australian Experience[M].Dordrecht,The Netherlands:Kluwer Academic Press,2000.
[44]Holechek J L,Pieper R D,Herbel C H.Range Management:Principles and Practices[M].Englewood Cliffs,New Jersey:Prentiss-Hall,Inc.,1989.
[45]Westoby M,Walker B,Noy-Meir I.Opportunistic management for rangelands not at equilibrium[J].Journal of Range Management,1989,42(4):266-274.
[46]Scarnecchia D L.Viewpoint:the range utilization concept,allocation arrays,and range management science[J].Journal of Range Management,1999,52(2):157-160.
[47]Mace R.Overgrazing overstated[J].Nature,1991,349:280-281.
[48]Macoon B,Sollenberger L E,Moore J E,etal.Comparison of three techniques for estimating the forage intake of lactating dairy cows on pasture[J].Journal of Animal Science,2003,81:2357-2366.
[49]Moore J E.Practical approaches to estimating pasture intake[A].In:Joost R E,Roberts C A.Nutrient Cycling in Forage Systems[M].Columbia,MO:Center for Pasture Management,Univ.of Missouri,1996.
[50]Hobbs N T,Carpenter L H.Viewpoint:animal-unit equivalents should be weighted by dietary differences[J].Journal of Range Management,1986,39(5):470.
[51]Scarnecchia D L.The animal-unit and animal-unit-equivalent concepts in range science[J].Journal of Range Management,1985,38(4):346-349.
[52]Scarnecchia D L.Viewpoint:animal-unit-equivalents cannot be meaningfully weighted by indices of dietary overlap[J].Journal of Range Management,1986,39(5):471.
[53]Scarnecchia D L.The relationship of stocking intensity and stocking pressure to other stocking variables[J].Journal of Range Management,1985,38(6):558-559.
[54]Scarnecchia D L,Kothmann M M.A dynamic approach to grazing management terminology[J].Journal of Range Management,1982,35(2):262-264.
[55]Sampson A W.Range and Pasture Management[M].New York:John Wiley and Sons,Inc.,1923.
[56]Pickford G D.Range survey methods in western United States[J].Herbage Reviews,1940,8:1-12.
[57]Stoddart L A,Smith A D.Range Management(1st ed)[M].New York:McGraw-Hill Book Co.,1943.
[58]Voisin A.Grass Productivity[M].New York:Philosophical Library,Inc.,1959.
[59]Vallentine J F.An improved AUM for range cattle[J].Journal of Range Management,1965,18(6):346-348.
[60]Robert K L,Richard V M.Stocking Rate:The Key Grazing Mmanagement Fecision.Texas Agricultural Extension Service,L-5400[S].The Texas A&M University System,2001.
[61]Meyer T L.Estimating Livestock Forage Demand:Defining the Animal Unit[D].Lincoln:University of Nebraska-Lincoln,Master of Science,2010.
[62]Waller S S,Moser L E,Anderson B.EC86-113:A Guide for Planning and Analyzing a Year-Round Forage Program[R].Historical Materials from University of Nebraska-Lincoln Extension,1986.
[63]Colin McLaren A.Agriculture Note,AG0590:Dry Sheep Equivalents for Comparing Different Classes of Livestock[R].State of Victoria,Department of Primary Industries,1997.
[64]Coop I E.A review of the ewe equivalent system[J].New Zealand Agricultural Science,1965,1(3):13-18.
[65]Fawcett E J,Paton W N.A review of livestock production in New Zealand during the past 26seasons(1901/2to 1926/27),based on standard values and units[A].The New Zealand Year-Book (Vol.37)[M].Wellington:Government Printer,1929.
[66]Parker W J.Standardisation between livestock classes:the use and misuse of the stock unit system[J].Proceedings of the New Zealand Grassland Association,1998,60:243-248.
[67]任繼周.草原調(diào)查與規(guī)劃[M].北京:中國農(nóng)業(yè)出版社,1985.
[68]Jahnke H E.Livestock Production Systems and Livestock Development in Tropical Africa[M].Kiel:Kieler Wissenschaftsverlag Vauk,1982.
[69]Scarnecchia D L,GasLins C T.Modeling animal-unit equivalents for beef cattle[J].Agricultural Systems,1987,23(1):19-26.
[70]Edwards P J.Terms describing aspects of vegetation and its management[A].In:Tainton N M.Veld and Pasture Management in South Africa[M].Pietermaritzburg:Shuter and Shooter(Pty)Ltd.,1981.
[71]Hutton J B.The effect of lactation on intake in the dairy cow[J].Proceedings of the New Zealand Society of Animal Production,1963,23:39-52.
[72]Manske L L.Animal Unit Equivalent for Beef Cattle Based on Metabolic Weight,Range Management Report DREC 98-1020[R].Dickinson,ND:NDSU Dickinson Research Extension Center,1998.
[73]Kleiber M.The Fire of Life:An Introduction to Animal Energetics[M].New York:R.E.Krieger.,1975.
[74]Cook C W.Common use of summer range by sheep and cattle[J].Journal of Range Management,1954,7(1):10-13.
[75]Smith A D.Determining common use grazing capacities by application of the key species concept[J].Journal of Range Management,1965,18(4):196-201.
[76]Ospina E.A proposal for a research agenda on the economics of multispecies grazing[A].In:Baker F H,Jones R K.Proceedings of A Conference on Multispecies Grazing[M].Morrilton,Ark.,USA:Winrock International Institute for Agricultural Development,1985.
[77]Holechek J L,Pieper R D,Herbel C H.Range Management:Principles and Practices[M].USA:Prentice-Hall,1995.
[78]Squires V R.Dietary overlap between sheep,cattle,and goats when grazing in common[J].Journal of Range Management,1982,35(1):116-119.
[79]MacArthur R H,Wilson E O.The Theory of Island Biogeography[M].Princeton:Princeton University Press,1967.
[80]Voeten M M,Prins H H T.Resource partitioning between sympatric wild and domestic herbivores in the Tarangire region of Tanzania[J].Oecologia,1999,120(2):287-294.
[81]Gordon I J,Illius A W.Resource partitioning by ungulates on the Isle of Rhum[J].Oecologia,1989,79(3):383-389.
[82]Hobbs N T,Swift D M.Estimates of habitat carrying capacity incorporating explicit nutritional constraints[J].Journal of Wildlife Management,1985,49(3):814-822.
[83]Horn H S.Measurement of“Overlap”in comparative ecological studies[J].The American Naturalist,1966,100:419-424.
[84]Flinders J T,Conde L.An Alternate Method of Calculating Animal-unit-equivalents Helps Ranchers and Wildlife(Abstract)[C].San Diego,California:The 33rd Annual Meeting,Society for Range Manage,1980.
[85]Celaya R,Oliván M,F(xiàn)erreira L M M,etal.Comparison of grazing behaviour,dietary overlap and performance in non-lactating domestic ruminants grazing on marginal heathland areas[J].Livestock Science,2007,106(2):271-281.
[86]Dunbar R I M.Competition and niche separation in a high altitude herbivore community in Ethiopia[J].African Journal of E-cology,1978,16(3):183-199.
[87]Taylor C A.The Botanical Composition of Cattle Diets on a 7-pasture High-intensity Low-frequency Grazing System[D].Texas:Texas A&M Univerisity,1973.
[88]Roux P W.The Autecology of Tetrachne Dregei Nees[D].Pietermaritzburg:University of Natal,1968.
[89]Schoener T W.The Anolis lizards of Bimini:resource partitioning in a complex fauna[J].Ecology,1968,49(4):704-726.
[90]Abrams P.Some comments on measuring niche overlap[J].Ecology,1980,61(1):44-49.
[91]Pianka E R.The structure of lizard communities[J].Annual Review of Ecology and Systematics,1973,4:53-74.
[92]Botha P,Blom C D,Sykes E,etal.A comparison between the diets of small and large stock on mixed Karoo veld[J].Proceedings of the Annual Congresses of the Grassland Society of Southern Africa,1983,18(1):101-105.
[93]Hanley T A.The nutritional basis for food selection by ungulates[J].Journal of Range Management,1982,35(2):146-151.
[94]Biot Y.How long can high stocking densities be sustained[A].In:Behnke R H,Scoones I,Kerven C.Range Ecology at Disequilibrium New Models of Natural Variability and Pastoral Adaption in African Savannas[M].London:Overseas Development Institute,1993.
[95]Phillipson J.Rainfall,primary production and‘carrying capacity’of Tsavo National Park(East),Kenya[J].African Journal of Ecology,1975,13(3-4):171-201.
[96]Coughenour M B,Singer F J.Yellowstone elk population responses to fire:A comparison of landscape carrying capacity and spatial-dynamic ecosystem modeling approaches[A].In:Greenlee J M.The Ecological Implications of Fire in Greater Yellowstone.Proceedings of The 2nd Biennial Conference on The Greater Yellowstone Ecosystem[C].Fairfield,Washington:International Association of Wildland Fire,1996.
[97]Coughenour M B.Elk carrying capacity on Yellowstone’s northern elk winter range:Preliminary modeling to integrate climate,landscape,and nutritional requirements[A].In:Despain D,Hampre R.Plants and Their Environments:Proceedings of the first Biennial Scientific Conference on the Greater Yellowstone Ecosystem[C].Denver,Colorado,USA:US Department of the Interior,National Park Service,Rocky Mountain Region,Yellowstone National Park,1994.
[98]趙彥光,洪瓊花,謝萍,等.云貴高原石漠化地區(qū)人工草場營養(yǎng)價值評價研究[J].草業(yè)學報,2012,21(1):1-9.
[99]Holechek J L,Gomez H,Molinar F,etal.Grazing studies:what we’ve learned[J].Rangelands,1999,21(2):12-16.
[100]McDaniel K C,Tiedeman J A.Sheep use on mountain winter range in New Mexico[J].Journal of Range Management,1981,34(2):102-104.
[101]Squires V R.Distance to water as a factor in performance of livestock on arid and semiarid rangelands[A].In:Mayland H F.Water-Animal Relations Proceedings[M].Kimberley,Idaho:Water-Animal Relations Symposium Committee,1973.
[102]Sheehy D P,Thorpe J,Kirychuk B.Rangeland,livestock and herders revisited in the northern pastoral region of China[A].In:Bedunah D J,McArthur E D,F(xiàn)ernandez-Gimenez M C.Rangelands of Central Asia:Proceedings of the Conference on Transformations,Issues,and Future Challenges(2004January 27;Salt Lake City,UT.Proceeding RMRS-P-39)[M].Fort Collins,CO U.S.:Department of Agriculture,F(xiàn)orest Service,Rocky Mountain Research Station,2006.
[103]Briske D D,F(xiàn)uhlendorf S D,Smeins F E.Vegetation dynamics on rangelands:a critique of the current paradigms[J].Journal of Applied Ecology,2003,40(4):601-614.
[104]Campbell B M,Dore’D,Luckert M,etal.Economic comparisons of livestock production in communal grazing lands in Zimbabwe[J].Ecological Economics,2000,33:413-438.
[105]Sandford S,Scoones I.Opportunistic and conservative pastoral strategies:Some economic arguments[J].Ecological Economics,2006,58(1):1-16.
[106]Subcommittee on Beef Cattle Nutrition,Committee on Animal Nutrition,National Research Council NRC.Nutrient Requirements of Beef Cattle,Seventh Revised Edition[M].Washington D.C.:National Academy Press,2000.
[107]楊博,吳建平,楊聯(lián),等.中國北方草原草畜代謝能平衡分析與對策研究[J].草業(yè)學報,2012,21(2):187-195.
[108]張英俊,楊高文,劉楠,等.草原碳匯管理對策[J].草業(yè)學報,2013,22(2):290-299.