王舒悰，顧　琪，陳霜霜，趙劉杰，焦月瀟，王福升

(南京林業(yè)大學(xué) 南方現(xiàn)代林業(yè)協(xié)同創(chuàng)新中心，南京林業(yè)大學(xué) 竹類研究所，江蘇 南京 210037)

地被竹的枝葉稠密，根系發(fā)達，對土壤的需求不高，具有抗鹽堿的優(yōu)點[1-3]，繁殖能力強且病蟲害少，因此在園林中常被用來作為地被植物。我國在園林的造景利用上，自80年代才開始將觀賞地被竹類從日本進行引種[4]，隨后在南京、蘇州、連云港、上海等地應(yīng)用于公園中。目前我國北方地區(qū)的植被綠化配置上層大型喬木及灌木林基本成型，下層地被類型的多元化配置亟待解決，因此對地被竹的引種到北方后的生長發(fā)育進行了大量研究[5]，對地被竹的耐寒性進行研究，為其在園林中作為林下植被提供了科學(xué)基礎(chǔ)，從而為引種工作提供依據(jù)[6-8]。翠竹(Pleioblastuspygmaea)屬于禾本科竹亞科苦竹屬植物，耐陰性矮小灌木狀地被竹類，竹根根系較淺，竹鞭單軸生長，混生型竹種。翠竹的葉片小且翠綠，水平排列成扇形，翠竹的地上部分高度約20～40 cm，稈低矮密生，修剪后翠竹生長周期延長，高度降低，高生長呈現(xiàn)“慢-快-慢”規(guī)律[9]，適宜人為控制植株高度進行修剪。翠竹性耐寒，適宜疏松的砂質(zhì)土，喜肥。翠竹的繁殖可通過播種途徑，亦可通過對竹鞭進行剪鞭埋鞭進行栽種，一般對翠竹的移栽是在春季或者梅雨季節(jié)帶宿土進行移植。通常將翠竹栽植于花壇中作為庭院布置植被，或者是作為公園的邊緣植物盆景欣賞，是城市園林綠化中優(yōu)良的地被植物種類。翠竹在園林的綠化中作為綠籬，具有隔斷效果，在應(yīng)用中具有不同的搭配形式，竹草搭配也是造園的一種綠化形式。地被竹不僅具有較強的園林綠化觀賞價值，因其鞭根發(fā)達且錯綜復(fù)雜地串鞭在一起，因此被用于固土護坡，對于水土保持具有較高的利用價值[10]。作為園林地被的一種新穎植物類群[11]，研究地被竹在園林綠化中的運用，對于將綠色生態(tài)引入城市，構(gòu)建多層帶的立體綠化，豐富城市植被類型具有重要作用。因此用地被竹類植物替代部分草坪植物，為垂直綠化尋找良好的材料配置，對于豐富生態(tài)多樣性、豐富地被植物品種具有重要作用[12]。

1　材料與方法

1.1　研究區(qū)概況

研究地設(shè)置在江蘇省南京市溧水區(qū)白馬鎮(zhèn)的白馬教學(xué)科研基地，白馬鎮(zhèn)位于溧水的東南部，屬于典型的寧鎮(zhèn)揚低丘陵區(qū)，地理坐標(biāo)介于31°23′～31°48′N，118°51′～119°14′E之間。屬于北亞熱帶的季風(fēng)氣候，年平均氣溫為15.5 ℃，年平均降水量為1 037 mm。四季分明，梅雨季節(jié)從6月中旬到7月上旬，平均約22 d，梅雨季節(jié)短期降水強度大，研究地土壤主要為下蜀系黃棕壤土。

1.2　研究材料

研究對象為翠竹的1年生實生苗，實生苗來源為白馬教學(xué)科研基地，地理坐標(biāo)為31°36′49″N，119°9′15″E。研究用的翠竹為南京林業(yè)大學(xué)2015年3月份開花的翠竹播種苗。

1.3　研究方法

研究通過人為控制光照強度進行研究，通過在花盆上方覆蓋不同處理的雙色板蓋板來控制光照，對照不覆蓋蓋板。通過在不同蓋板的正中央位置打不同的圓形孔孔徑來控制光照強度，雙色板圓形孔孔直徑大小分別為5、10、15、20、25 cm，分別標(biāo)記為A處理、B處理、C處理、D處理、E處理，對照標(biāo)為F處理，不同處理的孔徑面積占花盆的內(nèi)徑面積的照度百分比依次為1.29%、5.16%、11.62%、20.66%、32.28%、100%，每個處理重復(fù)36次。將2叢翠竹1年生實生苗移栽到塑料花盆正中位置，內(nèi)置營養(yǎng)土與黃土混合土厚度為15 cm；混合土為m(營養(yǎng)土)∶m(黃心土)=1∶2，塑料花盆的規(guī)格為外徑51 cm，內(nèi)徑44 cm，高度35 cm。蓋板材料為雙色板，大小為60 cm×60 cm，厚度為0.5 cm。

2016年3月上旬移栽成活，成活1周后蓋塑料蓋板，待生長穩(wěn)定后，從2016年4月15日起開始第1次取樣，之后隔25 d左右進行下次取樣，共取樣12次，記為12次調(diào)查時間。1-12次調(diào)查時間的具體日期分別為2016年4月15日、5月9日、6月3日、6月27日、7月21日、8月12日、9月5日、9月23日、10月17日、11月10日、12月2日和12月26日。每次取樣每個處理隨機抽取3盆，用水管沖洗土后取出樣本，分別裝自封袋，并用黑色記號筆標(biāo)記好序號，帶回實驗室處理。

將測定的樣品分別放入GZX-9146MBE電熱恒溫鼓風(fēng)干燥箱(上海博迅實業(yè)有限公司醫(yī)療設(shè)備廠出品)進行烘干。首先選擇在105 ℃條件下對翠竹殺青2 h，之后在75 ℃條件下烘48 h至恒重。烘干之后，分別區(qū)分開翠竹的葉片、地上部分和地下部分，用電子天平(奧豪斯-托利多儀器(上海)有限公司出品，測量范圍0～210 g，測量精度0.000 1 g)測定烘干后的翠竹的生物量干質(zhì)量：葉片的干質(zhì)量；地上部分干質(zhì)量；地下部分干質(zhì)量。

2　數(shù)據(jù)處理

研究的影響因素主要是6個不同的光照處理以及12次調(diào)查時間，因此，首先選用線性模型分別分析兩個主要因素與翠竹各生長指標(biāo)之間的關(guān)系，作出數(shù)據(jù)分布的箱線圖。為了避免雙因素之間存在著非線性關(guān)系的影響，因此還調(diào)用了廣義可加模型(GAM)，以光照處理和調(diào)查時間為自變量，以各部件生物量為因變量，研究雙因素對翠竹干質(zhì)量生物量的影響，并作出偏殘差圖[13-14]，對各影響指標(biāo)結(jié)果進行預(yù)測。

為了分析翠竹生物量的分配，研究翠竹葉片干質(zhì)量與地上部分干質(zhì)量和地下部分干質(zhì)量生物量之間的相關(guān)關(guān)系，加入不偏長軸法(RMA)指數(shù)模型[15]，對三者之間的相關(guān)系數(shù)進行檢驗。

試驗數(shù)據(jù)在Microsoft Excel 2007軟件中進行整理、簡單計算和部分圖表生成，在R 3.0.0 (R Development Core Team，2013)軟件中進行線性擬合回歸和部分圖像制作。

3　結(jié)果與分析

3.1　翠竹各構(gòu)件干質(zhì)量的研究

對不同光照強度下的翠竹的葉片、地上部分和地下部分的干質(zhì)量進行統(tǒng)計，翠竹各構(gòu)件的干質(zhì)量生物量分布的箱線圖見圖1。從圖1A、1B、1C中可以看出：隨著光照強度的增加，翠竹各部件干質(zhì)量隨著增加。由圖1A可知：A和B處理下的葉片干質(zhì)量變化不大，葉片干質(zhì)量生物量在0.232 1～0.240 6 g，A、B、C、D和 E處理下葉片干質(zhì)量均值分別是F對照下的11.75%、12.18%、21.6%、35.97%和49.5%。翠竹地上部分干質(zhì)量生物量增加，在F對照下地上部干質(zhì)量最大為2.998 3 g，A處理下地上部分干質(zhì)量最小為0.329 9 g，A、B、C、D和 E處理下地上部分干質(zhì)量均值分別是F對照下的11%、11.57%、21.83%、36.36%和48.08%。A和B光照度處理下，地下部分干質(zhì)量生物量均值相差不大，分別為0.333 6 g和0.351 5 g， A、B、C、D和 E處理下地下部分干質(zhì)量均值分別是F對照下均值的9.29%、9.78%、13.32%、24.75%和37.45%。隨著光照時間的增加，翠竹各構(gòu)件干質(zhì)量呈現(xiàn)先增加后降低的波動趨勢，不同光照強度下各構(gòu)件干質(zhì)量的均值在第9次調(diào)查時間時達到最大，干質(zhì)量積累最大。

從圖2A、2C、2E中可以看出光照強度對翠竹各構(gòu)件干質(zhì)量的影響呈現(xiàn)為正效果，即隨著光照強度的增加，各構(gòu)件的干質(zhì)量積累增加。從圖2B和2D可以看出，光照時間在研究前期對葉片及地上部分干質(zhì)量影響呈現(xiàn)正效果，在第10次調(diào)查時間各構(gòu)件的干質(zhì)量積累最大，之后呈現(xiàn)為局部負(fù)效果。模型預(yù)測結(jié)果與實際有所偏差，但總體趨勢與實際相符。即在10月中旬至11月中旬時期積累達到最大，而試驗過程中翠竹的葉片在11月就開始有枯黃現(xiàn)象產(chǎn)生，甚至大量葉片凋落，造成后期生物量積累有所下降。從圖2F可以看出研究過程中，光照時間對地下部分干質(zhì)量呈現(xiàn)為正效果。模型預(yù)測在第10次調(diào)查時間時葉片干質(zhì)量最大，因此地上部干質(zhì)量積累也最大。光照處理和調(diào)查時間對地下部分干質(zhì)量積累的預(yù)測結(jié)果與試驗結(jié)果有偏差，但整體上的變化趨勢與試驗結(jié)果是一致的。

3.2　翠竹葉片干質(zhì)量占翠竹總干質(zhì)量百分比的研究

圖1　不同處理和調(diào)查時間下翠竹各構(gòu)件的干質(zhì)量分布Fig.1　The dry weight of different parts of P. pygmaea under various treatments and investigation times說明：箱線圖中邊框中黑色粗線表示數(shù)據(jù)中的中位數(shù)，箱線圖中的上下邊緣分別代表數(shù)據(jù)的上四分位數(shù)和下四分位數(shù)。箱線圖中的上下邊緣延伸線即“觸須”，表示數(shù)據(jù)正常的分布范圍。如果數(shù)據(jù)無異常值，則箱體的上下邊緣分別代表數(shù)據(jù)的最大值和最小值；如出現(xiàn)異常值，則外端點處表示了箱體四分位數(shù)間距的1.5倍值。

圖2　不同光照處理和調(diào)查時間對翠竹各構(gòu)件干質(zhì)量影響的偏殘差圖Fig.2　Partial residuals of the dry weight of various parts of P. pygmaea under different light treatments and investigation time

圖3　不同光照處理和調(diào)查時間對葉片干質(zhì)量百分比影響的偏殘差圖Fig.3　Partial residuals of the proportion of leaf dry weight to total weight under different light treatments and investigation time

3.3　翠竹各構(gòu)件干質(zhì)量之間的相關(guān)關(guān)系

在研究翠竹各構(gòu)件生物量分配上，選用不偏長軸法(RMA)進行研究分析，調(diào)用了RMA指數(shù)模型，對地上部分干質(zhì)量和葉質(zhì)量進行線性擬合，擬合后的地上部分干質(zhì)量與葉質(zhì)量之間的RMA指數(shù)為0.996 2，擬合程度很高，表明在翠竹的幾種不同光照處理下，翠竹的葉片和地上部分干質(zhì)量之間線性相關(guān)。在對翠竹的地下部分干質(zhì)量生物量和葉片干質(zhì)量之間進行擬合時，翠竹地下部分干質(zhì)量與葉質(zhì)量之間的RAM指數(shù)為0.941 9，擬合優(yōu)度較高呈線性關(guān)系。翠竹的葉質(zhì)量與翠竹的地上部分和地下部分的干質(zhì)量影響線性相關(guān)。因此在對翠竹地上部分干質(zhì)量和地下部分干質(zhì)量進行RMA指數(shù)模型擬合時，表現(xiàn)出很強的線性相關(guān)，對應(yīng)的RMA指數(shù)為0.947 6。研究表明：翠竹的葉片干質(zhì)量生物量與地上部分干質(zhì)量生物量和地下生物量之間線性相關(guān)，翠竹的葉片生物量積累與地上部分生物量積累和地下部分生物量積累變化趨勢一致。因此在野外試驗中可以用翠竹的葉片干質(zhì)量生物量來研究翠竹的生物量積累大小。

4　討論

葉片的干質(zhì)量生物量可以反應(yīng)植物對環(huán)境中的光照強度的響應(yīng)，王德賢[16-17]通過對玉米進行了一系列研究，用葉質(zhì)量來速度估測單株葉面積。植物通過葉片的變形來適應(yīng)光照強度的變化，如通過降低單位葉面積干質(zhì)量[18]，或者通過降低葉片的厚度等，單位葉面積的生物量就相應(yīng)增加，有助于葉片來捕獲更多能量[19]，吳能表等[20]研究表明在低光照條件下，植物葉片的厚度隨著光照強度的降低而降低，通過增加單位葉片生物量來增加對光照的攔截，因此隨著光照強度的增加，翠竹的干質(zhì)量生物量增加。

綜合不同研究結(jié)果來看，光照強度的大小影響翠竹的干質(zhì)量生物量和生物量分配。由于光強影響著植物的光合作用，因此在弱光照處理下，植物葉片對光能的利用能力較弱[21]，進而影響植物積累光合產(chǎn)物，翠竹各部件的生物量隨光照處理強度的增加而增加，這與李濱勝等[22]的研究結(jié)果一致。對翠竹葉片干質(zhì)量生物量的研究，研究結(jié)果與孫曉玉等和吳能表等[19-20]的研究結(jié)果一致，翠竹在弱光照條件下，通過葉片對光照的適應(yīng)，來捕獲更多能量供給生命活動需求。隨著植物生命活動過程中對營養(yǎng)物質(zhì)的需求消耗植物的生物量積累，因此隨著光照時間的增加，葉片和地上部分干質(zhì)量生物量出現(xiàn)先增加后降低的趨勢。在對翠竹干質(zhì)量占總質(zhì)量百分比的研究中結(jié)果表明光照度影響著翠竹各部件的生物量分配，在較低的光照處理下，翠竹的生長量出現(xiàn)分化，翠竹的鞭根等生長量小，葉片干質(zhì)量生長量占據(jù)的比例就大。當(dāng)植物在遭受脅迫條件時，就會通過葉片的迅速凋落來降低死亡率，建立植物體系內(nèi)部的生理平衡[21]，因此，隨著調(diào)查時間的增加相當(dāng)于增加了光照的時間，葉片的干質(zhì)量百分比就下降。

綠色植物主要是通過葉片獲取光能，通過葉片來進行光合作用，積累產(chǎn)物，植物的葉片生物量直接影響著植物的總的生物量，因此翠竹的地上部分和地下部分及葉片干質(zhì)量生物量之間線性相關(guān)。林新春等[22]對苦竹(Pleioblastusamarus)的研究結(jié)果表明隨著竹子年齡的增加，在個體的生物量分配上出現(xiàn)先增加后降低的趨勢，翠竹的葉片干質(zhì)量在后期出現(xiàn)單峰型變化，與蘇文會等[23]對大木竹(Bambusawenchouensis)的研究結(jié)果一致，葉片生物量的比例也出現(xiàn)先增加后降低的趨勢。植物的生長過程中，為了對外界環(huán)境的適應(yīng)，不斷通過對個體的調(diào)整，優(yōu)化個體生物量的分配[24]。在生產(chǎn)上也可以用器官的生物量來代替整個植株的生物總量[25-26]，在生產(chǎn)中具有重要的實踐意義。通過對翠竹各部件生物量的研究，為后期建立模型定性從植物的葉片質(zhì)量來推斷出翠竹整體生物量提供理論依據(jù)。

研究竹類植物各構(gòu)件生物量，對于不同立地條件下植物生物量的估測具有重要意義[27]。同時在眾多的研究中，通過模型的構(gòu)建，研究了竹子各器官生物量的分配，對于竹類植物的集約化經(jīng)營提供了合理的理論依據(jù)[28-29]。探討光照強度對翠竹干質(zhì)量的影響，有利于為同一生態(tài)條件下進行植被配置提供參考以便對其在園林中的應(yīng)用進行合理安排設(shè)計，創(chuàng)造宜居的城市生態(tài)，發(fā)揮更佳的生態(tài)價值[30]。

[2] 李瑞昌, 曹幫華. 鹽堿脅迫對2種地被竹生長指標(biāo)的影響[J]. 北方園藝, 2011(13): 83-87.

[3] 王建武, 吳純?nèi)A, 趙亮, 等. 6種觀賞地被竹的抗逆性研究[J]. 山東林業(yè)科技, 2015, 45(5): 48-51. DOI: 10.3969/j.issn.1002-2724.2015.05.011.

[4] 陳天國. 觀賞地被竹類引種選育及其園林應(yīng)用[J]. 世界竹藤通訊, 2010, 8(1): 23-26. DOI: 10.3969/j.issn.1672-0431.2010.01.007.

[5] 張玲, 曹幫華, 賈波, 等. 山東地區(qū)地被竹引種試驗初步研究[J]. 世界竹藤通訊, 2009, 7(1): 4-7. DOI: 10.3969/j.issn.1672-0431.2009.01.002.

[6] 秦宇. 山東省引種觀賞竹的抗寒性研究[D]. 泰安：山東農(nóng)業(yè)大學(xué), 2014.

[7] 孫飛飛. 觀賞竹的耐寒性研究及園林應(yīng)用[D]. 杭州：浙江農(nóng)林大學(xué), 2015.

[8] 黃滔, 黃程前, 劉瑋, 等. 湖南長沙地區(qū)13種觀賞竹抗寒性研究[J]. 中國農(nóng)學(xué)通報, 2015, 31(25): 6-12.

[9] 劉國華, 王福升. 修剪對6種地被竹高生長的影響[J]. 南京林業(yè)大學(xué)學(xué)報(自然科學(xué)版), 2008, 32(2): 145-147. DOI: 10.3969/j.issn.1000-2006.2008.02.033.

[10] 陳雙林, 楊清平, 郭子武. 主要環(huán)境因素對小佛肚竹出筍、成竹和畸形稈率的影響[J]. 四川農(nóng)業(yè)大學(xué)學(xué)報, 2008, 26(1): 117-120. DOI: 10.3969/j.issn.1000-2650.2008.01.025.

[11] 史軍義, 易同培, 馬麗莎, 等. 園林地被竹及其開發(fā)利用[J]. 四川林業(yè)科技, 2006, 27(6): 95-100. DOI: 10.3969/j.issn.1003-5508.2006.06.020.

[12] 魏娜, 馮小虎. 5種地被竹在北京園林中的適應(yīng)性評價與應(yīng)用[J]. 世界竹藤通訊, 2015，13(2): 32-35. DOI: 10.13640/j.cnki.wbr.2015.02.009.

[13] HASTIE T J，TIBSHIRANI R J．Generalized additive models[J]．Statistical Science, 1990, 1(3): 297-328．

[14] WOOD S N. Generalized additive models: An introduction with R [M]. London: Chapman and Hall, 2006, 119-323.

[15] LEDUC D J. A comparative analysis of the reduced major axis technique of fitting [J]. Canadian Journal of Forest Research, 1987, 17(17): 654-659.

[16] 王得賢. 用單葉重速測青貯玉米單株總?cè)~面積[J]. 青海師范大學(xué)學(xué)報(自然科學(xué)版), 1996(2): 46-50. DOI: 10.16229/j.cnki.issn1001-7542.1996.02.010.

[17] 王得賢. 用單葉干重速測青貯玉米地上生物量[J]. 中國草食動物科學(xué), 1997(3): 16-19.

[18] 柳鳳娟, 向雙, 陽小成, 等. 兩種光照生境下4種常綠闊葉樹的單位葉面積干重、光合能力與化學(xué)防御物質(zhì)含量比較[J]. 應(yīng)用與環(huán)境生物學(xué)報, 2010, 16(4): 462-467. DOI: 10.3724/SP.J.1145.2010.00462.

[19] 孫曉玉, 陸兆華, 李鵬輝, 等. 紫莖澤蘭(Eupatoriumadenophorum)對光強的生態(tài)適應(yīng)性[J]. 林業(yè)研究(英文版), 2006, 17(2): 116-120.

[20] 吳能表, 談鋒, 肖文娟, 等. 光強因子對少花桂幼苗形態(tài)和生理指標(biāo)及精油含量的影響[J]. 生態(tài)學(xué)報, 2005, 25(5): 1159-1164. DOI:10.3321/j.issn:1000-0933.2005.05.031.

[21] 田媛, 唐立松, 喬瑞平. 梭梭幼苗個體生長規(guī)律與死亡率關(guān)系研究[J]. 植物研究, 2016, 36(1): 84-89. DOI: 10.7525/j.issn.1673-5102.2016.01.012.

[22] 林新春, 方偉, 俞建新, 等. 苦竹各器官生物量模型[J]. 浙江林學(xué)院學(xué)報, 2004, 21(2): 168-171. DOI:10.3969/j.issn.2095-0756.2004.02.011.

[23] 蘇文會, 顧小平, 官鳳英, 等. 大木竹種群生物量結(jié)構(gòu)及其回歸模型[J]. 南京林業(yè)大學(xué)學(xué)報(自然科學(xué)版), 2006, 30(5): 51-54. DOI:10.3969/j.issn.1000-2006.2006.05.012.

[24] FABBRO T, KORNER C. Altitudinal differences in flower traits and reproductive allocation [J]. Flora, 2004, 199(1): 70-81.

[25] 劉長成, 魏雅芬, 劉玉國, 等. 貴州普定喀斯特次生林喬灌層地上生物量[J]. 植物生態(tài)學(xué)報, 2009, 33(4): 698-705. DOI:10.3773/j.issn.1005-264x.2009.04.008.

[26] 楊前宇, 謝錦忠, 張瑋, 等. 椽竹各器官生物量模型[J]. 浙江農(nóng)林大學(xué)學(xué)報, 2011, 28(3): 519-526. DOI:10.3969/j.issn.2095-0756.2011.03.027.

[27] 鄭郁善, 梁鴻燊. 臺灣桂竹各器官生物量模型研究[J]. 竹子研究匯刊, 1998, 17(1): 37-41.

[28] 徐振國, 黃大勇, 郭起榮, 等. 麻竹地上器官生物量分配及其模型優(yōu)化[J]. 中南林業(yè)科技大學(xué)學(xué)報, 2016, 36(4): 80-83. DOI:10.14067/j.cnki.1673-923x.2016.04.015.

[29] 邱銀河. 撐麻7號竹地上部分生物量分配研究[J]. 世界竹藤通訊, 2007, 5(3): 29-31. DOI:10.3969/j.issn.1672-0431.2007.03.008.

[30] 曹洪霞. 淺談光照對植物生態(tài)作用的影響[J]. 安徽農(nóng)學(xué)通報, 2010, 16(3): 56-57. DOI:10.3969/j.issn.1007-7731.2010.03.025.