桉树Eucalyptus spp.在中国南方的种植已经超过367万hm^2[1]，因其显著的速生性和丰产性而被列为我国短周期工业原料林造林的首选树种之一，它用不到3%的林地面积解决了18%的工业用材需求，其良好的短期经济效益为其大面积的发展起到了积极的作用^[2]。然而，由于桉树引种和推广种植的地区大多为台风频发的东南沿海地区，又因杂交无性系的推出使其品种繁多，其适生性和抗风性又各不同，造成明显产量差异，因此有必要对广泛推广种植的品种作系统的多点对比，为其造林品种的选择提供科学的依据。目前关于桉树适生性和抗性方面的研究多集中于桉树的抗病性^[3-5]和抗寒性^[6-7]，虽然前人亦有相关桉树的适生性研究^[8]和抗风性对比试验^[9-13]，但多为单一地区的研究，缺乏对多个立地、多个桉树无性系的系统比较，且研究的品种和立地情况也各有不同。本文对10个目前主要推广种植的桉树无性系，在广东和广西的6个不同立地的生长情况进行对比研究，以期对其抗性和适生性作出科学评价，为“适地适树”的选择提供科学参考。

1 材料与方法 1.1 试验地概况

广东3个试验地位于雷州林业局的唐家林场（20°50′01″ N，109°51′76″ E）和遂溪林场（21°21′81″ N，110°01′66″ E），及湛江林业局属下的吴川林场（21°18′37″ N，110°31′13″ E）。唐家林场林地为玄武岩发育的厚有机质层赤红壤，土壤黏重，肥力较好；吴川林场林地为花岗岩发育的赤红壤，中壤土，肥力中等；遂溪林场林地为浅海沉积物发育的砂壤土，松散而少黏粒，肥力中下。广西3个试验地：钦廉林场（21°45′46″ N，109°03′08″E），浅海沉积物发育的中壤土，肥力中上；防城林地（21°48′02″N，108°25′16″E），沙页岩发育的赤红壤，土层较薄，中壤土，肥力中下；东门林场（22°28′N，107°51′E），石灰岩发育的赤红壤，土层深厚黏重，肥力较好，相对干旱。各试验地土壤分析结果见表 1。除广西东门林场试验地属于南亚热带湿润季风气候外，其余5个试验地均属于典型热带北缘海洋性季风气候，干湿季明显，夏季盛行东南风、南风及西南风，冬季多北风和东北风，每年6—11月为热带风暴和台风季节，这为各无性系的抗风性比较提供了条件。

1.2 供试材料

供试桉树无性系有来自广西东门林场的尾巨桉E. urophylla×E. grandis系列DH3226、DH3213、DH3327、DH3222和巨赤桉E. grandis× camaldulensis DH2012，尾细桉E.urophylla×E.tereticornis系列中有来自雷州林业局林业科学研究所的M1和来自热带林业研究所的LH9211及LH9224，天然杂交种有雷州林业局选育的石岭1号（SH1）和来自湛江的尾叶桉E. urophylla自然杂交种U6，合计10个无性系，全部使用扦插苗，苗高25~30 cm，挑选生长基本均匀一致的苗木，统一于2002年6月雨后造林。

1.3 试验设计和造林

6个试验地，全部采用随机区组设计，广东各试验地重复4次，广西各试验地重复3次。6月间定植，每小区10行每行20株，株行距2 m×3 m，小区面积1 200 m²。3种桉树基肥（S-芬兰的Kemira“星王”复合肥N₁₅P₁₅K₁₅，R-热林所Ritf桉树专用配方肥N₅P₁₈K₈B_0.1，L-雷林1号桉树专用配方肥N₁₀P₁₀K₁₀B_0.1）处理10个无性系(各肥料名称中元素N、P、K、B及下标分别表示N、P₂O₅、K₂O、B的质量百分比)，等质量施加3种桉树基肥，造林前统一每穴施桉树基肥500 g，造林后2个半月结合抚育除草，每株分别追施尿素50 g和氯化钾30 g。第2和第3年4—5月间抚育追施桉树专用肥N₁₆P₄K₈B_0.02 250 g·株^-1，广东各试验地机械开沟小穴种植，广西防城试验地人工垦带后挖穴种植，其余2个试验地采用机耕全垦小穴造林。

1.4 调查指标及方法

试验的调查指标包括成活率、树高、胸径、发病率和风害率等，每小区从第3行第3株开始，每木进行检测，观测3行每行10株，即每小区有效观测30株，3次重复。成活率调查在植后1个月内，各试验地各处理调查3个重复，每重复观测30株。风害调查在台风“伊布都”和“黑格比”登陆广东后进行，其中心风力分别为11级和14级，台风过后4 d内完成。

风害率=受害株数/调查总株数×100%；风折率=折断株数/调查总株数×100%^[9-11]。文中单株材积和蓄积量均取自当地采伐数据，由于年均出材量数据有所缺失，本文中年均出材量（MAI）=单株材积×种植密度（每公顷株数）/材积年份。

1.5 模糊判别的隶属函数

单一指标的排序结果无法恰当反映多指标的综合效应，所以采用模糊数学中隶属函数的方法对各桉树无性系的适生性进行综合评价^[14]。用于分析的隶属函数R(X_i)计算方程为：

$R(X_i)=(X_i-X_{\text{min}})/(X_{\text{max}}-X_{\text{min}})，$

(1)

$R(X_i)=1-(X_i-X_{\text{min}})/(X_{\text{max}}-X_{\text{min}})，$

(2)

式中，X_i为指标测定值，X_min和X_max分别为所有处理的某一指标的最小值和最大值。如果某一指标与生长表现为正相关，则采用式(1) 计算隶属值，反之则用式(2)。累加具体各指标的隶属值，并求出平均值进行比较。

1.6 数据统计与分析

所有试验数据采用Microsoft Excel 2010输入并整理绘图，使用SPSS19.0软件进行方差分析、相关分析和Duncan’s多重比较。

2 结果与分析 2.1 成活率分析

在同一生长环境条件及相同栽培管理的情况下，造林成活率、抗逆性（如抗风、旱、寒、病虫害等）能力的大小，基本可以反映品系（树种/种源/无性系）对环境的适应程度及遗传差异^[12]。造林成活率的高低是树种适生性的第一反应，调查结果（表 2）表明，10个桉树无性系的成活率平均为96.0%，LH9224的成活率(平均93.7%)相对较差，显著低于其余品系（95.2%~97.4%）；试验点间，黏重土壤（唐家）和降雨量相对较少的立地（东门）有相对较低的成活率，分别为93.9%和92.2%，低于平均值。

由表 2可知，桉树成活率的平均值在无性系间的差异小于4%，立地间的差异小于6%，说明绝大部分杂交桉树无性系仍然保持桉树适生性强的优良特性，不论在黏重的立地（唐家）还是沙质土地（遂溪），不论在降雨量大（防城年降雨量>2 000 mm）或小（东门年降雨<1 200 mm）的地方均有理想的成活率。

2.2 抗病性分析

对各试验地各无性系桉树的抗病性调查发现，当年生桉树以LH9224的发病率最高，主要为青枯病造成的死亡或黄株落叶，在唐家、吴川和钦廉其发病率均超过10%，6个试验地其平均发病率为6.7%；其次为M1无性系，其平均发病率为4.08%；以LH9211抗病性最强，在6个试验地1 620株有效调查株中，没有1株病例，其次为DH2012和SH1，发病率只为0.12%（图 1）。从立地类型上分析，桉树在黏重的土壤上发病率较高，同在广东雷州半岛地区的3个试验地，唐家土壤黏重，其平均发病率为3.63%，而中壤土的吴川林地发病率为1.62%，沙质土的遂溪林地发病率为0；东门林场发病率较低则可能是其降雨量明显减少所致（图 2）。

2.3 抗风性分析

由于桉树大多种植在台风多发的东南沿海地区，所以桉树抗风性强弱是其推广种植的重要指标^[6-7]。试验林营建后第2年受台风“伊布都”的影响，广东3个试验地的桉树都受到较大程度的台风危害，经过试验地的风力8~10级。对桉树试验林的抗风情况调查结果列于表 3，认为抗风无性系（风害率1%~10%）有LH9211、LH9224和SH1，较抗风无性系（风害率10%~20%）有U6和M1，抗风性较差的无性系（风害率20%~40%）有DH3213、DH2012和DH3222，抗风性差的无性系（风害率>50%）为DH3327和DH3226。调查还发现DH2012无性系易风折，其次为DH3226，而M1则易折顶（树顶1.0~1.5 m折断），在吴川受风害严重的小区DH2012风折率高达56%，主要在树干基部1~2 m处折断，且发生在生长最好的IIR5小区，在唐家IIL1小区DH3226风折率也高达52%，在吴川M1风折率高达46%。

表 3中结果同时表明，几乎任何无性系其风害率都随风力的增大而显著增大，吴川风力8级，其平均风害率为16.1%，遂溪风力9级，其平均风害率为25.8%，唐家风力10级，其平均风害率为38.2%。

桉树试验林6年生时，台风“黑格比”在广东阳江地区登陆，中心风力14级，主要对其邻近的吴川和遂溪2个试验地试验林造成了危害，吴川风力12级，遂溪风力11级。对所有观察株（每无性系360株）的调查结果（表 4）显示，6年生桉树试验林以M1的抗风性最差，平均风折率达到17.5%，其次为DH2012，平均风折率12.4%，以LH9211和DH3222最为抗风，风折率为1.8%，风倒情况少有发生。

表 3、表 4还表明，不同林龄桉树受台风危害的情况有所不同，幼林易受风害，当风力8级以上时，幼林平均风害率达到16.1%，而当风力11级时成林平均风害率仅为6.9%。幼林主要以风倒或倾斜为主，因其枝叶特别茂盛，有头重脚轻的情况，但无性系间亦有差别，DH2012风折情况同样严重就是例证；成林则主要以风折为主，少见风倒情况，不同林龄不同无性系其抗风性表现亦不相同。总体而言，以来自广东的LH9211和LH9224抗风性最强，其次为SH1和U6，以来自广西东门的DH系列无性系抗风性相对较差。

2.4 速生性分析

短轮伐期桉树的砍伐周期一般为6年。对6.5年生桉树无性系各立地的生长调查结果列于表 5和表 6，从树木的速生性来看，除受台风影响最为严重的唐家试验地外，其余各试验地无论树高和胸径均以DH3213表现最为良好，以U6表现为差，其他无性系在各试验地生长表现各异。

8.5年生时，保存最为完好的钦廉试验地为代表的单株材积见图 3，速生性最好的为DH3213，其单株材积达到0.222 m³·株^-1，年均出材量（MAI）为38.87 m³·hm^-2·a^-1，其次为DH3327和DH3222，单株材积分别为0.167和0.146 m³·株^-1，MAI为29.12和25.46 m³·hm^-2·a^-1，再次为抗风性强的LH9224，单株材积为0.140 m³·株^-1，MAI为24.52 m³·hm^-2·a^-1，表现最差的U6单株材积为0.08 m³·株^-1，MAI为17.19 m³·hm^-2·a^-1。

图 4结果则表明，立地对桉树产量的影响极为明显，好的立地比差的立地增产1倍以上，6.5年生时，遂溪的平均出材蓄积量为79.6 m³·hm^-2，MAI为12.24 m³·hm^-2·a^-1，而东门的平均出材蓄积量为175.3 m³·hm^-2，MAI为26.97 m³·hm^-2·a^-1。

2.5 综合评价

采用模糊数学中隶属函数的方法对各桉树无性系的适生性进行综合评价，隶属值越大，说明桉树的生长表现越好。由表 7结果可知，在广东和广西6个试验地条件下，无性系LH9211、DH3213、SH1、DH3222和U6的隶属值均不小于3.5，高于其他5个无性系，其中，最大隶属值与最小隶属值相差近1倍。表 8表明，由于较低的成活率和最高发病率，种植于唐家的桉树生长综合表现较差，隶属值最低，其他5个立地的隶属值相差不大。就具体某个地区而言：从表 9可以看出，唐家最适合M1和LH9211这2个无性系，隶属值相对较高；从表 10可以看出，吴川适宜种植DH3213、LH9211和SH1这3个无性系，这3个无性系的隶属值排在前三；从表 11可以看出，适宜在遂溪种植的桉树无性系为DH3213和DH3222。

3 讨论与结论

杂交桉树无性系保留了桉树适应性强的特性，各无性系造林成活率均能达到93%以上，平均为96.0%，对比桉树实生苗造林^[13]没有太大的差异，保留了桉树强大生命力的优点，在桉树的适宜种植地区，其成活率均不成为问题。

桉树无性系的选择一要考虑其速生性，二要考虑其抗风性，特别是台风频发的沿海地区。试验中DH系列有相对理想的速生性，但其抗风能力严重不足，在风力9级时即可造成65%以上的风害。因此，建议台风频发地区可多种植抗风性强的尾细桉系列的LH9224无性系，以减少因风倒伏以致重新种植的麻烦，产量亦能达中等以上水平，其在福建闽南地区亦有良好表现^[12]。在内陆少有风害的地区，DH3213和DH3327是较为理想的种植品系，其能使产量提高30%以上，并取得较为理想的经济效益，DH系列的速生性也被前人的试验所肯定^{[13, 15]}。

10个桉树无性系，在唐家和吴川的综合评价隶属值较差，因此，对于成活率差和发病率最高的唐家，只能选择成活率最高的M1和具抗病性的LH9211，而对于出材率较差的吴川和遂溪，适宜选择速生性最佳的DH3213无性系。

立地条件的好坏显著地影响桉树单位面积的产量，良好立地条件下的桉树产量可以比较差立地条件下桉树产量多1倍甚至3倍以上，因此，林地的选择极为重要，虽然桉树的适生性很强，但不同无性系产值有着极大的差别，因此，选好地、种好桉树具有极为现实的意义。

桉树各无性系的抗病性虽然有差别，但目前桉树种植更多地使用组培苗，并在造林地附近直接培育进行造林，一定程度上减少了病害的传播。桉树成林的病害也可能在桉树与相思Acacia confusa、桉树与西南桦Betula spp、桉树与红椎Castanopsis hystrix等的混交种植模式中得以减少，改变桉树大面积纯林种植的方式是减少桉树病虫害发生的有效途径。