基于構件特征的內蒙古典型草原植物羊草個體地上生物量估算

胡靜，侯向陽，薩茹拉，郭豐輝，丁勇*

(1.中國農業科學院草原研究所，內蒙古 呼和浩特 010010；2.中國農業科學院研究生院，北京 100081；

3.農業部牧草資源與利用重點實驗室，內蒙古 呼和浩特 010010)

基于構件特征的內蒙古典型草原植物羊草個體地上生物量估算

胡靜1,2，侯向陽1,3，薩茹拉1，郭豐輝1,2，丁勇1,3*

(1.中國農業科學院草原研究所，內蒙古 呼和浩特 010010；2.中國農業科學院研究生院，北京 100081；

3.農業部牧草資源與利用重點實驗室，內蒙古 呼和浩特 010010)

摘要：生物量是評價草原生態系統生產力的重要指標，草原主要植物個體地上生物量估算模型的建立，可為采用非刈割方法獲取較準確的草原地面數據資料提供新方法。通過分析內蒙古典型草原植物羊草的主要構件特征(葉片數、葉長寬積、莖粗、莖高和株高)與其個體地上生物量之間的關系，比較并篩選出能較好預測羊草個體生物量變化的構件指標，采用Minitab軟件的回歸分析方法，分別利用6月、7月和8月采集的90株和總體270株羊草植物樣本，建立羊草不同生長階段(生長初期、生長旺盛期和生長末期)以及總體樣本的個體地上生物量估算方程，并對實測值與模擬值進行t檢驗分析。結果表明，葉長寬積和株高相對于其他構件信息能夠更好地描述羊草個體地上生物量，方程Y=-0.193+0.009X2+0.011X5(X2為葉長寬積，X5為株高，P<0.01，R2=0.854)為最優羊草個體生物量估算模型。

關鍵詞：羊草；構件特征；個體地上生物量；典型草原

DOI：10.11686/cyxb2015053http://cyxb.lzu.edu.cn

胡靜，侯向陽，薩茹拉，郭豐輝，丁勇. 基于構件特征的內蒙古典型草原植物羊草個體地上生物量估算. 草業學報, 2015, 24(8): 211-217.

Hu J, Hou X Y, Sa R L, Guo F H, Ding Y. Estimating individual above-ground biomass ofLeymuschinensisplants. Acta Prataculturae Sinica, 2015, 24(8): 211-217.

收稿日期：2015-01-27；改回日期：2015-03-18

基金項目：國家重點基礎研究發展計劃(973計劃)項目(2014CB138805)，國家自然科學基金項目(71103185,71311120089)，中央公益性科研院所基本科研業務費項目(1610332015001)，國家科技支撐計劃項目(2012BAD13B07)和中國農業科學院科技創新工程(CAAS-ASTIP-IGR 2015-05)資助。

作者簡介：胡靜(1986-)，女，內蒙古呼和浩特人，在讀博士。E-mail：hu1129@sina.com

通訊作者*Corresponding author. E-mail: dingyong228@126.com

Estimating individual above-ground biomass ofLeymuschinensisplants

HU Jing1,2, HOU Xiang-Yang1,3, SA Ru-La1, GUO Feng-Hui1,2, DING Yong1,3*

1.InstituteofGrasslandResearchofChineseAcademyofAgriculturalScience,Hohhot010010,China; 2.GraduateSchoolofChineseAcademyofAgriculturalScience,Beijing100081,China; 3.KeyLaboratoryofGrasslandResourcesandUtilization,MinistryofAgriculture,Hohhot010010,China

Abstract:Biomass is an important component of grassland ecosystem productivity. Models able to describe the productivity of the main plant species in grassland offer the potential to estimate aboveground biomass without physical cutting. In this research the relationships among plant characteristics (leaf number, leaf length×leaf width, stem diameter, stem height and plant height) and individual above-ground biomass was analyzed, subsequently those characteristics able to best predict changes of individual above-ground biomass were selected for modelling. Ninety Leymus chinensis plants were sampled in June (early growing season), July (middle growing season) and August (late growing season), respectively, and total of 270 observations were used to build models able to estimate individual plant above-ground biomass. Equations were obtained using regression analysis using Minitab 16.0, and T-Test analysis was used to compare observed and predicted values. The results showed that leaf length×leaf width and plant height could describe individual above-ground biomass of L. chinensis plants; model Y=-0.193+0.009X2+0.011X5(X2indicated leaf length×leaf width，X5indicated plant height， P<0.01， R2=0.854).

Key words:Leymus chinensis; module characteristics; individual above-ground biomass; typical steppe

羊草(Leymuschinensis)草原是歐亞大陸草原的特有草地類型，集中分布于中國東北平原和內蒙古高原東部，是內蒙古典型草原植物群落的建群種或優勢種[1]。羊草為賴草屬根莖型多年生禾草，莖稈細嫩，葉量豐富，營養價值高，具有較好的耐鹽堿、耐旱、耐踐踏性，是一種優良的牧草資源[2]。生物量是評價草原生態系統生產力的重要指標，準確的測定草地植物生物量不僅是監測草原環境變化的依據，也是深入開展生態學、生理學及相關學科研究的基礎[3]。羊草生物量變化是羊草草原生態系統群落結構優劣和功能高低的最直接體現，在大量的放牧[4-5]、刈割[6-7]、氣候變化[8-9]、土壤營養元素變化[10-11]等相關研究中均有涉及。個體生物量是構成種群生物量和群落生產力的基礎，個體地上生物量常與其主要構件的形態特征密切相關。與20世紀80年代相比，草原生產力持續衰減，嚴重地區可以下降60%~80%[12]。安淵等[13]研究發現草地退化后植被發生明顯的矮化，同時伴隨著生物量的減少。2000年之后，隨著我國草原退化形勢日趨嚴峻，很多研究以過度放牧為背景，闡述了羊草構件主要形態指標的變化，并明確指出，植物形態變化是生物量降低的主要原因，如王煒等[14-15]、Gallacher和Hill[16]的研究均認為在過度放牧下，草原植物性狀整體表現為“小型化”(individual miniaturization)，植株高度、構件生物量較未退化樣地正常植株下降30%~80%，隨著放牧脅迫增強，植物首先采取高度和生物量降低的適應策略。鐘夢瑩等[17]研究發現，在過度放牧干擾下羊草發生矮化現象，具體表現為株高明顯降低，且葉片、根長、根節長減少，這也是導致羊草個體生物量明顯下降的主要原因。李西良等[18]在呼倫貝爾草原的研究也表明，羊草莖高、株高、葉面積等對放牧響應極為敏感，放牧干擾會造成羊草株高和個體地上生物量顯著降低。

可見，羊草植物形態特征與羊草地上生物量之間存在著明顯的相關關系，但這種關系多為定性描述，仍缺乏對植物個體地上生物量與其主要易測形態指標的定量關系的研究。在科學研究和生產實踐中，常以破壞性較大的樣方法來進行草原生物量或生產力的地面調查和動態監測。草原主要植物個體地上生物量估算模型的建立，可為采用非刈割方法獲取較準確的草原地面數據資料提供重要的方法和途徑。本研究作為探索性研究，以內蒙古典型草原優勢植物羊草作為研究對象，通過對植物葉片數、葉長寬積、莖高、莖粗、株高5個主要構件與個體地上生物量的分析，構建一系列多元線性回歸模型，并進行模型優化和驗證，希冀該研究能夠為草原植物個體生物量和群落生產力的模擬研究提供新的思路和啟示。

1材料與方法

1.1　研究區域概況

研究區域位于內蒙古錫林郭勒盟錫林浩特市朝克烏拉蘇木中國農業科學院草原研究所草原生態系統保護和可持續利用研究與示范基地，地理位置系43°26′-44°08′ N，116°04′-117°05′ E，平均海拔1100 m。該區域屬溫帶半干旱草原氣候，年平均溫度0.5~1℃，年平均降水350 mm左右，其中60%~80%的降水量集中在5-8月，年蒸發量1600~1800 mm，約為降水量的4~5倍。2014年生長季前半段較為濕潤，但溫度較低，對植物生長產生了一定影響，自7月下旬開始，發生較為嚴重的干旱(圖1)，對植物的生長產生較為嚴重的影響，致使部分植物提前枯黃。試驗區自2007年開始禁牧，只作為割草利用，采樣區植被類型為地帶性典型草原植被，優勢植物為羊草，主要植物包括大針茅(Stipagrandis)、糙隱子草(Cleistogenessquarrosa)、米氏冰草(Agropyronmichnoi)和黃囊苔草(Carexkorshinskyi)等，土壤為鈣栗土。

圖1　2014年研究區4月至10月降雨量與溫度月動態(錫林浩特氣象站)Fig.1　The dynamics of precipitation and temperature from April to October of study area in 2014

1.2　研究方法

試驗設置3個采樣區作為重復，分別在2014年6月12日(生長初期)、7月12日(生長旺盛期)和8月12日(生長末期)進行取樣，隨機選取單株羊草齊地刈割，每個采樣區每月采集羊草樣品40株，合計120株，全年共采集樣品360株。采集后的樣品帶回實驗室對株高、齊地刈割部位的莖粗、葉片數、葉長、葉寬和莖高等指標進行測量和計數，后將樣品放入烘箱，65℃烘48 h至恒重，稱重后獲取地上個體生物量數據。

隨機選取270株(6月、7月、8月各90株)羊草作為個體生物量估算模型的基礎數據，通過Minitab軟件建立個體生物量與各構件(葉片數、葉長寬積、莖粗、莖高、株高)數量特征的多元線性回歸模型，并對羊草構件特征進行篩選以確定最優模型。本文中葉長寬積為每株羊草總葉片的累加值。通過剩余90株(6月、7月、8月各30株)羊草對模型進行驗證，模型驗證采用Independent-Samples T-Text法對模擬值和實測值進行比較。

2結果與分析

2.1　羊草地上個體生物量與主要構件的相關性分析

選取常規易測且具有相似變化規律的羊草主要構件作為估算模型的基礎變量，研究假定各構件對羊草個體地上生物量的貢獻值相同，由此可建立羊草個體地上生物量的多元線性估算模型為：

Y=β0+β1X1+β2X2+β3X3+β4X4+β5X5+μ

式中，X1為葉片數，X2為葉長寬積，X3為莖粗，X4為莖高，X5為株高，μ為隨機變量，β為變量系數。

利用Minitab軟件的相關分析功能，對羊草地上個體生物量、葉長寬積、莖粗、莖高和株高進行相關性分析，結果表明，羊草葉長寬積與羊草個體地上生物量的相關系數最高，是較適合的預測變量，二者間的R2(判定系數)為0.935，表明羊草地上個體生物量的變異性有93.50%可以通過葉生物量的增加或減少來解釋。莖高和株高對羊草地上個體生物量的影響也不可忽視，相關系數也較高。通常習慣將相關系數較小的變量予以剔除，以確保模型的準確性。本研究中莖粗與羊草地上個體生物量的相關系數較低，而葉片數為離散變量，為了更好的擬合模型，將莖粗和葉片數這兩個變量剔除(圖2)。

圖2　羊草地上個體生物量與

因此羊草地上個體生物量估算模型修正為：

Y=β0+β2X2+β4X4+β5X5+μ

式中，X2為葉長寬積，X4為莖高，X5為株高，μ為隨機變量，β為變量系數。

2.2　參數估算

通過Minitab中的回歸分析進行參數估算，其中T是統計量的值，即描述樣本特征的值，T值的絕對值越大，表示羊草個體地上生物量和各構件的關系越密切。P值為t檢驗(假設系數為0時，自變量和因變量無關)的判定值，P值小于0.05時，該假設不成立。因此T值越大，P值越小，則方程中的參數值越好。S表示殘差的標準差，其值越小越好。R2是回歸模型誤差占總誤差的百分比，取值越接近1，表明回歸模型與數據吻合的越好。R2(adj)是經過Minitab回歸分析后的調整值，R2和R2(adj)越接近，表明回歸模型越可靠。

從表1可以看出，羊草的葉長寬積、莖高和株高能夠較好的描述個體地上生物量，特別是在樣本量較大(3月總樣本)的情況下，方程參數的P值均為顯著水平，S較小，R2較大，且R2和R2(adj)相接近，表明羊草個體生物量估算模型具有較好的擬合性(表1)。

表1　羊草個體地上生物量估算模型

2.3　模型優化

方程參數的估算結果表明，6月、7月、8月和總樣本的多元回歸模型中，莖高參數β4的P值均大于0.05，t檢驗結果不顯著，可以考慮剔除該項。模型中莖高和株高均是與高度有關的構件特征，相比較來說，羊草株高與個體地上生物量的相關關系要好于莖高，且在試驗操作時，株高的測量較方便，因此剔除莖高對羊草個體地上生物量估算模型進行進一步優化。

優化后的羊草地上個體生物量估算模型修正為：

Y=β0+β2X2+β5X5+μ

式中，X2為葉長寬積，X5為株高，μ為隨機變量，β為變量系數。

方程參數如表2所示。優化后各估算模型中所有參數的P值均達到顯著水平，S和R2未有明顯的波動，說明莖高的剔除有助于提高方程自變量系數的準確性。因此在本研究中，選擇羊草葉長寬積和株高這兩個構件特征值作為地上個體生物量的預測因子。6月至8月羊草個體地上生物量估算模型的R2呈逐漸增加趨勢，這是由于羊草各構件生長趨于穩定，數據變異性逐漸減小所致，8月的構件特征數據能夠更好的描述羊草個體地上生物量。同時加大樣本量也有助于提高模型的準確性(表2)。

表2　羊草個體地上生物量估算模型優化

2.4　模型驗證

通過剩余的90株羊草樣品(6月、7月、8月各30株)對以上模型進行驗證，以t檢驗分析對羊草個體地上生物量的實測值和模擬值進行比較，如表3所示，模擬值與實測值之間t檢驗結果均大于0.05，差異不顯著，表明羊草個體地上生物量估算模型較可靠，可以對實際情況進行模擬(表3)。

表3　羊草個體地上生物量估算模型驗證

3討論

本研究所獲取的羊草主要構件數量特征在不同的生長階段均呈現出了相同的變化趨勢，表現為6月生長初期<7月生長旺盛期>8月生長末期，這與植物體內含水量隨著氣候發生波動有關[19-20]，同時也可能是因為羊草在抽穗后，將更多的營養投入到生殖構件中，營養構件中的投入就會相應的減少所致[21]。因此在建立生物量估算模型時，要根據實際情況選取適當的取樣時間。對比優化后羊草地上個體生物量估算模型的判定系數R2，可知8月的樣品可以得到較好的地上個體生物量估算模型。如果為連續取樣，則樣品量較大時方程的擬合性好。通過本研究，我們認為應該進行多年的監測取樣，通過對大量樣本的研究，能夠更清晰地刻畫出羊草植物主要構件特征在生長季不同階段的變化趨勢和規律，以避免極端氣候事件對研究結果的影響。

對羊草構件特征的研究表明，植物構件與其生物量之間存在著一定的相關關系。從羊草地上個體生物量估算模型建立的過程中可以看出，并不是構件疊加越多，得到的方程擬合性越好，有的構件對于個體生物量的作用存在冗余現象，剔除不但不影響模型的準確性，反而會使模型變得簡化和優化。另外，從植物構件特征可以繼續對模型進行擴展，如羊草葉長寬積可以用來計算葉面積，李亞軍等[22]通過分析羊草葉面積和葉長寬積之間的相關關系，并與數字圖像法相比較，得出羊草葉面積的計算公式為：A=k(L×W)，k=0.6555，式中，A代表葉面積，k為葉長和葉寬乘積的修正系數，L為葉長，W為葉寬。在此次取樣中雖未對葉面積進行測定，但在未來的研究中可以將葉面積作為變量參與模型的構建，并結合光譜技術[23]，將單株的植物模型擴展到群落水平。本研究為多元線性回歸模型的建立過程提供了一種篩選適合因子方式，不僅適用于生物量模型，也適用于其他多元線性回歸模型的優化。

研究植物個體地上生物量模型模擬，具有重要的實踐價值。草原植物多樣性豐富，每種植物主要構件形態特征與個體生物量間是否存在普遍的相關性仍值得深入研究。目前植株的形態研究在農作物方面較多，并建立了多種作物生長模型[24]，而關于草原主要植物的相關研究還較少。植物模型的研究已趨于多維化和可視化，但在未知領域基礎數據的研究仍必不可少。在以往的草地監測和測產中，傳統的樣方法具有一定的破壞性，加之草原生態系統的異質性，往往使動態監測的結果存在著一定的誤差，即使通過多點重復，也難以使誤差消除。通過模型對植物生物量進行估算，一方面可以實現保護性測定，另一方面，通過定點、定樣和定株監測估產，可有效提高測定結果的準確性。

4結論

內蒙古典型草原植物羊草葉長寬積和株高相對于葉片數、莖高和莖粗能夠更好地描述其個體地上生物量的變化；通過Minitab軟件的建模和模型優化處理，方程Y=-0.193+0.009X2+0.011X5(X2為葉長寬積，X5為株高，P<0.01，R2=0.854)為最優羊草個體生物量估算模型。

References：

[1]Pan Q M, Bai Y F, Han X G,etal. Effects of nitrogen additionson aLeymuschinensispopulation in typical steppe of Inner Mongolia. Acta Phytoecologica Sinica, 2005, 29(2):311-317.

[2]Li L H, Wang Q B, Bai Y F,etal. Soil respiration of aLeymuschinensisgrassland stand in the Xilin river basin as affected by over-grazing and mate. Acta Phytoecologica Sinica, 2000, 24(6):680-686.

[3]Dong X Y, Fu H, Li X D,etal. Effects on plant biomass and CNP contents of plants in grazed and fenced steppe grasslands of the Loess Plateau. Acta Prataculturae Sinica, 2010, 19(2):175-182.

[4]Wang R Z, Ripley E A. Effects of grazing on aLeymuschinensisgrassland on the Songnen plain of north-eastern China. Journal of Arid Environments, 1997, 36(2):307-318.

[5]Wang Y H, He X Y, Zhou G S. Study on the responses ofLeymuschinensissteppe to grazing in Songnen plain. Acta Agrestia Sinica, 2002, 10(1):45-49.

[6]Baoyin T G T, Li Y H, Bao Q H,etal. Effects of mowing regimes and climate variability on hay production ofLeymuschinensisgrassland in northern China. The Rangeland Journal, 2014, 36(6):593-600.

[7]He F, Li X L, Wan L Q. Influence of rainfall in growing season and elipping intensity on aboveground biomass ofLeymuschinensissteppe. Pratacultural Science, 2009, 26(4):28-32.

[8]Xu Z Z, Zhou G S. Combined effects of water stress and high temperature on photosynthesis, nitrogen metabolism and lipid peroxidation of a perennial grassLeymuschinensis. Planta, 2006, 224(5):1080-1090.

[9]Wang R Z, Gao Q. Climate-driven changes in shoot density and shoot biomass inLeymuschinensis(Poaceae) on the North-east China Transect (NECT). Global Ecology and Biogeography, 2003, 12(3):249-259.

[10]Ma L, Yuan S, Guo C,etal. Carbon and nitrogen dynamics of nativeLeymuschinensisgrasslands along a 1000 km longitudinal precipitation gradient in northeastern China. Biogeosciences, 2014, 11(24):7097-7106.

[11]He N P, Yu Q, Wu L,etal. Carbon and nitrogen store and storage potential as affected by land-use in aLeymuschinensisgrassland of northern China. Soil Biology and Biochemistry, 2008, 40(12):2952-2959.

[12]Han G D, Jiao S Y, Ao D,etal. Effects of plant species diversity and productivity under different stocking rates in theStipabrevifloraGriseb. desert steppe. Acta Ecologica Sinica, 2007, 27(1):182-188.

[13]An Y, Li B, Yang C,etal. Influence of grazing rate on population structure ofStipagrandis. Acta Phytoecologica Sinica, 2002, 26(2):163-169.

[14]Wang W, Liang C Z, Liu Z L,etal. Analysis of the plant individual, behaviouring the degradation and restoring succession in steppe community. Acta Phytoecologica Sinica, 2000, 24(3):268-274.

[15]Wang W, Liang C Z. Mechanism of degradation succession inLeymuschinensis+Stipagrandissteppe community. Acta Phytoecologica Sinica, 2000, 24(4):468-472.

[16]Gallacher D, Hill J. Effects of camel grazing on the ecology of small perennial plants in the Dubai (UAE) inland desert. Journal of Arid Environments, 2006, 66(4):738-750.

[17]Zhong M Y, Fan Q L, Zhang Y J,etal. Morphological plasticity and the biomass allocation models of miniaturizedLeymuschinensis. Acta Agrestia Sinica, 2013, 21(2):260-264.

[18]Li X L, Hou X Y, Wu X H,etal. Plastic responses of stem and leaf functional traits inLeymuschinensisto long-term grazing in a meadow steppe. Acta Phytoecologica Sinica, 2014, 38(5):440-451.

[19]Wullschleger S D, Hanson P J, Tschaplinski T J. Whole-plant water flux in understory red maple exposed to altered precipitation regimes. Tree Physiology, 1998, 18(2):71-79.

[20]Savi T, Andri S, Nardini A. Impact of different green roof layering on plant water status and drought survival. Ecological Engineering, 2013, 57(8):188-196.

[21]Liu P Y, Zhang Q L, Yang Y F. Biomass structure and quantitative relationship models of modules in clonal population ofPuccinilliachinampoensisin Songnen plain. Chinese Journal of Applied Ecology, 2004, 15(4):543-548.

[22]Li Y J, Ji S Q, Zhang S W,etal. The determination method of leaf area onLeymuschinensis. Journal of Anhui Agricultural Sciences, 2009, 37(15):6819-6851.

[23]Lv X D, Wang J G, Sun Q Z,etal. Research on the hyperspectral remote sensing estimation models for the fresh yield of alfalfa grassland. Acta Prataculturae Sinica, 2014, 23(1):84-91.

[24]Lai Q, Li Q F, Mo R G,etal. A study on factors affecting moisture content measurement and ratio of forage fresh weight to dry weight. Chinese Journal of Grassland, 2008, 30(4):73-77.

參考文獻：

[1]潘慶民, 白永飛, 韓興國, 等. 氮素對內蒙古典型草原羊草種群的影響. 植物生態學報, 2005, 29(2):311-317.

[2]李凌浩, 王其兵, 白永飛, 等. 錫林河流域羊草草原群落土壤呼吸及其影響因子的研究. 植物生態學報, 2000, 24(6):680-686.

[3]董曉玉, 傅華, 李旭東, 等. 放牧與圍封對黃土高原典型草原植物生物量及其碳氮磷貯量的影響. 草業學報, 2010, 19(2):175-182.

[5]王玉輝, 何興元, 周廣勝. 放牧強度對羊草草原的影響. 草地學報, 2002, 10(1):45-49.

[7]何峰, 李向林, 萬里強. 生長季降水量和刈割強度對羊草群落地上生物量的影響. 草業科學, 2009, 26(4):28-32.

[12]韓國棟, 焦樹英, 敖登, 等. 短花針茅草原不同載畜率對植物多樣性和草地生產力的影響. 生態學報, 2007, 27(1):182-188.

[13]安淵, 李博, 楊持, 等. 不同放牧率對大針茅種群結構的影響. 植物生態學報, 2002, 26(2):163-169.

[14]王煒, 梁存柱, 劉鐘齡, 等. 草原群落退化與恢復演替中的植物個體行為分析. 植物生態學報, 2000, 24(3):268-274.

[15]王煒, 梁存柱. 羊草+大針茅草原群落退化演替機理的研究. 植物生態學報, 2000, 24(4):468-472.

[17]鐘夢瑩, 樊青麗, 張亞軍, 等. 矮化羊草形態可塑性及生物量分配模式研究. 草地學報, 2013, 21(2):260-264.

[18]李西良, 侯向陽, 吳新宏, 等. 草甸草原羊草莖葉功能性狀對長期過度放牧的可塑性響應. 植物生態學報, 2014, 38(5):440-451.

[21]劉佩勇, 張慶靈, 楊允菲. 松嫩平原朝鮮堿茅無性系種群構件生物量結構及相關模型分析. 應用生態學報, 2004, 15(4):543-548.

[22]李亞軍, 紀澍琴, 張守偉, 等. 羊草葉面積的測定方法. 安徽農業科學, 2009, 37(15):6819-6851.

[23]呂小東, 王建光, 孫啟忠, 等. 苜蓿人工草地高光譜遙感估產模型的研究. 草業學報, 2014, 23(1):84-91.

[24]來強, 李青豐, 莫日根, 等. 影響牧草含水量測定以及牧草干鮮比的主要因素. 中國草地學報, 2008, 30(4):73-77.