油松葉功能性狀分布特征及其控制因素

張 凱,侯繼華,*,何念鵬

1 北京林業(yè)大學(xué), 北京 100083 2 中國科學(xué)院地理科學(xué)與資源研究所, 北京 100101

油松葉功能性狀分布特征及其控制因素

張 凱1,侯繼華1,*,何念鵬2

1 北京林業(yè)大學(xué), 北京 100083 2 中國科學(xué)院地理科學(xué)與資源研究所, 北京 100101

植物葉功能性狀能夠直接或間接反應(yīng)植物對環(huán)境的適應(yīng)策略,這種適應(yīng)策略會間接影響植物的生長、繁殖和生存。目前已有大量關(guān)于植物葉功能性狀與環(huán)境間關(guān)系的研究,但這些研究多使用性狀平均值代替物種,忽略了性狀的種內(nèi)變化。油松是我國重要的造林樹種,對我國陸地生態(tài)系統(tǒng)起著重要的作用,研究其葉功能性狀與環(huán)境的關(guān)系有助于更好的的理解種內(nèi)功能性狀與環(huán)境間的關(guān)系。于2014年6—9月對遼寧、內(nèi)蒙、北京、山西、陜西、寧夏、青海等地天然油松進(jìn)行采樣并對其重要的8個葉功能性狀的分布特征及主控環(huán)境因子進(jìn)行了研究,結(jié)果表明：1)油松各葉性狀值存在較大的種內(nèi)變異系數(shù)(4.82%—25.85%),除1年生葉碳含量(LCC)油松各葉功能性狀值在不同研究地點(diǎn)間差異顯著(P<0.05)；2)油松葉長(LL)、葉厚(LT)、比葉面積(SLA)、氣孔密度(SD)、葉氮含量(LNC)存在較弱的經(jīng)度格局,LT、SD、LNC存在較弱的緯度格局(0.05

油松；葉功能性狀；分布特征；主控因素；通徑分析

植物功能性狀是植物個體可以測量的、能夠直接或間接反應(yīng)植物適應(yīng)策略的性狀[1],這種適應(yīng)策略會間接影響植物的生長、繁殖和生存[2]。國外研究表明,植物生長策略和植物對資源的利用能力與葉功能性狀具有緊密的聯(lián)系,同時植物對環(huán)境變化形成的生存對策也可以通過葉功能性狀反映出來[3]。隨著植物葉經(jīng)濟(jì)譜理論的提出[4],葉功能性狀及其與環(huán)境間關(guān)系的研究越發(fā)成為生態(tài)學(xué)研究的熱點(diǎn)問題。研究表明,常綠植物的比葉重隨著年平均氣溫和年降水量的增加而降低[5],植物葉片中的氮、磷含量極易受環(huán)境條件(如土壤的養(yǎng)分含量,生長季的溫度及年降水量)的影響,全球452個樣點(diǎn)的1280種植物葉氮和葉磷含量都隨著緯度升高和溫度降低而增加,氮磷比則降低[4]。我國對于葉功能性狀與環(huán)境間關(guān)系的研究起步較晚,多集中在小尺度海拔梯度上葉性狀與環(huán)境之間關(guān)系的研究[6- 7]或大尺度環(huán)境梯度下同一科屬植物葉功能性狀與環(huán)境因子之間的關(guān)系[8- 9]。然而在環(huán)境梯度上植物功能性狀的變化不僅受到環(huán)境的影響,還受系統(tǒng)發(fā)育的影響[10- 11],種內(nèi)植物性狀會沿著環(huán)境梯度而發(fā)生變化,這種變化能夠使植物種在新的環(huán)境下存活、生長和繁殖[12],用單一物種研究環(huán)境梯度下性狀間的關(guān)系可以更好的區(qū)分環(huán)境對性狀的影響[11]。戚德輝等[13]研究表明不同環(huán)境因素在不同程度上對鐵桿蒿(Artemisiasacrorum)功能性狀產(chǎn)生影響,對鐵桿蒿功能性狀影響大小依次是年均降雨量、年均溫和年均蒸發(fā)量。楊利民等[14]研究表明水分條件是影響中國東北樣帶羊草(Leymuschinensis)氣孔密度變化的主要因子。蘇文華等[15]人對滇石櫟(Lithocarpusdealbatus)的研究發(fā)現(xiàn)其功能性狀存在較大的種內(nèi)變異并且與生境緯度和溫度有關(guān)。張慧文等[7]的研究發(fā)現(xiàn)土壤含水量和土壤氮含量是天山云杉葉片功能特征沿海拔梯度變化的主要驅(qū)動因子。這些在種內(nèi)水平的研究可以弱化系統(tǒng)發(fā)育的影響,凸顯環(huán)境對功能性狀變化的影響。綜上所述,在研究植物功能性狀與環(huán)境間關(guān)系時(尤其在大尺度上),往往使用性狀平均值代替物種,忽略性狀的種內(nèi)變化。然而對于分布區(qū)較廣的物種而言,由于受到不同環(huán)境條件以及系統(tǒng)發(fā)育的影響,其種內(nèi)功能性狀在環(huán)境梯度上的變異不容忽視。本研究有助于更好的區(qū)分環(huán)境對油松葉性狀的影響進(jìn)一步了解植物種內(nèi)葉功能性狀與環(huán)境之間的關(guān)系,同時也是對植物功能性狀與環(huán)境間關(guān)系研究的重要補(bǔ)充。

油松(Pinustabuliformis)是我國特有種,具有耐低溫、干旱和瘠薄的生態(tài)學(xué)特性,是優(yōu)良的造林樹種,在中國遼寧、吉林、內(nèi)蒙古、河北、河南、山西、陜西、山東、甘肅、寧夏、青海以及四川北部等地皆有分布[16]。天然油松林對中國北方森林群落具有很強(qiáng)的代表性,同時對陸地生態(tài)系統(tǒng)碳循環(huán)和碳儲存也起到了重要的作用[17]。國內(nèi)對于油松的研究較早,但早期研究多集中在其遺傳特征[18-21]與形態(tài)特征[22]之上,隨著近年來植物功能性狀研究的開展,人們逐漸意識到功能性狀對于植物生長發(fā)育繁殖以及其與環(huán)境之間相互關(guān)系的重要性,對油松功能性狀的研究理應(yīng)受到重視。油松天然分布區(qū)緯度變化范圍在31°—44°N之間；經(jīng)度變化范圍在101°30′—124°45′E；海拔變化范圍在100—2800 m之間,分布區(qū)內(nèi)氣候條件,土壤狀況,地形因子等都存在較大差異[23]。Diaz等[24]在區(qū)域尺度上的研究發(fā)現(xiàn)植物功能性狀與氣候的關(guān)系比全球尺度更加緊密,在不同區(qū)域間和不同性狀間也存在差異。由此我們推斷天然油松各葉功能性狀在其分布區(qū)內(nèi)亦存在較大的差異,并且存在一定的分布格局,而且各性狀分布的主控因素不同。本文通過在油松天然分布區(qū)內(nèi)9個代表性樣點(diǎn)進(jìn)行取樣研究,試圖1)探明油松這一物種的各個葉功能性狀分布特征2)揭示油松各性狀分布的主控因素,進(jìn)一步理解種內(nèi)功能性狀與環(huán)境間的相互關(guān)系進(jìn)而對我們在全球變化的背景下更好的保護(hù)發(fā)展和利用油松這一物種提供理論指導(dǎo)。

圖1 油松分布圖及取樣點(diǎn)Fig.1 Distribution of Piuns tabuliformis and the sampling points

1 材料和方法

1.1 研究區(qū)概況

本研究在天然油松林分布區(qū)內(nèi)選擇9個代表性的樣點(diǎn),包括北京市延慶縣松山自然保護(hù)區(qū)、遼寧省鞍山市千山自然保護(hù)區(qū)、遼寧省綏中縣三山林場、內(nèi)蒙古黑里河自然保護(hù)區(qū)、山西省太岳山靈空山自然保護(hù)區(qū)、陜西省寧陜縣火地糖、陜西省黃龍縣蔡家川林場、青海省海南藏族自治區(qū)貴德縣、寧夏省銀川市賀蘭山自然保護(hù)區(qū)(圖1)。各研究區(qū)位于101°—124°E, 33°—42°N之間,海拔變化范圍在138—2804 m之間,年均溫4.8—10℃,年降水量254.2—1100 mm之間,土壤以棕壤和褐土為主。各研究地區(qū)環(huán)境條件差異較大,詳細(xì)情況見表1。

表1 研究地區(qū)自然狀況

1.2 樣地設(shè)置和樣品采集

本研究于2014年6月至9月按照從南到北的順序依次在各個研究地點(diǎn)進(jìn)行樣地設(shè)置和樣品采集。為減少林分條件及地形因子的影響在所有研究地點(diǎn)選取林齡相近(林齡通過咨詢當(dāng)?shù)亓謽I(yè)工作人員及油松胸徑綜合判斷)林分情況較為相似的天然油松成熟林,并在每一個研究地點(diǎn)選擇坡度較為平緩的陽坡或半陽坡設(shè)置3個20m×20m樣地,用GPS測定每一個樣地的經(jīng)緯度與海拔并記錄,用EM50儀器測定每一樣地內(nèi)土壤體積含水率(土壤中水分體積與土壤總體積的比值)與溫度。然后在每一個樣地內(nèi)挑選3株長勢良好無病蟲害的油松采集葉片,每株油松按其東南西北四個方向各取一枝生長良好的枝條,在每個枝條上選取無病蟲害的1年生與多年生松針各50束,將采集的松針置于兩片濕潤的濾紙之間,分開裝進(jìn)自封袋,放入保溫盒中儲存。最后在每個樣地對角線等距離用體積為100cm3的環(huán)刀取3份土樣裝入鋁盒中并做好標(biāo)記。

1.3 測定指標(biāo)及方法

1.3.1 植物樣品的測定

在每一個自封袋內(nèi)隨機(jī)選取20束松針,取出后迅速用濾紙吸干表面的水分。用1/10000的電子天平稱量其飽和鮮重并用精度為0.01mm的數(shù)顯游標(biāo)卡尺測定每一束松針的葉長(LL)、中間部位葉寬(LD)以及葉厚(LT),油松松針單面弧面面積采用公式S=π×LD×LL/2計算[25],油松單束松針面積=2S+2LD×LL。另隨機(jī)選取5—10束松針減掉其頭尾兩部分保留3—5cm置于福爾馬林溶液中保存?；氐綄?shí)驗(yàn)室后將松針放入85℃烘箱內(nèi)烘干48h至恒重并測定其干重,取出浸泡在福爾馬林溶液中的松針風(fēng)干24h,葉片氣孔個數(shù)通過SEM(Scanning Electron Microscope S- 3400N)拍攝照片計數(shù)。再據(jù)此計算其他指標(biāo)。計算公式如下：葉干物質(zhì)含量(LDMC)=葉片干重(mg)/葉片鮮重(mg)；比葉面積(SLA)=葉面積(cm2)/葉干重(g)；氣孔密度(SD)=氣孔個數(shù)(stomatal)/照片面積(mm2)。將烘干葉片磨粉過細(xì)篩后,采用重鉻酸鉀外加熱法測定有機(jī)碳含量(LCC)；經(jīng)H2SO4-H2O2法消煮后,用凱氏定氮法測定全氮含量(LNC),鉬銻抗比色法測定全磷含量(LPC)[26]。

1.3.2 土壤樣品的測定

將鋁盒中的土樣取出放于塑料薄膜上自然風(fēng)干并測定其干重。土壤容重(soil unit weight, SUW)=土壤干重(g)/100(cm3)。然后去除草根、石塊后研磨過0.15mm篩經(jīng)四分法取樣用于土壤C、N、P含量的測量。采用重鉻酸鉀外加熱法測定有機(jī)碳含量；經(jīng)H2SO4-H2O2法消煮后,用凱氏定氮法測定全氮含量；經(jīng)H2SO4-HCLO4法消煮后,用鉬銻抗比色法測定全磷含量[26]。土壤體積含水率由EM50直接讀出。

1.4 環(huán)境因子數(shù)據(jù)集建立

不同經(jīng)緯度環(huán)境條件對植物功能性狀的影響主要來源于氣候、土壤以及地形因子的差異[24]。通過數(shù)據(jù)的調(diào)查與整理,我們建立了油松葉功能性狀主要影響因素的數(shù)據(jù)集,其中氣候數(shù)據(jù)通過中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http://www.escience.gov.cn/)以及相關(guān)文獻(xiàn)[27-34]查詢得出,土壤數(shù)據(jù)與海拔由本實(shí)驗(yàn)測得：1)氣候因素：年均溫(MAT)、年降水量(MAP)、日照百分率(SP)、空氣相對濕度(ARH)；2)土壤因素：土壤體積含水率(SMCV)、土壤容重(SUW)、土壤氮含量(SNC)、土壤磷含量(SPC)、土壤碳含量(SCC)；3)地形因素：海拔(Alt)

1.5 數(shù)據(jù)分析

首先對油松各葉性狀值進(jìn)行描述性分析,得到各葉性狀值的分布特征及其變異程度,并采用線性回歸分析探討各葉性狀值與經(jīng)緯度之間的相關(guān)關(guān)系及其在經(jīng)緯度上的分布格局；然后通過單因素方差分析檢驗(yàn)油松各葉功能性狀值在不同地區(qū)是否存在顯著差異；最后通過通徑分析得出各環(huán)境因子對油松功能性狀的直接影響及其通過其他環(huán)境因子對油松功能性狀的間接影響大小,找出直接影響油松各功能性狀的主要環(huán)境因子。采用SPSS 17.0(SPSS, Chicago, USA)軟件進(jìn)行統(tǒng)計分析,采用Excel 2010(Microsoft Office 2013,New York,USA)和Minitab 16(State College, PA, USA)軟件作圖。

2 結(jié)果與分析

2.1 油松葉性狀分布特征

2.1.1 油松葉性狀統(tǒng)計特征

油松各葉性狀值基本滿足正態(tài)分布趨勢且都存在較大的標(biāo)準(zhǔn)差及變異系數(shù),1年生油松LL、LT、SLA、LDMC、SD、LNC、LPC、LCC平均值分別為9.65cm、1.38mm、127.59cm2/g、369.8mg/g、68.68個/mm2、1.21%、0.13%、49.77%。其中油松LPC含量變異系數(shù)最大為23.29%,LDMC變異系數(shù)最小為5.91%。多年生油松LL、LT、SLA、LDMC、SD、LNC、LPC、LCC平均值分別為11.16cm、1.41mm、102.85cm2/g、460mg/g、64.2個/mm2、1.17%、0.077%、51.35%。其中油松LPC含量變異系數(shù)最大為25.85%；LDMC變異系數(shù)最小為4.82%。其余各指標(biāo)詳細(xì)分布特征見圖2,圖3。

圖2 油松一年葉各性狀統(tǒng)計分析Fig.2 Statistical characteristics of annual leaf functional traitsLL： 葉長 Leaf length; LT： 葉厚 Leaf thickness; SLA： 比葉面積 Specific leaf area; LDMC： 葉干物質(zhì)含量 Leaf dry matter content; SD： 氣孔密度 Stomatal density; LNC： 葉氮含量 Leaf nitrogen concentration; LPC： 葉磷含量 Leaf phosphorus concentration; LCC： 葉碳含量 Leaf carbon concentration

圖3 油松多年葉各性狀頻度分布圖Fig.3 Statistical characteristics of perennial leaf functional traits

2.1.2 油松葉性狀分布空間格局

圖4 油松葉性狀經(jīng)度分布格局Fig.4 The longitude distribution pattern of the leaf functional traits of pinus tabuliformisLL: 葉長 Leaf length; LT: 葉厚 Leaf thickness; SLA: 比葉面積 Specific leaf area; SD: 氣孔密度 Stomatal density; LNC: 葉氮含量 Leaf nitrogen concentration;只將與經(jīng)度存在顯著回歸關(guān)系的性狀作圖列出

圖5 油松葉性狀緯度分布格局Fig.5 The latitude distribution pattern of the leaf functional traits of pinus tabuliformisLT: 葉厚 Leaf thickness; SD: 氣孔密度 Stomatal density; LNC: 葉氮含量 Leaf nitrogen concentration;只將與緯度存在顯著回歸關(guān)系的性狀作圖列出

油松1年生LL、LT、SLA、SD、LNC隨經(jīng)度的增加而增加；LT、SD、LNC隨緯度的增加而增加。油松多年生LT,SLA,SD隨經(jīng)度的增加而增加；LT、SD隨緯度的增加而增加。以上各性狀值與分布區(qū)經(jīng)緯度存在顯著相關(guān)關(guān)系(P<0.05),但各擬合方程的R2較低(0.05

2.2 油松葉性狀分布控制因素

2.2.1 葉功能性狀地理差異分析

除1年生LCC外油松各葉性狀值在不同研究地點(diǎn)之間皆存在顯著性差異(P<0.05),這種顯著性差異是油松對不同生境的適應(yīng)性變化產(chǎn)生的結(jié)果。各研究地之間環(huán)境因子存在較大的差異且各環(huán)境因子之間相互影響共同作用于油松,通過通徑分析有助于人們進(jìn)一步了解各環(huán)境因子對油松葉功能性狀的直接以及間接影響程度。

2.2.2 油松各葉性狀分布通徑分析

通過表3可知油松LL主要受到年平均降水量和海拔高度的影響,其中年降水量對其產(chǎn)生正效應(yīng)而海拔高度產(chǎn)生負(fù)效應(yīng)；LT主要受到年均溫和土壤體積含水率的影響,且都產(chǎn)生顯著的負(fù)效應(yīng)；SLA主要受到年平均降水量土壤氮含量的影響,海拔高度對油松的直接影響較小,但其通過其他環(huán)境因子對油松比葉面積起著顯著的負(fù)效應(yīng)；1年生LDMC主要受到土壤體積含水率與土壤氮含量的影響而多年生LDMC主要受到土壤磷含量和空氣相對濕度影響；SD主要受到海拔高度的影響,同時也在一定程度上受到年均溫和年降水量的影響,其中年降水量對氣孔密度產(chǎn)生正效應(yīng),年均溫度與海拔高度對氣孔密度產(chǎn)生負(fù)效應(yīng)；1年生LNC主要受到空氣相對濕度和土壤磷含量的影響而多年生LNC主要受到海拔和年降水量的影響；土壤磷含量對LPC產(chǎn)生顯著的正效應(yīng),年降水量對1年生LPC影響較小但對多年生LPC產(chǎn)生顯著的負(fù)效應(yīng)。

3 討論

3.1 油松各葉性狀分布特征

由于物種的分布區(qū)內(nèi)生境不是一致的,物種要對不同的生境做出適應(yīng)性變化。功能性狀的變化就是這種適應(yīng)性變化的表現(xiàn)或反應(yīng)形式之一,每個物種都有一定的適應(yīng)幅度,在功能性狀方面它們的變化也是有限的[15]。Hallik等[35]的研究表明功能性狀的種間變異系數(shù)可高達(dá)95%而種內(nèi)變異系數(shù)則不超過30%。本研究中油松各功能性狀的種內(nèi)變異系數(shù)在4.82%—25.85%之間,支持上述結(jié)論。雖然油松各性狀種內(nèi)變異系數(shù)沒有超過30%,但在不同研究地點(diǎn)間各功能性狀值都存在顯著差異(表2),由此可見種內(nèi)變異雖然較小卻也是生物對不同生境適應(yīng)性變化的直接體現(xiàn),因此在研究植物功能性狀與環(huán)境間相關(guān)關(guān)系時,不能簡單的使用性狀平均值代替物種而忽略種內(nèi)變異。

隨著全球環(huán)境變化問題的日益凸顯,對于能夠響應(yīng)和影響環(huán)境變化的重要植物葉功能性狀在全球分布格局上的研究受到了更大的關(guān)注[36],這些研究主要針對全球范圍內(nèi)的不同植物葉片、種子、樹干性狀[37-40]等,本研究發(fā)現(xiàn)油松1年生LL、LT、SLA、SD、LNC呈現(xiàn)出經(jīng)度格局；LT、SD及LNC呈現(xiàn)出緯度格局。油松多年生LT、SLA、SD呈現(xiàn)出經(jīng)度格局；LT、SD呈現(xiàn)緯度格局,并且都表現(xiàn)出顯著的相關(guān)性,其余各性狀則并未表現(xiàn)出與經(jīng)緯度的相關(guān)性。由此可知油松物理性狀相較與化學(xué)性狀而言體現(xiàn)出更明顯的經(jīng)緯度分布格局,這與不同經(jīng)緯度地區(qū)溫度以及降水的梯度性變化息息相關(guān),因此在探討植物物種功能性狀分布格局及其驅(qū)動力時,選擇物理性狀更為合適。根據(jù)經(jīng)典的植被、物種與氣候的關(guān)系,全球、洲際和區(qū)域尺度的氣候與植物功能性狀應(yīng)該有密切的聯(lián)系,經(jīng)向和緯向地形地貌的地理空間變異對植物葉功能性狀的影響主要來自于氣溫、降水和土壤特性的梯度變化[41]。在我國油松天然分布區(qū)內(nèi),隨著緯度增加,不僅氣溫降低并且降水量減少[23],這可能是油松葉性狀緯度分布格局形成的主要原因；自西向東由于我國地形地貌以及降水量的變化油松分布平均海拔高度下降[16],影響油松不同分布區(qū)內(nèi)的溫度與降水情況,這可能是油松葉性狀經(jīng)度分布格局形成的主要原因。然而由于同一地區(qū)不同個體之間各性狀存在很大的變異(表2)即存在非常大的局部效應(yīng)(Local effect),顯著降低了大尺度上葉性狀的經(jīng)緯度格局,導(dǎo)致各性狀值雖與經(jīng)緯度存在顯著的相關(guān)性(0.05

3.2 油松各葉性狀分布主控因素

大量研究表明,植物葉功能性狀與氣候條件、土壤養(yǎng)分、海拔高度等存在顯著的相關(guān)關(guān)系[4,42-45]。在區(qū)域尺度上,植物功能性狀與氣候的關(guān)系更加緊密,但不同區(qū)域間和不同性狀間也存在不少差異[4]。本研究發(fā)現(xiàn)油松各葉功能性狀受到多種因素的共同作用,但不同的性狀主控因素各不相同,1年生葉與多年生葉性狀間控制因素也存在一定差異(表3),因此在對油松等針葉樹種葉功能性狀進(jìn)行研究時,葉齡對其的影響不容忽視。

LL主要受到年平均降水量和海拔高度的影響(表3)。Jacobs[46]等研究表明葉片長寬與年平均降水量正相關(guān),這與本文研究結(jié)果一致。水分是植物細(xì)胞擴(kuò)張生長的動力,植物細(xì)胞在擴(kuò)張生長的過程中,需要充足的水分使細(xì)胞產(chǎn)生膨脹壓力,隨著降水量的增加油松葉片在生長發(fā)育過程中可利用的水分含量增加,充足的水分條件有利于油松葉細(xì)胞的擴(kuò)張生長,增加油松針葉長度。另外在油松天然分布區(qū)內(nèi)隨著海拔的升高,降水量與溫度降低(表1),溫度降低會使油松生長季縮短,進(jìn)而導(dǎo)致其葉長減小。年均溫與土壤體積含水率對油松LT產(chǎn)生顯著的負(fù)效應(yīng)(表3)。大量研究表明葉厚與溫度存在負(fù)相關(guān)關(guān)系,Higuchi等[47]認(rèn)為這是由于隨著溫度升高植物葉片柵欄組織與海綿組織細(xì)胞的層數(shù)及厚度減少導(dǎo)致的而劉全宏等[48]認(rèn)為這可能是低溫影響到葉原基的分化以及幼葉細(xì)胞的生長導(dǎo)致的,本研究支持上述結(jié)論但其機(jī)理還有待進(jìn)一步研究。土壤體積含水率體現(xiàn)了當(dāng)年生長季油松林土壤中可利用水分狀況,其數(shù)值能在一定程度上表征當(dāng)年氣候干旱程度。在干旱環(huán)境下,植物葉片傾向于變得小而厚,葉肉細(xì)胞排列緊密,細(xì)胞很少出現(xiàn)質(zhì)壁分離,這樣的組織結(jié)構(gòu)可以較好的吸收和保持水分,提高水分傳輸效率[49],是植物對環(huán)境的一種適應(yīng)機(jī)制。

表2 葉性狀方差分析表

LL： 葉長 Leaf length; LT： 葉厚 Leaf thickness; SLA： 比葉面積 Specific leaf area; LDMC： 葉干物質(zhì)含量 Leaf dry matter content; SD： 氣孔密度 Stomatal density; LNC： 葉氮含量 Leaf nitrogen concentration; LPC： 葉磷含量 Leaf phosphorus concentration; LCC： 葉碳含量 Leaf carbon concentration;*P<0.05; **P<0.01.

SLA主要受到年平均降水量、土壤氮含量的影響。海拔高度的直接通徑系數(shù)較小分別為-0.009和-0.004,但其通過年降水量以及土壤氮含量對油松SLA起著顯著的負(fù)效應(yīng)(表3)。Baigorri等(1999)[50]研究表明,比葉重與環(huán)境中水分狀況密切相關(guān),在氣候干燥、土壤水分短缺的環(huán)境中生長的植物葉片比葉重增加,而相對含水量降低,從而導(dǎo)致SLA減小。以往研究發(fā)現(xiàn)多數(shù)植物的SLA都是隨著海拔的上升減小,導(dǎo)致這種現(xiàn)象的根本原因是隨著海拔上升溫度的降低及生長季的縮短[51],但在本研究中油松SLA隨海拔升高而降低主要是由于降水量的減少。一方面降水減少會顯著增加油松葉片厚度(表3),導(dǎo)致油松SLA降低；另一方面本研究中海拔與緯度共同對溫度產(chǎn)生影響,各研究地區(qū)海拔較高之處緯度較低,從而削弱了海拔對溫度的影響作用。通過表3可以發(fā)現(xiàn)土壤氮含量對葉氮含量會產(chǎn)生微弱的正效應(yīng),但對油松比葉面積卻產(chǎn)生顯著的負(fù)效應(yīng)(表3),而大量研究表明葉氮含量與比葉面積存在正相關(guān)關(guān)系。由此可見環(huán)境條件雖然會對植物不同性狀產(chǎn)生不同的影響但植物功能性狀間的相互關(guān)系主要是由于植物對資源利用的權(quán)衡策略決定的,氣候及環(huán)境條件對其的影響是適中的[4]。

LDMC能夠在一定程度上反映植物對其生境資源的利用狀況,并且它不容易受到其它葉片性狀的限制,如葉片組分和葉片厚度,但是容易受到水因子的束縛[52]。油松LDMC控制因素較為復(fù)雜,其中1年生LDMC主控因素為土壤體積含水率與土壤氮含量；多年生LDMC主控因素為土壤磷含量和相對濕度(表3)?？傮w上LDMC主要受到水分條件以及土壤養(yǎng)分的控制,但在不同生長階段對其產(chǎn)生影響的關(guān)鍵因子并不相同。這可能是由于油松分布區(qū)較廣,具有豐富的遺傳多樣性[21],本研究未考慮遺傳因素的影響而導(dǎo)致的。由圖2、圖3可以看出油松葉干物質(zhì)含量總體變異系數(shù)最小,其在各樣點(diǎn)間的差異可能主要是由于地區(qū)間的遺傳差異導(dǎo)致的,環(huán)境因子對其解釋力較低。

SD主要受到海拔高度的影響(表3),與徐鵬彬等[53]對色季拉山東坡急尖長荀冷杉葉氣孔密度的研究結(jié)果一致,同時也受到年均溫(負(fù)效應(yīng))和年降水量(正效應(yīng))的影響。費(fèi)松林等[54]研究表明,溫度升高有利于植物氣孔密度增加,李芬蘭等[49]分析認(rèn)為溫度升高導(dǎo)致葉片蒸騰作用加強(qiáng),植物會面臨大量的失水,因此小而密的氣孔可以防止由于蒸騰而導(dǎo)致的水分散失,這與本研究得到的結(jié)果相反,分析可能原因有1):溫度、海拔、降水都對油松氣孔密度產(chǎn)生顯著的影響(表3),而海拔差異通常導(dǎo)致溫度、降水、光照和土壤等環(huán)境因子發(fā)生相應(yīng)的變化,造成熱量和水分重新分配,使植物生長的環(huán)境條件更為復(fù)雜,也使植物葉功能性狀對環(huán)境產(chǎn)生復(fù)雜的適應(yīng)性變化[7],它們之間通過復(fù)雜的相互關(guān)系對油松氣孔密度產(chǎn)生影響,不能單獨(dú)用溫度解釋其變化原因。2)研究物種的差異：研究發(fā)現(xiàn)不同物種氣孔密度與溫度降水海拔等的關(guān)系并不相同[55-57]。

油松1年生LNC主要受到空氣相對濕度(負(fù)效應(yīng))以及土壤磷含量的影響。多年生LNC主要受到海拔高度(正效應(yīng))和年平均降水量(負(fù)效應(yīng))的影響,而土壤氮含量對其影響較小(表3)。以往有大量針對植物葉片氮磷含量與溫度降水以及海拔之間相互關(guān)系的研究,發(fā)現(xiàn)氣候越是干旱,葉片具有越高的氮磷含量[42],植物葉片的氮磷含量隨海拔升高而增加[45],這與本研究得到的結(jié)論基本相同。而水分是植物生長發(fā)育所必需的,與其他因子相比,植物功能性狀對水分的響應(yīng)更為顯著,但研究結(jié)果發(fā)現(xiàn)油松葉片在不同生長發(fā)育階段對其葉氮含量起到關(guān)鍵性作用的水因子并不相同(表3)。LPC主要受到土壤磷含量的影響,其中年降水量還對多年生葉磷含量產(chǎn)生顯著的負(fù)效應(yīng)(表3)。磷元素作為構(gòu)成細(xì)胞的重要元素,植物對磷元素的獲取主要是通過根系從土壤中吸收,而油松作為演替后期地帶性頂級群落,由于生產(chǎn)力的下降、新鮮凋落物與腐殖質(zhì)中N∶P比的增大,會導(dǎo)致其受到P元素的限制[58],因此土壤中磷元素的含量極大的限制了油松葉磷含量。

3.3 磷元素對油松生長的限制性作用

國外研究表明,一個穩(wěn)定且長期未受到災(zāi)難性干擾的森林生態(tài)系統(tǒng),在其演替后期由于生產(chǎn)力下降、新鮮凋落物與腐殖質(zhì)中N:P比的增大,會導(dǎo)致該森林生態(tài)系統(tǒng)逐漸受到P元素的限制,并隨著時間的推移,這種限制性作用隨之增強(qiáng)[58]。而國內(nèi)研究也發(fā)現(xiàn)中國植被的氮磷比較全球平均值高[41],容易受到磷元素的限制。本研究發(fā)現(xiàn)油松1年生葉磷含量為1.3%(圖2)多年生葉葉磷含量為0.077%(圖3),明顯低于Reich等研究發(fā)現(xiàn)的全球平均水平0.15%[37]。由通徑分析的結(jié)果可知油松LT、LDMC、LNC、LPC都在一定程度上受到土壤磷含量的影響,由此我們推斷,在油松天然分布區(qū)內(nèi),磷元素含量可能是各油松成熟林生長的重要限制性因素。

4 結(jié)論

1) 油松各葉性狀值基本滿足正態(tài)分布趨勢且都存在較大的標(biāo)準(zhǔn)差及變異系數(shù),其中LPC含量變異系數(shù)最大,1年生葉為23.29%,多年生葉為25.85%；LDMC變異系數(shù)最小,一年葉為5.91%,多年葉為4.82%。部分葉功能性狀分布存在較弱經(jīng)緯度格局,水熱條件的變化以及較大的局部效應(yīng)是造成這種經(jīng)緯度格局的主要原因。

2) 除1年生LCC外油松各葉功能性狀在不同研究區(qū)域間差異極顯著且在同一區(qū)域不同個體間存在較大變異。

3) 油松各葉性狀分布的主控因素不同。其中LL主要受到年平均降水量和海拔高度的影響；LT主要受到年均溫和土壤體積含水率的影響；SLA主要受到年平均降水量和土壤氮含量的影響；LDMC主要受到水因子的影響；SD主要受到海拔高度的影響；LNC主要受到海拔和水因子的影響；LPC主要受到土壤磷含量的影響。

[1] Violle C, Navas M L, Vile D, Kazakou E, Fortunel C, Hummel I, Garnier E. Let the concept of trait be functional!. Oikos, 2007, 116(5): 882- 892.

[2] Bernhardt-R?mermann M, R?mermann C, Nuske R, Parth A, Klotz S, Schmidt W, Stadler J. On the identification of the most suitable traits for plant functional trait analyses. Oikos, 2008, 117(10): 1533- 1541.

[3] Westoby M. A leaf-height-seed (LHS) plant ecology strategy scheme. Plant and Soil, 1998, 199(2): 213- 227.

[4] Wright I J, Reich P B, Westoby M, Ackerly D D, Baruch Z, Bongers F, Cavender-Bares J, Chapin T, Cornelissen J H C, Diemer M, Flexas J, Garnier E, Groom P K, Gulias J, Hikosaka K, Lamont B B, Lee T, Lee W, Lusk C, Midgley J J, Navas M L, Niinemets ü, Oleksyn J, Osada N, Poorter H, Poot P, Prior L, Pyankov V I, Roumet C, Thomas S C, Tjoelker M G, Veneklaas E J, Villar R. The worldwide leaf economics spectrum. Nature, 2004, 428(6985): 821- 827.

[5] Wright I J, Reich P B, Cornelissen J H C, Falster D S, Groom P K, Hikosaka K, Lee W, Lusk C H, Niinemets ü, Oleksyn J, Osada N, Poorter H, Warton D I, Westoby M. Modulation of leaf economic traits and trait relationships by climate. Global Ecology and Biogeography, 2005, 14(5): 411- 421.

[6] 祁建, 馬克明, 張育新. 遼東櫟(Quercusliaotungensis)葉特性沿海拔梯度的變化及其環(huán)境解釋. 生態(tài)學(xué)報, 2007, 27(3): 930- 937.

[7] 張慧文, 馬劍英, 孫偉, 陳發(fā)虎. 不同海拔天山云杉葉功能性狀及其與土壤因子的關(guān)系. 生態(tài)學(xué)報, 2010, 30(21): 5747- 5758.

[8] 董莉莉, 劉世榮, 史作民, 馮秋紅. 中國南北樣帶上栲屬樹種葉功能性狀與環(huán)境因子的關(guān)系. 林業(yè)科學(xué)研究, 2009, 22(4): 463- 469.

[9] 馮秋紅, 史作民, 董莉莉, 劉世榮. 南北樣帶溫帶區(qū)櫟屬樹種功能性狀間的關(guān)系及其對氣象因子的響應(yīng). 植物生態(tài)學(xué)報, 2010, 34(6): 619- 627.

[10] Falsenstein J. Phylogenies and the comparative method. The American Naturalist, 1985, 125(1): 1- 15.

[11] Wright I J, Ackerly D D, Bongers F, Harms K E, Ibarra-Manriquez G, Martinez-Ramos M, Mazer S J, Muller-Landau H C, Paz H, Pitman N C A, Poorter L, Silman M R, Vriesendorp C F, Webb C O, Westoby M, Wright S J. Relationships among ecologically important dimensions of plant trait variation in seven Neotropical forests. Annals of Botany, 2007, 99(5): 1003- 1005.

[12] Byars S G, Papst W, Hoffmann A A. Local adaptation and cogradient selection in the alpine plant,Poahiemata, along a narrow altitudinal gradient. Evolution, 2007, 61(12): 2925- 2941.

[13] 戚德輝, 溫仲明, 楊士梭, 王紅霞, 郭茹. 基于功能性狀的鐵桿蒿對環(huán)境變化的響應(yīng)與適應(yīng). 應(yīng)用生態(tài)學(xué)報, 2015, 26(7): 1921- 1927.

[14] 楊利民, 韓梅, 周廣勝, 李建東. 中國東北樣帶關(guān)鍵種羊草水分利用效率與氣孔密度. 生態(tài)學(xué)報, 2007, 27(1): 16- 24.

[15] 蘇文華, 施展, 楊波, 楊建軍, 趙冠華, 周睿. 滇石櫟沿緯度梯度葉片功能性狀的種內(nèi)變化. 植物分類與資源學(xué)報, 2015, 37(3): 309- 317.

[16] 徐化成. 油松. 北京: 中國林業(yè)出版社, 1990: 18- 23.

[17] 楊洪曉, 吳波, 張金屯, 林德榮, 常順利. 森林生態(tài)系統(tǒng)的固碳功能和碳儲量研究進(jìn)展. 北京師范大學(xué)學(xué)報: 自然科學(xué)版, 2005, 41(2): 172- 177.

[18] 徐曉春, 榮文琛, 吳天林, 岳水林. 馬尾松種源、林分和單株遺傳變異的初步研究. 林業(yè)科技通訊, 1990, (8): 4- 7.

[19] 翁殿伊, 王同立, 楊井泉. 油松地理種源選擇試驗(yàn)初探. 河北林業(yè)科技, 1983, (3): 5- 8.

[20] 徐化成, 孫肇鳳, 郭廣榮, 馮林. 油松天然林的地理分布和種源區(qū)的劃分. 林業(yè)科學(xué), 1981, 17(3): 258- 270.

[21] 王磊, 樊軍鋒, 劉永紅, 楊培華, 王孟昌. 我國油松主要分布區(qū)種質(zhì)資源遺傳多樣性. 東北林業(yè)大學(xué)學(xué)報, 2009, 37(12): 3- 7, 13- 13.

[22] 徐化成, 唐謙, 張淑娟, 馬虹. 油松氣候生態(tài)型的研究. 林業(yè)科學(xué), 1986, 22(1): 10- 20.

[23] 徐化成. 油松地理變異和種源選擇. 北京: 中國林業(yè)出版社, 1992: 6- 16.

[24] Diaz S, Cabido M, Casanoves F. Plant functional traits and environmental filters at a regional scale. Journal of Vegetation Science, 1998, 9(1): 113- 122.

[25] 張林, 金冬梅. 針葉樹種葉面積的主要測定方法//第八屆全國生物多樣性保護(hù)與持續(xù)利用研討會論文集. 金華: 國際生物多樣性計劃中國委員會, 中科院生物多樣性委員會, 2008.

[26] 鮑士旦. 土壤農(nóng)化分析. 北京: 中國農(nóng)業(yè)出版社, 2000: 135- 155.

[27] 張赟, 趙亞洲, 張春雨, 趙秀海. 北京松山油松種群結(jié)構(gòu)及空間分布格局. 應(yīng)用與環(huán)境生物學(xué)報, 2009, 15(2): 175- 179.

[28] 張赟, 張春雨, 趙秀海. 內(nèi)蒙古黑里河天然油松林主要樹種的空間分布格局. 西北植物學(xué)報, 2009, 29(1): 167- 173.

[29] 奇凱, 張春雨, 侯繼華, 吳蕾, 鄒璐, 趙秀海. 赤峰市黑里河天然油松林空間分布格局研究. 水土保持研究, 2010, 17(4): 239- 242.

[30] 陳振舉, 孫雨, 何興元, 陳瑋, 邵雪梅, 張華宇, 王忠鈺, 劉筱彧. 千山油松年輪寬度年表的建立及其與氣候的關(guān)系. 應(yīng)用生態(tài)學(xué)報, 2007, 18(10): 2191- 2201.

[31] 趙鵬祥, 趙臻, 郝紅科. 黃龍山蔡家川林場森林類型碳密度及其變化研究. 西北農(nóng)林科技大學(xué)學(xué)報: 自然科學(xué)版, 2011, 39(7): 77- 82.

[32] 許浩, 張源潤, 胡天華, 季波, 何建龍, 蔡進(jìn)軍, 王繼飛. 寧夏賀蘭山自然保護(hù)區(qū)油松林碳儲量及分布格局. 生態(tài)環(huán)境學(xué)報, 2013, 32(11): 1785- 1789.

[33] 吳昊, 張明霞, 王得祥. 秦嶺南坡油松-銳齒槲櫟混交林群落不同層次多樣性特征及環(huán)境解釋. 西北植物學(xué)報, 2013, 33(10): 2086- 2094.

[34] 李永華, 羅天祥, 盧琦, 田曉婭, 吳波, 楊恒華. 青海省沙珠玉治沙站17種主要植物葉性因子的比較. 生態(tài)學(xué)報, 2005, 25(5): 994- 999.

[35] Hallik L, Niinemets ü, Wright I J. Are species shade and drought tolerance reflected in leaf-level structural and functional differentiation in Northern Hemisphere temperate woody flora?. New Phytologist, 2009, 184(1): 257- 274.

[36] Wright I J, Reich P B, Westoby M. Strategy shifts in leaf physiology, structure and nutrient content between species of high-and low-rainfall and high-and low-nutrient habitats. Functional Ecology, 2001, 15(4): 423- 434.

[37] Reich P B, Oleksyn J. Global patterns of plant leaf N and P in relation to temperature and latitude. Proceedings of the National Academy of Sciences of the United States of America, 2004, 101(30): 11001- 11006.

[38] Moles A T, Ackerly D D, Tweddle J C, Dickie J B, Smith R, Leishman M R, Mayfield M M, Pitman A, Wood J T, Westoby M. Global patterns in seed size. Global Ecology and Biogeography, 2007, 16(1): 109- 116.

[39] Moles A T, Warton D I, Warman L, Swenson N G, Laffan S W, Zanne A E, Pitman A, Hemmings F A, Leishman M R. Global patterns in plant height. Journal of Ecology, 2009, 97(5): 923- 932.

[40] Chave J, Coomes D, Jansen S, Lewis S L, Swenson N G, Zansen A E. Towards a worldwide wood economics spectrum. Ecology Letters, 2009, 12(4): 351- 366.

[41] 孟婷婷, 倪健, 王國宏. 植物功能性狀與環(huán)境和生態(tài)系統(tǒng)功能. 植物生態(tài)學(xué)報, 2007, 31(1): 150- 165.

[42] Prior L D, Bowman D M J S, Eamus D. Seasonal differences in leaf attributes in Australian tropical tree species: family and habitat comparisons. Functional Ecology, 2004, 18(5): 707- 718.

[43] Craine J M, Lee W G. Covariation in leaf and root traits for native and non-native grasses along an altitudinal gradient in New Zealand. Oecologia, 2003, 134(4): 471- 478.

[44] H?lscher D, Schmitt S, Kupfer K. Growth and leaf traits of four broad-leaved tree species along a hillside gradient. Forstwissenschaftliches Centralblatt, 2002, 121(5): 229- 239.

[45] Han W X, Fang J Y, Guo D L, Zhang Y. Leaf nitrogen and phosphorus stoichiometry across 753 terrestrial plant species in China. New Phytologist, 2005, 168(2): 377- 385.

[46] Jacobs B F. Estimation of rainfall variables from leaf characters in tropical Africa. Paliogeography, Palioclimatology, Palioecology, 1999, 145(1/3): 231- 250.

[47] Higuchi H, Sakuratani T, Utsunomiya N. Photosynthesis, leaf morphology, and shoot growth as affected by temperatures in cherimoya (AnnonacherimolaMill.) trees. Scientia Horticulturae, 1999, 80(1/2): 91- 104.

[48] 劉全宏, 王孝安, 田先華, 肖婭萍. 太白紅杉(Larixchinensis)葉的形態(tài)解剖學(xué)特征與環(huán)境因子的關(guān)系. 西北植物學(xué)報, 2001, 21(5): 885- 893.

[49] 李芳蘭, 包維楷. 植物葉片形態(tài)解剖結(jié)構(gòu)對環(huán)境變化的響應(yīng)與適應(yīng). 植物學(xué)通報, 2005, 22(S): 118- 127.

[50] Baigorri H, Antolín M C, Sánchez-Díaz M. Reproductive response of two morphologically different pea cultivars to drought. European Journal of Agronomy, 1999, 10(2): 119- 128.

[51] 楊冬梅, 章佳佳, 周丹, 錢敏杰, 鄭瑤, 金靈妙. 木本植物莖葉功能性狀及其關(guān)系隨環(huán)境變化的研究進(jìn)展. 生態(tài)學(xué)雜志, 2012, 31(3): 702- 713.

[52] Wilson P J, Thompson K, Hodgson J G. Specific leaf area and leaf dry matter content as alternative predictors of plant strategies. New Phytologist, 1999, 143(1): 155- 162.

[53] 徐鵬彬. 中國溫帶地區(qū)4種針葉林樹種葉屬性對海拔梯度的響應(yīng)研究[D]. 蘭州: 蘭州大學(xué), 2012.

[54] 費(fèi)松林, 方精云, 樊擁軍, 趙坤, 劉雪皎, 崔克明. 貴州梵凈山亮葉水青岡葉片和木材的解剖學(xué)特征及其與生態(tài)因子的關(guān)系. 植物學(xué)報, 1999, 41(9): 1002- 1009.

[55] Sam O, Jeréz E, Dell′Amico J, Ruiz-Sanchez M C. Water stress induced changes in anatomy of tomato leaf epidermes. Biologia Plantarum, 2000, 43(2): 275- 277.

[56] 王瑞麗, 于貴瑞, 何念鵬, 王秋鳳, 趙寧, 徐志偉. 氣孔特征與葉片功能性狀之間關(guān)聯(lián)性沿海拔梯度的變化規(guī)律: 以長白山為例. 生態(tài)學(xué)報, 2016, 36(8): 1- 10.

[57] 張曉艷, 楊惠敏, 侯宗東, 王根軒. 土壤水分和種植密度對春小麥葉片氣孔的影響. 植物生態(tài)學(xué)報, 2003, 27(1): 133- 136.

[58] Wardle D A, Walker L R, Bardgett R D. Ecosystem properties and forest decline in contrasting long-term chronosequences. Science, 2004, 305(5683): 509- 513.

Leaf functional trait distribution and controlling factors ofPinustabuliformis

ZHANG Kai1, HOU Jihua1,*, HE Nianpeng2

1BeijingForestryUniversity,Beijing100083,China2InstituteOfGeographicSciencesAndNaturalResourcesResearch,ChineseAcademyofScience,Beijing100101,China

Leaf functional traits can directly or indirectly reflect the adaptation strategy of plants to the environment, influencing their survival, growth, and reproduction. Until recently, there have been many studies on the relationship between leaf functional traits and environment, but most have used the single or mean value of the traits to represent one species, ignoring considerable intraspecific variations.Pinustabuliformisis one of the main afforestation tree species in China, playing an important role in terrestrial ecosystems. Thus, the relationship between the leaf functional traits ofP.tabuliformisand environmental factors could elucidate the relationship between intraspecific variations and environmental factors. The results of this study could provide a theoretical basis for better protection, development, and utilization ofP.tabuliformisforests under global climate change. In the present paper, we analyzed the environment factors and distribution of eight leaf functional traits ofP.tabuliformisin its natural distribution range, including leaf length (LL), leaf thickness (LT), specific leaf area (SLA), leaf dry matter content (LDMC), stomatal density (SD), leaf nitrogen concentration (LNC), leaf phosphorus concentration (LPC), and leaf carbon concentration (LCC), from leaves collected in Inner Mongolia, Liaoning, Beijing, Shanxi, Shaanxi, Ningxia and Qinghai Provinces, from June to September 2014. The results showed that 1) The leaf functional traits ofP.tabuliformisshowed significant intraspecific variations, and the variation coefficient was between 4.82% and 25.85%. The traits, except the LCC of annual leaves, showed significant differences between different research sites (P<0.05), the intraspecific variation of LPC was the highest (23.29% for annual leaves and 25.85% for perennial leaves), and the intraspecific variation of LDMC was the lowest (5.91% for annual leaves and 4.82% for perennial leaves). 2) There were significant correlations between the leaf trait and the longitude and latitude (P<0.05). LL, LT, SLA, SD, and LNC ofP.tabuliformisleaves exhibited a weak longitude distribution pattern, and LT, SD, and LNC exhibited a weak latitude distribution pattern. The combined effect of changing hydrothermal conditions along longitude and latitude and strong local effects could contribute to the formation of these patterns. There was also significant variation among different individuals at the same site, suggesting large local effects, which could weaken the longitude and latitude distribution pattern on the large scale, leading to a smallRvalue (0.05

Pinustabuliformis; leaf functional traits; distribution characteristics; main controlling factor; path analysis

國家自然科學(xué)基金青年基金項(xiàng)目(31000263)；教育部科學(xué)技術(shù)重點(diǎn)項(xiàng)目(107014)； 生物多樣性與生態(tài)工程教育部重點(diǎn)實(shí)驗(yàn)室開放課題資助；中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金(TD2011-07)；國家自然科學(xué)基金面上項(xiàng)目(31470506)；中國科學(xué)院地理研究所“可楨杰出青年學(xué)者”項(xiàng)目(2013RC102)資助

2015- 08- 29;

日期:2016- 06- 13

10.5846/stxb201508291793

*通訊作者Corresponding author.E-mail: houjihua@bjfu.edu.cn

張凱,侯繼華,何念鵬.油松葉功能性狀分布特征及其控制因素.生態(tài)學(xué)報,2017,37(3):736- 749.

Zhang K, Hou J H, He N P.Leaf functional trait distribution and controlling factors ofPinustabuliformis.Acta Ecologica Sinica,2017,37(3):736- 749.