杉木(Cunninghamia lanceolata)是我国亚热带地区主要的速生用材树种之一，目前，其人工林面积达1 096万hm²，在我国林业生产中占有重要地位。为了追求较高的经济效益，长期以来，杉木采取纯林、连栽、短轮伐期和皆伐等经营方式，导致林分生产力和土壤质量不断降低^[1-2]，严重影响森林的可持续经营。研究表明，将杉木纯林改造为针阔混交林，可以改良土壤，提高土壤有机质、氮和磷等养分质量分数，进而提高林分生产力^[3]。目前，研究一般是基于少数树种混交得出的结论^[3]，实际的森林经营则采用多树种混交，而不同树种的凋落物分解和根系周转也存在很大差异^[4]；因此，林分改造对土壤的影响存在一定的不确定性。

土壤质量是维持生物生产力、维持环境质量以及促进生命健康的能力^[5]，被认为是土壤物理、化学和生物学特征的联合，对外界环境的变化做出整体响应的体现^[6]；因此，较高的土壤质量是确保林分改造成功的重要条件。20世纪90年代以来，土壤质量评价广泛用于比较不同植被类型对土壤的影响^[7-9]，然而，不同树种组成的同一植被类型的土壤质量评价，尚缺乏深入研究，这既不能加深对森林—土壤相互作用的认识，也不能为低效林分改造提供指导。筛选适宜的树种配置，尤其是针阔叶树种的组合，以加快土壤质量的提高，从而确保林分改造的成功率和效率，是当前我国林业生产的迫切要求。笔者选取南亚热带不同树种配置的杉阔混交林作为研究对象，以杉木纯林的土壤化学性质为对照，分析比较不同林分的土壤物理、化学和微生物特征，并基于主成分分析，对土壤质量进行综合评价，筛选土壤改良效果较好的树种配置，以期为本地区的森林土壤改良和可持续经营，提供科学依据。

研究区位于广东省佛山市高明区(E 112°38′26″~112°42′25″，N 22°41′54″~22°46′50″)。属南亚热带湿润季风气候，年平均气温22.0 ℃，年降水量平均达2 000 mm，有明显的干湿季之分，4—9月的降水量占全年的80%。土壤为花岗岩发育的酸性赤红壤，土层深厚。2002年起，对杉木纯林进行皆伐，保留单株杉木萌芽条，引入阔叶树种，进行随机块状混交造林。造林前整地，当年开展铲草、松土、扩穴和培土等抚育2次，次年和第3年各抚育1次。

2014年12月—2015年1月，选择不同树种组成的11~12年杉阔混交林作为研究对象，建立样地，另选取13年杉木纯林(仅测定其土壤化学性质)作为对照，样地地形和朝向等因素基本一致，坡度为20~30°，坡向朝南，林分概况见表 1。选择林分内代表性部位，在上、中和下坡位，分别挖取3个土壤剖面，用100 cm³环刀和土钻分别对0~25 cm土壤取样，每个坡位3个重复，取约200 g土样装入密封袋，带回实验室分析。

表 1(Table 1)

表 1 试验林概况 Table 1 Characteristics of different stands 林分 Stand 树种组成 Species composition 平均胸径 Average DBH/cm 平均树高 Average height/m 郁闭度 Canopy closure A RC+SC+MG+EG+CL 9.5 7.7 0.70 B CB+ML+CH+ES+CL 8.1 7.7 0.75 C CB+GH+SW+EF+CL 11.1 9.0 0.75 D CB+ML+SW+CH+CL 10.4 8.6 0.75 E SS+LS+SW+SH+CL 11.7 9.1 0.75 F SH+LS+SW+EV+CL 9.7 8.1 0.75 G SH+LS+EG+BB+CL 9.7 8.7 0.70 H SS+LS+SH+CH+BB+CL 10.9 8.7 0.70 I ML+MG+WC+CB+CL 10.8 9.1 0.70 J ML+CC+CB+MM+CL 8.1 7.3 0.70 K ML+LF+CC+WC+CL 9.4 9.7 0.70 L ML+LF+CC+ES+CL 11.7 11.7 0.70 CK CL 11.2 11.4 0.75  注：RC：红花荷 Rhodoleia championii；SC：火焰木 Spathodea campanulata；MG：灰木莲 Mangletia glauca；EG：楝叶吴茱萸 Evodia glabrifolia；CB：阴香 Cinnamomum burmanni；ML：米老排 Mytilaria laosensis；CH：红锥 Castanopsis hystrix；ES：山杜英 Elaeocarpus sylvestris；GH：海南石梓 Gmelina hainanensis；SW：红荷 Schima wallichii；EF：格木 Erythrophleum fordii；SS：木荷 Schima superba；LS：大叶紫薇 Lagerstroemia speciosa；SH：海南蒲桃 Syzygium hainanense；EV：刺桐 Erythrina variegata；BB：红花羊蹄甲 Bauhinia blakeana；WC：盆架子 Winchia calophylla；MM：火力楠 Michelia macclurei；LF：枫香 Liquidambar formosana；CC：香樟 Cinnamomum camphora；CL：杉木 Cunninghamia lanceolata。

表 1 试验林概况 Table 1 Characteristics of different stands

采用环刀法测定杉阔混交林的土壤密度、自然含水量、总孔隙度、毛管孔隙度和非毛管空隙度，土壤粒径组成用筛分法与比重计法测定^[10]。

将杉阔混交林和杉木纯林的土壤样品风干、碾磨和过筛后，用玻璃电极法测定pH值，用K₂Cr₂O₇氧化-外加热法测定土壤有机质，用半微量凯氏法测定全氮，用NaOH碱熔-钼锑抗比色法测定全磷，用碱解-扩散法测定速效氮，用1 mol/L KCl浸提-连续流动分析仪法测定硝态氮和铵态氮质量分数，用盐酸-氟化铵提取-钼锑抗比色法测定速效磷^[10]。

分别采用氯仿熏蒸-K₂SO₄浸提法和氯仿熏蒸-K₂SO₄提取-氮自动分析仪法，测定杉阔混交林的土壤微生物生物量碳和氮^[11]。

第1步，采用土壤质量指数，对不同树种组成的林分土壤质量进行评价。以土壤pH值、土壤有机质、全氮、全磷、速效氮、铵态氮、硝态氮、速效磷、微生物生物量碳、微生物生物量氮、密度、土壤含水量、总孔隙度、砂粒、粉粒和黏粒质量分数作为土壤质量评价的指标。因各指标的重要性不同，量纲也存在很大差异，且各指标的变化具有连续性；故综合评价时，先采用连续性质的隶属度函数进行标准化^[12]，以消除量纲的影响，并确定其权重。其中，土壤密度和黏粒质量分数用降型函数(式1)，其他指标用升型函数(式2) 表示。

式中：F(X_i)为各指标的隶属度值，X_ij为各指标观测值，X_imax和X_imin分别为第i项因子中的最大值和最小值。

第2步，对隶属度值进行主成分分析，得到各指标的公因子方差，再计算各指标的权重

式中：w_i为第i项因子的权重，C_i为第i项因子的公因子方差。

第3步，应用模糊数学中的模糊集加权综合方法，计算土壤质量指数SQI

利用SPSS 16.0进行单因素方差分析和主成分分析，不同林分之间各指标的差异，采用Duncan法进行多重比较，显著水平α=0.05。结果中各指标的数值为平均值±标准误。

由表 2可知，不同林分土壤物理性质的差异较大。其中，林分B、D、I、J、K和L的土壤密度相对较小，土壤含水量、总孔隙度、毛管孔隙度和非毛管孔隙度较大；林分A、C、E、F、G和H的密度相对较大，土壤含水量、总孔隙度、毛管孔隙度和非毛管孔隙度相对较小。土壤颗粒组成在林分间的差异相对较小，以砂粒最高，粉粒其次，黏粒最低。除非毛管孔隙度和砂粒百分比外，其余指标在不同林分间的差异显著(P < 0.05)。

表 2(Table 2)

表 2 不同林分的土壤物理性质 Table 2 Soil physical characteristics of different stands 林分 Stand 密度 Bulk density/ (g·cm-3) 土壤含水量 Soil water content/% 总孔隙度 Total porosity/% 毛管孔隙度 Capillary porosity/% 非毛管孔隙度 Non-capillary porosity/% 砂粒比例 Sand ratio/% 粉粒比例 Silt ratio/% 黏粘粒比例 Clay ratio/% A 1.25±0.03bc 30.0±2.4abcd 49.6±4.7abcd 42.3±3.0abc 7.3±3.2a 52.7±7.4a 28.3±4.0c 19.0±3.5abc B 1.13±0.13cd 39.9±9.9a 54.9±8.0a 48.2±5.3a 6.7±2.8a 53.3±2.5a 34.3±0.6ab 12.3±2.1c C 1.36±0.11ab 26.0±7.1bcde 45.2±4.6abcde 40.2±2.3abc 5.0±2.3a 53.0±4.0a 32.3±2.5abc 14.7±5.5bc D 1.05±0.12d 32.0±2.9abc 52.5±4.0ab 43.2±3.1abc 9.3±4.8a 48.7±0.6a 29.3±1.5bc 22.0±1.7a E 1.44±0.05a 20.8±5.0de 40.9±1.6cde 38.7±4.0abc 2.2±2.4a 52.7±6.4a 28.3±3.5c 19.0±4.6abc F 1.42±0.05a 21.4±6.1cde 39.4±7.8de 36.7±6.5bc 2.8±1.4a 47.7±4.0a 31.7±0.6abc 20.7±3.5ab G 1.47±0.09a 18.6±5.9e 35.7±7.8e 33.3±8.4c 2.4±1.1a 50.0±3.6a 28.7±6.7c 21.3±5.0ab H 1.37±0.04ab 23.9±5.1bcde 42.3±5.2bcde 39.8±3.6abc 2.5±1.6a 46.3±4.5a 35.3±2.5a 18.3±4.0abc I 1.15±0.12cd 31.0±2.4abcd 50.6±3.8abc 46.7±2.6ab 4.0±2.5a 47.7±2.5a 33.3±2.1abc 19.0±1.0abc J 0.98±0.18d 34.5±2.4ab 55.8±5.6a 49.0±3.0a 6.7±6.4a 49.7±2.1a 32.7±0.6abc 17.7±1.5ab K 1.04±0.04d 26.7±7.9bcde 53.2±7.4a 43.3±7.7abc 9.9±4.7a 50.3±3.5a 30.0±1.7abc 19.7±4.0ab L 0.99±0.04d 31.6±5.5abcd 50.9±1.8abc 44.9±9.3ab 6.0±7.8a 50.7±2.3a 31.0±1.0abc 18.3±3.2abc   注：表中同一列数据后不同字母表示差异显著(P < 0.05)，下同。Note: The different letters at the same column of data in the table indicate significant differences (P < 0.05). The same below.

表 2 不同林分的土壤物理性质 Table 2 Soil physical characteristics of different stands

由表 3可知：杉木林的土壤pH值(4.06) 低于混交林(4.15~4.71)，除林分D外，差异均达到显著水平(P < 0.05)；杉木林的速效氮质量分数(30.8 mg/kg)显著小于所有混交林(81.4~152.8 mg/kg)(P < 0.05)，速效磷质量分数(0.8 mg/kg)也小于所有混交林(0.9~2.4 mg/kg)；杉木林土壤全磷质量分数(0.34 g/kg)高于除林分C(0.45 g/kg)外的其他混交林(0.14~0.32 g/kg)。混交林土壤有机质、全氮、铵态氮和硝态氮的变化，其相对于杉木林没有一致的规律。

表 3(Table 3)

表 3 不同林分的土壤化学性质 Table 3 Soil chemical properties in different stands 林分 Stand pH值 pH value 土壤有机质 Soil organic matter/ (g·cm-3) 全氮 Total N/ (g·kg-1) 全磷 Total P/ (g·kg-1) 速效氮 Available N/ (mg·kg-1) 铵态氮 NH4+-N/ (mg·kg-1) 硝态氮 NO3--N/ (mg·kg-1) 速效磷 Available P/ (mg·kg-1) A 4.71±0.14 a 31.7±4.6 cd 1.16±0.27 cde 0.15±0.01 e 115.3±21.1 abc 4.0±1.8 c 2.8±0.5 c 1.2±0.4 cde B 4.36±0.07 de 30.6±4.8 d 1.05±0.20 de 0.14±0.01 e 97.9±7.7 bc 4.2±0.7 c 1.7±1.1 c 1.4±0.4 bcde C 4.64±0.13 ab 16.4±3.4 f 0.66±0.16 f 0.45±0.26 a 87.4±15.9 c 7.0±2.3 bc 0.9±0.0 c 1.2±0.5 cde D 4.15±0.07 fg 37.2±1.8 c 1.34±0.06 bc 0.15±0.01 e 110.1±5.2 bc 6.1±3.7 c 1.2±0.8 c 0.9±0.1 de E 4.41±0.05 cd 25.5±4.1 de 1.15±0.14 cde 0.21±0.02 bcde 93.1±10.8 c 11.0±2.6 ab 6.1±1.9 ab 2.1±0.7 abc F 4.40±0.04 cd 24.1±2.2 e 1.00±0.070 e 0.19±0.03 cde 92.1±3.7 c 15.7±2.7 a 3.4±2.3 bc 1.6±0.4 abcde G 4.53±0.10 bc 21.4±3.6 ef 0.99±0.17 e 0.15±0.02 e 109.5±58.5 bc 11.9±1.9 a 1.0±0.8 c 2.4±0.9 a H 4.37±0.09 de 21.3±1.9 ef 0.96±0.09 e 0.18±0.03 de 81.4±7.2 c 15.4±2.4 a 2.3±0.9 c 2.4±0.8 a I 4.38±0.07 d 45.1±0.4 b 1.45±0.12 b 0.32±0.01 abc 119.5±15.4 abc 11.9±3.5 a 2.6±2.3 c 1.9±0.1 abcd J 4.27±0.05 def 52.7±1.8 a 1.78±0.00 a 0.32±0.03 abc 133.4±7.4 ab 12.4±3.1 a 1.3±1.0 c 1.7±0.4 abcde K 4.22±0.08 ef 49.8±7.1 ab 1.85±0.16 a 0.31±0.01 abcd 152.8±25.4 a 11.7±0.6 a 6.6±3.8 a 2.3±0.2 ab L 4.26±0.02 def 51.3±4.6 ab 1.72±0.10 a 0.29±0.01 bcd 136.4±18.8 ab 11.6±3.1 a 2.8±2.2 c 2.1±0.6 abc CK 4.06±0.07 g 23.1±1.1 e 1.28±0.19 bc 0.34±0.02 ab 30.8±4.4 d 11.6±1.6 a 1.7±0.5 c 0.8±0.2 e

表 3 不同林分的土壤化学性质 Table 3 Soil chemical properties in different stands

由图 1可知，土壤微生物生物量碳和氮质量分数，在不同林分间存在显著差异(P < 0.05)，分别介于346.8~821.7与32.5~65.4 mg/kg之间，二者在不同林分间的大小顺序不一致。

图 1(Figure 1)

不同小写字母表示不同林分间差异显著P < 0.05。Values with different letters refer to significant differences at 0.05 level among different stands. 图 1 不同林分土壤微生物生物量碳和氮质量分数 Figure 1 Soil microbial biomass carbon (a) and nitrogen (b) content of different stands

按照特征值>1的原则，提取出4个主成分(表 4)，累积方差贡献率为85.1%，说明这4个相互独立的主成分，可以反映土壤各指标分量总变异的80%以上。第1主成分方差贡献率为35.3%，主要综合了土壤有机质、全氮、密度、土壤含水量和总孔隙度的信息，因子载荷都大于0.8；第2主成分方差贡献率为24.2%，主要综合了铵态氮和速效磷的信息，因子载荷分别为0.839和0.764；第3、第4主成分因子中，粉粒、黏粒和砂粒质量分数的因子载荷相对较大。

表 4(Table 4)

表 4 土壤质量指标的主成分因子荷载、公因子方差及权重 Table 4 Principal factor loading of soil quality indicators, and communalities and weights 指标 Indicator 主成分 Principal factor 公因子方差 Communalities 权重 Weight 1 2 3 4 pH值 pH value -0.507 -0.572 -0.255 0.253 0.713 0.052 4 土壤有机质 Soil organic matter 0.889 0.417 0.015 0.052 0.968 0.071 1 全氮 Total N 0.800 0.571 -0.084 0.100 0.984 0.072 3 全磷 Total P 0.211 0.104 0.384 0.486 0.439 0.032 2 速效氮 Available N 0.771 0.440 -0.230 0.235 0.896 0.065 8 铵态氮 NH4+-N -0.351 0.839 0.287 -0.103 0.921 0.067 6 硝态氮 NO3--N 0.035 0.493 -0.284 0.542 0.619 0.045 5 速效磷 Available P -0.290 0.764 0.123 0.280 0.762 0.056 0 土壤微生物生物量碳 Soil microbial biomass C 0.647 -0.564 -0.343 -0.153 0.877 0.064 4 土壤微生物生物量氮 Soil microbial biomass N 0.352 -0.108 -0.738 -0.268 0.752 0.055 2 密度 Bulk density 0.974 0.099 0.130 -0.066 0.979 0.071 9 土壤含水量 Soil water content 0.819 -0.421 0.315 -0.136 0.966 0.070 9 总孔隙度 Total porosity 0.956 -0.183 0.169 -0.007 0.975 0.071 6 砂粒比例 Sand ratio 0.112 -0.622 -0.166 0.708 0.929 0.068 2 粉粒比例 Silt ratio 0.040 -0.077 0.898 -0.266 0.884 0.064 9 黏粒比例 Clay ratio 0.133 -0.608 0.650 0.379 0.953 0.070 0 特征值 Eigenvalue 5.641 3.878 2.514 1.585 方差贡献率 Contribution rate of variance/% 35.254 24.236 15.710 9.905 累积贡献率 Cumulative contribution rate/% 35.254 59.490 75.200 85.105

表 4 土壤质量指标的主成分因子荷载、公因子方差及权重 Table 4 Principal factor loading of soil quality indicators, and communalities and weights

不同杉阔混交林中，土壤质量综合指数的变化范围是0.25~0.64(图 2)。其中，林分J、K和L的土壤质量指数均大于0.6，林分A、B和I介于0.4~0.6之间，其余林分均小于0.4。

图 2(Figure 2)

图 2 不同林分的土壤质量指数 Figure 2 Soil quality indices of different stands

笔者研究发现，土壤理化性质在不同林分间存在显著差异。由于这些林分在改造前都长期连续种植杉木，土地利用历史一致，而林龄和气候等其他条件也一致，因此，不同林分间的土壤理化性质差异，主要来源于树种的影响。不同树种的凋落物和细根周转动态差异较大^[4]，如生长较快的米老排可以向土壤输入更多的有机质^[13]。较高质量分数的土壤有机质，意味着为土壤动物和微生物提供更多的食物；同时，植物根系更为活跃，周转更快。土壤动物、微生物和根系的物理活动的提高，都有利于土壤变得疏松透气，形成更多的毛管，以储存更多的水^[14]。因此，林分I、J、K和L的土壤有机质质量分数更高，全氮、速效氮质量分数也相应更高，土壤密度更小，含水量和孔隙度更大；而土壤颗粒组成，主要取决于漫长的母岩风化过程^[15]，相比之下，生物在短时间内的作用影响显得微乎其微。因此，笔者研究中，土壤颗粒组成差异不明显。

研究中多数混交林的土壤全磷质量分数低于杉木林，而pH值、有效氮和有效磷质量分数均高于杉木林，表明林分改造缓解了土壤酸化的情况；同时，促进氮和磷等养分元素的释放，与前人的研究结果一致^{[3, 16]}。不同树种的凋落物混合分解时，可能存在促进或抑制效应^[17-18]，如桤木(Alnus cremastogyne)、火力楠可以促进杉木分解，而红锥、樟树和木荷等则无效^[4]，这可能是不同混交林土壤化学性质差异的主要原因。此外，杉木凋落物较难分解，导致土壤中有机酸积累^[19]，而磷在强酸性的土壤中，易与土壤中的Al³⁺、Fe³⁺等金属离子形成螯合物，沉淀下来^[15]，最终表现为土壤pH值和有效磷质量分数等降低。

土壤微生物生物量碳和氮在不同混交林间差异显著，但规律性不明显。大量研究表明，植物物种对土壤微生物的数量和组成有影响^[20-21]，但其机制尚不清楚。G.Berg等^[20]认为，主要是每种植物特定的根系分泌物，决定了植物物种与微生物之间特定的互利关系。D. R. Zak等^[21]发现，微生物生物量取决于植物生产力和多样性的共同作用，而非单纯的植物多样性。

笔者研究不同杉阔混交林的土壤质量指数介于0.25~0.64之间，变化范围相对较大。在选取土壤质量评价指标时，土壤含水量、密度和有机质质量分数等理化性质是关键指标^[7-9]，笔者研究也证实了这一点。此外，一些研究认为，土壤生物学特征是影响土壤质量的主要指标之一^[7-8]，笔者研究未有类似发现，这可能是受采样时间的影响。如罗达等^[22]发现，最大微生物生物量在旱季的马尾松林中出现，而雨季则在格木林中出现。

笔者研究表明，将杉木纯林改造为混交林时，土壤质量将因树种组成而异，这与广州南沙地区的研究结论一致^[23]。可见树种的选择，对于土壤改良效果至关重要。米老排和香樟是林分J、K和L的共有造林树种，现有研究也表明，这2个树种生物量和凋落物量较大，土壤改良效果较好^[13]，林分J、K和L的树种配置模式，适宜在南亚热带的林分改造中推广应用；而在土壤质量较低的林分E、F、G和H中，海南蒲桃和大叶紫薇是共有造林树种，前者属于小乔木，后者是园林中常用的落叶树种，在自然生境中，生长相对较慢，生物量和凋落物量较低，可能是林分土壤质量较差的主要原因。因此，在林分改造中，这些树种可能只适合作为伴生树种，不宜作为建群种。

南亚热带杉木纯林改造为杉阔混交林后，土壤pH值、速效氮和速效磷质量分数均升高，全磷质量分数降低，不同杉阔混交林的土壤质量指数介于0.25~0.64之间。其中，米老排+枫香+香樟+盆架子+杉木、米老排+香樟+阴香+火力楠+杉木、米老排+枫香+香樟+尖叶杜英+杉木这3种树种配置方式的土壤改良作用明显。在进行多树种混交的林分改造时，树种组成的不同，导致土壤质量产生较大差异。笔者筛选的树种配置方式，只能作为整体应用，还不能量化每一树种的具体作用，无疑将限制其实际应用范围；因此，优良树种的种间关系、混交比例和造林密度等对土壤质量的影响，是下一步研究中应当关注的重点。