柚木(Tectona grandis)为马鞭草科(Verbenaceae)柚木属落叶半落叶大乔木，树高35～45 m，胸径0.9～2.5 m，材质优良，是国际上最重要的热带珍贵用材树种之一，可用做高级家具及装饰材料、乐器以及车船等。我国无柚木天然分布，引种历史可追溯到1820年，1960年后开始大量发展，现已引种至南方十省(区)60多个县(市)。近十多年来柚木天然林急剧减少，传统出口国相继减少甚至停止原木出口，国际市场上柚木原木与锯材价格逐年攀升；随着国内经济发展和人民生活水平提高，珍贵实木家具和室内装饰的需求量越来越大，国内柚木人工林种植拥有广阔前景。在印度、缅甸和尼日利亚等国，柚木人工林发展时间长，研究历史悠久，尤其是在立地质量评价和土壤方面。国外学者通常采用地位指数法(Malende et al.，1990； Akindele，1991； Nanang et al.，1999； Sajjaduzzaman et al.，2005； Upadhyay et al.，2005； Jerez-Rico et al.，2011； Torres et al.，2012)评价柚木立地质量，此外在柚木 林土壤方面研究颇多(Zech et al.，1991； Drechsel et al.， 1994； Amonsah et al.，2000； Rugmini et al.，2007；Watanabe et al.，2010； Samndi et al.，2012； Ziblim et al.，2012)。我国引种柚木接近两百年，在柚木立地方面也有不少，卢俊培(1994)以海南32块样地的地形因子和9个土壤理化性质对海南柚木人工林进行立地分类；陶国祥(2005)采用地形因子、土壤厚度和质地但缺乏土壤理化性质来划分云南柚木人工林立地类型；张树芬等(2005)从海拔、地形因子、土层厚度和植被等分析对柚木生长的影响，此外一些学者(热带林业研究所柚木调查组，1977； 武勇，2005； 郭华北，2009； 顾永顺等，2011； 施荣达，2013)也对柚木人工林生长与立地的关系做了初步分析。本研究以我国柚木主要栽培区的云南西双版纳12～26年生、具一定规模的柚木人工林为对象，通过51块样地的植被生长、地形因子和19个土壤理化性质调查，采用主成分分析和因子分析划分土壤养分等级，尝试以土壤质地、土壤养分及土壤密度划分立地类型，旨在为柚木人工林种植的立地选择及合理经营等提供科学依据，以期促进我国柚木人工林栽培和推广。

研究区位于云南省南端的西双版纳傣族自治州景洪市勐罕镇，勐腊县的勐伴、勐仑和勐棒镇(21°24′29.1″—21°54′46.3″E，101°05′05.1″—101°36′46.2″N)，属北回归线以南的热带北缘区，受印度洋季风气候控制。年降水量1 100～1 900 mm，5—10月为雨季；年均气温18～22 ℃，最冷月均气温8.8～15.6 ℃，≥10 ℃年积温5 062～8 000 ℃，海拔800 m以下地区2013年年积温皆在7 500 ℃以上；气候温暖、湿润，长夏无冬，夏热多雨，秋春相连且为期较短。海拔1 500 m以下的低山丘陵和低中山地以砖红壤和赤红壤为主，林下植被以华南毛蕨(Cyclosorus parasiticus)、蔓生莠竹(Microstegium vagans)、海金沙(Lygodium japonicum)和飞机草(Eupatorium odoratum)居多。

2 014年11—12月对区域内柚木人工林设置样地，要求是同一起源、未间伐，疏密度在0.4以上，面积2 hm²以上的人工林，且为样地的坡向、坡度、坡位一致的同龄纯林，共设置51块20 m×20 m样地，其中景洪市勐罕镇勐波100 hm²柚木人工林设置11块样地，勐腊县的勐伴镇300 hm²设置36块样地和其他镇共设置4块样地。记录每块样地的地貌、坡位、坡度和坡向等地形因子，使用GPS仪确定每块样地的经纬度和海拔高度；对样地内柚木每木调查树高、胸径和枝下高，选取5株最高柚木树高平均值为优势高；面向样地沿样地左下角向右上角的对角线方向，选取样地左下角、中间和右上角3个点挖土壤剖面，按0～20 cm和20～40 cm分层取3点混合土样，进行土壤化学性质分析，并用环刀在每点每层取原状土测定土壤密度、孔隙度、田间持水量、机械组成等土壤物理性质。

电位法测定土壤pH值(水溶液浸提)；烘干法、环刀法测定土壤密度；比重法测定土壤机械组成。重铬酸钾氧化-外加热法测定有机质含量；凯氏消煮法测全N含量；硫酸-高氯酸消煮-钼锑抗分光光度计法测定全P含量；高氯酸-氢氟酸消煮-火焰光度法测定全K含量；盐酸和硫酸溶液浸提法测定有效P含量；0.1 mol·L^-1盐酸浸提原子吸收光谱法测定有效Zn²⁺和Fe(Fe²⁺，Fe³⁺)含量；1 mol·L^-1乙酸铵浸提-火焰光度法测定有效K含量；乙酸铵交换-原子吸收分光光度法测定交换性Ca²⁺和Mg²⁺含量；1 mol·L^-1氯化钾-中和滴定法测定交换性酸、交换性Al³⁺含量；依据中华人民共和国林业行业标准LY/T1225-1999，测定阳离子交换量和盐基饱和度。

数据处理利用Excel 2010和SPSS 19。

采用主成分分析、相关性分析、回归分析和数量化方法Ⅰ(董文泉，1979)对立地因子进行处理。数量化方法Ⅰ原理的核心是通过赋值(0，1)使定性因子定量后建立数学模型，确定变量对基准变量的影响程度并进行预测，从而发现内在联系和数量规律。

各样地调查主要立地因子与优势木年平均高生长量如表 1。根据立地因子的变化幅度及对林木生长发育的影响程度，并结合相关林业调查技术标准，对各因子分级处理，其类目划分如表 2。

表 1 各样地的立地因子与优势木年平均高生长量 Tab.1 Site factors and mean annual increments of dominant tree height of different plots

表 1 各样地的立地因子与优势木年平均高生长量 Tab.1 Site factors and mean annual increments of dominant tree height of different plots 样地 Plot 坡度 Slope degree 坡向 Aspect 坡位 Slope position 土壤密度 Soil density/ (g·cm-3) 质地 Texture 海拔 Elevation/m 优势木年 平均高生 长量 Mean annual in crement of dominant tree height/m 1 13.8 半阴坡Half-shady slope 中坡Middle slope 1.15 粉砂质壤土Silt loam 741 1.95 2 13.9 半阴坡Half-shady slope 上坡Upper slope 1.12 重壤土Heavy loam 761 1.90 3 13.8 半阳坡 Half-sunny slope 上坡Upper slope 1.15 壤土Loam 785 1.50 4 16.8 阴坡Shady 中坡Middle slope 1.14 粉砂质壤土Silt loam 763 1.83 5 30.0 阳坡Sunny slope 下坡Lower slope 1.03 重壤土Heavy loam 725 1.41 6 21.0 半阴坡Half-shady slope 上坡Upper slope 0.76 重壤土Heavy loam 863 0.43 7 20.6 阳坡Sunny slope 上坡Upper slope 0.92 重壤土Heavy loam 898 0.60 8 19.9 阳坡Sunny slope 下坡Lower slope 1.12 壤土Loam 894 1.32 9 26.0 阳坡Sunny slope 中坡Middle slope 1.01 壤土Loam 896 1.33 10 35.5 阳坡Sunny slope 上坡Upper slope 1.00 重壤土Heavy loam 909 1.12 11 28.6 半阳坡Half-sunny slope 中坡Middle slope 0.90 壤土Loam 954 1.26 12 28.7 阳坡Sunny slope 下坡Lower slope 1.05 重壤土Heavy loam 932 1.36 13 20.6 半阴坡Half-shady slope 中坡Middle slope 0.90 壤土Loam 934 1.24 14 26.5 半阳坡Half-sunny slope 中坡Middle slope 0.99 壤土Loam 911 1.54 15 22.8 半阴坡Half-shady slope 下坡Lower slope 1.04 壤土Loam 882 1.29 16 26.2 阳坡Sunny slope 下坡Lower slope 0.94 壤土Loam 854 1.27 17 32.6 半阳坡Half-sunny slope 下坡Lower slope 0.90 重壤土Heavy loam 857 1.33 18 23.7 半阳坡Half-sunny slope 上坡Upper slope 0.97 粉砂质壤土Silt loam 950 1.41 19 29.2 阳坡Sunny slope 中坡Middle slope 1.09 壤土Loam 954 1.45 20 14.9 半阳坡Half-sunny slope 上坡Upper slope 1.16 壤土Loam 932 1.10 21 13.4 半阳坡Half-sunny slope 中坡Middle slope 0.89 重壤土Heavy loam 880 0.50 22 12.4 半阳坡Half-sunny slope 上坡Upper slope 0.88 壤土Loam 970 0.97 23 26.0 半阳坡Half-sunny slope 下坡Lower slope 1.08 粉砂质壤土Silt loam 919 1.66 24 25.8 半阳坡Half-sunny slope 上坡Upper slope 0.98 重壤土Heavy loam 940 1.16 25 26.0 半阴坡Half-shady slope 中坡Middle slope 0.88 壤土Loam 975 0.95 26 27.9 阳坡Sunny slope 中坡Middle slope 1.10 壤土Loam 915 1.33 27 14.1 半阴坡Half-shady slope 中坡Middle slope 0.84 壤土Loam 877 1.09 28 30.8 阳坡Sunny slope 下坡Lower slope 0.91 重壤土Heavy loam 900 1.31 29 31.0 阳坡Sunny slope 上坡Upper slope 0.91 重壤土Heavy loam 910 1.06 30 21.4 阳坡Sunny slope 上坡Upper slope 0.97 壤土Loam 999 0.82 31 20.3 阳坡Sunny slope 上坡Upper slope 1.02 砂质壤土Sandy loam 994 1.08 32 8.8 半阴坡Half-shady slope 中坡Middle slope 1.03 壤土Loam 770 1.35 33 23.5 阳坡Sunny slope 上坡Upper slope 1.00 壤土Loam 820 1.24 34 29.0 阳坡Sunny slope 下坡Lower slope 0.99 重壤土Heavy loam 844 1.50 35 18.5 半阴坡Half-shady slope 中坡Middle slope 0.97 壤土Loam 853 1.39 36 27.5 阳坡Sunny slope 中坡Middle slope 1.15 重壤土Heavy loam 836 1.30 37 0.0 平坡Flat slope 平地Flat ground 1.17 壤土Loam 582 1.17 38 23.7 半阳坡Half-sunny slope 下坡Lower slope 1.10 重壤土Heavy loam 638 0.96 39 20.0 半阳坡Half-sunny slope 下坡Lower slope 1.25 砂质壤土Sandy loam 650 1.31 40 30.9 半阳坡Half-sunny slope 中坡Middle slope 1.15 砂质壤土Sandy loam 663 1.26 41 24.2 半阴坡Half-shady slope 中坡Middle slope 0.94 壤土Loam 956 1.16 42 18.0 半阳坡Half-sunny slope 上坡Upper slope 0.88 壤土Loam 930 1.12 43 18.0 半阳坡Half-sunny slope 上坡Upper slope 0.88 壤土Loam 938 1.16 44 19.7 半阳坡Half-sunny slope 中坡Middle slope 1.00 重壤土Heavy loam 908 1.12 45 22.6 半阳坡Half-sunny slope 下坡Lower slope 0.97 壤土Loam 898 1.01 46 34.0 半阴坡Half-shady slope 中坡Middle slope 0.83 重壤土Heavy loam 917 0.99 47 0.0 平坡Flat slope 平地Flat ground 0.98 壤土Loam 952 1.19 48 0.0 平坡Flat slope 平地Flat ground 0.93 重壤土Heavy loam 928 1.15 49 18.4 半阳坡Half-sunny slope 下坡Lower slope 0.99 壤土Loam 832 0.98 50 37.0 半阳坡Half-sunny slope 下坡Lower slope 0.86 重壤土Heavy loam 738 1.10 51 25.0 半阳坡Half-sunny slope 下坡Lower slope 0.94 壤土Loam 788 1.11

表 2 立地因子类目划分 Tab.2 Category division of site factors

表 2 立地因子类目划分 Tab.2 Category division of site factors 项目Item 类目 Category 质地 Texture 砂壤土 Sandy loam 粉砂壤土 Silt loam 壤土 Loam 重壤土 Heavy loam 坡向 Aspect 阳坡 Sunny slope 半阳坡 Half-sunny slope 半阴坡 Half-shady slope 阴坡 Shady slope 平坡 Flat slope 坡位 Slope position 上坡 Upper slope 中坡 Middle slope 下坡Lower slope 平地 Flat ground 坡度 Slope degree 平缓坡 Conservative slope(0～5°) 缓坡Gentle slope(5～15°) 斜坡Slope(15～25°) 陡坡Steep slope (25～35°) 土壤密度Soil density ＞1.10 g·cm-3 1.10～1.00 g·cm-3 1.00～0.90 g·cm-3 ＜0.90 g·cm-3 海拔 Elevation ＜700 m 700～800 m 800～900 m ＞900 m

对土壤化学性质与柚木优势木年均高生长量进行相关性分析，结果表明(表 3)：土壤pH值、全K含量、盐基饱和度、交换性Mg²⁺含量以及上层土壤有效P含量与柚木优势木年均高生长量极显著正相关(P＜0.01)；土壤交换性Al³⁺含量、交换性酸含量和有效Fe含量与柚木优势木年均高生长量极显著负相关(P＜0.01)。盐基饱和度是衡量土壤肥力的重要因素，其值越高说明土壤越肥沃；在一定范围内，土壤pH值越高，柚木生长越好；全K和交换性Mg²⁺含量越高，交换性Al³⁺、交换性酸和有效Fe含量越低，柚木生长也越好。

表 3 土壤化学性质与柚木优势木年均高生长量的相关性^① Tab.3 Correlation analysis of soil chemical properties and mean annual increment of dominant tree height

表 3 土壤化学性质与柚木优势木年均高生长量的相关性① Tab.3 Correlation analysis of soil chemical properties and mean annual increment of dominant tree height 土层 Soil layer pH 有机质 含量 Organic matter content 全N 含量 Total N content 全P 含量 Total P content 全K 含量Total K content 盐基 饱和度 Base satur- ation 速效K 含量 Availa- ble K content 交换性Exchangeable 阳离子 交换量 Cation exchange capacity 有效Available Ca2+ Mg2+ Al3+ Acid P Zn2+ Fe2+/Fe3+ 0～20 cm 0.541** -0.225 -0.014 -0.019 0.407** 0.468** 0.023 0.110 0.389** -0.521** -0.542** -0.231 0.368** 0.090 -0.405** 20～40 cm 0.419** 0.201 -0.051 -0.060 0.429** 0.521** -0.028 0.205 0.485** -0.536** -0.546** -0.165 0.049 -0.022 -0.414** ①**：α=0 .01

将土壤化学性质因子与优势木年均高生长量做逐步回归分析，选择偏相关系数较大的pH值、交换酸含量、全K含量、全P含量、有效P含量、交换性Ca^2＋含量、交换性Al^3＋含量和有效Fe含量8项化学性质指标，采用主成分分析法划分土壤养分等级。

根据所测得8项化学性质，对其进行因子分析，结果如表 4。前3个主成分的累积贡献率接近80%(表 5)，故可选择前3个主成分进行分析。根据因子分析选择的3个主成分的特征根λ及各成分矩阵，计算标准化特征向量。在对原始数据变量进行标准化后，计算各主成分，以各标准化后的变量矩阵乘以标准化后的正交特征向量矩阵，矩阵的乘法在Excel中进行。计算出的各主成分分别乘以各方差分析比，最终得出综合得分，并以此为依据51个样地划分为4个土壤养分等级(表 6)。

表 4 化学性质指标分析 Tab.4 Factor analysis of chemical properties

表 4 化学性质指标分析 Tab.4 Factor analysis of chemical properties 成份 Component 初始特征值 Initial eigenvalues 合计 Total 方差 Variance(%) 累积 Accumulation(%) 1 3.822 47.775 47.775 2 1.476 18.444 66.219 3 1.094 13.676 79.895 4 0.649 8.116 88.011 5 0.455 5.682 93.692 6 0.317 3.959 97.651 7 0.183 2.285 99.936 8 0.005 0.064 100.000

表 5 因子分析成分矩阵 Tab.5 Factor analysis of component matrix

表 5 因子分析成分矩阵 Tab.5 Factor analysis of component matrix 化学性质 Chemical properties 成份 Component a1 a2 a3 pH -0.811 0.167 0.085 交换性酸含量 Exchangeable acid content 0.939 0.001 0.139 全P 含量 Total P content -0.519 0.421 -0.371 全K 含量 Total K content -0.299 0.325 0.849 有效P 含量 Available P content -0.141 0.911 -0.085 交换性Ca2+含量 Exchangeable Ca2+ content -0.708 -0.439 -0.199 交换性Al3+含量 Exchangeable Al3+ content 0.950 0.002 0.124 有效Fe含量 Available Fe content 0.707 0.377 -0.383

表 6 土壤养分等级分组 Tab.6 Classification of sample plots based on soil nutrient grades

表 6 土壤养分等级分组 Tab.6 Classification of sample plots based on soil nutrient grades 土壤养分等级 Soil nutrient grades 样地 Sample plot 贫养分 Poverty nutrient 6，7，11，16，17，21，24，27，28，29，30 低养分 Low nutrient 9，10，12，13，15，20，22，25，32，33，35，40 中养分 Middle nutrient 3，8，18，19，23，26，31，34，36，37，38，39， 42，43，44，45，46，47，48，49，51 高养分 High nutrient 1，2，4，5，14，41，50

根据调查样地情况，结合以上类目划分表(表 2)，列出[0，1]反应表，按照数量化理论Ⅰ模型计算各类目与柚木优势木年均高生长量的回归系数及有关参数。回归分析计算结果如表 7所示。

表 7 立地因子数量化回归分析结果 Tab.7 Quantity regression result of site factors

表 7 立地因子数量化回归分析结果 Tab.7 Quantity regression result of site factors 项目 Item 类目 Category 代码 Code 得分 Score 得分范围及所 占比例 Score range and Proportion 偏相关系数 Coefficient of partial correlation t检验 t-test 常量项constant 1.268 贫养分 Poverty nutrient X11 -0.176 养分等级 低养分 Low nutrient X12 0.003 0.209(11.741%) 0.435** 3.168 Nutrient gradeX1 中养分 Middle nutrient X13 0 高养分 High nutrient X14 0.033 ＞1.10 X21 0 土壤密度 1.10~1.00 X22 -0.181 0.450(25.280%) 0.665** 5.843 Soil density X2 1.00~0.90 X23 -0.224 ＜0.90 X24 -0.450 陡坡 Steep slope X31 0.245 坡度 斜坡 Slope X32 0 0.245(13.764%) 0.584** 4.712 Slope degree X3 缓坡 Gentle slope X33 0.145 平缓坡 Conservative slope X34 0 ＜700 X41 -0.144 海拔 700~800 X42 0.099 0.243(13.652%) 0.347* 2.423 Elevation X4 800~900 X43 0 ＞900 X44 0.053 半阳坡 Half-sunny slope X51 0.008 坡向 半阴坡 Half-shady slope X52 0.142 0.142(7.978%) 0.349* 2.446 Aspect X5 阴坡 Shady slope X53 0.048 平坡 Flat slope X54 0 阳坡 Sunny slope X55 0 上坡 Upper slope X61 0.006 坡位 下坡 Lower slope X62 0 0.113(6.348%) 0.312* 2.153 Slope position X6 平地 Flat ground X63 0.110 中坡 Middle slope X64 0.113 粉砂质壤土 Silt loam X71 0.237 质地 砂质壤土 Sandy loam X72 -0.043 0.378(21.236%) 0.534** 4.138 Texture X7 重壤土 Heavy loam X73 -0.141 壤土 Loam X74 0

数量化回归模型如下：

从立地因子回归分析结果可知，所选立地因子中，除海拔外，与优势木年均高生长量的偏相关系数均达显著水平(P＜0.05)，其中土壤因子、坡度因子均达极显著水平(P＜0.01)。从得分范围来看，土壤性质(土壤养分、土壤密度和土壤质地)在所选立地因子中得分占比58%，其次是地势和地形条件(坡度、海拔、坡向和坡位)。回归模型的复相关系数R=0.865，达到极显著水平。将回归模型中7个立地因子(养分等级X₁，土壤密度X₂，坡度X₃，海拔X₄，坡向X₅，坡位X₆和土壤质地X₇)与优势木年均高生长模型进行适用性检验，得出F_(7，43)=18.236＞F_0.01(7，43)=3.087，即7个立地因子与优势木年均高生长相关性是极紧密的，用该方程来评价西双版纳柚木人工林立地质量是可靠的。

根据样地立地因子回归模型，各样地回归预测得分，以近似均分为原则划分4个生产力等级组，每生产力等级根据样地内土壤质地、土壤养分及土壤密度特征划分立地类型，共11个立地类型(表 8)。根据所分各立地生产力等级立地特点，对其进行评述。各主要养分含量比较见图 1。

表 8 立地类型划分结果 Tab.8 Classification result of site type

表 8 立地类型划分结果 Tab.8 Classification result of site type 生产力等级 Site productivity groups 立地类型Site type 优势高年均 生长量预测 Mean annual increment of dominant height predicted/m 优势高年均 生长量实测 Mean annual increment of dominant height measured/m 样地 Sample plot 高产组Ⅰ High production groupⅠ Ⅰ1 粉砂壤Silt loam， 高中等养分High & middle nutrient 1.832 1.814 1，4，23 Ⅰ2 中等养分Middle nutrient， 土壤密度Soil density >1.10 1.526 1.701 2，3 Ⅱ1 壤土、重壤土Loam & heavy Loam， 土壤密度Soil density>1.00 1.367 1.313 5，8，9，12，19，20，26， 32，36，37 中产组Ⅱ Middle production groupⅡ Ⅱ2 粉砂壤Silt loam， 1.00>土壤密度Soil density>0.90 1.351 1.406 18 Ⅱ3 砂壤Sandy loam 土壤密度Soil density ＞1.10 1.270 1.283 39，40 Ⅱ4 壤土、重壤土Loam & heavy loamg 1. 001＞土壤密度Soil density＞0.90 1.251 1.218 10，13，14，16，17，24， 25，29，34，41，47，51 Ⅲ1 壤土Loam， 1.10＞土壤密度Soil densit＞0.90 1.116 1.120 15，30，33，35，45，49 低产组Ⅲ Low production groupⅢ Ⅲ2 砂壤土Sandy loam， 1.10＞土壤密度 Soil density>1.00 1.106 1.078 31 Ⅲ3 重壤土Heavy loam， 土壤密度Soil density ＜1.00 1.072 1.107 28，38，44，46，48，50 Ⅲ4 壤土Loam， 土壤密度Soil density ＜0.90 0.949 1.120 11，22，27，42，43 劣产组Ⅳ Poor production groupⅣ Ⅳ1 重壤土Loam， 贫养分Poverty nutrient， 土壤密度Soil density ＜1.00 0.679 0.511 6，7，21

图 1 各生产力等级主要土壤理化性质 Fig.1 Main chemical and physical properties for each productivity group

高产组立地条件优越、土壤肥力高，物理性质优良、生产潜力最大，其中潜力最高的立地类型为Ⅰ₁，土壤质地为粉砂质壤土，砂黏比例适中，通气透水性好，适合柚木生长。本组立地土壤密度值较高，上层均值0.97 g·cm^-3，下层均值1.12 g·cm^-3，孔隙状况好。土壤pH最高为6.00，交换性酸含量22.25 cmol(H⁺)·kg^-1最低，交换性Al³⁺和有效Fe含量均最低，符合柚木在微酸至中性条件立地生长的特性。盐基饱和度78%，表明土壤肥力水平高。从养分元素来看，全K、有效P、交换性Ca^2＋、交换性Mg^2＋含量在所有生产力等级中最高。林木生长旺盛、整齐，年均高生长1.43 m、年均胸径生长1.51 cm、蓄积量年均增长达20.949 m³·hm^-2。本组因个体长势旺，须适时进行间伐，调整林分密度，并加强林分管理，立地生产潜力可发挥至最大，是培育柚木大径材最理想的立地类型。

中产组立地条件良好，土壤肥力较高，有良好的生产潜力。土壤质地大部分为壤土，少部分或为砂壤土、或为重壤土。上层和下层土壤密度均为1.01 g·cm^-3，孔隙状况适中，通气透水性良好。土壤pH均值5.32，交换性酸含量低，交换性Al^3＋和有效Fe含量低。土壤养分含量较高，全K、交换性Ca^2＋、交换性Mg^2＋和盐基饱和度高。本组立地条件较优，年均高生长1.02 m、年均胸径生长1.25 cm、年均蓄积量增长13.847 m³·hm^-2。本组在适宜经营管理下，可培育柚木中大径材。

低产组立地条件一般，肥力低，生产潜力一般。土壤质地多为壤土或重壤土。上层土壤密度为0.95 g·cm^-3，下层为0.94 g·cm^-3，孔隙状况较差，保水保肥能力较差。土壤酸性大，有效肥力受到抑制。养分含量较低，肥力较差，全K含量、有效P含量、盐基饱和度以及交换性Ca^2＋，Mg^2＋含量较低。林木生长中下，年均高生长0.94 m、年均胸径生长1.19 cm、年均蓄积量12.265 m³·hm^-2。本组立地若加强经营管理，仍可培育中小径材的柚木。

劣产组立地条件最差，土壤肥力最差。土壤质地类型均为重壤土，土壤密度低，总孔隙度大，土壤各养分含量均较低，肥力贫瘠，且土壤强酸性。林下植被白茅(Imperata cylindric)几乎覆盖整个样地，根系可深达60 cm，与柚木竞争严重影响柚木根系吸收养分。林木长势极差，尖削度大，部分成“小老头树”，年均高生长0.38 m、年均胸径生长0.53 cm、年均蓄积量增长0.930 m³·hm^-2。本组立地保留柚木没有发展前景，应及时进行改造，减少损失。

柚木对土壤的水分情况非常敏感，在通气透水性较好且排水性能优良的立地生长良好(Watanabe et al.，2010)。砂质土壤粒径较大，土壤骨架松散，粒间孔隙大，尽管通气透水性好，但保肥保水能力较差，养分含量少，不耐旱、不耐涝(Christian，2008)。而黏重土壤虽然大孔隙较少，但由于其小孔隙众多，单位体积的干土质量比较小，也可能表现出较小密度(Ziblim et al.，2012)。 本研究中，柚木人工林的土壤密度最大1.26 g·cm^-3，说明当初种植柚木选择立地时，已注意到土壤疏松度的问题，但密度偏低的立地对柚木生长影响较大(Rugmini et al.，2007)。

本研究所分各生产力等级中，劣产组的土壤密度均较低，但其质地却较黏重。黏重土壤的黏粒含量高，土壤骨架紧实黏结，粒间孔隙多而狭小，养分含量丰富，且黏质土本身能够吸引带正电荷的游离态养分，保肥保水性能优良，但是其孔隙往往被水占据，通气透水性能差，排水困难而雨水难以下渗，容易产生滞水或干旱情况(Anwar et al.，2007)。在泰国和利比里亚的柚木人工林研究中发现，湿低洼地和排水状况较差的立地容易造成顶梢枯死和叶黄化(Zech et al.，1991；Tanaka et al.，1998)。

柚木喜好在微酸至中性土壤中生长，在pH值6.5～7.5的土壤能自然生长并成为优势种(Kaosaard，1981)。当土壤的pH值处于中性附近范围时，pH值变化对柚木生长影响不显著(Watanabe et al.，2010； Amonsah et al.，2000)。但是当土壤酸度较高时，则会限制柚木生长(Rugmini et al.，2007)，如本研究中，高产组土壤pH值高，交换性酸低，而劣产组土壤显著酸许多，且相关性分析结果显示土壤pH值和交换酸的含量分别与柚木生长呈极显著正相关和极显著负相关。土壤酸度与土壤中的交换性Al^3＋含量关系很大(Sanchez，1977)。当土壤酸性较强时，无机离子态铝会对植物有毒害，减少植物Mg²⁺，Fe²⁺或Ca²⁺等阳离子吸收量，或与一些阴离子如PO₄^3-形成复杂络合物而使P变为无效P(Hartemink，2003)。在碱性或中性环境下，铝所形成的一系列聚合羟基形态或氢氧化铝固体对植物的毒性很小，但是当酸性较强时，铝将以Al³⁺，Al(OH)²⁺或Al(OH)₂⁺形式存在，而这些形态的铝对植物有很大毒害，影响植物根尖细胞伸长和细胞分裂。铝与磷酸根离子、羟基等极性基团结合，与质膜的脂质或脂蛋白结合影响膜的结构和功能，降低Mg²⁺-ATPase和K⁺-Mg²⁺-ATPase活性，破坏质膜结构，从而影响Mg²⁺，Fe²⁺或Ca²⁺等阳离子的吸收(孔繁翔等，2000)。本研究中，上层土壤pH值为4.76～6.47，下层为4.32～6.44，属于强酸性至弱酸性土壤，上下层土壤pH值与柚木生长极显著正相关，说明在一定程度下土壤pH值提高，有助于提高柚木林产量。因此，对强酸土壤可施石灰等碱性肥料，中和土壤酸性，提高交换性Mg²⁺和速效P的含量，促进柚木生长(潘一峰等，1997)。

土壤中除了C和N外，K是吸收最多的大量元素。K元素在植物体中有重要作用，它在渗透压、离子调节、氨基酸和蛋白质的合成以及抵消体内过量N毒害中必不可少。本研究区以砖红壤为主，其供K能力极弱，K含量的高低制约着柚木生长(Fernández-Moya et al.，2015)。从相关分析结果可知，上下2层土壤全K含量与柚木优势木年高生长量均达极显著正相关，因此，柚木造林适当追施钾肥，可提高生长量。

Mg^2＋是植物进行光合作用的重要元素，此外它对P也很重要，Mg²⁺对植物体内P的运输和叶片及茎中C的移动起重要作用(McAlister et al.，1998)。本区土壤上下层交换性Mg²⁺与柚木的生长呈极显著正相关，说明Mg²⁺对柚木生长的重要性。

Attiwill et al.，1993 土壤中P元素对柚木生长至关重要(Imoro et al.，2012； Fernández-Moya et al.，2015)。在热带森林中，P可能是主要的养分限制因子(Attiwill et al.，1993； Swaine，1996)，土壤酸性较大，影响土壤中P的有效性，进而影响植物根系生长。Drechsel等(1994)研究发现，P与柚木的生长呈正相关关系。本研究区土壤有效P含量普遍较低，尽管土壤中全P与柚木生长的相关性不显著，但是上层土壤的有效P含量与柚木生长极显著正相关。梁坤南等(2002)的柚木幼林施肥试验表明，高磷水平的柚木生长是低磷水平的2.24倍。

盐基饱和度是判断土壤肥力情况的重要指标，且其值与pH值关系紧密，通常盐基饱和度大于80%的土壤肥力高，50%～80%的土壤肥力中等，而小于50%的土壤不肥沃(黄昌勇等，2010)。盐基饱和的土壤，交换性Ca^2＋，Mg^2＋等盐基离子含量较高，而交换性Al^3＋、交换性酸的含量较低。在本区所分4个生产力等级中，高产组盐基饱和度最高，上层均值80%，下层均值70%，而劣产组盐基饱和度远低于50%，土壤肥力贫瘠。

粉砂质土壤的砂黏适中，大小孔隙比例适当，柚木生长最优，而在黏重土壤和砂质土壤生长较差。土壤密度值越小对柚木生长越不利，本研究区土壤密度1.10～1.26 g·cm^-3的立地柚木生长最佳。

柚木在本区易因土壤酸性过强而导致生长不良，强酸性易造成Al毒害和影响Mg，Fe，Ca，P等元素的有效性及吸收。在选择种植柚木立地时，应选择微酸至中性的土壤，避免因土壤过酸造成生长不良。土壤pH值低很可能是本区柚木生长的限制因子，这也有待进一步研究。 本研究区土壤的全K、交换性Mg^2＋、有效P等的含量对柚木生长至关重要，在柚木生长过程中应根据立地实际情况适当补充。此外，盐基饱和度也是本区立地条件优良的重要指标。

按照以上所述立地要求，选择合适的造林地，是柚木人工林丰产的基础。如不能合理选择立地而盲目种植，则其经济和生态效益将会大打折扣。虽然柚木经营时间较长，但选择通气透水性较好、高盐基、富K，低Fe、低Al的微酸至中性立地，采用完善的经营管理措施，可培育出柚木大径材，同时合理利用小径材，会使柚木在云南西双版纳地区的发展前景更广阔。