吳虎平，孫海，金橋，阮音音，錢佳奇，梁浩，嚴(yán)光玉，張亞玉1,※

[1.中國(guó)農(nóng)業(yè)科學(xué)院特產(chǎn)研究所，吉林 長(zhǎng)春 130112；2.成都大學(xué)藥學(xué)與生物工程學(xué)院，四川 成都 610106；3.吉林省中藥材種植（養(yǎng)殖）重點(diǎn)實(shí)驗(yàn)室，吉林 長(zhǎng)春 130112；4.瓦屋山藥業(yè)有限責(zé)任公司，四川 眉山 620365]

雅連為毛莨科植物三角葉黃連（Coptis deltoidea C.Y. Cheng et Hsiao），屬多年草本植物，以根入藥，是四川著名道地藥材[1]。雅連味苦，性寒，入心、肝、胃、大腸經(jīng)，具有清熱燥濕、瀉火解毒的功效[2]，曾被古代醫(yī)者推崇，視為黃連中的珍品[3-5]。《本草綱目》中記載“黃連今雖吳、蜀皆有，惟以雅州、眉州為良”[1]，其中所述的黃連正是如今的雅連。目前，面對(duì)野生雅連瀕臨滅絕的現(xiàn)狀，人們采用仿生栽培模式種植雅連，在實(shí)際生產(chǎn)中，選林不當(dāng)會(huì)造成雅連病害嚴(yán)重，影響雅連的存苗及產(chǎn)量，而關(guān)于林地環(huán)境條件對(duì)雅連生長(zhǎng)發(fā)育的研究尚屬于空白。依據(jù)仿生栽培建立的林藥復(fù)合體系中，已經(jīng)明確了林下參的存苗率、病蟲害發(fā)病率、根部的生長(zhǎng)發(fā)育及品質(zhì)的形成與其生長(zhǎng)的森林生態(tài)環(huán)境條件密切相關(guān)[6-8]，通過雅連實(shí)際種植地的調(diào)查研究也發(fā)現(xiàn)了類似的問題。因此，選擇適宜的林種是雅連仿生栽培技術(shù)關(guān)鍵步驟之一。

雅連生長(zhǎng)發(fā)育及其次生代謝產(chǎn)物積累與其生境條件密切相關(guān)，生境條件中主要包括土壤特性、伴生樹種、海拔高度、坡度坡向、光照強(qiáng)度和遮陰度等因素[9]，其中土壤特性是影響雅連最直接因素之一。而伴生樹種是影響土壤質(zhì)量的關(guān)鍵性因素。其原因如下：一是凋落物形成層是森林生態(tài)系統(tǒng)中植物和土壤兩個(gè)亞系統(tǒng)之間物質(zhì)循環(huán)和能量轉(zhuǎn)化的中間環(huán)節(jié)，是土壤養(yǎng)分的主要來源，同時(shí)影響植物和土壤間碳氮循環(huán)[10,11]；二是森林生態(tài)系統(tǒng)中的凋落物如枯枝落葉等的腐爛分解、樹木根系在土壤中的穿插和死亡留下的孔隙都會(huì)使得土壤變得疏松，從而改變土壤松緊度和保持水分的能力等[12]。在森林生態(tài)系統(tǒng)內(nèi)，植物從土壤中吸收的超90%的氮和磷以及超過60%的其他礦質(zhì)元素均來自于樹種等凋落物的分解[11]。不同樹種凋落物的輸入量顯著提高表層土壤的堿解氮和有機(jī)質(zhì)含量，對(duì)于提高土壤養(yǎng)分具有重要的意義[13]。

土壤肥力是保障中藥材生產(chǎn)的根本，其高低直接影響著中藥材的生長(zhǎng)發(fā)育和次生代謝產(chǎn)物積累[14]?？茖W(xué)評(píng)價(jià)土壤肥力對(duì)于指導(dǎo)農(nóng)業(yè)生產(chǎn)至關(guān)重要，由于土壤肥力指標(biāo)較多，如何對(duì)指標(biāo)進(jìn)行去冗余的同時(shí)能最大程度地反映土壤質(zhì)量特征是土壤質(zhì)量評(píng)價(jià)中首要考慮的問題。鑒于此，本文以雅連主產(chǎn)區(qū)土壤為研究對(duì)象，利用典范對(duì)應(yīng)分析（Canonical Correlation Analysis,CCA）和因子分析相結(jié)合的方式，以期全面合理評(píng)價(jià)不同樹種下雅連土壤肥力質(zhì)量，為雅連選林和仿生栽培提供技術(shù)支撐。

1 試驗(yàn)材料及方法

1.1 采樣區(qū)概況

雅連土壤樣品采自四川省眉山市洪雅縣瓦屋山鎮(zhèn)，該區(qū)域是雅連的主產(chǎn)區(qū)。多年平均氣溫13.2 ℃，1 月平均氣溫8.5 ℃，7 月平均氣溫23.1 ℃，生長(zhǎng)期多年平均200 d左右，無霜期較長(zhǎng)，平均為240 d。年平均日照時(shí)數(shù)2 100 h。年平均降水量2 500 mm，年平均降水日數(shù)230 d。土壤類型以山地黃棕壤為主，雅連主產(chǎn)區(qū)樹種信息見表1。

表1 雅連主產(chǎn)區(qū)土壤樣點(diǎn)地理信息Table 1 Geographical information of soil samples in C.deltoidea main producing area

1.2 試驗(yàn)材料與方法

1.2.1 試驗(yàn)材料 同一基地，采集不同樹種少花桂（Cinnamomum pauciflorum Nees）、珙桐（Davidia involucrate Baill）、核桃樹（Juglans）、貓兒屎（Decaisnea insignis）、楓楊（Pterocarya stenoptera）、漆樹（Rhus verniciflua Stokes）和對(duì)照組（無樹種植雅連地）的雅連及其土壤，采樣時(shí)間為2020 年9 月20 日，采集雅連時(shí)，從雅連主產(chǎn)區(qū)選擇了具有代表性的4 個(gè)采樣基地，每一基地選擇相同的樹種，樹種間距大于50 m，在每一樹種下選擇3~5 個(gè)采樣點(diǎn)，每個(gè)樣點(diǎn)隨機(jī)采3~5 株雅連。采集土壤時(shí)，首先去除枯枝落葉層，取雅連根部周圍的土壤，去除土壤中石礫碎屑，再將土壤帶回實(shí)驗(yàn)室風(fēng)干，利用0.85 mm 和0.150 mm 網(wǎng)格篩分均勻，裝入自封袋，備用。

1.2.2 測(cè)試方法 土壤pH 按照土水比1: 2.5 比例，采用梅特勒SK220 pH 計(jì)測(cè)定；土壤全氮、全碳利用元素分析Vario EL Ⅲ（德國(guó)產(chǎn)）測(cè)定；土壤有機(jī)質(zhì)含量為全碳含量1.724；土壤全磷用HNO3-HCIO4消化-鉬銻抗比色法；土壤全鉀用NaOH 熔融 火焰光度法；土壤銨態(tài)氮和硝態(tài)氮用1mol/L KCl 浸提，AA-3 流動(dòng)分析儀測(cè)定（Brawn Luyy，德國(guó)）；速效磷采用碳酸氫鈉浸提 鉬銻抗比色測(cè)定；速效鉀采用乙酸銨浸提 火焰光度法測(cè)定；中微量元素鐵、錳、銅和鋅等采用HNO3-HCIO4消解ICP710ES 法測(cè)定，土壤含水量利用烘干法測(cè)定[15]。

1.3 分析方法

1.3.1 數(shù)據(jù)分析方法 運(yùn)用Excel 2019 和SPASS 對(duì)得到的數(shù)據(jù)進(jìn)行基本統(tǒng)計(jì)分析，單因素方差分析（Oneway ANOVA）和LSD 差異顯著性檢驗(yàn)?zāi)J(rèn)設(shè)定水平為P ＜0.05。不同樹種下雅連土壤指標(biāo)對(duì)雅連單支重量和折干率的影響利用Canoco5.0 進(jìn)行典范對(duì)應(yīng)分析并排序，篩選出影響不同樹種下雅連單支重量和折干率的主要土壤質(zhì)量指標(biāo)。

1.3.2 土壤肥力質(zhì)量的綜合評(píng)價(jià)方法 經(jīng)CCA 排序后確定的土壤肥力質(zhì)量指標(biāo)，去掉不需要的指標(biāo)，在SPASS 中經(jīng)過數(shù)據(jù)標(biāo)準(zhǔn)化，繼而進(jìn)行Bartlett 球形度檢驗(yàn)。檢驗(yàn)后的數(shù)據(jù)經(jīng)因子分析得到相關(guān)因子特征值和方差貢獻(xiàn)率、因子載荷矩陣、各群落因子得分和綜合評(píng)價(jià)值及其排名，由公式（1）計(jì)算各群落類型下土壤肥力質(zhì)量評(píng)價(jià)分值[16,17]。

式中，SFQAV 為土壤肥力質(zhì)量評(píng)價(jià)分值；ai為各因子的方差貢獻(xiàn)率；zi為因子得分，zi=∑wijxij，wij為第i 個(gè)變量在第j 個(gè)因子處的因子得分系數(shù)，xij為第i 個(gè)變量在第j 個(gè)因子處的標(biāo)準(zhǔn)化值。

2 結(jié)果與分析

2.1 不同樹種下雅連土壤理化性質(zhì)比較

結(jié)合全國(guó)第2 次土壤普查土壤養(yǎng)分分級(jí)標(biāo)準(zhǔn)[18]，對(duì)洪雅縣雅連主產(chǎn)區(qū)不同樹種下土壤理化性質(zhì)統(tǒng)計(jì)分析表明（表2）：6 種不同樹種和對(duì)照組的雅連土壤理化性質(zhì)均存在不同程度的差異。在珙桐、少花桂、貓兒屎和楓楊下的土壤含水率較高，分別為56%、57%、53%和55%，顯著高于核桃樹、漆樹和對(duì)照土壤含水率，其中對(duì)照組土壤含水量最低為37%。表明林地土壤保持水分的能力明顯得到提高，但是不同類型的樹種保持水分的能力不相同。

表2 不同樹種下雅連土壤pH 和理化性質(zhì)Table 2 Soil pH and physical and chemical properties of C.deltoidea under different tree species

不同樹種下土壤pH不同，除漆樹下土壤pH（4.43）小于對(duì)照（4.88）外，其余樹種下土壤pH 總體大于對(duì)照，其中貓兒屎、漆樹和對(duì)照組土壤pH 小于5，酸性較強(qiáng)，珙桐下雅連土壤pH（6.22）最高，除珙桐外，其余6 種不同樹種下土壤的平均pH范圍為4.7~5.35，不同樹種下土壤pH 存在顯著性差異。有機(jī)質(zhì)含量在貓兒屎和珙桐樹下的土壤中含量最高，在少花桂下的土壤中含量低，且其含量遠(yuǎn)遠(yuǎn)高于國(guó)家一級(jí)水平標(biāo)準(zhǔn)（40 g/kg）。

不同樹種下土壤有機(jī)質(zhì)以貓兒屎最高（196.37g/kg），少花桂最低（79.361 g/kg），二者呈顯著差異水平（P ＜0.05），其他樹種下土壤有機(jī)質(zhì)含量差異不顯著。不同樹種下土壤全氮含量差異不顯著，但是全磷和全鉀存在顯著性差異且均為少花桂下土壤中含量最高，其中土壤全磷在漆樹下含量最低（3.3 g/kg），且高于國(guó)家一級(jí)標(biāo)準(zhǔn)（1g/kg）；土壤全鉀在核桃樹下含量最低（7.2g/kg），處于國(guó)家五級(jí)水平，在少花桂下含量最高（13.56 g/kg），處于國(guó)家四級(jí)水平，不同樹種下土壤全鉀含量總體均低；速效磷在少花桂下含量最高（312.5 mg/kg），在天竺桂下含量最低（81.15 mg/kg），其含量高于國(guó)家一級(jí)水平（40 mg/kg），速效鉀在珙桐下含量最高（780 g/kg），在漆樹下含量最低（247.04 mg/kg），其含量高于國(guó)家一級(jí)水平（200 mg/kg）。

通過對(duì)不同樹種下雅連土壤中微量元素的分析，見圖1，除鐵元素外，土壤中的中微量元素鈉、鎂、鋁、鈣等均存在差異性，其中，鋁、鈣和鋅元素在珙桐和少花桂下總體含量最高，在楓楊和漆樹下土壤中總體含量最低；鈉和銅元素在貓兒屎下土壤中含量最高，在少花桂和楓楊下土壤中含量最低。結(jié)果表明，少花桂和珙桐樹下土壤中的各種礦質(zhì)元素含量較高，而漆樹下土壤中的各種礦質(zhì)元素含量較低。

圖1 不同樹種下雅連土壤的礦質(zhì)元素含量Fig.1 Contents of mineral elements in C.deltoidea soil under different tree species

2.2 不同樹種與其下雅連土壤相關(guān)性分析

不同樹種下雅連的生長(zhǎng)發(fā)育受各種環(huán)境因子的綜合影響，CCA 排序的結(jié)果表明，不同的樹種類型和土壤理化性質(zhì)有極顯著的相關(guān)性（Monte Carlo Test，P =0.002 ＜0.01），圖2 中第一排序軸和第二排序軸中的貢獻(xiàn)率分別為51.52%和20.69%，二者累計(jì)貢獻(xiàn)率達(dá)到72.21%，說明CCA 結(jié)果可用于篩選土壤因子。在影響不同樹種下雅連的土壤理化性質(zhì)中（圖2），AK（速效鉀）、pH 和AP（速效磷）的箭頭較長(zhǎng)且與所在坐標(biāo)軸的夾角較小，說明這3 個(gè)因子對(duì)不同樹種下土壤理化性質(zhì)影響大且相關(guān)性高。由CCA 分析可知AK、pH和AP 在不同樹種下土壤中存在顯著性，其中珙桐和少花桂下土壤中pH 高于其他樹種和對(duì)照組，且珙桐下土壤中pH 與其他樹種下土壤pH 存在顯著性（P ＜0.01），少花桂下土壤中AP 高于其他樹種下土壤AP且在不同樹種下土壤中存在顯著性差異（P ＜0.05），珙桐下土壤中AK 高于其他樹種下土壤AK 且在不同樹種下土壤中存在極顯著性差異（P＜0.001）。結(jié)合圖2 結(jié)果可知這3 個(gè)土壤因子在一定程度可以作為不同樹種下土壤的理化性質(zhì)的關(guān)鍵因子。

圖2 不同樹種與土壤因子的典范對(duì)應(yīng)分析Fig.2 Canonical correspondence analysis between different tree species and soil factors

Cu、Mg、Na 和OM（土壤有機(jī)質(zhì)）等4 個(gè)因子的箭頭連線較短，與所在坐標(biāo)軸的夾角較大，且與樹種無顯著的相關(guān)性（monte carlo test，P ＜0.05），結(jié)果表明，這4 種因子對(duì)不同樹種下土壤理化性質(zhì)的影響不明顯。為了使所得結(jié)果更加準(zhǔn)確，選擇CCA 和箱型圖相結(jié)合的分析方式，進(jìn)一步對(duì)篩選所得的14 個(gè)土壤理化指標(biāo)進(jìn)行土壤質(zhì)量綜合評(píng)價(jià)。

2.3 不同樹種下雅連土壤綜合評(píng)價(jià)

因子分析是從多個(gè)變量里面篩選部分綜合變量的一種降維多元統(tǒng)計(jì)方法，用以達(dá)到數(shù)據(jù)簡(jiǎn)化的目的[14]。本研究對(duì)除去Cu、Mg、OM 和Na 等因子后的14 個(gè)土壤理化性質(zhì)進(jìn)行的Bartlett 球形檢驗(yàn)的顯著系數(shù)為P=0.00 ＜0.05，表明數(shù)據(jù)適合做因子分析。經(jīng)因子分析得到因子特征值和方差貢獻(xiàn)率、因子載荷矩陣、不同樹種下土壤的因子得分及綜合評(píng)價(jià)值，通過各個(gè)土壤因子量化的方法得到最后的綜合評(píng)價(jià)，用以篩選出不同樹種下土壤綜合質(zhì)量排名。

按照特征值＞1 的原則，抽取了4 個(gè)公因子（表3），其特征值分別為5.12、2.95、2.17 和1.40，累計(jì)方差貢獻(xiàn)率達(dá)83.13%，表明因子分析用于評(píng)價(jià)土壤質(zhì)量是可靠的。由因子載荷矩陣（表4）及因子特征值和方差貢獻(xiàn)率（表3）可知，因子1 包含了原始信息的36.57%，主要支配TK、TP、Ca、Fe、Mn 和Zn 等指標(biāo)（因子載荷值的絕對(duì)值＞0.65）；因子2包含了原始信息的21.09%，主要支配指標(biāo)為AP、NH4+-N 和NO3--N；因子3 包含了原始信息的15.47%，主要支配pH、AK；因子4 包含了原始信息的10.00%，主要支配指標(biāo)為Al。

表3 因子特征值和方差貢獻(xiàn)率Table 3 Factor eigenvalue and variance contribution rate

表4 因子載荷矩陣Table 4 Factor load matrix

經(jīng)因子分析進(jìn)行處理的14 個(gè)土壤理化指標(biāo)，其得分系數(shù)矩陣和標(biāo)準(zhǔn)化數(shù)據(jù)值加權(quán)，4 個(gè)公因子可獲得在不同樹種下土壤水平中的得分，所以由各因子得分和方差貢獻(xiàn)率加權(quán)得到不同樹種下土壤質(zhì)量綜合評(píng)價(jià)分為（表5）：少花桂（0.337）＞珙桐（0.240）＞對(duì)照（0.037）＞核桃樹（-0.067）＞貓兒屎（-0.070）＞楓楊（-0.096）＞漆樹（-0.381）。這一結(jié)果表明相比對(duì)照組，少花桂和珙桐下土壤綜合質(zhì)量高且與對(duì)照差異明顯，核桃樹、貓兒屎、楓楊和漆樹下土壤質(zhì)量較低，其中核桃樹和貓兒屎綜合評(píng)分幾乎相同，漆樹下土壤質(zhì)量評(píng)分最低且與其他樹種下土壤評(píng)分差異明顯。

3 討論

土壤中重要的化學(xué)因子pH 是影響仿生栽培類植物選地的關(guān)鍵性指標(biāo)之一，其對(duì)有機(jī)質(zhì)的積累和分解以及營(yíng)養(yǎng)元素的釋放、轉(zhuǎn)化、遷移及其有效性都有直接的影響[19,20]。本研究發(fā)現(xiàn)，珙桐、少花桂、貓兒屎、核桃樹和楓楊下土壤pH 高于對(duì)照組，漆樹下土壤pH 低于對(duì)照組，這與雷澤勇等[21]研究相類似。分析原因可能是樹種類型、樹齡和所取土壤的深度（本研究取土為雅連根部土壤，范圍為15~30cm）和含水率等會(huì)影響凋落物在土壤中的分解、分解后的養(yǎng)分組成類型和比例以及土壤微生物群落等，從而影響土壤酸堿度。另外，已有研究表明適宜雅連生長(zhǎng)的土壤pH為5.5~6.5[22]，但隨著雅連種植年限的增加，土壤趨向于酸化[23]，本研究也發(fā)現(xiàn)珙桐和少花桂下的土壤pH 接近最適pH的范圍，但漆樹下土壤pH 較低，可能會(huì)出現(xiàn)酸化現(xiàn)象。由于每種藥材的生長(zhǎng)都有其適宜的pH，土壤過酸或過堿都會(huì)抑制藥材的生長(zhǎng)[24]，可見選擇適宜的樹種后土壤pH 改變，這可能有助于提高雅連生長(zhǎng)過程中抵抗環(huán)境脅迫的能力，提升雅連品質(zhì)和存苗率等。

不同樹種對(duì)雅連生長(zhǎng)土壤中的有機(jī)質(zhì)、鉀元素和磷元素均有影響。在森林生態(tài)系統(tǒng)中，凋落物是土壤養(yǎng)分主要來源，樹種和落葉量更是直接影響土壤養(yǎng)分含量[25]，但是不同樹種對(duì)土壤養(yǎng)分產(chǎn)生不同的影響，如張璐等[26]研究發(fā)現(xiàn)樟樹和唐竹有利于提高土壤堿解氮，但不利于提高磷含量。本研究中珙桐樹下的土壤速效鉀和有機(jī)質(zhì)的含量高于其他處理，可能是珙桐在該地區(qū)的生長(zhǎng)年限較長(zhǎng)，產(chǎn)生了較厚的枯落層，經(jīng)過微生物分解所產(chǎn)生的腐殖質(zhì)較多，因此珙桐樹下土壤中的有機(jī)質(zhì)含量也較多。同時(shí)，作為植物生長(zhǎng)所需的主要肥力之一，鉀元素也是影響雅連生長(zhǎng)的主要因素，雖然速效鉀是能夠被植物直接吸收的類型，但容易發(fā)生淋溶與洗脫，因此在山地環(huán)境中，速效鉀會(huì)因降雨而發(fā)生沿坡向流失或積累，這將導(dǎo)致其碎片化隨機(jī)分布[27]。本文采樣地坡度范圍為55~60°，坡度較陡，在6 種不同樹種中珙桐的樹齡較高，枝葉覆蓋面積較廣，能夠有效地緩解由降雨帶來的速效鉀流失的現(xiàn)象，而漆樹樹齡較低且分布不廣，則使速效鉀含量較低，所以這可能是造成珙桐和漆樹下土壤速效鉀差異較大的原因之一。另外，本研究還發(fā)現(xiàn)少花桂下雅連土壤中速效磷和全磷的含量高于其他處理，由于磷是一種長(zhǎng)期沉積性礦物質(zhì)元素，主要來源于巖石的風(fēng)化和凋落物的礦化[28-30]，此過程相對(duì)漫長(zhǎng)，而該地區(qū)存在較多人工種植的少花桂，其凋落物層較厚，所以這些凋落物可能通過影響巖石分化，進(jìn)一步影響磷的含量，最終使少花桂下土壤的磷含量總體較高，這一結(jié)果也與許彥崟[31]的研究相符。總體而言，本研究發(fā)現(xiàn)少花桂和珙桐下土壤養(yǎng)分供應(yīng)能力總體大于其他樹種和對(duì)照組，而氮、磷、鉀是植物生長(zhǎng)發(fā)育的三要素，全量和速效養(yǎng)分含量分別代表土壤養(yǎng)分潛在供應(yīng)能力和實(shí)時(shí)供應(yīng)能力[32]，說明這兩種樹種有利于雅連的生長(zhǎng)發(fā)育。

某些無機(jī)元素對(duì)中藥材的生長(zhǎng)、藥效物質(zhì)的積累及功效的發(fā)揮具有重要的作用，如鉀、硼、鉬、銅、鋅等元素是植物體內(nèi)酶或輔酶的重要組成部分，在植物生長(zhǎng)發(fā)育過程中起著不可缺少的作用[33]，周利和郭蘭萍等[34]研究表明不同產(chǎn)地的黃連樣品中無機(jī)元素的含量差異較大，其中K、Mg、Mn、Zn、Fe 含量最高，還有研究認(rèn)為土壤中的中微量元素含量會(huì)在一定程度上影響植物藥材的微量元素的含量，并對(duì)植物的不同部位產(chǎn)生不同的影響[35,36]。在本研究中，土壤中的Fe、K、Mn、Zn等在不同樹種和對(duì)照存在差異（圖1 和圖2），在少花桂和珙桐樹下的總體含量高，可能的原因是珙桐樹和少花桂的凋落層多，影響土壤中的中微量元素的積累。另外，這些元素的積累可能影響雅連有效成分和品質(zhì)，所以在后續(xù)的研究中，我們將針對(duì)土壤中的中微量元素、雅連中的中微量元素及其有效成分和獨(dú)特的功效展開研究。

4 結(jié)論

本文通過典型對(duì)應(yīng)分析和因子分析，對(duì)四川省眉山市洪雅縣瓦屋山鎮(zhèn)6 種不同樹種下的土壤質(zhì)量進(jìn)行了綜合評(píng)價(jià)，發(fā)現(xiàn)土壤pH、速效鉀和速效磷是影響最大的土壤因子；土壤無機(jī)元素（Ca、K、Mn 和Zn 等）在不同樹種均存在差異性，且在少花桂和珙桐樹下土壤中含量最高。本研究表明少花桂和珙桐樹下雅連土壤質(zhì)量綜合評(píng)價(jià)最佳，適宜于雅連的仿生栽培；漆樹和楓楊下綜合評(píng)價(jià)最低，不宜于雅連的仿生栽培。綜上所述，土壤pH、速效磷、有效鉀及鈣、錳、鋅可以作為雅連產(chǎn)區(qū)土壤質(zhì)量評(píng)價(jià)監(jiān)視因子，并建議優(yōu)先選擇少花桂和珙桐作為選林的樹種。