李蓓蕾，宋照亮,2,3，姜培坤,2，周國(guó)模,2，李自民,4

(1.浙江農(nóng)林大學(xué) 環(huán)境與資源學(xué)院，浙江 臨安 311300;2.浙江農(nóng)林大學(xué) 浙江省森林生態(tài)系統(tǒng)碳循環(huán)與固碳減排重點(diǎn)實(shí)驗(yàn)室,浙江 臨安311300;3.中國(guó)科學(xué)院 地球化學(xué)研究所 環(huán)境地球化學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,貴州 貴陽(yáng) 550002;4.中國(guó)農(nóng)業(yè)科學(xué)院 茶葉研究所/農(nóng)業(yè)部茶樹(shù)生物學(xué)與資源利用重點(diǎn)實(shí)驗(yàn)室，浙江 杭州310008)

植硅體，又稱 “植硅石”[1]，存在于很多植物細(xì)胞中，禾本科Poaceae植物細(xì)胞中其含量特豐富[2]。植硅體是高等植物在生長(zhǎng)過(guò)程中沉淀在細(xì)胞壁、細(xì)胞腔內(nèi)以及細(xì)胞間隙中的非晶質(zhì)二氧化硅礦物顆粒[3]。植硅體結(jié)構(gòu)復(fù)雜，其主要成分是二氧化硅(670～950 g°kg－1)，水(10～120 g°kg－1)，碳(1～60 g°kg－1)以及一些微量元素鈉、鉀、鈣、鐵、鋁、鈦等[4]。植物死亡或凋落以后，大量的植硅體被釋放到土壤中并很好的保存起來(lái)[5]，對(duì)全球碳匯起到重要作用[2]。竹子作為典型的硅超富集禾本科植物[6－7]，在全球范圍內(nèi)的覆蓋面積約為22×106hm2[8]。中國(guó)是世界上竹子分布最廣的國(guó)家之一，現(xiàn)存的竹林面積約有5×106hm2[9]，其中約有2/3是毛竹Phyllostachys edulis[10]。毛竹是中國(guó)經(jīng)營(yíng)歷史悠久的竹種[11]，具有巨大的生物量，它的生態(tài)學(xué)功能在陸地硅和碳循環(huán)中起到重要的作用[12]。在植物的生長(zhǎng)過(guò)程中，受蒸騰作用控制，硅元素主要以植硅體的形式富集在竹子等植物的葉片中[2,13－14]。植硅體主要通過(guò)竹葉凋落物歸還于土壤中[2]。本研究選取浙江省不同地區(qū)的毛竹作為研究對(duì)象，通過(guò)對(duì)不同巖性土壤上發(fā)育的同一竹齡毛竹竹葉和同一巖性土壤上發(fā)育的不同竹齡的毛竹竹葉中植硅體含量的分布特征進(jìn)行分析，來(lái)闡明毛竹中植硅體的分布規(guī)律，為了解毛竹林植硅體碳匯調(diào)控以及植硅體在竹林生態(tài)系統(tǒng)硅和碳在生物地球化學(xué)循環(huán)中的作用提供科學(xué)參考。

1 材料與方法

本研究所選取的毛竹竹葉樣品的采集點(diǎn)分別位于浙江省臨安市青山(30°13′N，119°46′E)，安吉縣山川鄉(xiāng)船壩(30°27′N，119°41′E)，紹興市新昌縣巧英(29°23′N，121°11′E)和紹興市新昌縣大市聚(29°28′N，120°59′E)。所處地區(qū)巖性分別為頁(yè)巖、花崗閃長(zhǎng)巖、花崗巖和玄武巖。不同地區(qū)毛竹林均已去除林下灌木、雜草。樣品采集時(shí)間為5月竹葉成熟期。頁(yè)巖地區(qū)選取1年生、3年生和5年生的毛竹，分為上(距離毛竹冠頂0～1 m)，中(距離毛竹冠頂2～3 m)和下(距離毛竹冠頂3～4 m)等3個(gè)部位，分別采集竹葉150 g左右，編號(hào)裝袋；按照同樣的方法，在花崗閃長(zhǎng)巖、花崗巖和玄武巖地區(qū)分別采集5年生毛竹3個(gè)部位的竹葉(其中1年生的竹子實(shí)際為1年2個(gè)月生的竹子)。同時(shí)，采集各點(diǎn)附近0～20 cm的土壤500 g左右，編號(hào)裝袋。將采集回來(lái)的竹葉樣品用自來(lái)水反復(fù)沖洗，再用去離子水沖洗干凈，放置鼓風(fēng)烘箱中105℃殺青后，75℃左右烘干至恒量。將每袋植物樣品充分混合，分成2份，一份用剪刀將樣品其剪至0.5 cm左右供植硅體提取分析，另一份用植物粉碎機(jī)粉碎，供總硅測(cè)試分析。將采集回來(lái)的土壤樣品自然風(fēng)干，挑根，研磨并過(guò)篩后供土壤pH、有機(jī)碳和有效硅的分析測(cè)試。

植物樣品中總硅的測(cè)定采用偏硼酸鋰熔融-分光光度法[15]，并利用植物灌木枝葉標(biāo)樣CBW07602進(jìn)行監(jiān)測(cè)。竹葉樣品中植硅體采用改進(jìn)的微波消解法提取[16]，同時(shí)運(yùn)用Walkley-Black方法檢驗(yàn)并除去植硅體外緣可能存在的有機(jī)碳[17]。將已提取的植硅體放入鼓風(fēng)烘箱中70℃烘干，然后對(duì)植硅體進(jìn)行稱量和記錄，計(jì)算樣品中植硅體質(zhì)量分?jǐn)?shù)。土壤pH值、有機(jī)碳和有效硅質(zhì)量分?jǐn)?shù)等測(cè)定參照文獻(xiàn)[15]方法。

本研究的結(jié)果數(shù)據(jù)為3次重復(fù)的平均值，運(yùn)用Duncan新復(fù)極差法進(jìn)行多重比較分析。數(shù)據(jù)的統(tǒng)計(jì)分析均在Microsoft，Excel和SPSS等軟件上進(jìn)行。

2 結(jié)果與分析

2.1 毛竹土壤的基本理化參數(shù)

表1給出了不同巖性上發(fā)育的同一竹齡毛竹土壤的基本理化參數(shù)。不同巖性下毛竹土壤pH值差異顯著，變化范圍為pH 4.33～5.15，呈現(xiàn)為酸性。4種不同巖性下的土壤中有機(jī)碳的質(zhì)量分?jǐn)?shù)差異顯著，花崗巖和花崗閃長(zhǎng)巖上發(fā)育的土壤中有機(jī)碳分別為32.9和29.3 g°kg－1，要遠(yuǎn)高于頁(yè)巖和玄武巖。玄武巖上發(fā)育的土壤中有效硅的含量為237.5 mg°kg－1，顯著高于其他3種巖性(表1)。

表1 不同巖性類型毛竹土壤的基本理化參數(shù)Table 1 Basic physicochemical parameters of soil from different lithologies

2.2 不同竹齡以及不同巖性下毛竹中二氧化硅和植硅體的平均質(zhì)量分?jǐn)?shù)

同一巖性(頁(yè)巖)下不同竹齡的毛竹竹葉中二氧化硅平均質(zhì)量分?jǐn)?shù)變化范圍是51.8～62.6 g°kg－1，植硅體平均質(zhì)量分?jǐn)?shù)為50.8～57.6 g°kg－1，差異不顯著(表2)。不同巖性下同一竹齡的毛竹竹葉中二氧化碳以及植硅體平均質(zhì)量分?jǐn)?shù)變化見(jiàn)表3。由表3可知:不同巖性下所發(fā)育的毛竹竹葉中二氧化硅質(zhì)量分?jǐn)?shù)以及植硅體質(zhì)量分?jǐn)?shù)的變化顯著，變化范圍分別為 56.1～103.7 g°kg－1，50.8～99.1 g°kg－1。

表2 頁(yè)巖土壤上不同竹齡毛竹竹葉中二氧化硅及植硅體平均質(zhì)量分?jǐn)?shù)Table 2 Effects of plant age based on the same lithology on the average contents of SiO2and phytolith in moso bamboo leaves

表3 不同巖性土壤上5年生毛竹竹葉中二氧化硅及植硅體平均質(zhì)量分?jǐn)?shù)Table3 Effectsoflithologywiththesameplantageontheaverage contents of SiO2and phytolith in moso bamboo leaves

2.3 不同竹齡和不同部位毛竹竹葉中植硅體質(zhì)量分?jǐn)?shù)

如圖1中所示:通過(guò)對(duì)青山地區(qū)不同竹齡同一部位竹葉中植硅體質(zhì)量分?jǐn)?shù)對(duì)比，不同竹齡的上部和中部葉子中植硅體沒(méi)有明顯差異，但是和下部葉子中植硅體質(zhì)量分?jǐn)?shù)差異明顯。青山湖地區(qū)同一竹齡竹葉中植硅體的質(zhì)量分?jǐn)?shù)基本上是由下部到上部遞增的，除了青山3年生毛竹中不同部位的葉子中植硅體沒(méi)有顯著性差異，青山1年生毛竹和青山5年生毛竹中具有顯著差異。

2.4 不同巖性和不同部位毛竹竹葉中植硅體質(zhì)量分?jǐn)?shù)

不同巖性對(duì)其上發(fā)育的毛竹中植硅體質(zhì)量分?jǐn)?shù)的影響是比較大的，如圖2所示。不同巖性上毛竹的同一部位竹葉的對(duì)比，植硅體質(zhì)量分?jǐn)?shù)的變化都是比較大的，并且都有顯著差異。花崗巖以及玄武巖上發(fā)育的毛竹，其不同部位竹葉中植硅體的質(zhì)量分?jǐn)?shù)差異不大；在頁(yè)巖和花崗閃長(zhǎng)巖上發(fā)育的毛竹中上部、中部和下部葉子中植硅體質(zhì)量分?jǐn)?shù)變化比較大，差異顯著。

3 討論

3.1 毛竹竹葉中植硅體質(zhì)量分?jǐn)?shù)及其影響因素

毛竹竹葉中植硅體質(zhì)量分?jǐn)?shù)變化范圍為50.8～99.1 g°kg－1，同一巖性土壤上發(fā)育的不同竹齡的毛竹竹葉中二氧化硅和植硅體的平均質(zhì)量分?jǐn)?shù)的不顯著變化可能是與竹葉平均2 a為周期的更新有關(guān)。此外，毛竹林的凋落量還受大、小年(大年是指出筍多的年份，小年是指出筍少的年份)的影響。一般來(lái)講，小年毛竹林的凋落量要大于大年，凋落量的變化也可能影響竹葉中二氧化硅和植硅體的質(zhì)量分?jǐn)?shù)[18]。不同巖性條件下發(fā)育的毛竹林竹葉中二氧化硅以及植硅體的質(zhì)量分?jǐn)?shù)差異都比較大，這與母巖及其形成土壤的理化性質(zhì)有關(guān)。不同母巖發(fā)育的土壤的pH值、有機(jī)碳和有效硅的質(zhì)量分?jǐn)?shù)差異比較大，從而影響其上發(fā)育的毛竹林對(duì)土壤溶液中硅的吸收以及竹葉中植硅體的產(chǎn)生。如相對(duì)于玄武巖和頁(yè)巖地區(qū)的土壤，以花崗巖和花崗閃長(zhǎng)巖為母巖的土壤中有效硅的質(zhì)量分?jǐn)?shù)雖然不高，但是在這2種巖性土壤上發(fā)育的毛竹竹葉中植硅體的質(zhì)量分?jǐn)?shù)較高，這可能與這2種土壤的pH值較低、土壤中有機(jī)碳的質(zhì)量分?jǐn)?shù)相對(duì)較高，土壤較疏松，能夠促進(jìn)植物對(duì)硅的吸收，有利于植物中植硅體的積累有關(guān)[13,19]。然而到底是哪種原因占主要地位，或是幾種原因同時(shí)作用，亦或是還有其他的可能原因的影響，都還需要更進(jìn)一步的調(diào)查研究。

圖1 頁(yè)巖土壤上不同竹齡和不同部位毛竹葉中植硅體的質(zhì)量分?jǐn)?shù)Figure 1 Contents of phytolith in different parts of the leaves within Moso bamboo among different ages

圖2 不同巖性和不同部位5年生毛竹葉中植硅體的質(zhì)量分?jǐn)?shù)Figure 2 Contents of phytolith in different parts of the leaves within Moso bamboo among different lithologies

植硅體質(zhì)量分?jǐn)?shù)與二氧化硅質(zhì)量分?jǐn)?shù)間極顯著的相關(guān)性(R2=0.852 1，P＜0.01，圖3)說(shuō)明:植硅體與硅的吸收與積累量密切相關(guān)。毛竹竹葉中植硅體質(zhì)量分?jǐn)?shù)基本呈現(xiàn)從上部到下部遞減的趨勢(shì)且差異性顯著(圖1，圖2)，主要與不同部位對(duì)硅吸收與積累不同有關(guān)。竹子對(duì)硅的吸收主要受蒸騰作用控制，上部的蒸騰作用要強(qiáng)于中部和下部，從而積累更多的硅[13,20]。

3.2 毛竹竹葉中植硅體產(chǎn)生通量及其影響因素

據(jù)報(bào)道，毛竹竹葉初級(jí)凈生產(chǎn)力約為3 049～6 114 kg°hm－2°a－1[21]。據(jù)此，我們估算毛竹中植硅體產(chǎn)生的平均通量為209.5～420.2 kg°hm－2°a－1。不同巖性上毛竹竹葉植硅體產(chǎn)生通量從大到小的順序依次為: 花崗巖(605.9 kg°hm－2°a－1)＞花崗閃長(zhǎng)巖(520.3 kg°hm－2°a－1)＞玄武巖(411.5 kg°hm－2°a－1)＞頁(yè)巖(310.6 kg°hm－2°a－1)。相對(duì)于玄武巖和頁(yè)巖地區(qū)的毛竹林，花崗巖和花崗閃長(zhǎng)巖上發(fā)育的毛竹中植硅體的產(chǎn)生通量相對(duì)較高，這主要與花崗巖和花崗閃長(zhǎng)巖上發(fā)育的毛竹葉片中植硅體含量較高有關(guān)。

除了不同母質(zhì)發(fā)育的土壤可通過(guò)改變植硅體質(zhì)量分?jǐn)?shù)對(duì)植硅體產(chǎn)生通量具有影響外，竹葉初級(jí)凈生產(chǎn)力的變化也是不可忽略的因素。毛竹林的葉凋落全年都會(huì)發(fā)生，但其數(shù)量隨著季節(jié)和竹子自身的生長(zhǎng)特點(diǎn)而變化[21]。此外，經(jīng)營(yíng)模式的不同也會(huì)對(duì)竹葉初級(jí)凈生產(chǎn)力產(chǎn)生一定的影響。一般來(lái)講，粗放經(jīng)營(yíng)要大于集約經(jīng)營(yíng)[21]。

圖3 植硅體質(zhì)量分?jǐn)?shù)與二氧化硅質(zhì)量分?jǐn)?shù)的相關(guān)性Figure 3 Correlation of phytolith content and SiO2content

3.3 對(duì)中國(guó)毛竹林植硅體碳匯調(diào)控的啟示

竹子(主要是毛竹)主要分布在熱帶和亞熱帶地區(qū)，約占世界森林總面積的1.5%～2.0%。中國(guó)現(xiàn)存的竹林面積約為5.0×106hm2[9]，毛竹面積約為3.3×106hm2[10]，集中分布在浙江、江西和湖南等中亞熱帶紅壤地區(qū)[22－23]。若按植硅體產(chǎn)生通量為 209.5～420.2 kg°hm－2°a－1計(jì)算，那么全國(guó)毛竹林的植硅體產(chǎn)生量約為 0.7～1.4 Tg°a－1(1 Tg=1012g)。在植硅體形成的過(guò)程中，植硅體內(nèi)可以包裹 1%～6%的有機(jī)碳[2,24－25]，稱為植硅體碳，其在全球碳循環(huán)和氣候變化中發(fā)揮著重要的作用[26－27]。若以植硅體中植硅體碳含量(3±1)%[14]和植硅體平均產(chǎn)生通量計(jì)算，那么中國(guó)毛竹每年通過(guò)植硅體大約可以固定(76.1～152.5)×106kg二氧化碳。竹子具有生長(zhǎng)快、繁殖快、易再生等優(yōu)點(diǎn)[28]。如果未來(lái)二三十年竹子面積繼續(xù)按目前每年2%～3%的速率增長(zhǎng)[8]，那么，到2050年竹子的面積將增加1倍以上[29]。若以毛竹植硅體產(chǎn)生平均通量計(jì)算，全國(guó)毛竹林植硅體及其包裹有機(jī)碳的含量都將增加1倍以上。

相對(duì)于玄武巖和頁(yè)巖等其他巖性上發(fā)育的毛竹林，花崗巖和花崗閃長(zhǎng)巖上毛竹植硅體產(chǎn)生通量較高(約為前者的1.5倍)，植硅體在土壤中的穩(wěn)定性相對(duì)較高(因土壤具有較低pH值)[13]，碳匯潛力較大。今后開(kāi)展竹林造林—再造林時(shí)，可以盡量選擇花崗巖和花崗閃長(zhǎng)巖等地區(qū)。此外，施用硅肥[30－34]或高供硅能力的巖粉，來(lái)提高毛竹竹葉中植硅體的產(chǎn)生通量都可以增加毛竹林植硅體產(chǎn)生量，從而提高植硅體內(nèi)所包裹的有機(jī)碳含量，固定更多大氣中的二氧化碳，對(duì)增加全球碳匯起到重要作用。

4 結(jié)論

不同巖性上發(fā)育的毛竹竹葉中植硅體質(zhì)量分?jǐn)?shù)為花崗巖(99.1 g°kg－1)＞花崗閃長(zhǎng)巖(85.1 g°kg－1)＞玄武巖(67.3 g°kg－1)＞頁(yè)巖(56.1 g°kg－1)。不同的巖性對(duì)其上發(fā)育的毛竹中植硅體的產(chǎn)生通量影響比較大，花崗巖(606.0 kg°hm－2°a－1)＞花崗閃長(zhǎng)巖(520.4 kg°hm－2°a－1)＞玄武巖(411.6 kg°hm－2°a－1) ＞頁(yè)巖(310.6 kg°hm－2°a－1)。

按毛竹林面積的3.3×106hm2和植硅體產(chǎn)生通量為209.5～420.2 kg°hm－2°a－1計(jì)算，全國(guó)毛竹林通過(guò)竹葉凋落物形式產(chǎn)生的植硅體量為0.7～1.4 Tg°a－1。若以植硅體中含碳量(3±1)%計(jì)算，那么全國(guó)毛竹大約可以固定二氧化碳(76.1～152.5)×106kg°a－1。

選擇花崗巖和花崗閃長(zhǎng)巖等毛竹植硅體產(chǎn)生通量高的巖性上開(kāi)展竹林造林—再造林，或施用硅肥或高供硅能力的巖粉，或增加竹林的面積可顯著提高毛竹竹葉中植硅體及其包裹碳的產(chǎn)生通量，增強(qiáng)植硅體碳匯能力。

[1]呂厚遠(yuǎn)，賈繼偉，王偉銘，等.“植硅體”含義和禾本科植硅體的分類[J].微體古生物學(xué)報(bào)，2002，19(4):389－396.Lü Houyuan，JIA Jiwei，WANG Weiming，et al.On the meaning of phytolith and its classification in gramineac[J].Acta Micropalaeontol Sin，2002，19(4): 389－396.

[2]PARR J，SULLIVAN L，CHEN Bihua，et al.Carbon bio-sequestration within the phytoliths of economic bamboo species [J].Global Change Biol，2010，16(10): 2661－2667.

[3]PIPERNO D R.Phytolith Analysis: An Archaeological and Geological Perspective [M].San Diego: Academic Press，1988.

[4]呂后東，李榮華，徐慶梅.植物硅酸體分析與中草藥鑒定[J].時(shí)珍國(guó)藥研究，1998，8(5):28－30.Lü Houdong，LI Ronghua，XU Qingmei.Analysis of opal phytoliths and identification of Chinese herbal medicine[J].Shizhen J Tradit Chin Med Res，1998，8(5): 28－30.

[5]ALEXANDRE A，MEUNIER J D，COLIN F，et al.Plant impact on the biogeochemical cycle of silicon and related weathering processes [J].Geochim Cosmochim Acta，1997，61(3): 677－682.

[6]LI Zhenji，LIN Peng，HE Jianyuan，et al.Silicon’s organic pool and biological cycle in moso bamboo community of Wuyishan Biosphere Reserve [J].J Zhejiang Univ Sci B，2006，7(11): 849－857.

[7]黃張婷，姜培坤，宋照亮，等.不同竹齡雷竹中硅及其他營(yíng)養(yǎng)元素吸收和積累特征[J].應(yīng)用生態(tài)學(xué)報(bào)，2013，24(5): 1347－1353.HUANG Zhangting，JIANG Peikun，SONG Zhaoliang，et al.Uptake and accumulation characteristics of silicon and other nutritional elements in different age Phyllostachys praecox plants [J].Chin J Appl Ecol，2013,24(5):1347－1353.

[8]郭起榮，楊光耀，杜天真，等.中國(guó)竹林的碳素特征[J].世界竹藤通訊，2005，3(3):25－28.GUO Qirong，YANG Guangyao，DU Tianzhen，et al.Carbon character of Chinese bamboo forest[J].World Bamboo Rattan，2005，3(3): 25－28.

[9]GUI Renyi，LENG Huanan，ZHUANG Shunyao，et al.Aluminum tolerance in moso bamboo (Phyllostachys pubescens) [J].Bot Rev，2011，77(3): 214－222.

[10]JIANG Peikun，MENG Cifu，ZHOU Guomo，et al.Comparative study of carbon storage in different forest stands in subtropical China [J].Bot Rev，2011，77(3): 242－251.

[11]應(yīng)葉青，郭璟，魏建芬，等.水分脅迫下毛竹幼苗光合及葉綠素?zé)晒馓匦缘捻憫?yīng)[J].北京林業(yè)大學(xué)學(xué)報(bào)，2009，31(6): 128－133.YING Yeqing，GUO Jing，WEI Jianfen，et al.Photosynthetic and chlorophyll fluorescent responses of Phyllostachys pubescens seeding to water deficiency stress [J].J Beijing For Univ，2009，31(6): 128－133.

[12]DU Huaqiang，ZHOU Guomo，F(xiàn)AN Wenyi，et al.Spatial heterogeneity and carbon contribution of aboveground biomass of moso bamboo by using geostatistical theory [J].Plant Ecol，2010，207(1): 131－139.

[13]SONG Zhaoliang，WANG Hailong，STRONG P J，et al.Plant impact on the coupled terrestrial biogeochemical cycles of silicon and carbon: implications for biogeochemical carbon sequestration [J].Earth-Sci Rev，2012，115(4): 319－331.

[14]SONG Zhaoliang，LIU Hongyan，LI Beilei，et al.The production of phytolith-occluded carbon in China’s forests:implications to biogeochemical carbon sequestration [J].Global Change Biol，2013，19(9): 2907－2915.

[15]魯如坤.土壤農(nóng)業(yè)化學(xué)分析方法[M].北京:中國(guó)農(nóng)業(yè)科技出版社，2000:12－202.

[16]PARR J F，DOLIC V，LANCASTER G，et al.A microwave digestion method for the extraction of phytoliths from herbarium specimens [J].Rev Palaeobot Palynol，2001，116(3/4): 203－212.

[17]WALKLEY A，BLACK I A.An examination of the degtjareff method for determining soil organic matter，and a proposed modification of the chromic acid titration method [J].Soil Sci，1934，37(1): 29－38.

[18]陳靈芝，黃建輝，嚴(yán)昌榮.中國(guó)森林生態(tài)系統(tǒng)養(yǎng)分循環(huán)[M].北京:氣象出版社，1997:3－10..

[19]李建武.海南島玄武巖發(fā)育土壤的痕量元素和Sr-Nd同位素特征與大氣物源貢獻(xiàn)[D].南京:中國(guó)科學(xué)院南京土壤研究所，2012.LI Jianwu.Trace Elemental and Sr-Nd Isotopic Characteristics in Relation to Dust Accretion of Basalt-derived Soils in Hainan Island [D].Nanjing: Chinese Academy of Sciences.Institure of Soil Science，2012.

[20]DING Tiping，ZHOU Jianxiong，WAN Defang，et al.Silicon isotope fractionation in bamboo and its significance to the biogeochemical cycle of silicon [J].Geochim Cosmochim Acta，2008，72(5): 1381－1395.

[21]周國(guó)模，吳家森，姜培坤.不同管理模式對(duì)毛竹林碳儲(chǔ)量的影響[J].北京林業(yè)大學(xué)學(xué)報(bào)，2006，28(6):51－55.ZHOU Guomo，WU Jiasen，JIANG Peikun.The impacts of different management modes on the carbon storage within moso bamboo [J].J Beijing For Univ，2006，28(6): 51－55.

[22]周國(guó)模，姜培坤，徐秋芳.竹林生態(tài)系統(tǒng)中碳的固定與轉(zhuǎn)化[M].北京:科學(xué)出版社，2010:1－5.

[23]汪慶華，唐根年，李睿.浙江省特色農(nóng)產(chǎn)品立地地質(zhì)背景研究[M].北京:地質(zhì)出版社，2007:10－180.

[24]PARR J F，SULLIVAN L A，QUIRK R.Sugarcane phytoliths: encapsulation and sequestration of a long-lived carbon fraction [J].Sugar Technol，2009，11(1): 17－21.

[25]PARR J F，SULLIVAN L A.Phytolith occluded carbon and silica variability in wheat cultivars [J].Plant Soil，2011，342(1－2): 165－171.

[26]孟賜福，姜培坤，徐秋芳，等.植物生態(tài)系統(tǒng)中的植硅體閉蓄有機(jī)碳及其在全球土壤碳匯中的重要作用[J].浙江農(nóng)林大學(xué)學(xué)報(bào)，2013，30(6):921－929.MENG Cifu，JIANG Pengkun，XU Qiufang，et al.PhytOC in plant ecologica system and its important roles in the global soil carbon sink [J].J Zhejiang A ＆ F Univ，2013，30(6): 921－929.

[27]李自民，宋照亮，李蓓蕾，等.杭州西溪濕地植物植硅體產(chǎn)生及其影響因素[J].浙江農(nóng)林大學(xué)學(xué)報(bào)，2013，30(4): 470－476.LI Zimin，SONG Zhaoliang，LI Beilei，et al.Phytolith production in wetland plants of the Hangzhou Xixi Wetlands ecosystem [J].J Zhejiang A ＆ F Univ，2013，30(4): 470－476.

[28]XU Yong，WONG Minghung，YANG Jianli，et al.Dynamics of carbon accumulation during the fast growth period of bamboo plant[J].Bot Rev，2011，77(3): 287－295.

[29]CHEN Xiangang，ZHANG Xiaoquan，ZHANG Yiping，et al.Changes of carbon stocks in bamboo stands in China during 100 years [J].For Ecol Manage，2009，258(7): 1489－1496.

[30]MATICHENKOV V，CALVERT D，SNYDER G.Silicon fertilizers for Citrus in Florida [J].Proc Fla State Hort Soc，1999，112: 5－8.

[31]ALVAREZ J，DATNOFF L E.The economic potential of silicon for integrated management and sustainable rice production [J].Crop Prot，2001，20(1): 43－48.

[32]MA Jianfeng，TAKAHASHI E.Soil，F(xiàn)ertilizer，and Plant Silicon Research in Japan [M].Amsterdam: Elsevier 2002:9-26.

[33]LIANG Yongchao，HUA Haixia，ZHU Yongguan，et al.Importance of plant species and external silicon concentration to active silicon uptake and transport[J].New Phytol，2006，172(1): 63－72.

[34]MECFEL J，HINKE S，GOEDEL W A，et al.Effect of silicon fertilizers on silicon accumulation in wheat [J].J Plant Nutr Soil Sci，2007，170(6): 769－772.