戴耀良,何國強(qiáng),藺萬煌,尹婷輝,
(1.深圳市仙湖植物園,廣東 深圳 518001; 2.深圳大學(xué)有害生物防控中心,廣東 深圳 518060;3.湖南農(nóng)業(yè)大學(xué)生物科學(xué)技術(shù)學(xué)院,湖南 長沙 410128)
道路綠地是城市綠地系統(tǒng)的有機(jī)構(gòu)成要素,是城市的交通空間,它與建筑空間的地位同等重要。一般在城市干道上,主要設(shè)置行道樹和各種綠帶(包括人行道綠帶、分車綠帶和路側(cè)綠帶),并按照道路綠地設(shè)計(jì)的基本原則進(jìn)行有規(guī)律的布局,使道路綠化在變化和統(tǒng)一中完成,從而形成優(yōu)良的城市道路綠地景觀效果[1]。分車帶作為城市道路車行道及人行道的分隔帶,貫穿整個(gè)城市,其綠化植物的生長狀況直接影響了分車帶景觀及城市生態(tài)效果。
深圳地處南亞熱帶,屬亞熱帶海洋性氣候,依山傍海,氣候溫暖宜人,良好的氣候和地理環(huán)境,為植物的生長提供了得天獨(dú)厚的條件[2]。但部分道路存在忽視分車帶寬度片面追求復(fù)合式多層次的植物配置模式、群落結(jié)構(gòu)不合理、植物搭配不當(dāng)?shù)葍A向,直接影響城市道路的景觀效果和生態(tài)效益。目前,主成分分析法廣泛用于農(nóng)林業(yè)生產(chǎn)、園林綠化及其他植物栽培等方面的研究[3-12],本研究基于主成分分析法探討影響城市道路分車帶綠地群落的變量因子,了解其制約因素,以期為城市綠化實(shí)踐提供相應(yīng)的技術(shù)支持。
根據(jù)深圳市道路分車帶寬度及綠化配置模式類型,選擇寬度為10、6、8、4 m代表性強(qiáng)的各類型分車綠化帶,共設(shè)置樣地31個(gè),面積合計(jì)5825 m2,中央分車帶樣地11個(gè),快慢分車帶樣地14個(gè),慢車道與人行道分車帶樣地6個(gè),樣地長度統(tǒng)一為10 m。涵蓋寶安大道、北環(huán)大道、濱河大道、龍城大道、龍翔大道、深南大道、濱海大道、沙河?xùn)|路、裕安一路、新湖路等,其中南北走向的道路13條,東西走向的道路18條。
調(diào)查內(nèi)容包括樣地面積、位置、分車帶寬度、群落結(jié)構(gòu)等相關(guān)信息。將群落因子調(diào)查匯總為喬木物種數(shù)、喬木株數(shù)密度、灌木物種數(shù)、灌木株數(shù)密度、地被物種數(shù)、地被面積、草坪物種數(shù)、草坪面積、群落層次、上層郁閉度、總郁閉度、重疊郁閉度、總物種數(shù)、喬灌總株數(shù)等14個(gè)因子。因本項(xiàng)目研究的樣地立足于土壤、氣候條件及養(yǎng)護(hù)措施相同,立地因子僅對(duì)分車帶寬度進(jìn)行調(diào)查分析。
采用Excel 2007和SPSS 23.0軟件對(duì)樣地進(jìn)行數(shù)據(jù)處理和運(yùn)算,其中SPSS 23.0自動(dòng)對(duì)原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理。
主成分降維分析的條件是要求KMO檢驗(yàn)滿足KMO值>0.5且bartlett球形度檢驗(yàn)的P值<0.05,才有統(tǒng)計(jì)學(xué)意義[13]。因此,將14個(gè)群落調(diào)查因子及1個(gè)立地因子(分車帶寬度)結(jié)合植物群落總體狀況進(jìn)行分析,通過對(duì)不同因子逐一刪減組合分別進(jìn)行多次演算比對(duì),篩選確定灌木物種數(shù)(X1)、喬木物種數(shù)(X2)、喬木株數(shù)密度(X3)、群落層次(X4)、上層郁閉度(X5)、總郁閉度(X6)、總物種數(shù)(X7)、喬灌總株數(shù)(X8)、分車帶寬度(X9)及地被面積(X10)為滿足主成分降維分析檢驗(yàn)要求的最佳10個(gè)因子組合。
本項(xiàng)目以SPSS 23.0軟件進(jìn)行操作運(yùn)算。其中軟件自動(dòng)運(yùn)算出初始因子載荷矩陣,每個(gè)載荷量表示主成分與對(duì)應(yīng)變量的相關(guān)關(guān)系,但不可直接使用于主成分表達(dá)式中[11-13],主成分表達(dá)式特征向量計(jì)算式為:Zij=aij/SQRT。
由于不同的指標(biāo)存在不同的量綱,而不同量綱的指標(biāo)是不能進(jìn)行相互比較的[14]。因此,為了消除由于量綱和數(shù)量不同可能帶來的影響,對(duì)統(tǒng)計(jì)的原始數(shù)據(jù)必須按主成分分析的方法進(jìn)行標(biāo)準(zhǔn)化處理,調(diào)查因子的數(shù)據(jù)經(jīng)標(biāo)準(zhǔn)化處理后的相關(guān)系數(shù)矩陣詳見表1。
主成分分析是從眾多的原始指標(biāo)中綜合出少數(shù)具有代表性的因子,因而要求原有指標(biāo)之間具有相關(guān)性[15-16]。在對(duì)影響因子相關(guān)性關(guān)系的45個(gè)比對(duì)中,有35個(gè)呈正相關(guān)關(guān)系,占總比對(duì)數(shù)的77.78%;10個(gè)呈負(fù)相關(guān)關(guān)系,占總比對(duì)數(shù)的22.22%。表1相關(guān)性矩陣數(shù)據(jù)說明,各因子的選擇具有充分的合理性。
標(biāo)準(zhǔn)化處理后的數(shù)據(jù)經(jīng)KMO和 Bartlett 球度檢驗(yàn)的結(jié)果詳見表2。結(jié)果表明,KMO值為0.527>0.5,且Bartlett球度檢驗(yàn)相應(yīng)的顯著水平為0.000,適合做主成分分析;10個(gè)變量因子的公因子方差提取均大于0.8,其中6個(gè)變量超過0.9。這說明提取的公因子能夠很好地反映原始變量的主要信息。
表1 標(biāo)準(zhǔn)化數(shù)據(jù)的相關(guān)性矩陣
經(jīng)主成分降維分析得到相關(guān)矩陣的特征值及貢獻(xiàn)率詳見表3。結(jié)果表明前5個(gè)主成分的累積方差達(dá)88.365%。實(shí)際中通常取累計(jì)貢獻(xiàn)率達(dá)到85%的m個(gè)主成分[17-18]。從第1~第4主成分的初始特征值均大于1,從第1~第5個(gè)主成分累計(jì)可以解釋>85%總方差。
表2 KMO和巴特利特檢驗(yàn)
因此,將影響道路分車帶綠地的10個(gè)指標(biāo)的全部信息歸結(jié)為若干個(gè)成分,總方差累計(jì)為88.365%的前5個(gè)主成分能充分反映道路分車帶綠地的大部分信息,即前5個(gè)主成分較好地保留了原來指標(biāo)的絕大部分信息量,因此只需選擇前5個(gè)主成分進(jìn)行分析。
表3 總方差解釋
設(shè)Y1為第1個(gè)主成分,Y2為第2個(gè)主成分…Y5為第5個(gè)主成分,第1、2、3、4、5主成分的因子載荷量aij見表4。從表4主成分aij值可見:在10個(gè)因子中具有較高正載荷量的因子依次是喬木株數(shù)密度(0.87)、喬木物種數(shù)(0.84)、總物種數(shù)(0.77)、喬灌總株數(shù)(0.73)、群落層次(0.70)、分車帶寬度(0.68)、上層郁閉度(0.66)、總郁閉度(0.64)、灌木物種數(shù)(0.61),但也表現(xiàn)出有較高負(fù)載荷量的是群落層次(-0.54)、地被面積(-0.53)、灌木物種數(shù)(-0.49)、總物種數(shù)(-0.45)、喬木株數(shù)密度(-0.31)。說明群落株數(shù)密度(特別是喬木株數(shù)密度)、總物種數(shù)、分車帶寬度、群落層次、上層郁閉度這幾個(gè)因子可較大地制約分車帶群落的整體狀況。
表4 因子載荷量矩陣
根據(jù)表3中第1、2、3、4、5主成分提取的載荷平方和分別為3.565、2.164、1.281、1.068、0.758,各主成分的特征向量根據(jù)公式Z1j=a1j/SQRT(3.565)、Z2j=a2j/SQRT(2.164)、Z3j=a3j/SQRT(1.281)、Z4j=a4j/SQRT(1.068)、Z5j=a5j/SQRT(0.758),計(jì)算可得出各主成分的表達(dá)式:
Y1=0.08X1+0.44X2+0.46X3+0.1X4-0.04X5+0.22X6+0.41X7+0.39X8+0.36X9+0.29X10
Y2=0.42X1-0.09X2+0.09X3+0.47X4+0.39X5+0.27X6-0.3X7+0.35X8-0.09X9-0.36X10
Y3=-0.43X1+0.12X2+0.07X3-0.23X4+0.58X5+0.57X6-0.07X7-0.27X8-0.02X9+0.04X10
Y4=0.56X1-0.06X2-0.3X3-0.53X4+0.11X5+0.17X6-0.18X7+0.18X8+0.01X9+0.47X10
Y5=-0.09X1-0.46X2+0.07X3+0.1X4+0.27X5-0.25X6-0.17X7-0.04X8+0.75X9+0.17X10
綜合表3、表4及主成分表達(dá)式,可見在第1主成分中,貢獻(xiàn)率是35.648%,以株數(shù)密度(X3、X8)影響最大,其次是喬木物種數(shù)(X2)及總物種數(shù)(X7),再次是分車帶寬度(X9)的影響,無負(fù)數(shù)的載荷量指標(biāo)數(shù),較高的載荷量為X3>X2>X7>X8>X9。在第2主成分中,貢獻(xiàn)率是21.642%,群落層次(X4)、灌木物種數(shù)(X1)及上層郁閉度(X5)對(duì)分車帶群落影響最大,較高的正載荷量為X4>X1>X5;總物種數(shù)(X7)、地被面積(X10)對(duì)綠化植物生長存在較大負(fù)相關(guān)影響,負(fù)載荷量絕對(duì)值X7>X10。在第3主成分中,貢獻(xiàn)率是12.815%,上層郁閉度(X5)及總郁閉度(X6)對(duì)分車帶群落影響最大,較高的正載荷量為X5>X6;灌木物種數(shù)(X1)、群落層次(X4)對(duì)分車帶群落結(jié)構(gòu)也有較大的負(fù)面影響,負(fù)載荷量絕對(duì)值X1>X4。在第4主成分中,貢獻(xiàn)率是10.675%,以灌木物種數(shù)(X1)、地被面積(X10)影響大,載荷量為X1>X10;有較高負(fù)載荷量的為群落層次。在第5主成分中,貢獻(xiàn)率是7.584%,有較高正載荷量的為X9、有較高負(fù)載荷量的為X2,即分車帶寬度對(duì)植物生長影響較大且成正相關(guān),而喬木物種數(shù)過多,也不利于綠化植物整體健康狀況。
城市道路分車帶綠地對(duì)于城市生態(tài)環(huán)境、城市景觀等具有舉足輕重的作用,因此認(rèn)為分車帶綠地主要的評(píng)價(jià)標(biāo)準(zhǔn)是生態(tài)效益及景觀效果。分車帶綠化植物的種類構(gòu)成、數(shù)量結(jié)構(gòu)、空間結(jié)構(gòu)及年齡結(jié)構(gòu)等直接影響群落發(fā)育及生態(tài)系統(tǒng)的持續(xù)穩(wěn)定性,從而影響景觀效果。植物種類的構(gòu)成與結(jié)構(gòu)狀況可通過樣地調(diào)查獲得,在專業(yè)調(diào)查中通常有林分調(diào)查及立地調(diào)查。同一城市氣候條件相同,分車帶綠地的土壤通常改良為種植土,因此,本項(xiàng)目研究立足于認(rèn)為土壤、氣候條件及養(yǎng)護(hù)措施相同。此外,不以林木蓄積量為追求目標(biāo)是城市綠化與林業(yè)生產(chǎn)的最大區(qū)別。基于這些因素,本研究選擇及篩出的喬木物種數(shù)及株數(shù)密度、灌木物種數(shù)及株數(shù)密度、地被物種數(shù)及面積、群落層次、總郁閉度、上層郁閉度、地被蓋度、總物種數(shù)、喬灌總株數(shù)及分車帶寬度等變量因子均與群落生態(tài)效益及景觀效果密切相關(guān)。
本研究表明具有較高正載荷量的因子依次為喬木株數(shù)密度>喬木物種數(shù)>總物種數(shù)>喬灌總株數(shù)>分車帶寬度>灌木物種數(shù)>群落層次>上層郁閉度>總郁閉度。負(fù)載荷量絕對(duì)值較大的因子是群落層次>地被面積>灌木物種數(shù)>總物種數(shù)>喬木株數(shù)密度。其中,株數(shù)密度(喬木株數(shù)密度及喬灌總株數(shù))和郁閉度(上層郁閉度、總郁閉度)為群落密度指標(biāo),均表現(xiàn)了較高的載荷量。這一分析結(jié)果符合林學(xué)理論“群落密度在很大程度上決定了群落的內(nèi)部結(jié)構(gòu)”[19-20]。而喬木是綠地造景的基本元素,不僅可美化城市、凈化空氣、遮蔭降溫,更可以減輕車輛行駛造成的噪聲污染和防止眩光提高行車安全等,在景觀功能、生態(tài)效益、藝術(shù)風(fēng)格等方面能起主導(dǎo)作用。也有研究表明,大樹的遮蔭效果以及通過樹葉光合作用吸收二氧化碳呼出氧氣,疏適空氣,減輕由車輛帶來的噪聲、灰土和炫光等污染,調(diào)節(jié)氣候等生態(tài)功能是草地和灌木的數(shù)倍乃至數(shù)十倍之多[21],這說明喬木在分車帶綠化中的重要性。但在主成分分析中,喬木株數(shù)密度也體現(xiàn)了較大負(fù)載荷量,因?yàn)橹陻?shù)密度過大,個(gè)體實(shí)際占有空間低于其健康生長所需的生活空間時(shí),樹木生長必然會(huì)受到影響,并且隨著生長空間的減小,光照、通風(fēng)條件不足,影響愈明顯。樹木根系得不到伸展,樹冠發(fā)育不良,大量樹木歪斜倒伏,風(fēng)害、病蟲害發(fā)生率大增等問題。當(dāng)生長空間的減小達(dá)到致死空間時(shí),樹木生長停止直至死亡,自疏現(xiàn)象嚴(yán)重。
總物種數(shù)既表現(xiàn)出較高載荷量,也表現(xiàn)有較高負(fù)載荷量。這符合種群生長的邏輯斯蒂(Logistic)理論,即某個(gè)地區(qū)(空間)可以維持種群(物種)最大個(gè)體數(shù)目是由平均氣候條件、生境的構(gòu)造、食物供給量等決定的[22]。因此,在分車帶這種有限的綠地空間內(nèi)不能片面追求物種多樣性。
同時(shí)表現(xiàn)出較高載荷量及高負(fù)載荷量的群落層次可直接體現(xiàn)植物配置狀況。植物配置合理能充分利用空間,結(jié)構(gòu)優(yōu)化能發(fā)揮最大的生態(tài)效益,若群落層次結(jié)構(gòu)不當(dāng),不僅無景觀可言,而且群落的植物種類組成及數(shù)量也不穩(wěn)定,植物個(gè)體健康狀況差,從而無法發(fā)揮最大的生態(tài)效益。
群落郁閉度也是密度指標(biāo),體現(xiàn)群落的樹冠垂直投影狀況,郁閉度受株數(shù)密度及樹冠生長的制約,主成分分析結(jié)果表明上層郁閉度對(duì)群落的狀況影響較大。因?yàn)橛糸]度因子直接影響群落內(nèi)的光照強(qiáng)度從而影響群落狀況,從而表現(xiàn)出較高的因子載荷量。
作為立地因子的分車帶寬度,在分車帶綠地中與綠化植物的生長關(guān)系往往未引起人們的關(guān)注,主成分分析結(jié)果揭示分車帶寬度有較高正載荷量,與群落狀況密切相關(guān)。這是因?yàn)榉周噹挾炔粌H直接關(guān)系到地下營養(yǎng)空間的大小,而且當(dāng)太陽輻射從分車帶兩側(cè)照射時(shí),分車帶寬度的大小也直接影響群落內(nèi)的光照強(qiáng)度,這對(duì)群落的植物種類組成、植物生長狀況及相對(duì)穩(wěn)定都有極大的關(guān)系。
因?yàn)橹鞒煞址治鍪且韵嚓P(guān)系數(shù)矩陣為出發(fā)點(diǎn)進(jìn)行因子分析,所以主成分分析表達(dá)式中的X1~X10是經(jīng)過標(biāo)準(zhǔn)化變換后的標(biāo)準(zhǔn)變量。Y1、Y2、Y3、Y4、Y5分別代表了分車帶綠地群落中不同主導(dǎo)因子影響的5個(gè)主成分得分,通過主成分表達(dá)式可分別算出目標(biāo)樣地各主成分的得分,可為計(jì)算綜合得分提供數(shù)據(jù)支持。在評(píng)價(jià)某一分車帶綠地群落總體狀況時(shí),綜合得分Y=Y1+Y2+Y3+Y4+Y5,因此,主成分表達(dá)式可用于分車帶綠地群落的評(píng)價(jià)。
根椐深圳城市道路分車帶綠地群落狀況的主成分分析結(jié)果,株數(shù)密度(特別是喬木株數(shù)密度)、總物種數(shù)、分車帶寬度、群落層次及上層郁閉度等因子較大地制約分車帶群落的整體狀況。因此,建議在分車帶綠地群落構(gòu)建工作中可采取以下措施:①應(yīng)以分車帶寬度為基礎(chǔ),分車帶寬度越大越適宜配置喬木、中下層則適宜選用耐蔭物種,分車帶寬度偏小時(shí),則應(yīng)以灌木及地被植物為主。②在喬木株數(shù)密度控制方面,行列式配置原則上應(yīng)執(zhí)行國家相關(guān)標(biāo)準(zhǔn)規(guī)范,但采用自然式配置的分車帶綠地,目前尚無嚴(yán)格的株距標(biāo)準(zhǔn)規(guī)范,在種植設(shè)計(jì)方面多數(shù)設(shè)計(jì)師往往憑經(jīng)驗(yàn)或直覺來確定冠幅尺度大小[23],需要根據(jù)不同的徑級(jí)生長階段及郁閉度要求,參考相應(yīng)自由樹的冠徑比指標(biāo)計(jì)算出合理的株距來控制株數(shù)密度,才更利于植物健康生長并提高空間利用程度,發(fā)揮最好的生態(tài)效益[24]。③群落層次及上層郁閉度直接影響群落內(nèi)的光照強(qiáng)度及空間結(jié)構(gòu)與利用,應(yīng)避免上層郁閉度過大,以保證中下層有適度的光照強(qiáng)度。群落中不同層次樹種的選擇必須綜合考慮樹種對(duì)光照的要求、樹冠形狀及生態(tài)位等因素;此外,發(fā)生樹冠重疊時(shí),應(yīng)控制陽性樹種的樹冠重疊率[25],不同樹種在城市綠地群落中的樹冠可重疊系數(shù)有待研究。④在分車帶綠地群落中,應(yīng)遵循生態(tài)學(xué)原理合理配置植物,不可片面追求物種多樣性及無序堆砌植物。因?yàn)閱挝幻娣e內(nèi)物種數(shù)及每物種個(gè)體數(shù)量的多少與環(huán)境容量直接相關(guān),至于不同分車帶適宜的物種數(shù)則有待進(jìn)一步研究。
*:夏德美,梁國華等人參與外業(yè)調(diào)查,特此致謝!
[1]李一倫.杭州市城區(qū)主干道路綠地植物景觀調(diào)查與分析[D].杭州:浙江農(nóng)林大學(xué),2012.
[2]深圳市中國科學(xué)院仙湖植物園.深圳植物志:第1-4卷 [M].北京:中國林業(yè)出版社,2015.
[3]孫麗靜,蔣益成.嶺南不同園林植物根區(qū)土壤微生物功能多樣性季節(jié)動(dòng)態(tài)變化[J].江蘇農(nóng)業(yè)科學(xué),2016,44(10):478-487.
[4]王思?jí)?,翁沁玉,代知廣,等.基于主成分分析的武漢市中山公園水體污染研究[J].湖北大學(xué)學(xué)報(bào):自然科學(xué)版,2017,39(2):206-210.
[5]張永江,鄧茂,黃曉容,等.典型生態(tài)旅游城市黔江區(qū)大氣污染物分析及健康風(fēng)險(xiǎn)評(píng)估[J].西南師范大學(xué)學(xué)報(bào):自然科學(xué)版,2017,42(4):81-87.
[6]張亞,楊亞麗,沈濤,等.云南省食用薯蕷中無機(jī)元素含量主成分分析與聚類分析[J].湖北農(nóng)業(yè)科學(xué),2015,54(6):1400-1401.
[7]李文敏,魏虹,李昌曉,等.基于高光譜參數(shù)的楓楊葉綠素含量估算模型優(yōu)化[J].林業(yè)科學(xué),2014,50(4):56-59.
[8]李世瑤,蔡煥杰,陳新明.基于主成分分析的畦灌質(zhì)量評(píng)價(jià)[J].農(nóng)業(yè)工程學(xué)報(bào),2013,29(24),86-93.
[9]吳玉紅,田霄鴻,同延安,等.基于主成分分析的土壤肥力綜合指數(shù)評(píng)價(jià)[J].生態(tài)學(xué)雜志,2010,29(1):173-180.
[10]董瑩瑩,王紀(jì)華,李存軍,等.基于主成分分析的葉面積指數(shù)尺度效應(yīng)[J].農(nóng)業(yè)工程學(xué)報(bào),2012,28(11):164-169.
[11]公麗艷,孟憲軍,劉乃僑,等.基于主成分與聚類分析的蘋果加工品質(zhì)評(píng)價(jià)[J].農(nóng)業(yè)工程學(xué)報(bào),2014,30(13):276-283.
[12]高吉喜,段飛舟,香寶.主成分分析在農(nóng)田土壤環(huán)境評(píng)價(jià)中的應(yīng)用[J].地理研究,2006,25(5):836-842.
[13]劉大海,李寧,晁陽.統(tǒng)計(jì)分析從入門到精通[M].北京:清華大學(xué)出版社,2008:115-146.
[14]彭惠,王嵐,徐芬,等.基于主成分分析的園林植物耐蔭性指標(biāo)研究[J].草原與草坪,2016,36(3):90-95.
[15]梁曼,黃富榮,何學(xué)佳,等.熒光光譜成像技術(shù)結(jié)合聚類分析及主成分分析的藻類鑒別研究[J].光譜學(xué)與光譜分析,2014,134(8):132-136.
[16]崔雪梅,郭海如,方嬙,等.基于主成分分析的油菜鹽害生理反應(yīng)規(guī)律[J].江蘇農(nóng)業(yè)科學(xué),2016,44(2):105.
[17]宋之杰,高曉紅.一種多指標(biāo)綜合評(píng)價(jià)中確定指標(biāo)權(quán)重的方法[J].燕山大學(xué)學(xué)報(bào),2002,26(1):20-23.
[18]張子龍,王文全,繆作清,等.主成分分析在三七連作土壤質(zhì)量綜合評(píng)價(jià)中的應(yīng)用[J].生態(tài)學(xué)雜志,2013,32(6):1636-1644.
[19]汪瑛.北京市行道樹結(jié)構(gòu)分析與健康評(píng)價(jià)[D].北京:中國林業(yè)科學(xué)研究院,2011.
[20]孟憲宇.測(cè)樹學(xué)[M].北京:中國林業(yè)出版社,1995:112-121.
[21]王虎彪.城市道路綠化的幾點(diǎn)思考[J].山西建筑,2009,35(16):353.
[22]溫國勝,楊京平,陳秋夏.園林生態(tài)學(xué)[M].北京:化學(xué)工業(yè)出版社,2007.
[23]王旭東,楊秋生,張慶費(fèi).常見園林樹種樹冠尺度定量化研究[J].中國園林,2016(10):73-77.
[24]何國強(qiáng),謝良生,戴耀良,等.100種園林樹木冠徑比的初步研究[J].廣東農(nóng)業(yè)科學(xué),2017,44(5):57-64.
[25]許安芳,吳隆高,胡中成,等.杉木地理種源胸徑與冠幅相關(guān)檢驗(yàn)及其應(yīng)用[J].浙江林學(xué)院學(xué)報(bào),1998,15(2):131-137.