楊 濤 邱勇斌 沈 漢 鄭成忠 張 振 王文月 金國(guó)慶 周志春

（1. 中國(guó)林業(yè)科學(xué)研究院亞熱帶林業(yè)研究所 浙江省林木育種重點(diǎn)實(shí)驗(yàn)室 杭州 311400；2. 華中農(nóng)業(yè)大學(xué) 武漢 430070；3. 浙江省開化縣林場(chǎng) 開化 324300）

全球森林儲(chǔ)存著近3千億t碳，林木是森林生態(tài)系統(tǒng)中最重要的碳匯，在維持全球碳平衡、減緩CO2濃度上升等方面起著重要作用（郭金堂等，2016）。森林的適應(yīng)性取決于林木基因型對(duì)未來(lái)氣候變化的反應(yīng)，對(duì)林木生產(chǎn)力的形成和固碳能力具有重要影響，國(guó)內(nèi)外學(xué)者研究輻射松（Pinus radiata）、火炬松（Pinus taeda）、落葉松（Larix gmelinii）、冷杉（Abies fabri）和白樺（Betula platyphylla）等樹種的生物量累積與分配、含碳率及含碳量表明，固碳能力在不同樹種、同一樹種不同品種、樹齡或器官間、不同立地間均存在顯著差異（Eliaset al., 2003；Zhanget al., 2009；賈慶彬等，2014；Wanget al., 2015；Donget al., 2015）。同一樹種含碳量差異主要是因種源、家系、無(wú)性系和個(gè)體間生長(zhǎng)速度、生物量與含碳率的差異引起，選育林木良種造林可以顯著提高碳匯價(jià)值和經(jīng)濟(jì)效益。目前，較為精確的含碳量測(cè)量方法是在測(cè)定各器官生物量和含碳率的基礎(chǔ)上推算出來(lái)的，由于該方法耗時(shí)耗力，且對(duì)林地具破壞性（林瑋等，2020），難以應(yīng)用到大群體樣本，而利用異速生長(zhǎng)方程估算含碳量是常用的替代方法（Widagdoet al., 2020；張悅等，2022）。

家系和無(wú)性系是當(dāng)前主要造林樹種的良種推廣方式。無(wú)性系綜合利用加性與非加性遺傳效應(yīng)，具有較高的遺傳增益，而全同胞和半同胞家系僅分別包含1/2和1/4的加性方差（Wenget al.，2008；康向陽(yáng)，2019）；同時(shí)，通過(guò)無(wú)性系育種可獲得目標(biāo)性狀表現(xiàn)突出的林木新品種，是利用群體中遺傳增益最為顯著的個(gè)體而不是群體平均值，具有家系林業(yè)無(wú)法比擬的優(yōu)勢(shì)（Wenget al., 2008）；并且，無(wú)性系林木的生長(zhǎng)性狀表現(xiàn)穩(wěn)定、一致，便于集約化經(jīng)營(yíng)，可大幅提高林業(yè)管理效率（Wu, 2019）。但從育種研究的角度，同一樹種無(wú)性系與家系的固碳能力是否存在差異還有待揭示。

柏木（Cupressus funebris）是栽培歷史悠久的生態(tài)和用材樹種，主要分布在我國(guó)長(zhǎng)江流域及以南地區(qū)，其材質(zhì)優(yōu)良，耐瘠薄，喜鈣質(zhì)土，適生性和涵養(yǎng)水源能力強(qiáng)，是石漠化地區(qū)和石質(zhì)山地生態(tài)修復(fù)的先鋒樹種，近年來(lái)在我國(guó)南方廣泛推廣栽植，為速生優(yōu)質(zhì)高固碳品種選育提供了巨大空間。本研究項(xiàng)目組自20世紀(jì)90年代開始收集全國(guó)主產(chǎn)區(qū)的柏木優(yōu)樹資源，系統(tǒng)評(píng)價(jià)育種群體的遺傳多樣性，選育出一批適生性強(qiáng)的優(yōu)良種質(zhì)，從種子園豐產(chǎn)、育苗技術(shù)以及“采穗圃營(yíng)建-截頂促萌-幼化-扦插”等方面實(shí)現(xiàn)了柏木苗木的規(guī)模化繁育（樓君等，2014；Yanget al., 2016；鄭一等，2020）。本研究以營(yíng)建在浙江省開化縣林場(chǎng)9年生柏木無(wú)性系和家系試驗(yàn)測(cè)定林為對(duì)象，探究其生長(zhǎng)、生物量、含碳量及其分配情況，篩選并驗(yàn)證通用異速生長(zhǎng)方程的準(zhǔn)確性，嘗試評(píng)價(jià)不同基因型的固碳潛力，以期為我國(guó)南方高固碳柏木優(yōu)良品系選擇及其推廣應(yīng)用奠定理論基礎(chǔ)。

1 材料與方法

1.1 試驗(yàn)材料與栽植

柏木無(wú)性系和家系試驗(yàn)測(cè)定林來(lái)自浙江省開化縣林場(chǎng)，地理坐標(biāo)為118°24′24″E，29°8′17″N。屬亞熱帶季風(fēng)性氣候，年均氣溫16.4 ℃，晝夜溫差平均10.5 ℃，年均降雨量1 814 mm，年均日照時(shí)數(shù)1 679.8 h，無(wú)霜期252天。試驗(yàn)點(diǎn)為柏木在該地區(qū)造林的代表性立地，其中，家系測(cè)定林材料共34個(gè)系號(hào)，包括33個(gè)家系（1代無(wú)性系種子園中分單系采集的種子）和1個(gè)對(duì)照（CK1，1代無(wú)性系種子園中的混合種子），于2012年3月在苗圃內(nèi)播種育苗；無(wú)性系測(cè)定林材料共36個(gè)系號(hào)，包括35個(gè)無(wú)性系（1代無(wú)性系種子園中采集種子育苗后選擇的優(yōu)良單株）和1個(gè)對(duì)照（CK2，1代無(wú)性系種子園中的混合種子），于2012年5月選取10 cm長(zhǎng)半木質(zhì)化插條在苗床上扦插育苗。2013年?duì)I建試驗(yàn)測(cè)定林，采用隨機(jī)完全區(qū)組設(shè)計(jì)，設(shè)置5個(gè)試驗(yàn)小區(qū)，10株小區(qū)雙行排列，株行距2 m×2 m。

1.2 生長(zhǎng)調(diào)查與生物量取樣

1.2.1 生長(zhǎng)調(diào)查與材積計(jì)算 9年生時(shí)，全林調(diào)查樹高（tree height,H, m）、胸徑（diameter at breast height,DBH, cm）、枝下高（ height under branches, UBH, cm）和冠幅（crown width CW, m），單株材積（volume,V, m3）利用公式V=0.000 057 173 591×DBH1.31.8813305×H0.9956884計(jì)算（楊育林等，2014）。

1.2.2 取樣與生物量測(cè)定 選取樹高、胸徑、材積、枝下高和冠幅作為評(píng)價(jià)指標(biāo)，設(shè)置各性狀權(quán)重依次為1.5、1.5、1.5、1.0和1.0，采用TOPSIS綜合評(píng)價(jià)法（Meshramet al., 2020）分別對(duì)無(wú)性系和家系試驗(yàn)林各系號(hào)的生長(zhǎng)性狀優(yōu)劣進(jìn)行排序，綜合篩選出優(yōu)、中、差3個(gè)等級(jí)，在3個(gè)等級(jí)類型中各選擇3個(gè)無(wú)性系（家系），每個(gè)無(wú)性系（家系）有3個(gè)區(qū)組，在每個(gè)區(qū)組中選擇3個(gè)單株，無(wú)性系和家系分別伐取樣木81株。

生物量測(cè)定分為地上和地下2部分。樹木伐倒后，樹干部分每1 m切斷并稱重記錄，之后每段裁取5 cm圓盤，分別測(cè)量每個(gè)圓盤（帶皮）的鮮質(zhì)量，標(biāo)記裝袋；樹冠分上、中、下3部分，分別測(cè)量每部分枝和葉的總質(zhì)量，每部分選取3條有代表性的枝條，分離枝葉并稱重、裝袋、標(biāo)記；地下生物量測(cè)定中，按大小分為粗根（直徑＞2 cm，包含根樁）和細(xì)根（≤2 cm），稱重并取樣（Zhuet al.，2013）。將所有樣品帶回實(shí)驗(yàn)室，于80 ℃烘箱烘至恒重，測(cè)量各部位樣品干質(zhì)量，計(jì)算含水率（%），然后根據(jù)各器官鮮質(zhì)量乘以含水率分別得出樹干干生物量（trunk biomass, TB, kg）、樹枝干生物量（branch biomass, BB, kg）、樹葉干生物量（leaf biomass, LB, kg）和根干生物量（root biomass, RB, kg），最后得出單株總干生物量（total biomass of a single tree,TBS, kg）：TBS=TB+BB+LB+TRB。

1.2.3 含碳率與含碳量測(cè)定 利用Elementar Vario EL（Germany）元素分析儀測(cè)定各器官含碳率（C, %），采用木材總碳（total carbon, TC）含量測(cè)定方法，溫度設(shè)置為1 100 ℃，進(jìn)樣量為20 mg（納曉瑩，2016）。根據(jù)各器官含碳率和干生物量分別得出樹干含碳量（trunk carbon stock, TCS, kg）、樹枝含碳量（branch carbon stock, BCS, kg）、樹葉含碳量（leaf carbon stock, LCS, kg）以及樹根含碳量（root carbon stock, RCS, kg），則單株總含碳量（total carbon stock of a single tree, TCSS, kg）計(jì)算公式為：TCSS=TCS+BCS+LCS+RCS。

1.3 數(shù)據(jù)處理與分析

1.3.1 方差分析和遺傳參數(shù) 采用以下線性混合模型對(duì)各性狀進(jìn)行生長(zhǎng)差異分析：

式中：μ為總體平均值；Fi為無(wú)性系或家系間效應(yīng)；Bj為區(qū)組間效應(yīng)；eijk為隨機(jī)誤差。

式中：δp為標(biāo)準(zhǔn)差；為性狀平均值。

1.3.2 異速生長(zhǎng)方程構(gòu)建 異速生長(zhǎng)方程采用冪函數(shù)（Y=aXb）的對(duì)數(shù)變換形式（lnY=a+blnX），由于模型存在異方差，使用對(duì)數(shù)形式模型是消除異方差現(xiàn)象的方法之一（張悅等，2022）。在多元回歸模型中，自變量間可能高度相關(guān)，導(dǎo)致回歸系數(shù)不穩(wěn)定，使用對(duì)數(shù)形式模型進(jìn)行線性回歸，可通過(guò)方差膨脹系數(shù)（variance inflation factor, VIF）衡量多重共線性嚴(yán)重程度，從而對(duì)自變量做出取舍（Linet al., 2017; Moussaet al., 2018）。本研究選用以下5個(gè)方程:

式中：TCSS為各器官的含碳量；k0、k1、k2、k3為擬合系數(shù)；D為胸徑；H為樹高；UBH為枝下高。式（2）、（3）和（5）通過(guò)VIF衡量多重共線性程度（Dutcet al.，2018）；采用決定系數(shù)（R2）和均方根誤差（root mean square error, RMSE）進(jìn)行模型效果評(píng)價(jià)和選優(yōu)。

統(tǒng)計(jì)分析采用SPSS 26.0軟件進(jìn)行，TOPSIS綜合評(píng)價(jià)應(yīng)用DPS 14.0軟件。

2 結(jié)果與分析

2.1 生長(zhǎng)性狀差異分析和遺傳參數(shù)估算

無(wú)性系的平均樹高、胸徑和單株材積分別為6.09 m、6.15 cm和0.010 6 m3，分別高出家系水平47.81%、28.39%和120.83%，家系水平的表型變異系數(shù)高于無(wú)性系，說(shuō)明家系具有較高的遺傳選擇響應(yīng)。樹高、胸徑、單株材積、枝下高和冠幅的無(wú)性系重復(fù)力幅度為0.651～0.756，家系遺傳力為0.132～0.322（表1）。

表1 柏木無(wú)性系和家系生長(zhǎng)性狀的差異分析和遺傳參數(shù)估算①Tab. 1 Difference analysis of growth traits and estimation of genetic parameters of cypress clones and families

方差分析結(jié)果表明，柏木無(wú)性系間與家系間各生長(zhǎng)性狀均存在極顯著差異（表1）。根據(jù)TOPSIS綜合評(píng)價(jià)結(jié)果（表2），無(wú)性系中排名1～12為優(yōu)等級(jí)、排名13～24為中等級(jí)、排名25～36為差等級(jí)，其中優(yōu)等級(jí)的Ci均值是中等級(jí)的159.49%、是差等級(jí)的298.15%；家系中排名1～11為優(yōu)等級(jí)、排名12～22為中等級(jí)、排名23～34為差等級(jí)，其中優(yōu)等級(jí)的Ci均值是中等級(jí)的153.46%、是差等級(jí)的317.89%。結(jié)合保存狀況，進(jìn)一步選擇出優(yōu)等無(wú)性系為C6、C9、C14，中等無(wú)性系為C1、C5、C11，差等無(wú)性系為C10、C20、C34，優(yōu)等無(wú)性系的胸徑均值和樹高均值分別是中等無(wú)性系的136.90%和113.84%、是差等無(wú)性系的182.89%和131.43%；選擇優(yōu)等家系為F8、F18和F24，中等家系為F5、F10和F22，差等家系為F11、F26和F32，其中，優(yōu)等家系的胸徑均值和樹高均值分別是中等家系的123.73%和116.32%、是差等家系的143.14%和126.99%。

表2 TOPSIS法無(wú)性系和家系各系號(hào)優(yōu)劣順序評(píng)價(jià)結(jié)果①Tab. 2 Obtaining the superior and inferior sequence of clone and family numbers by TOPSIS

2.2 生長(zhǎng)性狀與生物量、含碳量、含碳率的關(guān)系

由圖1、圖2可知，無(wú)論家系還是無(wú)性系水平，柏木胸徑和樹高均與單株總生物量存在顯著線性關(guān)系（P＜0.01），與單株總含碳量也存在顯著線性關(guān)系（P＜0.01），可見(jiàn)，生長(zhǎng)較好的基因型其生物量和含碳量也較高。胸徑和樹高與含碳率均無(wú)顯著線性關(guān)系，說(shuō)明含碳率受生長(zhǎng)狀況影響小。

圖1 柏木無(wú)性系樹高和胸徑與生物量、含碳量、含碳率的線性關(guān)系Fig. 1 Linear relationship between height and DBH of cypress clones and biomass, carbon storage and carbon content

圖2 柏木家系樹高和胸徑與生物量、含碳量、含碳率的線性關(guān)系Fig. 2 Linear relationship between height and DBH of cypress families and biomass, carbon storage and carbon content

2.3 生物量和含碳量的分配模式

由表3可知，柏木無(wú)性系地上部分生物量分配相比家系高2.12%，其中，地上部分樹干生物量分配比例最高（44.36%），但較家系低0.60%，而無(wú)性系枝生物量分配比例為18.45%，較家系高5.55%，葉生物量分配比例為23.22%，較家系低2.84%；無(wú)性系地下部分生物量分配比例相比家系低2.12%，其中，粗根生物量分配比例低出家系3.45%，而細(xì)根生物量分配比例高出家系1.34%，說(shuō)明無(wú)性系主要是以犧牲地下粗根的生物量分配、增加地上部分樹枝生物量和細(xì)根生物量來(lái)維持生長(zhǎng)優(yōu)勢(shì)。

表3 柏木各器官生物量、含碳率、含碳量分配特征Tab. 3 Biomass, carbon content rate and carbon storage distribution of cypress organs

無(wú)論無(wú)性系還是家系，柏木各器官含碳率排序均為葉＞枝＞粗根＞樹干＞細(xì)根。無(wú)性系各器官含碳量占比排序?yàn)闃涓桑緲淙~＞樹枝＞粗根＞細(xì)根，家系各器官碳含量占比排序均為樹干＞樹葉＞粗根＞樹枝＞細(xì)根，說(shuō)明樹干和樹葉是柏木含碳量的主要部分。整體上，無(wú)性系樹干、樹枝和細(xì)根的含碳量分配比例高于家系，而樹葉和粗根的含碳量分配比例低于家系。

2.4 含碳量異速生長(zhǎng)方程的建立

構(gòu)建單株總含碳量及其各器官含碳量與生長(zhǎng)性狀之間的異速生長(zhǎng)方程發(fā)現(xiàn)，以胸徑為單一變量擬合效果最差（式1），當(dāng)以樹高為第二自變量加入模型（式2），模型精度有所提高，但幅度有限（表4）。當(dāng)以樹高、胸徑、枝下高為獨(dú)立變量時(shí)，估算的含碳量（除無(wú)性系樹根含碳量外）精度較高，即模型（3）的R2達(dá)到最大，RMSE達(dá)到最小，且方差膨脹系數(shù)（VIF）較小，在2.977～3.982之間，不存在共線性問(wèn)題。而當(dāng)樹高與胸徑以D2H形式整合成一個(gè)參數(shù)作為自變量時(shí)（式4），與模型（3）相比，無(wú)性系和家系的單株總生物量、地上和地下部分生物量擬合度均有所下降，將枝下高作為獨(dú)立變量加入模型（4）后，R2均有提高，說(shuō)明枝下高作為獨(dú)立變量加入模型對(duì)含碳量模型的擬合度影響較大。此外，利用模型（2）和模型（3）估算的無(wú)性系樹根含碳量精度相近，R2均為0.880，模型（2）的RMSE和VIF較模型（3）分別下降3.06%和33.31%。綜合R2、RMSE以及VIF表現(xiàn)，篩選模型（3）為柏木生物量的最優(yōu)擬合模型。

表4 柏木各器官含碳量回歸模型①Tab. 4 Regression model of carbon storage of cypress organs

2.5 優(yōu)良品系選擇

利用篩選出的含碳量最佳擬合模型，計(jì)算無(wú)性系和家系試驗(yàn)林的單株總含碳量。由表5可知，含碳量的無(wú)性系重復(fù)力為0.647，家系遺傳力為0.161，按照20%入選率，無(wú)性系試驗(yàn)林材積和含碳量的遺傳增益分別為39.74%和51.48%，家系試驗(yàn)林分別為10.79%和16.14%。以單株材積和含碳量為指標(biāo)，采用獨(dú)立淘汰法篩選出生長(zhǎng)與固碳兼優(yōu)無(wú)性系6個(gè)，分別為C6、C9、C12、C16、C13和C2，其單株材積平均值高出CK1的191.49%，單株含碳量平均值高出CK1的107.11%。篩選出生長(zhǎng)與固碳兼優(yōu)家系6個(gè)，分別為F18、F9、F24、F8、F25和F13, 其單株材積平均值高出CK2的32.45%，單株含碳量平均值高出CK2的33.39%（圖3）。

圖3 獨(dú)立淘汰法篩選生長(zhǎng)和高固碳兼優(yōu)的無(wú)性系與家系Fig. 3 Screening clones and families with high carbon fixation and growth by independent elimination method

表5 20%入選率估算無(wú)性系與家系試驗(yàn)林的單株材積和含碳量的遺傳增益①Tab. 5 Estimation of genetic gain of individual volume and carbon storage in clonal and family experimental forests with 20% selection rate

3 討論

已有研究表明，通過(guò)測(cè)定后選擇優(yōu)良無(wú)性系造林，針葉樹可獲得額外5%～25%的遺傳增益，達(dá)到同世代家系試驗(yàn)林所獲遺傳增益的2倍（Wu, 2019）。柏木無(wú)性系生長(zhǎng)優(yōu)勢(shì)明顯，樹高、胸徑和材積分別高出家系水平47.81%、28.39%和120.83%，按照20%入選率，估算無(wú)性系試驗(yàn)林材積和含碳量的遺傳增益分別為39.74%和51.48%，家系試驗(yàn)林分別為10.79%和16.14%，無(wú)性系試驗(yàn)林相比于家系試驗(yàn)林的遺傳增益高。Mullin 等（1994）對(duì)10年生黑云杉（Picea mariana）無(wú)性系測(cè)定表明，通過(guò)顯性效應(yīng)和上位效應(yīng)獲得的額外遺傳增益幾乎與全同胞和半同胞選擇獲得的遺傳增益相等。Wu 等（2008）對(duì)新西蘭輻射松的研究同樣表明，選擇前5%無(wú)性系，材積遺傳增益可達(dá)24%，在相同選擇強(qiáng)度下，比家系測(cè)定林提高1倍以上。在火炬松無(wú)性系試驗(yàn)中，從每個(gè)家系選擇最佳單株繁育的無(wú)性系，前10名的遺傳增益幾乎是家系水平的2倍（Baltuniset al., 2007）。上述報(bào)道與本研究得出的結(jié)果類似，說(shuō)明柏木無(wú)性系不僅更好地利用了加性遺傳方差，而且獲得了更多非加性遺傳方差，保持了優(yōu)良性狀的遺傳值（Wenget al., 2008）；但是，額外遺傳增益差異主要由群體大小、家系內(nèi)測(cè)定的無(wú)性系數(shù)量、子代測(cè)定的準(zhǔn)確性等因素造成，因此，為了更精準(zhǔn)評(píng)價(jià)無(wú)性系和家系的增益水平，應(yīng)在多個(gè)氣候區(qū)評(píng)估家系/無(wú)性系基因型與環(huán)境的互作效應(yīng)，估計(jì)其加性和非加性遺傳方差，以預(yù)測(cè)無(wú)性系林業(yè)的穩(wěn)定性和遺傳增益（Wu, 2019）。

本研究表明，生物量、含碳量與樹高和胸徑均存在極顯著線性關(guān)系，說(shuō)明柏木生長(zhǎng)快慢對(duì)生物量和含碳量具有顯著影響（郭金堂等，2016），而含碳率與樹高和胸徑不存在顯著線性關(guān)系，與張悅等（2022）對(duì)長(zhǎng)白落葉松（Larix olgensis）各器官含碳量的測(cè)定結(jié)果類似，說(shuō)明柏木生長(zhǎng)快慢對(duì)含碳率影響較小，因此選擇生長(zhǎng)性狀優(yōu)良、生物量較大的優(yōu)良品系，其固碳能力也相應(yīng)較強(qiáng)。與家系不同，無(wú)性系的樹枝和樹葉分配比例均高于樹根，枝和葉的生物量占總生物量的41.67%，而家系的樹根分配比例高于樹枝，枝和葉的生物量占總生物量的38.96%，說(shuō)明無(wú)性系具有較大的樹冠，無(wú)性系的冠幅比家系長(zhǎng)33.33%（表1）。進(jìn)一步分析發(fā)現(xiàn)，無(wú)性系和家系的地上與地下生物量之比分別為6.16∶1和5.22∶1，樹干與枝的生物量之比分別為2.40∶1和3.49∶1，樹干與葉的生物量之比分別為1.91∶1和1.73∶1，生物量分配策略表明，無(wú)性系是以犧牲地下粗根生物量為代價(jià)增強(qiáng)細(xì)根和枝葉的生長(zhǎng)，無(wú)性系的枝具有更強(qiáng)的競(jìng)爭(zhēng)能力，充足的枝和葉是無(wú)性系表現(xiàn)出快速生長(zhǎng)的重要因素，生長(zhǎng)快速的無(wú)性系優(yōu)先分配更多的有機(jī)物給樹枝。Zhu等（2013）研究表明，在生理學(xué)上，生長(zhǎng)快速的樹木吸收更多的陽(yáng)光進(jìn)行光合作用，并在地上積累更多有機(jī)物質(zhì)。柏木喜鈣質(zhì)土，生長(zhǎng)在土壤瘠薄處，平衡生長(zhǎng)假說(shuō)認(rèn)為，為滿足快速生長(zhǎng)需求，需要更有效獲取環(huán)境中的資源，植株會(huì)將資源分配到最需要的器官中，這可能使柏木無(wú)性系的枝和細(xì)根器官分配更多的有機(jī)物（Stovallet al.,2013）?？梢灶A(yù)期，與家系相比，無(wú)性系可以實(shí)現(xiàn)更高的生產(chǎn)力和碳儲(chǔ)存能力（Wu, 2019）。

本研究對(duì)無(wú)性系和家系分別擬合含碳量異速生長(zhǎng)方程。胸徑作為最常用和最容易被精確測(cè)量的指標(biāo)，在生物量和含碳量模型中應(yīng)用最為廣泛，可獲得較高的擬合精度，但是本研究以胸徑為單一自變量時(shí)，模型R2均較低，這可能與樹齡有關(guān)，柏木試驗(yàn)林僅9年生，胸徑較小，表型變異系數(shù)高；當(dāng)以樹高和胸徑作為變量時(shí)，估算精度得到提高，有研究表明，生物量模型中加入樹高可抵消區(qū)域效應(yīng)（Feldpauschet al.,2012）；當(dāng)將枝下高作為自變量包含在模型中則能顯著提高估算精度，同時(shí)不存在多重共線性問(wèn)題，說(shuō)明柏木樹高、胸徑和枝下高與含碳量存在極強(qiáng)的內(nèi)在關(guān)聯(lián)，是含碳量模型的關(guān)鍵參數(shù)。柏木早期以高生長(zhǎng)為主，分枝多，枝下高決定樹冠高度，直接影響植株的光合作用以及生長(zhǎng)過(guò)程（李想等，2018；廖國(guó)莉等，2020），而樹冠特征也是影響枝葉生物量預(yù)測(cè)精度的重要變量（汪金松等，2011）。當(dāng)前，優(yōu)良林木基因型選擇是未來(lái)人工林生產(chǎn)和碳儲(chǔ)存能力關(guān)注的方向，從育種角度看，應(yīng)用遺傳選擇可提高柏木生長(zhǎng)潛力，表明栽植無(wú)性系可提高生產(chǎn)力。本研究按照20%入選率，評(píng)估無(wú)性系試驗(yàn)林材積和含碳量的遺傳增益分別為39.74%和51.48%，家系試驗(yàn)林分別為10.79%和16.14%，篩選出生長(zhǎng)與固碳兼優(yōu)的無(wú)性系和家系6個(gè)，為生產(chǎn)提供了重要材料。因此，在柏木早期生長(zhǎng)階段，考慮胸徑、樹高、枝下高以及分枝等性狀來(lái)選擇優(yōu)良品系，評(píng)估其固碳潛力更具有現(xiàn)實(shí)意義。

4 結(jié)論

柏木無(wú)性系生長(zhǎng)優(yōu)勢(shì)明顯，其樹高、胸徑和材積分別高出家系47.81%、28.39%和120.83%。柏木樹干生物量的分配比例最高，生物量分配模式差異是影響無(wú)性系和家系間生長(zhǎng)差異的重要因素，無(wú)性系主要是以犧牲粗根的生物量分配（較家系低3.45%）、增加樹枝生物量（較家系高5.55%）和細(xì)根生物量（較家系高1.34%）來(lái)維持生長(zhǎng)優(yōu)勢(shì)。生物量、含碳量均與生長(zhǎng)性狀呈極顯著線性關(guān)系，說(shuō)明選擇速生、生物量較大的優(yōu)良品系，可提高其固碳潛力，且無(wú)性系具有比家系更高的固碳能力。選擇將樹高、胸徑和枝下高作為特定參數(shù)，構(gòu)建生長(zhǎng)早期的含碳量模型精度高，預(yù)測(cè)效果好；進(jìn)一步結(jié)合材積性狀篩選出速生與固碳兼優(yōu)的無(wú)性系和家系各6個(gè)，將有助于提高柏木人工林的生產(chǎn)力和含碳量。