全 文 ：中亚热带南酸枣落叶阔叶林土壤磷素
空间变异及其影响因素
胡瑞彬１ 方　 晰１，２，３∗ 项文化１，２，３ 蒋　 芳１ 雷丕锋１，２，３ 赵丽娟１，２ 朱文娟１ 邓湘雯１，２，３
（ １中南林业科技大学生命科学与技术学院， 长沙 ４１０００４； ２南方林业生态应用技术国家工程实验室， 长沙 ４１０００４； ３湖南会同
杉木林生态系统国家野外科学观测研究站， 湖南会同 ４３８１０７）
摘　 要　 在湖南省长沙县大山冲国有林场南酸枣落叶阔叶林建立 １ ｈｍ２固定样地，基于植物
群落学调查数据和等距离网格布点取样的土壤养分测定数据，采用地统计学和 ＧＩＳ 相结合的
方法，研究土壤 Ｐ 含量的空间变异特征及其影响因子．结果表明： 腐殖质层及 ０ ～ １０、１０ ～ ２０、
２０～３０ ｃｍ土层的全 Ｐ、有效 Ｐ 含量的均为中等强度变异，同一土层有效 Ｐ 含量的变异程度较
全 Ｐ 含量高，随着土壤深度增加，全 Ｐ、有效 Ｐ 含量下降，变异程度增大．土壤全 Ｐ、有效 Ｐ 含量
具有较高的空间自相关性，主要受结构性因素影响．有效 Ｐ 含量的空间异质性较全 Ｐ 强，各土
层全 Ｐ、有效 Ｐ 含量空间自相关范围分别为 ９２．８０ ～ １６８．９０ ｍ、７９．４３ ～ １０６．２０ ｍ．同一土层全 Ｐ
含量分维数高于有效 Ｐ，空间格局较有效 Ｐ 复杂，有效 Ｐ 含量空间依赖性更强，具有更好的结
构性．腐殖质层及 ０～１０、１０～２０、２０～３０ ｃｍ土层的全 Ｐ、有效 Ｐ 含量呈现条带状和斑块状梯度
性分布，且全 Ｐ、有效 Ｐ 含量的空间变化趋势基本一致，高值出现在洼地，低值则出现在山脊．
全 Ｐ、有效 Ｐ 含量与高程、地表凋落物现存量呈显著负相关，与凹凸度、树种种数、株数、土壤
ｐＨ相关性不显著，与土壤有机 Ｃ、全 Ｎ含量呈显著正相关．这反映出土壤 Ｐ 的淋溶特性，其空
间变异性受到多种因子相互作用的影响．
关键词　 南酸枣落叶阔叶林； 土壤磷； 空间异质性； 地统计学； ＧＩＳ技术
本文由国家林业公益性行业科研专项项目（２０１３０４３１７）、科技部国际合作项目（２０１３ＤＦＡ３２１９０）和教育部高等学校博士学科点专项科研基金项
目（２０１２４３２１１１０００６）资助 Ｔｈｉｓ ｗｏｒｋ ｗａｓ ｓｕｐｐｏｒｔｅｄ ｂｙ ｔｈｅ Ｐｒｏｇｒａｍ ｏｆ Ｓｔａｔｅ Ｆｏｒｅｓｔｒｙ Ｓｐｅｃｉａｌ Ｆｕｎｄ ｆｏｒ Ｐｕｂｌｉｃ Ｗｅｌｆａｒｅ Ｉｎｄｕｓｔｒｙ ｏｆ Ｃｈｉｎａ （２０１３０４３１７），
ｔｈｅ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｏｐｅｒａｔｉｏｎ Ｐｒｏｊｅｃｔ ｏｆ Ｍｉｎｉｓｔｒｙ ｏｆ Ｓｃｉｅｎｃｅ ａｎｄ Ｔｅｃｈｎｏｌｏｇｙ ｏｆ Ｃｈｉｎａ （２０１３ＤＦＡ３２１９０）， ａｎｄ ｔｈｅ Ｒｅｓｅａｒｃｈ Ｆｕｎｄ ｆｏｒ ｔｈｅ Ｄｏｃｔｏｒａｌ Ｐｒｏｇｒａｍ
ｏｆ Ｈｉｇｈｅｒ Ｅｄｕｃａｔｉｏｎ ｏｆ Ｃｈｉｎａ （２０１２４３２１１１０００６）．
２０１５⁃０８⁃２４ Ｒｅｃｅｉｖｅｄ， ２０１５⁃１２⁃２４ Ａｃｃｅｐｔｅｄ．
∗通讯作者 Ｃｏｒｒｅｓｐｏｎｄｉｎｇ ａｕｔｈｏｒ． Ｅ⁃ｍａｉｌ： ｆａｎｇｘｉｚｈａｎｇ＠ ｓｉｎａ．ｃｏｍ
Ｓｐａｔｉａｌ ｈｅｔｅｒｏｇｅｎｅｉｔｙ ａｎｄ ｉｎｆｌｕｅｎｃｉｎｇ ｆａｃｔｏｒｓ ｏｆ ｓｏｉｌ ｐｈｏｓｐｈｏｒｕｓ ｃｏｎｃｅｎｔｒａｔｉｏｎ ｉｎ ａ ｍｉｄ⁃ｓｕｂ⁃
ｔｒｏｐｉｃａｌ Ｃｈｏｅｒｏｓｐｏｎｄｉａｓ ａｘｉｌｌａｒｉｓ ｄｅｃｉｄｕｏｕｓ ｂｒｏａｄ⁃ｌｅａｖｅｄ ｆｏｒｅｓｔ， Ｃｈｉｎａ． ＨＵ Ｒｕｉ⁃ｂｉｎ１， ＦＡＮＧ
Ｘｉ１，２，３∗， ＸＩＡＮＧ Ｗｅｎ⁃ｈｕａ１，２，３， ＪＩＡＮＧ Ｆａｎｇ１， ＬＥＩ Ｐｉ⁃ｆｅｎｇ１，２，３， ＺＨＡＯ Ｌｉ⁃ｊｕａｎ１，２， ＺＨＵ Ｗｅｎ⁃
ｊｕａｎ１， ＤＥＮＧ Ｘｉａｎｇ⁃ｗｅｎ１，２，３ （ １Ｃｏｌｌｅｇｅ ｏｆ Ｌｉｆｅ Ｓｃｉｅｎｃｅ ａｎｄ Ｔｅｃｈｎｏｌｏｇｙ， Ｃｅｎｔｒａｌ Ｓｏｕｔｈ Ｕｎｉｖｅｒｓｉｔｙ ｏｆ
Ｆｏｒｅｓｔｒｙ ａｎｄ Ｔｅｃｈｎｏｌｏｇｙ， Ｃｈａｎｇｓｈａ ４１０００４， Ｃｈｉｎａ； ２Ｎａｔｉｏｎａｌ Ｅｎｇｉｎｅｅｒｉｎｇ Ｌａｂｏｒａｔｏｒｙ ｆｏｒ Ａｐｐｌｉｅｄ
Ｆｏｒｅｓｔ Ｅｃｏｌｏｇｉｃａｌ Ｔｅｃｈｎｏｌｏｇｙ ｉｎ Ｓｏｕｔｈｅｒｎ Ｃｈｉｎａ， Ｃｈａｎｇｓｈａ ４１０００４， Ｃｈｉｎａ； ３Ｈｕｉｔｏｎｇ Ｎａｔｉｏｎａｌ
Ｆｉｅｌｄ Ｓｔａｔｉｏｎ ｆｏｒ Ｓｃｉｅｎｔｉｆｉｃ Ｏｂｓｅｒｖａｔｉｏｎ ａｎｄ Ｒｅｓｅａｒｃｈ ｏｆ Ｃｈｉｎｅｓｅ Ｆｉｒ Ｐｌａｎｔａｔｉｏｎ Ｅｃｏｓｙｓｔｅｍ ｉｎ Ｈｕｎａｎ
Ｐｒｏｖｉｎｃｅ， Ｈｕｉｔｏｎｇ ４３８１０７， Ｃｈｉｎａ） ．
Ａｂｓｔｒａｃｔ： Ｉｎ ｏｒｄｅｒ ｔｏ ｉｎｖｅｓｔｉｇａｔｅ ｓｐａｔｉａｌ ｖａｒｉａｔｉｏｎｓ ｉｎ ｓｏｉｌ ｐｈｏｓｐｈｏｒｕｓ （Ｐ） ｃｏｎｃｅｎｔｒａｔｉｏｎ ａｎｄ ｔｈｅ ｉｎ⁃
ｆｌｕｅｎｃｉｎｇ ｆａｃｔｏｒｓ， ｏｎｅ ｐｅｒｍａｎｅｎｔ ｐｌｏｔ ｏｆ １ ｈｍ２ ｗａｓ ｅｓｔａｂｌｉｓｈｅｄ ａｎｄ ｓｔａｎｄ ｓｔｒｕｃｔｕｒｅ ｗａｓ ｓｕｒｖｅｙｅｄ ｉｎ
Ｃｈｏｅｒｏｓｐｏｎｄｉａｓ ａｘｉｌｌａｒｉｅｓ ｄｅｃｉｄｕｏｕｓ ｂｒｏａｄｌｅａｖｅｄ ｆｏｒｅｓｔ ｉｎ Ｄａｓｈａｎｃｈｏｎｇ Ｆｏｒｅｓｔ Ｐａｒｋ ｉｎ Ｃｈａｎｇｓｈａ
Ｃｏｕｎｔｙ， Ｈｕｎａｎ Ｐｒｏｖｉｎｃｅ， Ｃｈｉｎａ． Ｓｏｉｌ ｓａｍｐｌｅｓ ｗｅｒｅ ｃｏｌｌｅｃｔｅｄ ｗｉｔｈ ｅｑｕｉｄｉｓｔａｎｔ ｇｒｉｄ ｐｏｉｎｔ ｓａｍｐｌｉｎｇ
ｍｅｔｈｏｄ ａｎｄ ｓｏｉｌ Ｐ ｃｏｎｃｅｎｔｒａｔｉｏｎ ａｎｄ ｉｔｓ ｓｐａｔｉａｌ ｖａｒｉａｔｉｏｎ ｗｅｒｅ ａｎａｌｙｚｅｄ ｂｙ ｕｓｉｎｇ ｇｅｏ⁃ｓｔａｔｉｓｔｉｃｓ ａｎｄ
ｇｅｏｇｒａｐｈｉｃａｌ ｉｎｆｏｒｍａｔｉｏｎ ｓｙｓｔｅｍ （ＧＩＳ） ｔｅｃｈｎｉｑｕｅｓ． Ｔｈｅ ｒｅｓｕｌｔｓ ｓｈｏｗｅｄ ｔｈａｔ ｔｈｅ ｖａｒｉａｔｉｏｎｓ ｏｆ ｔｏｔａｌ Ｐ
ａｎｄ ａｖａｉｌａｂｌｅ Ｐ ｃｏｎｃｅｎｔｒａｔｉｏｎｓ ｉｎ ｈｕｍｕｓ ｌａｙｅｒ ａｎｄ ｉｎ ｔｈｅ ｓｏｉｌ ｐｒｏｆｉｌｅ ａｔ ｄｅｐｔｈ ｏｆ ０－１０， １０－２０ ａｎｄ
２０－３０ ｃｍ ｗｅｒｅ ｍｏｄｅｒａｔｅ ａｎｄ ｔｈｅ ａｖａｉｌａｂｌｅ Ｐ ｓｈｏｗｅｄ ｈｉｇｈｅｒ ｖａｒｉａｂｉｌｉｔｙ ｉｎ ａ ｓｐｅｃｉｆｉｃ ｓｏｉｌ ｌａｙｅｒ ｃｏｍ⁃
ｐａｒｅｄ ｗｉｔｈ ｔｏｔａｌ Ｐ． Ｃｏｎｃｅｎｔｒａｔｉｏｎｓ ｏｆ ｔｏｔａｌ Ｐ ａｎｄ ａｖａｉｌａｂｌｅ Ｐ ｉｎ ｓｏｉｌ ｄｅｃｒｅａｓｅｄ， ｗｈｉｌｅ ｔｈｅ ｖａｒｉａｔｉｏｎｓ
ｉｎｃｒｅａｓｅｄ ｗｉｔｈ ｔｈｅ ｉｎｃｒｅａｓｅ ｉｎ ｓｏｉｌ ｄｅｐｔｈ． Ｔｈｅ ｔｏｔａｌ Ｐ ａｎｄ ａｖａｉｌａｂｌｅ Ｐ ｓｈｏｗｅｄ ｈｉｇｈ ｓｐａｔｉａｌ ａｕｔｏｃｏｒｒｅ⁃
应 用 生 态 学 报　 ２０１６年 ３月　 第 ２７卷　 第 ３期　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 ｈｔｔｐ： ／ ／ ｗｗｗ．ｃｊａｅ．ｎｅｔ
Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ， Ｍａｒ． ２０１６， ２７（３）： ７０５－７１５　 　 　 　 　 　 　 　 　 　 ＤＯＩ： １０．１３２８７ ／ ｊ．１００１－９３３２．２０１６０３．０１３
ｌａｔｉｏｎ， ｐｒｉｍａｒｉｌｙ ｒｅｓｕｌｔｅｄ ｆｒｏｍ ｔｈｅ ｓｔｒｕｃｔｕｒａｌ ｆａｃｔｏｒｓ． Ｔｈｅ ｓｐａｔｉａｌ ｈｅｔｅｒｏｇｅｎｅｉｔｙ ｏｆ ａｖａｉｌａｂｌｅ Ｐ ｗａｓ
ｓｔｒｏｎｇｅｒ ｔｈａｎ ｔｈａｔ ｏｆ ｔｏｔａｌ Ｐ， ａｎｄ ｔｈｅ ｓｐａｔｉａｌ ａｕｔｏｃｏｒｒｅｌａｔｉｏｎ ｒａｎｇｅｓ ｏｆ ｔｏｔａｌ Ｐ ａｎｄ ａｖａｉｌａｂｌｅ Ｐ ｖａｒｉｅｄ
ｆｒｏｍ ９２．８０ ｔｏ １６８．９０ ｍ ａｎｄ ｆｒｏｍ ７９．４３ ｔｏ １０６．２０ ｍ ｉｎ ｄｉｆｆｅｒｅｎｔ ｓｏｉｌ ｌａｙｅｒｓ， ｒｅｓｐｅｃｔｉｖｅｌｙ． Ａｔ ｔｈｅ
ｓａｍｅ ｓｏｉｌ ｄｅｐｔｈ， ｆｒａｃｔａｌ ｄｉｍｅｎｓｉｏｎｓ ｏｆ ｔｏｔａｌ Ｐ ｗｅｒｅ ｈｉｇｈｅｒ ｔｈａｎ ｔｈａｔ ｏｆ ａｖａｉｌａｂｌｅ Ｐ， ｗｉｔｈ ｍｏｒｅ ｃｏｍ⁃
ｐｌｅｘ ｓｐａｔｉａｌ ｐａｔｔｅｒｎ， ｗｈｉｌｅ ａｖａｉｌａｂｌｅ Ｐ ｓｈｏｗｅｄ ｓｔｒｏｎｇｅｒ ｓｐａｔｉａｌ ｃｏｒｒｅｌａｔｉｏｎ ｗｉｔｈ ｓｔｒｏｎｇｅｒ ｓｐａｔｉａｌ ｓｔｒｕｃ⁃
ｔｕｒｅ． Ｉｎ ｈｕｍｕｓ ｌａｙｅｒ ａｎｄ ｓｏｉｌ ｄｅｐｔｈｓ ｏｆ ０－１０， １０－２０ ａｎｄ ２０－３０ ｃｍ， ｔｈｅ ｓｐａｔｉａｌ ｖａｒｉａｔｉｏｎ ｐａｔｔｅｒｎ ｏｆ
ｔｏｔａｌ Ｐ ａｎｄ ａｖａｉｌａｂｌｅ Ｐ ｃｏｎｃｅｎｔｒａｔｉｏｎｓ ｓｈｏｗｅｄ ａｎ ａｐｐａｒｅｎｔ ｂｅｌｔ⁃ｓｈａｐｅｄ ａｎｄ ｓｐｏｔ ｍａｓｓｉｖｅ ｇｒａｄｉｅｎｔ
ｃｈａｎｇｅ． Ｔｈｅ ｈｉｇｈ ｖａｌｕｅ ａｐｐｅａｒｅｄ ａｔ ｌｏｗ ｅｌｅｖａｔｉｏｎ ａｎｄ ｖａｌｌｅｙ ｐｏｓｉｔｉｏｎ， ａｎｄ ｔｈｅ ｌｏｗ ｖａｌｕｅ ａｐｐｅａｒｅｄ ｉｎ
ｔｈｅ ｈｉｇｈ ｅｌｅｖａｔｉｏｎ ａｎｄ ｒｉｄｇｅ ａｒｅａ． Ｔｈｅ ｔｏｔａｌ Ｐ ａｎｄ ａｖａｉｌａｂｌｅ Ｐ ｃｏｎｃｅｎｔｒａｔｉｏｎｓ ｓｈｏｗｅｄ ｓｉｇｎｉｆｉｃａｎｔｌｙ
ｎｅｇａｔｉｖｅ ｃｏｒｒｅｌａｔｉｏｎ ｗｉｔｈ ｅｌｅｖａｔｉｏｎ ａｎｄ ｌｉｔｔｅｒ， ｂｕｔ ｔｈｅ ｒｅｌａｔｉｏｎｓｈｉｐ ｗｉｔｈ ｃｏｎｖｅｘｉｔｙ， ｓｐｅｃｉｅｓ， ｎｕｍｂｅｒｓ
ａｎｄ ｓｏｉｌ ｐＨ ｗａｓ ｎｏｔ ｓｉｇｎｉｆｉｃａｎｔ． Ｔｈｅ ｔｏｔａｌ Ｐ ａｎｄ ａｖａｉｌａｂｌｅ Ｐ ｅｘｈｉｂｉｔｅｄ ｓｉｇｎｉｆｉｃａｎｔ ｐｏｓｉｔｉｖｅ ｃｏｒｒｅｌａｔｉｏｎｓ
ｗｉｔｈ ｓｏｉｌ ｏｒｇａｎｉｃ ｃａｒｂｏｎ （ＳＯＣ）， ｔｏｔａｌ ｎｉｔｒｏｇｅｎ ｃｏｎｃｅｎｔｒａｔｉｏｎ， ｉｎｄｉｃａｔｉｎｇ ｔｈｅ ｌｅａｃｈｉｎｇ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ
ｏｆ ｓｏｉｌ Ｐ． Ｉｔｓ ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｗａｓ ａｆｆｅｃｔｅｄ ｂｙ ｍａｎｙ ｉｎｔｅｒａｃｔｉｖｅ ｆａｃｔｏｒｓ．
Ｋｅｙ ｗｏｒｄｓ： Ｃｈｏｅｒｏｓｐｏｎｄｉａｓ ａｘｉｌｌａｒｉｓ ｄｅｃｉｄｕｏｕｓ ｂｒｏａｄｌｅａｖｅｄ ｆｏｒｅｓｔ； ｓｏｉｌ ｐｈｏｓｐｈｏｒｕｓ； ｓｐａｔｉａｌ ｈｅｔｅｒｏ⁃
ｇｅｎｅｉｔｙ； ｇｅｏ⁃ｓａｔｉｓｔｉｃｓ； ＧＩＳ ｔｅｃｈｎｉｑｕｅ．
　 　 土壤养分空间异质性直接调控植物群落组成、
植被分布及其生物量格局，对生态系统演替有重要
作用，探讨森林土壤养分空间异质性及其影响因素
不仅有助于对土壤发育格局的认识［１］，对揭示土壤
与植物空间分布的关系、生物多样性维持和共存机
制具有重要的理论意义．因此，森林土壤养分空间变
异特征成为林学、生态学、土壤学的研究热点［２］ ．有
研究表明，森林土壤养分空间变异受地形、生物和气
候等结构性因素以及取样、测量误差等随机性因素
的共同作用［３－５］ ．但至今对小尺度森林群落土壤养
分、地形因子与植物分布三者之间复杂关系的认识
仍非常有限［２，６］ ．
磷（Ｐ）是植物生长发育不可缺少的营养元素，
以多种方式参与植物体内各种生理生化过程，是地
质时代尺度上植物生产力的限制性养分元素．土壤
Ｐ 成为了限制当前农林业产量的主要因素［７－８］，控
制着陆地生态系统的生态过程［９］ ．研究表明，土壤
Ｎ、Ｐ 供应存在明显的耦合过程［１０］，土壤缺 Ｐ 不仅限
制植物生长发育和生产力，对生态系统的其他组分
如微生物和 Ｃ、Ｎ矿化等生态过程也起限制作用，土
壤 Ｐ 的变化影响着包括 Ｃ、Ｎ 循环在内的其他元素
的生物地球化学循环［１１］ ．为此，研究森林土壤 Ｐ 的
空间分布格局及其影响因素，对认识森林土壤结构
与功能及其对植被的影响，调控和管理森林生态系
统 Ｐ 的生物地球化学循环具有重要意义［１２］ ．长期以
来，有关土壤 Ｐ 的研究多集中在农田、湿地、水域生
态系统，主要侧重于土壤 Ｐ 循环、转化，土壤 Ｐ 分级
及其有效性［１３］，土壤对 Ｐ 吸附与释放特性以及磷肥
施用效应等方面．２０ 世纪 ９０ 年代初，人们开始注意
到森林土壤 Ｐ 的活化、林木对 Ｐ 吸收利用以及林地
土壤 Ｐ 状况与农田土壤明显不同，陆续开展了低 Ｐ
胁迫对落叶松（Ｌａｒｉｘ ｇｍｅｌｉｎｉｉ）、马尾松（Ｐｉｎｕｓ ｍａｓｓｏ⁃
ｎａｎａ）等林木生长影响的研究［１４－１５］以及一些施肥试
验．结果表明，在许多情况下，Ｐ 似乎比 Ｎ、Ｋ 更为重
要，土壤 Ｐ 可能是制约我国南方人工林速生丰产，
影响地力、生态系统平衡和稳定的一个重要因
素［１６］ ．然而，目前对森林土壤 Ｐ 空间分布特征的研
究仍较少，小尺度森林群落土壤 Ｐ 空间异质性的研
究更少．土壤 Ｐ 的空间分布格局及其主导影响因素
不仅随着研究尺度的不同而异［１７］，即使同一研究尺
度的研究结果也不尽一致．张世熔等［１８］研究表明，
土壤 Ｐ 的空间分异主要受到随机因素（施肥）影响，
结构性因素的影响较小；而苏松锦等［５］和张娜等［１９］
研究发现，土壤 Ｐ 表现出强烈的空间自相关性，其
空间异质性主要受到结构性因素影响．秦松等［２０］认
为，土壤 Ｐ 空间分布与海拔呈显著正相关，但张娜
等［１９］和樊纲惟等［２１］研究得出，土壤 Ｐ 空间分布与
海拔、凹凸度呈显著负相关．
我国中亚热带地区森林类型多样，树种组成复
杂．南酸枣（Ｃｈｏｅｒｏｓｐｏｎｄｉａ ａｘｉｌｌａｒｉｓ）落叶阔叶林是我
国中亚热带低山丘陵区典型的次生林类型之一．国
内对南酸枣的研究主要集中在经济、药用价值上，对
其次生林生态系统结构和功能的研究很少．湖南省
长沙县大山冲林场保存着完整的南酸枣落叶阔叶次
生林，树种组成丰富，分布地形复杂，人为干扰
少［２２］，为土壤 Ｐ 的空间异质性研究提供了良好的场
所．本文在森林群落尺度上研究南酸枣落叶阔叶林
土壤 Ｐ 含量的空间变异特征，以及其与地形因子、
生物因子、土壤因子的关系，以期揭示亚热带森林土
壤 Ｐ 空间分布的调控机制，为亚热带森林健康经营
６０７ 应　 用　 生　 态　 学　 报　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 ２７卷
与管理提供科学依据．
１　 研究地区与研究方法
１􀆰 １　 研究区概况
研究区为湖南省长沙县大山冲国有林场
（２８°２３′—２８°２４′ Ｎ，１１３°１７′—１１３°１９′ Ｅ），为幕阜山
余脉西部边缘的湘中丘陵地区，海拔 ２４５～３２１ ｍ，属
中亚热带东南季风湿润气候，年均气温 １６．６ ～ １７．６
℃，极端最低气温－１１ ℃ （１ 月），极端最高气温 ４０
℃（７月），年降雨量 １４１２ ～ １５５９ ｍｍ，降雨主要集中
在 ５—７ 月．土壤为由板页岩和页岩发育而成的红
壤，地带性植被为亚热带常绿阔叶林，由于人为干
扰，原生植被已被破坏，经过 ５０多年的封山育林，现
保存有杉木（Ｃｕｎｎｉｎｇｈａｍｉａ ｌａｎｃｅｏｌａｔａ）、马尾松、毛
竹（Ｐｈｙｌｌｏｓｔａｃｈｙｓ ｈｅｔｅｒｏｃｙｃｌａ）等人工林群落，以及由
不同演替阶段树种组成的多种次生林．林场内次生
林人为干扰少．
１􀆰 ２　 样地设置
２０１３年 ３—４月，在南酸枣落叶阔叶林内选取
代表性地段，设置面积为 １ ｈｍ２（１００ ｍ×１００ ｍ）的固
定样地，涵盖了整个南酸枣落叶阔叶林，样地海拔
２５１～３１７ ｍ，为西北坡向，坡度为 ２２°，土壤 ｐＨ ４．１８～
５．０３．样地再分成 １００ 个 １０ ｍ×１０ ｍ 小样地，调查样
地植物组成特征，对胸径＞１ ｃｍ 的所有乔木植物挂
牌编号，记录植物种类，测定树高、冠幅、胸径等，计
算南酸枣落叶阔叶林主要树种密度、平均胸径、平均
树高、相对密度和重要值等林分特征指标［２３］，南酸
枣落叶阔叶林群落基本特征见表 １．
１􀆰 ３　 土壤样品的采集
２０１４年 ５—６月，分别在南酸枣落叶阔叶林 １００
个小样地的中心部位设置 １个面积为 ５０ ｃｍ×５０ ｃｍ
的小样方，测定地表凋落物现存量和采集土壤样品．
采样时，收集小样方内地表全部凋落物，采集小样方
内地表腐殖质层，最后挖土壤剖面，按 ０～ １０、１０ ～ ２０
和 ２０～３０ ｃｍ分层从下至上采集土壤样品（每个土
壤样品约 ３ ｋｇ）．采集地表凋落物样品 １００ 个、腐殖
质层样品 １００ 个、土壤样品 ３００ 个．土壤、腐殖质层
样品在室内风干，清除石头、根系等碎屑后分成两
份：一份研磨，过 ６０ 目土壤筛，用于化学分析；一份
用于测定土壤机械组成．凋落物样品于 ８５ ℃下烘干
至恒量，称量，计算出每个小样地单位面积地表的凋
落物现存量．
１􀆰 ４　 测定项目与方法
样地高程用罗盘仪测定，参照 ＣＴＦＳ 样地常用
的算法［２４］，首先从 １ ｈｍ２固定样地的 ４ 个顶点选出
最低点为基点，用罗盘仪测量 １００ 个 １０ ｍ×１０ ｍ 的
小样地两两角桩之间的倾角（即坡度），通过三角函
数计算出每个小样地 ４ 个顶点的高程，每个小样地
（１０ ｍ×１０ ｍ）的高程取该小样地 ４个顶点处的高程
平均值．
土壤含水率用 １０５ ℃烘干法测定，土壤、腐殖质
层样品全 Ｐ 含量用高氯酸⁃硫酸酸溶⁃钼锑抗比色法
测定，有效 Ｐ 含量用盐酸⁃氟化铵（Ｂｒａｙ⁃１）比色法测
定，ｐＨ用 ｐＨ 计法测定（水土比 ２ ∶ １），有机碳含量
用重铬酸钾⁃浓硫酸水合加热法测定，全 Ｎ含量用半
微量凯氏定氮法测定［２５］ ．
１􀆰 ５　 数据处理
用 Ｅｘｃｅｌ统计软件分析 １ ｈｍ２固定样地土壤、腐
殖质层全 Ｐ、有效 Ｐ 含量的最大值、最小值、平均值、
标准差和变异系数，用 ＳＰＳＳ 软件包中的单因素方
差分析（ｏｎｅ⁃ｗａｙ ＡＮＯＶＡ）分析不同土层全 Ｐ、有效
Ｐ含量的差异显著性（ α ＝ ０．０５） ，用Ｋ ⁃Ｓ检验对土
表 １　 南酸枣落叶阔叶林群落基本特征
Ｔａｂｌｅ １　 Ｃｏｍｍｕｎｉｔｙ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｔｈｅ Ｃｈｏｅｒｏｓｐｏｎｄｉａｓ ａｘｉｌｌａｒｉｓ ｄｅｃｉｄｕｏｕｓ ｂｒｏａｄｌｅａｖｅｄ ｆｏｒｅｓｔ
树种
Ｔｒｅｅ ｓｐｅｃｉｅｓ
密度
Ｄｅｎｓｉｔｙ
（ ｉｎｄ·ｈｍ－２）
平均胸径
Ａｖｅｒａｇｅ ＤＢＨ
（ｃｍ）
平均树高
Ａｖｅｒａｇｅ ｈｅｉｇｈｔ
（ｍ）
相对密度
Ｒｅｌａｔｉｖｅ
ｄｅｎｓｉｔｙ
重要值
Ｉｍｐｏｒｔａｎｃｅ
ｖａｌｕｅ
南酸枣 Ｃ． ａｘｉｌｌａｒｉｓ ２１４ ２０．４（１．０～５４．１） １４．１（２．０～３４．５） ４．０５ １９．９８
檵木 Ｌｏｒｏｐｅｔａｌｕｍ ｃｈｉｎｅｎｓｉｓ １３２５ ３．４（１．０～１４．８） ３．９（１．４～２１．５） ２５．０８ １３．４８
四川山矾 Ｓｙｍｐｌｏｃｏｓ ｓｅｔｃｈｕｅｎｓｉｓ ６１３ ３．１（１．０～２０．３） ３．５（１．４～８．９） １１．６０ ７．０７
满山红 Ｒｈｏｄｏｄｅｎｄｒｏｎ ｍａｒｉｅｓｉｉ ５６８ ２．４（１．０～１４．５） ３．２（１．３～７．５） １０．７５ ５．７３
格药柃 Ｅｕｒｙａ ｍｕｒｉｃａｔａ ４５１ ２．６（１．０～１４．７） ２．９（１．３～１２．０） ８．５４ ５．３５
木油桐 Ｖｅｒｎｉｃｉａ ｍｏｎｔａｎａ １３７ ５．７（１．０～４１．８） ６．５（１．５～２３．６） ２．５９ ４．１４
毛豹皮樟 Ｌｉｔｓｅａ ｃｏｒｅａｎａ ｖａｒ． ｓｉｎｅｎｓｉｓ １３６ ５．０（１．０～２９．３） ４．９（１．５～２１．５） ２．５７ ３．０７
其他 Ｏｔｈｅｒｓ １８３９ ０．４（１．０～３７．４） ９．３（１．３～２６．３） ３４．８２ ４１．１８
合计 Ｔｏｔａｌ ５２８３ ４．２（１．０～５４．１） ４．３（１．３～３４．５） １００ １００
括号内数据为胸径或树高的变化范围 Ｖａｌｕｅｓ ｉｎ ｐａｒｅｎｔｈｅｓｉｓ ｗｅｒｅ ｔｈｅ ｒａｎｇｅ ｏｆ ｔｒｅｅ ｄｉａｍｅｔｅｒ ａｔ ｂｒｅａｓｔ ｈｅｉｇｈｔ ｏｒ ｈｅｉｇｈｔ．
７０７３期　 　 　 　 　 　 　 　 　 　 胡瑞彬等： 中亚热带南酸枣落叶阔叶林土壤磷素空间变异及其影响因素　 　 　 　 　 　
壤、腐殖质层全 Ｐ、有效 Ｐ 含量进行正态性分布检
验，若 Ｐ＞０．０５，认为数据服从正态分布，对不能满足
正态分布的数据进行数据转换，使数据呈正态分布．
用阈值法检验特异值，即将处于样本均值加减
３倍标准差区间外的样点数据作为特异值［２６］ ．若数
值大于样本均值加上 ３ 倍标准差，则该特异值用同
一土层的正常最大值来代替；若数值小于样本减去
３ 倍标准差，则该特异值用同一土层的正常最小值
来代替．
运用半方差函数分析各土层全 Ｐ、有效 Ｐ 含量
的空间变异特征，通过半方差函数得到散点图，对散
点图采用球状模型、高斯模型、指数模型和线性模型
等理论模型进行拟合［２７］ ．模型拟合得到块金值
（Ｃ０）、基台值（Ｃ０＋Ｃ）和变程（Ａ０）３个评价空间变异
程度的重要参数．块金值与基台值之比（即结构比）
表示随机性因素引起的空间异质性占系统总变异的
比例［２８］ ．
半方差函数的计算公式：
γ（ｈ）＝ １
２Ｎ（ｈ）
∑
ｎ
ｉ＝１
［Ｚ（ｘｉ）－Ｚ（ｘｉ＋ｈ）］ ２
式中：ｒ（ｈ）为半方差函数；ｈ 为样本的间隔距离；Ｎ
（ｈ）是间距为 ｈ的观测样点的成对数；Ｚ（ｘｉ）为系统
某属性 Ｚ 在空间位置 ｘｉ处的值，Ｚ（ｘｉ＋ｈ）是（ｘｉ ＋ｈ）
处值的一个区域化变量［２９］ ．由半方差函数中 ｈ 和 ｒ
（ｈ）的回归曲线确定土壤养分空间分维数（Ｄ），表
征土壤养分空间异质性程度．
用 ＧＳ＋软件中的空间模块进行空间分布分析，
用 Ｋｒｉｇｉｎｇ插值法对固定样地的土壤全 Ｐ、有效 Ｐ 含
量进行空间插值，得到南酸枣落叶阔叶林固定样地
土壤全 Ｐ、有效 Ｐ 含量的空间分布特征图．
以高程、凹凸度为地形因子，分析地形因子对土
壤 Ｐ 空间分布的影响．凹凸度为小样地高程减去该
样地相邻的 ８个小样地高程的平均值，处于样地边
缘的小样地凹凸度为小样地中心的高程减去 ４个顶
点高程的平均值．若凹凸度为正值，说明该小样地比
周围小样地高；若为负值，说明该小样地高程比周围
小样地高程低［３０］ ．用 ＳＰＳＳ 软件包中 Ｐｅａｒｓｏｎ 相关分
析计算土壤层全 Ｐ、有效 Ｐ 含量与地形因子（高程、
凹凸度）、生物因子（地表凋落物现存量、物种数、株
数）、土壤因子（ｐＨ、有机碳、全 Ｎ）的相关系数．
２　 结果与分析
２􀆰 １　 土壤全 Ｐ、有效 Ｐ 含量的描述性统计特征
表 ２表明，腐殖质层及 ０～１０、１０～２０、２０～３０ ｃｍ
土层全 Ｐ 含量平均值分别为 ０．３６、０．２９、０．２７、０．２３
ｇ·ｋｇ－１，变化范围分别为 ０． ２２ ～ ０． ５４、０． １４ ～ ０． ５４、
０􀆰 １２～０．４６、０． １３ ～ ０． ４８ ｇ·ｋｇ－１，变异系数分别为
１９􀆰 １％、２３．４％、２８．２％、２８．８％，随着土层深度增加，
全 Ｐ 含量下降，变异程度增大，腐殖质层与其他各
土层之间差异显著，但 ０ ～ １０、１０ ～ ２０、２０ ～ ３０ ｃｍ 土
层间差异不显著．腐殖质层及 ０ ～ １０、１０ ～ ２０、２０ ～ ３０
ｃｍ土层有效 Ｐ 含量平均值分别为 １２．１３、２．４６、２．６３、
１．９０ ｍｇ·ｋｇ－１，变化范围分别在 ４．８８ ～ ２２．２９、０．１４ ～
０．５４、０．１２～０．４６、０．１３ ～ ０．４８ ｍｇ·ｋｇ－１，变异系数分
别为 ３１􀆰 ４％、３５．４％、４０．７％、４５．０％，随着土层深度
增加，有效 Ｐ 含量下降，变异程度增大，腐殖质层与
其他各土层差异显著，０～１０、１０～２０、２０～３０ ｃｍ土层
间差异不显著．腐殖质层及 ０～１０、１０～２０、２０～３０ ｃｍ
土层全 Ｐ、有效 Ｐ 含量的变异系数分别为 １９．１％ ～
２８．８％和 ３１．４％～４５．０％，属于中等程度变异，同一土
层有效 Ｐ 含量变异幅度高于全 Ｐ 含量．
阈值法检验结果表明，各土层全 Ｐ、有效 Ｐ 含量
的特异值均大于样本均值加上 ３倍标准差，其中，腐
殖质层全 Ｐ、有效 Ｐ 含量没有特异值；０～１０ ｃｍ土层
表 ２　 土壤全 Ｐ和有效 Ｐ含量的描述性统计
Ｔａｂｌｅ ２　 Ｄｅｓｃｒｉｐｔｉｖｅ ｓｔａｔｉｓｔｉｃｓ ｏｆ ｔｏｔａｌ ｐｈｏｓｐｈｏｒｕｓ ａｎｄ ａｖａｉｌａｂｌｅ ｐｈｏｓｐｈｏｒｕｓ ｃｏｎｃｅｎｔｒａｔｉｏｎ ｉｎ ｓｏｉｌ
项目
Ｉｔｅｍ
土层
Ｓｏｉｌ ｄｅｐｔｈ
（ｃｍ）
最小值
Ｍｉｎｉｍｕｍ
最大值
Ｍａｘｉｍｕｍ
平均值
Ａｖｅｒａｇｅ
标准差
ＳＤ
变异系数
ＣＶ
（％）
Ｋ⁃Ｓ检验 Ｐ值
Ｐ ｖａｌｕｅ ｏｆ
Ｋ⁃Ｓ ｔｅｓｔ
转换后的 Ｐ值
Ｐ ｖａｌｕｅ ａｆｔｅｒ
ｔｒａｎｓｆｏｒｍａｔｉｏｎ
ＴＰ Ｈ ０．２２ ０．５４ ０．３６ａ ０．０７ １９．１ ０．７０
（ｇ·ｋｇ－１） ０～１０ ０．１４ ０．５４ ０．２９ｂ ０．０７ ２４．０ ０．４７
１０～２０ ０．１２ ０．４６ ０．２７ｂ ０．０８ ２８．２ ０．８９
２０～３０ ０．１３ ０．４８ ０．２３ｂ ０．０７ ２８．８ ０．３０
ＡＰ Ｈ ４．８８ ２２．２９ １２．１３ａ ３．８０ ３１．４ ０．７５
（ｍｇ·ｋｇ－１） ０～１０ ０．９１ ５．５０ ２．４６ｂ ０．８７ ３５．４ ０．５９
１０～２０ １．２８ ６．６８ ２．６３ｂ １．０７ ４０．７ ０．０５
２０～３０ １．０２ ５．３４ １．９０ｂ ０．８６ ４５．０ ０．００ ０．４０５
ＴＰ： 全 Ｐ Ｔｏｔａｌ Ｐ； ＡＰ：有效 Ｐ Ａｖａｉｌａｂｌｅ Ｐ； Ｈ：腐殖质层Ｈｕｍｕｓ ｈｏｒｉｚｏｎ．下同 Ｔｈｅ ｓａｍｅ ｂｅｌｏｗ．不同字母表示不同土层间差异显著（Ｐ＜０．０５） Ｄｉｆｆｅｒ⁃
ｅｎｔ ｌｅｔｔｅｒｓ ｍｅａｎｔ ｓｉｇｎｉｆｉｃａｎｔ ｄｉｆｆｅｒｅｎｃｅ ａｍｏｎｇ ｄｉｆｆｅｒｅｎｔ ｓｏｉｌ ｌａｙｅｒｓ ａｔ ０．０５ ｌｅｖｅｌ．
８０７ 应　 用　 生　 态　 学　 报　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 ２７卷
全 Ｐ 含量有一个特异值，用正常最大值 ０．４８ ｇ·ｋｇ－１
代替，有效 Ｐ 含量有 ２个特异值，用正常最大值 ４．３０
ｍｇ·ｋｇ－１代替；１０ ～ ２０ ｃｍ 土层全 Ｐ 含量没有特异
值，有效 Ｐ 含量有 １ 个特异值，用正常最大值 ５．５８
ｍｇ·ｋｇ－１代替；２０～３０ ｃｍ土层全 Ｐ 含量有 １个特异
值，用正常最大值 ０．４０ ｇ·ｋｇ－１代替，有效 Ｐ 含量有
２个特异值，用正常最大值 ４．２８ ｍｇ·ｋｇ－１代替．Ｋ⁃Ｓ
检验结果（表 ２）表明，各土层全 Ｐ 含量均服从正态
分布，腐殖质层及 ０～１０、１０～２０ ｃｍ土层有效 Ｐ 含量
服从正态分布，２０ ～ ３０ ｃｍ 土层有效 Ｐ 含量不服从
正态分布，经倒数转换后服从正态分布，能满足地统
计学分析要求．
２􀆰 ２　 土壤全 Ｐ、有效 Ｐ 含量的空间变异及结构特征
样本间距 ｈ 的变化范围取最大间距的 １ ／ ２，为
６３．６４ ｍ，步长取最小间距，为 １０ ｍ．从图 １ 和表 ３ 可
以看出，腐殖质层、０ ～ １０ ｃｍ 土层全 Ｐ 含量的最佳
拟合模型为指数模型，１０ ～ ２０ ｃｍ 土层为球状模型，
２０～３０ ｃｍ土层为高斯模型，Ｒ２ ＝ ０．８８ ～ ０．９３，均达到
显著水平．腐殖质层、２０ ～ ３０ ｃｍ 土层有效 Ｐ 含量的
最佳拟合模型为球状模型，０～１０、１０～２０ ｃｍ 土层为
指数模型，Ｒ２ ＝ ０．８７ ～ ０．９９，均达到显著水平．表明用
半方差函数理论模型能较好地反映南酸枣落叶阔叶
图 １　 土壤全 Ｐ（ＴＰ）和有效 Ｐ（ＡＰ）含量的半方差函数理论模型
Ｆｉｇ．１　 Ｓｅｍｉｖａｒｉｏｇｒａｍ ｔｈｅｏｒｅｔｉｃａｌ ｍｏｄｅｌｓ ｏｆ ｔｏｔａｌ ｐｈｏｓｐｈｏｒｕｓ （ＴＰ） ａｎｄ ａｖａｉｌａｂｌｅ ｐｈｏｓｐｈｏｒｕｓ （ＡＰ） ｃｏｎｃｅｎｔｒａｔｉｏｎ ｉｎ ｓｏｉｌ．
ａ） 腐殖质层 Ｈｕｍｕｓ ｈｏｒｉｚｏｎ； ｂ） ０～１０ ｃｍ； ｃ） １０～２０ ｃｍ； ｄ） ２０～３０ ｃｍ．下同 Ｔｈｅ ｓａｍｅ ｂｅｌｏｗ．
９０７３期　 　 　 　 　 　 　 　 　 　 胡瑞彬等： 中亚热带南酸枣落叶阔叶林土壤磷素空间变异及其影响因素　 　 　 　 　 　
表 ３　 土壤全 Ｐ和有效 Ｐ半方差函数的模型类型及参数
Ｔａｂｌｅ ３　 Ｓｅｍｉｖａｒｉｏｇｒａｍ ｔｈｅｏｒｅｔｉｃａｌ ｍｏｄｅｌｓ ａｎｄ ｐａｒａｍｅｔｅｒｓ ｆｏｒ ｔｏｔａｌ ｐｈｏｓｐｈｏｒｕｓ ａｎｄ ａｖａｉｌａｂｌｅ ｐｈｏｓｐｈｏｒｕｓ ｉｎ ｓｏｉｌ
项目
Ｉｔｅｍ
土层
Ｓｏｉｌ ｄｅｐｔｈ
（ｃｍ）
模型
Ｍｏｄｅｌ
Ｃ０ Ｃ０＋Ｃ 结构比
ＮＳＲ
（％）
Ａ０
（ｍ）
Ｒ２ ＲＳＳ Ｄ
ＴＰ Ｈ Ｅ ０．００１９ ０．００４６ ４１．３ １６８．９０ ０．９３∗∗ ３．３１８×１０－８ １．９０３
０～１０ Ｅ ０．００２５ ０．００７２ ３４．７ １５２．９０ ０．８８∗∗ ２．１７５×１０－７ １．８６７
１０～２０ Ｓ ０．００３１ ０．００６４ ４８．４ ９２．８０ ０．９３∗∗ １．９８４×１０－７ １．９７８
２０～３０ Ｇ ０．００１８ ０．００５１ ３５．３ １０６．３７ ０．８９∗∗ ３．４７６×１０－７ １．８８２
ＡＰ Ｈ Ｓ ０．０１７３ １３．７９０ ０．１ ７９．４３ ０．８７∗∗ １．３２５×１０－３ １．８０３
０～１０ Ｅ ０．０６９０ ０．６７６０ １０．２ ８４．３０ ０．９７∗∗ ３．７０４×１０－４ １．８２１
１０～２０ Ｅ ０．４６４０ １．０８４０ ４２．８ １０６．２０ ０．９４∗∗ ３．５５２×１０－３ １．８６９
２０～３０ Ｓ ０．０１３７ ０．０３３４ ４１．０ ９３．８０ ０．９９∗∗ ７．８１３×１０－７ １．８２５
∗Ｐ＜０．０５； ∗∗Ｐ＜０．０１． 下同 Ｔｈｅ ｓａｍｅ ｂｅｌｏｗ． Ｅ： 指数模型 Ｅｘｐｏｎｅｎｔｉａｌ； Ｓ： 球状模型 Ｓｐｈｅｒｉｃａｌ； Ｇ： 高斯模型 Ｇａｕｓｓｉａｎ．
林样地土壤全 Ｐ、有效 Ｐ 含量的空间结构特征，且它
们的半方差函数曲线变化较为平稳，在整个研究尺
度上各种生态过程同等重要．
从表 ３可以看出，各土层全 Ｐ、有效 Ｐ 的块金效
应均为正值，同一土层有效 Ｐ 的块金值均高于全 Ｐ，
但均接近于 ０，表明由试验误差和小于最小取样尺
度引起的随机变异小．各土层全 Ｐ 含量结构比为
３４􀆰 ７％～４８．４％，表明土壤全 Ｐ 含量由随机因素引起
的空间异质性占其总空间异质性的 ３４．７％ ～４８．４％，
由结构性因素引起的空间异质性占其总空间异质性
的 ５１．６％～６５．３％，南酸枣落叶阔叶林样地土壤全 Ｐ
含量具有中等强度的空间自相关性，主要是结构性
因素引起的．同样，各土层有效 Ｐ 的结构比为
０．１％～４２􀆰 ８％，表明土壤有效 Ｐ 含量由随机因素引
起的空间异质性占其总空间异质性的 ０． １％ ～
４２􀆰 ８％，由结构性因素引起的空间异质性占其总空
间异质性的 ５７．２％ ～ ９９．９％．其中，腐殖质层、０ ～ １０
ｃｍ土层有效 Ｐ 含量具有强烈的空间自相关性，１０～
２０、２０～３０ ｃｍ土层有效 Ｐ 含量具有中等强度的空间
自相关性，主要都是结构性因素引起的．
从表 ３ 可以看出，南酸枣落叶阔叶林样地各土
层全 Ｐ 含量变程为 ９２．８０～１６８．９０ ｍ，有效 Ｐ 含量变
程为 ７９．４３～１０６．２０ ｍ，全 Ｐ 含量变程随土层深度增
加而下降，而有效 Ｐ 含量变程随土层深度增加而增
加，表明全 Ｐ 含量空间自相关性尺度随土层深度增
加而下降，有效 Ｐ 则随土层深度增加而增加，全 Ｐ
含量较有效 Ｐ 含量有较大的空间自相关性尺度，空
间连续性高于有效 Ｐ，其生态过程在较大尺度上起
作用．各土层全 Ｐ、有效 Ｐ 含量的分维数分别为
１􀆰 ８６７～１．９７８、１．８０３～１．８６９，同一土层（除 ２０～３０ ｃｍ
土层外）有效 Ｐ 分维数均低于全 Ｐ 分维数，表明土
壤全 Ｐ、有效 Ｐ 含量具有良好的分形特征，有效 Ｐ 具
有更好的结构性，空间分布简单，全 Ｐ 的空间格局
相对复杂．
２􀆰 ３　 土壤全 Ｐ、有效 Ｐ 含量的空间分布
如图 ２所示，各土层全 Ｐ、有效 Ｐ 含量呈明显的
条带状和斑块状梯度性分布，在一些小样地内为相
似的空间分布格局，全 Ｐ、有效 Ｐ 含量最高值基本出
表 ４　 土壤全 Ｐ、有效 Ｐ含量与地形因子、生物因子、土壤因子的相关系数
Ｔａｂｌｅ ４　 Ｃｏｒｒｅｌａｔｉｏｎ ｃｏｅｆｆｉｃｉｅｎｔｓ ｏｆ ｓｏｉｌ ｔｏｔａｌ ｐｈｏｓｐｈｏｒｕｓ ａｎｄ ａｖａｉｌａｂｌｅ ｐｈｏｓｐｈｏｒｕｓ ｃｏｎｃｅｎｔｒａｔｉｏｎｓ ｗｉｔｈ ｔｏｐｏｇｒａｐｈｉｃ ｆａｃｔｏｒｓ，
ｂｉｏｌｏｇｉｃａｌ ｆａｃｔｏｒｓ ａｎｄ ｓｏｉｌ ｆａｃｔｏｒｓ
项目
Ｉｔｅｍ
土层
Ｓｏｉｌ
ｄｅｐｔｈ
（ｃｍ）
地形因子
Ｔｏｐｏｇｒａｐｈｉｃ ｆａｃｔｏｒ
高程
Ｅｌｅｖａｔｉｏｎ
凹凸度
Ｃｏｎｖｅｘｉｔｙ
生物因子
Ｂｉｏｌｏｇｉｃａｌ ｆａｃｔｏｒ
凋落物
现存量
Ｌｉｔｔｅｒ
物种数
Ｓｐｅｃｉｅｓ
株数
Ｎｕｍｂｅｒｓ
土壤因子
Ｓｏｉｌ ｆａｃｔｏｒ
含水率
Ｍｏｉｓｔｕｒｅ
ｃｏｎｔｅｎｔ
ｐＨ 有机碳
Ｏｒｇａｎｉｃ Ｃ
全氮
Ｔｏｔａｌ Ｎ
ＴＰ Ｈ －０．３４１∗∗ －０．１８１ －０．２４０∗ ０．０１９ ０．１１４ － ０．２８９∗∗ ０．２６６∗∗ ０．３４２∗∗
０～１０ －０．２５１∗ －０．０６３ －０．２１７∗ －０．００７ ０．０９８ －０．２７０∗∗ ０．２２８∗ ０．２８７∗∗ ０．３２４∗∗
１０～２０ －０．２４０∗ －０．１１３ －０．２３１∗ ０．００３ ０．００６ －０．０４９ ０．０２０ ０．３８８∗∗ ０．３０８∗∗
２０～３０ －０．４１４∗∗ －０．２９０∗∗ －０．２１２∗ －０．２２０∗ －０．１７０ －０．０７１ ０．１２８ ０．５１７∗∗ ０．５７４∗∗
ＡＰ Ｈ －０．２１６∗ ０．０４６ －０．２０１∗ ０．０５９ ０．０９６ － ０．０５７ ０．３４１∗∗ ０．２８９∗∗
０～１０ －０．３１２∗∗ －０．１２５ －０．２５０∗ －０．０６５ －０．０１０ ０．２２５∗ ０．０２６ ０．４０４∗∗ ０．４４１∗∗
１０～２０ －０．２１８∗ －０．０７６ －０．２０９∗ －０．００２ ０．０７６ ０．０５８ ０．１９２ ０．４８２∗∗ ０．３２９∗∗
２０～３０ －０．３８３∗∗ －０．１０９ －０．２３０∗ －０．０９４ －０．０４１ ０．０６９ ０．２０４∗ ０．５３３∗∗ ０．３８６∗∗
０１７ 应　 用　 生　 态　 学　 报　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 ２７卷
图 ２　 土壤全 Ｐ（ＴＰ， ｇ·ｋｇ－１）和有效 Ｐ（ＡＰ， ｍｇ·ｋｇ－１）含量的空间分布格局
Ｆｉｇ．２　 Ｓｐａｔｉａｌ ｐａｔｔｅｒｎ ｏｆ ｔｏｔａｌ ｐｈｏｓｐｈｏｒｕｓ （ＴＰ， ｇ·ｋｇ－１） ａｎｄ ａｖａｉｌａｂｌｅ ｐｈｏｓｐｈｏｒｕｓ （ＡＰ， ｍｇ·ｋｇ－１） ｃｏｎｃｅｎｔｒａｔｉｏｎ ｉｎ ｓｏｉｌ．
现在高程较低的山谷洼地，最低值则出现在高程高
的山脊地带，全 Ｐ 含量随地形变化的规律不如有效
Ｐ 含量明显．
２􀆰 ４　 土壤全 Ｐ 和有效 Ｐ 含量变异的影响因子
从表 ４可以看出，土壤全 Ｐ、有效 Ｐ 含量与高程
呈显著负相关，与凹凸度也呈负相关，但除 ２０ ～ ３０
ｃｍ土层全 Ｐ 含量外，其他土层均未达到显著水平．
这表明高程、凹凸度对土壤 Ｐ 含量影响具有一致
性，高程对土壤 Ｐ 含量的影响显著高于凹凸度．
土壤全 Ｐ、有效 Ｐ 含量与地表凋落物现存量呈
显著负相关，与样地植物种数呈负相关，但除 ２０～３０
ｃｍ土层全 Ｐ 含量外，其他土层均未达到显著水平，
与样地植物株数不存在相关性．土壤全 Ｐ 含量与土
壤含水率呈负相关，除 ０ ～ １０ ｃｍ 土层外，其他土层
均未达到显著水平；土壤有效 Ｐ 含量与土壤含水率
呈正相关，除 ０ ～ １０ ｃｍ 土层外，其他土层均未达到
１１７３期　 　 　 　 　 　 　 　 　 　 胡瑞彬等： 中亚热带南酸枣落叶阔叶林土壤磷素空间变异及其影响因素　 　 　 　 　 　
显著水平．土壤全 Ｐ、有效 Ｐ 含量与土壤 ｐＨ 呈正相
关，除腐殖质层全 Ｐ、０ ～ １０ ｃｍ 土层全 Ｐ、２０ ～ ３０ ｃｍ
土层有效 Ｐ 外，其他土层与土壤 ｐＨ 相关性均未达
到显著水平．土壤全 Ｐ、有效 Ｐ 与土壤有机 Ｃ、全 Ｎ
均呈显著正相关．
３　 讨　 　 论
无论大尺度还是小尺度，土壤性质的空间异质
性普遍存在［３１－３２］ ．本研究中，各土层全 Ｐ、有效 Ｐ 含
量均为中等强度变异，且有效 Ｐ 含量的变异程度大
于全 Ｐ，可能与研究样地地形复杂多变、微生境丰
富［２２］、土壤 Ｐ 的淋溶，以及地上植物分布及其选择
性吸收利用特征有关．南酸枣落叶阔叶林 ０ ～ ３０ ｃｍ
土层全 Ｐ 算术平均含量为 ０．２６ ｇ·ｋｇ－１，接近贵州安
顺地区喀斯特小流域乔木林（０．２８ ｇ·ｋｇ－１） ［７］和浙
江天童阔叶林土壤全 Ｐ 含量（０．２６ ｇ·ｋｇ－１） ［１９］，低
于湖南会同杉木人工林土壤全 Ｐ 含量 （ ０􀆰 ３８
ｇ·ｋｇ－１） ［３３］，明显低于黑龙江西北部克山县土壤
（０􀆰 ７４～ ０． ８６ ｇ· ｋｇ－１ ） ［３４］ 和全国土壤全 Ｐ 含量
（０􀆰 ６０ ｇ·ｋｇ－１） ［３５］ ．南酸枣落叶阔叶林 ０ ～ ３０ ｃｍ 土
层有效 Ｐ 含量均值为 ２．３３ ｍｇ·ｋｇ－１，低于贵州安顺
地区喀斯特小流域乔木林（３．９６ ｍｇ·ｋｇ－１） ［７］、福建
省三明市格氏栲自然保护区（５．５８ ｍｇ·ｋｇ－１） ［５］和
全国土壤有效 Ｐ 含量（３．４０ ｍｇ·ｋｇ－１） ［３５］ ．各土层全
Ｐ、有效 Ｐ（除腐殖质层外）含量均处于极缺水平［３６］，
表明南酸枣落叶阔叶林土壤极缺 Ｐ．土壤养分有限
或缺乏影响林下植物对土壤有效养分的吸收与利
用，进而影响更新格局和过程［３７］，意味着土壤 Ｐ 含
量成为了南酸枣落叶阔叶林生长的限制因子．此外，
南酸枣落叶阔叶林土壤 Ｐ 具有较大的空间变异性，
可能对地上植被组成、空间分布及相关的生理生态
过程产生深远的影响［３０］，也可能是制约南酸枣落叶
阔叶林幼苗更新的重要因子．由于地表凋落物、根系
分解及其分泌物形成的有机质和养分首先进入土壤
表层，使得土壤养分（Ｎ、Ｐ）含量随土层深度增加而
下降．本研究中，随着土壤深度增加，土壤全 Ｐ、有效
Ｐ 含量呈递减趋势，变异程度增大，与一些研究［５，３８］
基本一致，表明南酸枣落叶阔叶林土壤 Ｐ 含量具有
明显的层次性，表层富集现象较明显，因此，应注意
保护表层植被，维持土壤养分库的稳定．
土壤养分空间异质性是结构性因素和随机性因
素综合作用的结果［３９］ ．通常，结构性因素强化其空
间相关性，而随机性因素促使其空间相关性减
弱［５，４０］ ．本研究中，各土层全 Ｐ、有效 Ｐ 含量的结构
比存在差异，同一土层全 Ｐ 含量的结构比高于有效
Ｐ（除 ２０～３０ ｃｍ 土层外），表明土壤有效 Ｐ 含量的
空间自相关性较全 Ｐ 含量强，明显受到结构性因素
的影响，可能与全 Ｐ、有效 Ｐ 淋溶特征和植物选择性
吸收利用有关．各土层全 Ｐ、有效 Ｐ 含量由结构性因
素引起的空间异质性分别占系统总变异比例的
５１．６％～６５．３％、５７．２％ ～ ９９．９％，表明土壤全 Ｐ、有效
Ｐ 含量在研究尺度上主要受结构性因素影响，即成
土母质、气候条件较为一致情况下，地形、林分结构、
土壤有机 Ｃ、Ｎ是土壤全 Ｐ、有效 Ｐ 空间异质性产生
的重要原因，可用相关系数、回归模型和分形维数等
定量表达，而随机性因素（采样、测定等人为干扰）
影响相对较小．有效变程（Ａ０）反映了区域变量空间
自相关性的范围，与观测尺度以及在取样尺度上影
响土壤养分的各种生态过程相互作用有关．在变程
内，变量具有空间自相关性，反之不存在．因此，Ａ０提
供了研究某种属性相似范围的一种测度，对土壤养
分取样设计有效性及区域变量的空间内插和制图具
有指导作用［５］ ．本研究中，各土层全 Ｐ、有效 Ｐ 含量
的空间自相关范围差异较大，如为传统统计分析构
建空间独立的数据，腐殖质层、０ ～ １０ ｃｍ 土层全 Ｐ
的取样尺度必须大于 １５０ ｍ，其他土层要大于 ９２ ｍ，
有效 Ｐ 必须大于 ７９ ｍ．若为土壤 Ｐ 的景观制图服
务，取样网格应小于 Ａ０，变程以外的取样对任何内
插和制图均无效［４１］ ．
研究表明，分维数高、养分空间格局复杂、随机
因素引起的异质性占系统变异的比值大，随机性较
强；相反，养分空间依赖性更强，具有更好的结构性，
空间分布简单，由结构性因素引起的空间异质性占
系统变异的比值大，与结构比的变化一致［４０］ ．本研
究中，各土层全 Ｐ、有效 Ｐ 含量的分维数分别为
１􀆰 ８６７～１．９７８、１．８０３ ～ １．８６９，尽管分维数差异不大，
但同一土层全 Ｐ 含量的分形维数比有效 Ｐ 高
０􀆰 ０６４～０􀆰 １０９，表明两者空间格局的局部变化，随机
性和结构性因素引起的空间异质性程度的不同，有
效 Ｐ 的空间依赖性更强，具有更好的结构性，空间
分布简单，而全 Ｐ 的空间格局较有效 Ｐ 复杂，随机
因素引起的空间异质性占系统变异性的比值大，随
机性较强，与结构比的变化一致（表 ３），腐殖质层、
０～１０ ｃｍ土层的有效 Ｐ 含量结构比低，由结构性因
素引起的空间异质性占系统变异的比值大，而各土
层全 Ｐ 含量结构比大，由随机因素引起的异质性程
度也较大．Ｋｒｉｇｉｎｇ 插值图可直观地反映各土壤养分
的空间分布特征．研究样地各土层全 Ｐ、有效 Ｐ 含量
２１７ 应　 用　 生　 态　 学　 报　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 ２７卷
具有高度的空间异质性，并决定了其空间分布格局
的存在．各土层全 Ｐ、有效 Ｐ 含量随着地形因子的变
化呈现条带状和斑块状梯度性分布，且全 Ｐ、有效 Ｐ
的变化趋势基本一致，全 Ｐ、有效 Ｐ 含量最高值基本
出现在高程较低的山谷洼地，最低值则出现在高程
高的山脊地带，有效 Ｐ 随地形因子的变化趋势比全
Ｐ 的变化趋势明显．可能是由于样地地形复杂、土壤
Ｐ 含量具有“洼积效应”，易在海拔低的坡谷汇集．表
明在某些小尺度下，全 Ｐ、有效 Ｐ 有较好的空间关联
性，土壤全 Ｐ、有效 Ｐ 含量均明显受到地形因子以及
多种影响因子的作用．此外，腐殖质层、０ ～ １０ ｃｍ 土
层全 Ｐ、有效 Ｐ 的 Ｋｒｉｇｉｎｇ 插值图是破碎的，空间异
质性明显，可能与 Ｐ 的水平迁移能力及稳定性、植
物选择吸收利用、林窗微生境等有关［５］ ．
地形是重要的土壤成土因子之一，调控太阳辐
射、降水的空间再分配，影响局部生境的小气候条件
及土壤厚度和养分的空间差异［４２］ ．研究表明，海拔、
凹凸度、坡度等地形因子对中小尺度土壤养分空间
变异有明显影响［１９，３２，４３］，制约着植物群落分布、生物
多样性及森林生态系统的健康发育［４３］ ．海拔梯度包
含了温度、湿度、光照和土壤属性等直接生境因子的
多尺度变化，是影响山地生境差异性的主导因
子［４４］ ．本研究中，土壤全 Ｐ、有效 Ｐ 含量与高程呈显
著负相关，也与凹凸度呈负相关，但不显著（除 ２０ ～
３０ ｃｍ土层全 Ｐ 含量外），这与张娜等［１９］、樊纲惟
等［２１］、张忠华等［３２］的研究结果基本一致，表明高
程、凹凸度对土壤 Ｐ 含量的影响具有一致性，显著
影响着土壤 Ｐ 的空间分布，低高程、洼地土壤 Ｐ 含
量较高，与 Ｋｒｉｇｉｎｇ 插值的空间分布格局一致，也反
映了土壤 Ｐ 的淋溶特征．
研究表明，土壤养分（Ｎ、Ｐ、Ｋ）与地表凋落物现
存量呈显著正相关［４５］，与乔木树种多样性呈正相
关，但不显著［４６］ ．也有研究发现，环境异质性与物种
多样性呈负相关性［４７］ ．本研究中，土壤全 Ｐ、有效 Ｐ
含量与凋落物现存量呈显著负相关，与样地树种种
数、株数相关性不显著．在调查过程中发现，样地的
凹凸度变化不十分明显，高海拔山脊处（凸地形）的
树种组成较为丰富，地表凋落物现存量较高，而凹地
形处地表凋落物现存量较少，可能是凹地形微环境
条件有利于凋落物的分解，加上凸地形土壤 Ｐ 向凹
地形样地淋溶汇集，使得凹地形土壤 Ｐ 含量较高，
导致土壤 Ｐ 含量与地表凋落物现存量呈显著负相
关．土壤全 Ｐ、有效 Ｐ 含量与土壤 ｐＨ呈正相关，表明
土壤 ｐＨ增高，有利于土壤 Ｐ 积累［４８］，与土壤 Ｐ 的
分布和可利用性首先取决于土壤 ｐＨ 值有关［３０］ ．水
分是影响土壤养分空间分布的一个重要因子，土壤
含水率与土壤全 Ｐ 含量呈负相关，与土壤有效 Ｐ 含
量呈正相关，但仅 ０ ～ １０ ｃｍ 土层达到显著水平，表
明土壤水分增加，减少 Ｐ 的吸附，易溶于水而随其
流动，有利于提高土壤 Ｐ 的有效性，不利于土壤全 Ｐ
的积累．土壤全 Ｐ、有效 Ｐ 含量与土壤有机 Ｃ、Ｎ含量
均呈显著正相关，表明土壤 Ｐ 与 Ｃ、Ｎ含量的变化趋
势基本一致，三者之间有着密切的耦合关系，影响着
土壤的生产力．研究表明，土壤有机 Ｃ可以显著活化
土壤 Ｐ，增加土壤 Ｐ 的有效性，减少 Ｐ 的吸附［４８］ ．由
于多种环境因子的综合作用，研究样地土壤 Ｐ 的空
间分布特征比较复杂，今后应加强对森林土壤养分
空间格局机理的研究．
参考文献
［１］　 Ｚｈａｎｇ Ｗ （张　 伟）， Ｌｉｕ Ｓ⁃Ｊ （刘淑娟）， Ｙｅ Ｙ⁃Ｙ （叶
莹莹）， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ａｎｄ ｉｎ⁃
ｆｌｕｅｎｃｉｎｇ ｆａｃｔｏｒｓ ｉｎ ｔｙｐｉｃａｌ Ｋａｒｓｔ ｖｉｒｇｉｎ ｆｏｒｅｓｔ． Ｔｒａｎｓａｃ⁃
ｔｉｏｎｓ ｏｆ ｔｈｅ Ｃｈｉｎｅｓｅ Ｓｏｃｉｅｔｙ ｏｆ Ａｇｒｉｃｕｌｔｕｒａｌ Ｅｎｇｉｎｅｅｒｉｎｇ
（农业工程学报）， ２０１３， ２９（１）： ９３－ １０１ （ ｉｎ Ｃｈｉ⁃
ｎｅｓｅ）
［２］　 Ｄｒａｈｏｒａｄ Ｓ， Ｆｅｌｉｘ⁃Ｈｅｎｎｉｎｇｓｅｎ Ｐ， Ｅｃｋｈａｒｄｔ ＫＵ， ｅｔ ａｌ．
Ｓｐａｔｉａｌ ｃａｒｂｏｎ ａｎｄ ｎｉｔｒｏｇｅｎ ｄｉｓｔｒｉｂｕｔｉｏｎ ａｎｄ ｏｒｇａｎｉｃ ｍａｔ⁃
ｔｅｒ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｂｉｏｌｏｇｉｃａｌ ｓｏｉｌ ｃｒｕｓｔｓ ｉｎ ｔｈｅ Ｎｅｇｅｖ
Ｄｅｓｅｒｔ （Ｉｓｒａｅｌ） ａｌｏｎｇ ａ ｒａｉｎｆａｌｌ ｇｒａｄｉｅｎｔ． Ｊｏｕｒｎａｌ ｏｆ Ａｒｉｄ
Ｅｎｖｉｒｏｎｍｅｎｔｓ， ２０１３， ９４： １８－２６
［３］　 Ｐｅｎｇ Ｗ⁃Ｘ （彭晚霞）， Ｓｏｎｇ Ｔ⁃Ｑ （宋同清）， Ｚｅｎｇ Ｆ⁃Ｐ
（曾馥平）， ｅｔ ａｌ． Ｔｈｅ ｃｏｕｐｌｉｎｇ ｒｅｌａｔｉｏｎｓｈｉｐｓ ｂｅｔｗｅｅｎ
ｖｅｇｅｔａｔｉｏｎ， ｓｏｉｌ， ａｎｄ ｔｏｐｏｇｒａｐｈｙ ｆａｃｔｏｒｓ ｉｎ Ｋａｒｓｔ ｍｉｘｅｄ
ｅｖｅｒｇｒｅｅｎ ａｎｄ ｄｅｃｉｄｕｏｕｓ ｂｒｏａｄｌｅａｆ ｆｏｒｅｓｔ． Ａｃｔａ Ｅｃｏｌｏｇｉｃａ
Ｓｉｎｉｃａ （生态学报）， ２０１０， ３０（１３）： ３４７２－３４８１ （ ｉｎ
Ｃｈｉｎｅｓｅ）
［４］　 Ｗａｎ Ｙ⁃Ｚ （王永壮）， Ｃｈｅｎ Ｘ （陈　 欣）， Ｓｈｉ Ｙ （史　
奕）． Ｐｈｏｓｐｈｏｒｕｓ ａｖａｉｌａｂｉｌｉｔｙ ｉｎ ｃｒｏｐｌａｎｄ ｓｏｉｌｓ ｏｆ Ｃｈｉｎａ
ａｎｄ ｒｅｌａｔｅｄ ａｆｆｅｃｔｉｎｇ ｆａｃｔｏｒｓ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ
Ｅｃｏｌｏｇｙ （应用生态学报）， ２０１３， ２４（１）： ２６０－ ２６８
（ｉｎ Ｃｈｉｎｅｓｅ）
［５］　 Ｓｕ Ｓ⁃Ｊ （苏松锦）， Ｌｉｕ Ｊ⁃Ｆ （刘金福）， Ｈｅ Ｚ⁃Ｓ （何中
声）， ｅｔ ａｌ． Ｔｈｅ ｓｐａｔｉａｌ ｈｅｔｅｒｏｇｅｎｅｉｔｙ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｉｎ
ａ ｍｉｄ⁃ｓｕｂｔｒｏｐｉｃａｌ Ｃａｓｔａｎｏｐｓｉｓ ｋａｗａｋａｍｉｉ ｎａｔｕｒａｌ ｆｏｒｅｓｔ．
Ａｃｔａ Ｅｃｏｌｏｇｉｃａ Ｓｉｎｉｃａ （生态学报）， ２０１２， ３２（１８）：
５６７３－５６８２ （ｉｎ Ｃｈｉｎｅｓｅ）
［６］　 Ｙａｖｉｔｔ ＪＢ， Ｈａｒｍｓ ＫＥ， Ｇａｒｃｉａ ＭＮ， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｈｅｔｅｒｏ⁃
ｇｅｎｅｉｔｙ ｏｆ ｓｏｉｌ ｃｈｅｍｉｃａｌ ｐｒｏｐｅｒｔｉｅｓ ｉｎ ａ ｌｏｗｌａｎｄ ｔｒｏｐｉｃａｌ
ｍｏｉｓｔ ｆｏｒｅｓｔ， Ｐａｎａｍａ． Ａｕｓｔｒａｌｉａｎ Ｊｏｕｒｎａｌ ｏｆ Ｓｏｉｌ Ｒｅ⁃
ｓｅａｒｃｈ， ２００９， ４７： ６７４－６８７
［７］　 Ｈｕ Ｚ⁃Ｌ （胡忠良）， Ｐａｎ Ｇ⁃Ｘ （潘根兴）， Ｌｉ Ｌ⁃Ｑ （李
恋卿）， ｅｔ ａｌ． Ｃｈａｎｇｅｓ ｉｎ ｐｏｏｌｓ ａｎｄ ｈｅｔｅｒｏｇｅｎｅｉｔｙ ｏｆ ｓｏｉｌ
ｏｒｇａｎｉｃ ｃａｒｂｏｎ， ｎｉｔｒｏｇｅｎ ａｎｄ Ｐｈｏｓｐｈｏｒｕｓ ｕｎｄｅｒ ｄｉｆｆｅｒｅｎｔ
ｖｅｇｅｔａｔｉｏｎ ｔｙｐｅｓ ｉｎ Ｋａｒｓｔ ｍｏｕｎｔａｉｎｏｕｓ ａｒｅａ ｏｆ ｃｅｎｔｒａｌ
Ｇｕｉｚｈｏｕ Ｐｒｏｖｉｎｃｅ， Ｃｈｉｎａ． Ａｃｔａ Ｅｃｏｌｏｇｉｃａ Ｓｉｎｉｃａ （生态
学报）， ２００９， ２９（８）： ４１８７－４１９５ （ｉｎ Ｃｈｉｎｅｓｅ）
［８］　 Ｗａｎｇ ＹＰ， Ｌａｗ ＲＭ， Ｐａｋ Ｂ． Ａ ｇｌｏｂａｌ ｍｏｄｅｌ ｏｆ ｃａｒｂｏｎ，
３１７３期　 　 　 　 　 　 　 　 　 　 胡瑞彬等： 中亚热带南酸枣落叶阔叶林土壤磷素空间变异及其影响因素　 　 　 　 　 　
ｎｉｔｒｏｇｅｎ ａｎｄ ｐｈｏｓｐｈｏｒｕｓ ｃｙｃｌｅｓ ｆｏｒ ｔｈｅ ｔｅｒｒｅｓｔｒｉａｌ ｂｉｏ⁃
ｓｐｈｅｒｅ． Ｂｉｏｇｅｏｓｃｉｅｎｃｅｓ， ２０１０， ７： ２２６１－２２８２
［９］　 Ｚｈａｎｇ Ｄ⁃Ｈ （张鼎华）， Ｌｉｎ Ｋ⁃Ｍ （林开淼）， Ｌｉ Ｂ⁃Ｆ
（李宝福）． Ｐｈｏｓｐｈｏｒｕｓ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｉｎ ｒｈｉｚｏｓｐｈｅｒｅ
ｓｏｉｌ ｏｆ Ｃｕｎｎｉｎｇｈａｍｉａ ｌａｎｃｅｏｌａｔａ， Ｐｉｎｕｓ ｍａｓｓｏｎｉａｎａ ａｎｄ
ｔｈｅｉｒ ｍｉｘｅｄ ｐｌａｎｔａｔｉｏｎｓ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏ⁃
ｌｏｇｙ （应用生态学报）， ２０１１， ２２ （ １１）： ２８１５ － ２８２１
（ｉｎ Ｃｈｉｎｅｓｅ）
［１０］　 Ｚｈａｎ Ｓ⁃Ｘ （詹书侠）， Ｃｈｅｎ Ｆ⁃Ｓ （陈伏生）， Ｈｕ Ｘ⁃Ｆ
（胡小飞）， ｅｔ ａｌ． Ｓｏｉｌ ｎｉｔｒｏｇｅｎ ａｎｄ ｐｈｏｓｐｈｏｒｕｓ ａｖａｉｌａ⁃
ｂｉｌｉｔｙ ｉｎ ｆｏｒｅｓｔ ｅｃｏｓｙｓｔｅｍｓ ａｔ ｄｉｆｆｅｒｅｎｔ ｓｔａｇｅｓ ｏｆ ｓｕｃｃｅｓ⁃
ｓｉｏｎ ｉｎ ｔｈｅ ｃｅｎｔｒａｌ ｓｕｂｔｒｏｐｉｃａｌ ｒｅｇｉｏｎ． Ａｃｔａ Ｅｃｏｌｏｇｉｃａ
Ｓｉｎｉｃａ （生态学报）， ２００９， ２９（９）： ４６７３－ ４６８０ （ ｉｎ
Ｃｈｉｎｅｓｅ）
［１１］　 Ｂｒｏｎｓｏｎ ＫＦ， Ｚｏｂｅｃｋ ＴＭ， Ｃｈｕａ ＴＴ， ｅｔ ａｌ． Ｃａｒｂｏｎ ａｎｄ
ｎｉｔｒｏｇｅｎ ｐｏｏｌｓ ｏｆ ｓｏｕｔｈｅｒｎ ｈｉｇｈ ｐｌａｉｎｓ ｃｒｏｐｌａｎｄ ａｎｄ
ｇｒａｓｓｌａｎｄ ｓｏｉｌｓ． Ｓｏｉｌ Ｓｃｉｅｎｃｅ Ｓｏｃｉｅｔｙ ｏｆ Ａｍｅｒｉｃａ Ｊｏｕｒｎａｌ，
２００４， ６８： １６９５－１７０４
［１２］　 Ｌｕｏ Ｘ⁃Ｘ （罗先香）， Ｄｕｎ Ｍ （郭　 萌）， Ｙａｎ Ｑ （闫　
琴）． Ｄｙｎａｍｉｃ ｄｉｓｔｒｉｂｕｔｉｏｎ ａｎｄ ｉｎｆｌｕｅｎｃｅ ｆａｃｔｏｒｓ ｏｆ ｓｏｉｌ
ｐｈｏｓｐｈｏｒｕｓ ｉｎ Ｙｅｌｌｏｗ Ｒｉｖｅｒ ｅｓｔｕａｒｙ ｗｅｔｌａｎｄ． Ｊｏｕｒｎａｌ ｏｆ
Ｓｏｉｌ ａｎｄ Ｗａｔｅｒ Ｃｏｎｓｅｒｖａｔｉｏｎ （水土保持学报）， ２０１１，
２５（５）： １５４－１６０ （ｉｎ Ｃｈｉｎｅｓｅ）
［１３］　 Ｑｉｎ Ｓ⁃Ｊ （秦胜金）， Ｌｉｕ Ｊ⁃Ｓ （刘景双）， Ｗａｎｇ Ｇ⁃Ｐ （王
国平）． Ｍｅｃｈａｎｉｓｍ ｏｆ ｐｈｏｓｐｈｏｒｕｓ ａｖａｉｌａｂｉｌｉｔｙ ｃｈａｎｇｉｎｇ
ｉｎ ｓｏｉｌ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ｓｏｉｌ Ｓｃｉｅｎｃｅ （土壤通报），
２００６， ３７（５）： １０１２－１０１６ （ｉｎ Ｃｈｉｎｅｓｅ）
［１４］　 Ｃｈｅｎ Ｌ⁃Ｘ （陈立新）， Ｙａｎｇ Ｃ⁃Ｄ （杨承栋）． Ｔｈｅ ｓｕｃ⁃
ｃｅｓｓｉｏｎ ｏｆ ｖａｒｉｏｕｓ ｔｙｐｅｓ ｏｆ ｐｈｏｓｐｈｏｒｕｓ， ｐｈｏｓｐｈａｔａｓｅ ａｃ⁃
ｔｉｖｉｔｙ， ａｎｄ ｔｈｅ ｒｅｌａｔｉｏｎｓｈｉｐ ｗｉｔｈ ｔｈｅ ｔｒｅｅ ｇｒｏｗｔｈ ｉｎ ｌａｒｃｈ
ｐｌａｎｔａｔｉｏｎｓ． Ｓｃｉｅｎｔｉａ Ｓｉｌｖａｅ Ｓｉｎｉｃａｅ （林业科学）， ２００４，
４０（３）： １２－１８ （ｉｎ Ｃｈｉｎｅｓｅ）
［１５］　 Ｙａｎｇ Ｑ （杨　 青）， Ｚｈａｎｇ Ｙ （张　 一）， Ｚｈｏｕ Ｚ⁃Ｃ （周
志春）， ｅｔ ａｌ． Ｇｅｎｅｔｉｃ ｖａｒｉａｔｉｏｎ ｉｎ ｒｏｏｔ ａｒｃｈｉｔｅｃｔｕｒｅ ａｎｄ
ｐｈｏｓｐｈｏｒｕｓ ｅｆｆｉｃｉｅｎｃｙ ｉｎ ｒｅｓｐｏｎｓｅ ｔｏ ｈｅｔｅｒｏｇｅｎｅｏｕｓ ｐｈｏｓ⁃
ｐｈｏｒｕｓ ｄｅｆｉｃｉｅｎｃｙ ｉｎ Ｐｉｎｕｓ ｍａｓｓｏｎｉａｎａ ｆａｍｉｌｉｅｓ． Ｃｈｉｎｅｓｅ
Ｊｏｕｒｎａｌ ｏｆ Ｐｌａｎｔ Ｅｃｏｌｏｇｙ （植物生态学报）， ２０１１， ３５
（１２）： １２２６－１２３５ （ｉｎ Ｃｈｉｎｅｓｅ）
［１６］　 Ｙｕ Ｊ （余　 健）． Ｕｎｄｅｒ Ｐｈｏｓｐｈｏｒｕｓ Ｓｔｒｅｓｓ Ｔｒｅｅｓ Ｓｅｃｒｅ⁃
ｔｉｏｎ ｏｆ Ｏｒｇａｎｉｃ Ａｃｉｄｓ ａｎｄ ｔｈｅ Ａｃｔｉｖａｔｉｏｎ ｏｆ Ｓｏｉｌ Ｐｈｏｓ⁃
ｐｈｏｒｕｓ． ＰｈＤ ｄｉｓｓｅｒｔａｔｉｏｎ． Ｎａｎｊｉｎｇ： Ｎａｎｊｉｎｇ Ｆｏｒｅｓｔｒｙ
Ｕｎｉｖｅｒｓｉｔｙ， ２００５ （ｉｎ Ｃｈｉｎｅｓｅ）
［１７］　 Ａｃｈｉｍ Ｄ， Ｔｈｏｍａｓ Ｇ， Ｎｉｅｌｓ Ｔ． Ｐｈｏｓｐｈｏｒｕｓ ｆｅｒｔｉｌｉｚｅｒ
ｅｆｆｅｃｔｓ ｏｎ ｓｏｉｌ ｐｈｏｓｐｈｏｒｕｓ ｐｏｏｌｓ ｉｎ ａｃｉｄ ｕｐｌａｎｄ ｓｏｉｌｓ． Ｓｏｉｌ
Ｓｃｉｅｎｃｅ， ２００２， ６６： ６５２－６６０
［１８］　 Ｚｈａｎｇ Ｓ⁃Ｒ （张世熔）， Ｈｕａｎｇ Ｙ⁃Ｆ （黄元仿）， Ｌｉ Ｂ⁃Ｇ
（李保国）， ｅｔ ａｌ． Ｔｅｍｐｏｒａｌ⁃ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｏｆ ｓｏｉｌ
ａｖａｉｌａｂｌｅ ｐｈｏｓｐｈｏｒｕｓ ａｎｄ ｐｏｔａｓｓｉｕｍ ｉｎ ｔｈｅ ａｌｌｕｖｉａｌ ｒｅｇｉｏｎ
ｏｆ ｔｈｅ Ｈｕａｎｇ⁃Ｈｕａｉ⁃Ｈａｉ Ｐｌａｉｎ． Ｐｌａｎｔ Ｎｕｔｒｉｔｉｏｎ ａｎｄ Ｆｅｒｔｉ⁃
ｌｉｚｅｒ Ｓｃｉｅｎｃｅ （植物营养与肥料学报）， ２００３， ９（１）：
３－８ （ｉｎ Ｃｈｉｎｅｓｅ）
［１９］　 Ｚｈａｎｇ Ｎ （张　 娜）， Ｗａｎｇ Ｘ⁃Ｈ （王希华）， Ｚｈｅｎｇ Ｚ⁃Ｍ
（郑泽梅）， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｈｅｔｅｒｏｇｅｎｅｉｔｙ ｏｆ ｓｏｉｌ ｐｒｏｐｅｒｔｉｅｓ
ａｎｄ ｉｔｓ ｒｅｌａｔｉｏｎｓｈｉｐｓ ｗｉｔｈ ｔｅｒｒａｉｎ ｆａｃｔｏｒｓ ｉｎ ｂｒｏａｄｌｅａｖｅｄ
ｆｏｒｅｓｔ ｉｎ Ｔｉａｎｔｏｎｇ ｏｆ Ｚｈｅｊｉａｎｇ Ｐｒｏｖｉｎｃｅ， Ｅａｓｔ Ｃｈｉｎａ．
Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应用生态学报），
２０１２， ２３（９）： ２３６１－２３６９ （ｉｎ Ｃｈｉｎｅｓｅ）
［２０］　 Ｑｉｎ Ｓ （秦　 松）， Ｆａｎ Ｙ （樊　 燕）， Ｌｉｕ Ｈ⁃Ｂ （刘洪
斌）， ｅｔ ａｌ． Ｓｔｕｄｙ ｏｎ ｔｈｅ ｒｅｌａｔｉｏｎｓ ｂｅｔｗｅｅｎ ｔｏｐｏｇｒａｐｈｉｃａｌ
ｆａｃｔｏｒｓ ａｎｄ ｔｈｅ ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎｓ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ．
Ｊｏｕｒｎａｌ ｏｆ Ｓｏｉｌ ａｎｄ Ｗａｔｅｒ Ｃｏｎｓｅｒｖａｔｉｏｎ （水土保持学
报）， ２００８， １５（１）： ４６－５２ （ｉｎ Ｃｈｉｎｅｓｅ）
［２１］　 Ｆａｎ Ｇ⁃Ｗ （樊纲惟）， Ｘｉａｎｇ Ｗ⁃Ｈ （项文化）， Ｌｅｉ Ｐ⁃Ｆ
（雷丕峰）， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ ａｎｄ ｄｒｉｖｉｎｇ ｆａｃｔｏｒｓ
ｏｆ ｓｏｉｌ ｐｈｏｓｐｈｏｒｕｓ ｉｎ ａ ｓｕｂｔｒｏｐｉｃａｌ ｅｖｅｒｇｒｅｅｎ ｂｒｏａｄｌｅａｆ
ｆｏｒｅｓｔ． Ｒｅｓｅａｒｃｈ ｏｆ Ａｇｒｉｃｕｌｔｕｒａｌ Ｍｏｄｅｒｎｉｚａｔｉｏｎ （农业现
代化研究）， ２０１４， ３５（３）： ３６７－３７０ （ｉｎ Ｃｈｉｎｅｓｅ）
［２２］　 Ｙｉ Ｈ （易 　 好）， Ｄｅｎｇ Ｘ⁃Ｗ （邓湘雯）， Ｘｉａｎｇ Ｗ⁃Ｈ
（项文化）， ｅｔ ａｌ． Ｔｈｅ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ａｎｄ ｒｅｇｅｎｅｒａｔｉｏｎ ｏｆ
ｔｈｅ Ｃｈｏｅｒｏｓｐｏｎｄｉａｓ ａｘｉｌｌａｒｉｅｓ ｂｒｏａｄ⁃ｌｅａｖｅｄ ｃｏｍｍｕｎｉｔｙ ｉｎ
ｔｈｅ ｈｉｌｌｙ ｒｅｇｉｏｎ ｏｆ ｃｅｎｔｒａｌ Ｈｕｎａｎ Ｐｒｏｖｉｎｃｅ， Ｃｈｉｎａ． Ａｃｔａ
Ｅｃｏｌｏｇｉｃａ Ｓｉｎｉｃａ （生态学报）， ２０１４， ３４（１２）： ３４６３－
３４７１ （ｉｎ Ｃｈｉｎｅｓｅ）
［２３］　 Ｚｈａｏ Ｌ⁃Ｊ （赵丽娟）， Ｘｉａｎｇ Ｗ⁃Ｈ （项文化）， Ｌｉ Ｊ⁃Ｘ
（李家湘）， ｅｔ ａｌ． Ｆｌｏｒｉｓｔｉｃ ｃｏｍｐｏｓｉｔｉｏｎ， ｓｔｒｕｃｔｕｒｅ ａｎｄ
ｐｈｙｔｏｇｅｏｇｒａｐｈｉｃ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｉｎ ａ Ｌｉｔｈｏｃａｒｐｕｓ ｇｌａｂｅｒ－
Ｃｙｃｌｏｂａｌａｎｏｐｓｉｓ ｇｌａｕｃａ ｆｏｒｅｓｔ ｃｏｍｍｕｎｉｔｙ ｉｎ ｔｈｅ ｓｕｂｔｒｏｐｉ⁃
ｃａｌ ｒｅｇｉｏｎ． Ｓｃｉｅｎｔｉａ Ｓｉｌｖａｅ Ｓｉｎｉｃａｅ （林业科学）， ２０１３，
４９（１２）： １０－１７ （ｉｎ Ｃｈｉｎｅｓｅ）
［２４］　 Ｈａｒｍｓ ＫＥ， Ｃｏｎｄｉｔ Ｒ， Ｈｕｂｂｅｌｌ ＳＰ， ｅｔ ａｌ． Ｈａｂｉｔａｔ ａｓｓｏ⁃
ｃｉａｔｉｏｎｓ ｏｆ ｔｒｅｅｓ ａｎｄ ｓｈｒｕｂｓ ｉｎ ａ ５０⁃ｈａ ｎｅｏｔｒｏｐｉｃａｌ ｆｏｒｅｓｔ
ｐｌｏｔ． Ｊｏｕｒｎａｌ ｏｆ Ｅｃｏｌｏｇｙ， ２００１， ８９： ９４７－９５９
［２５］　 Ｔｈｅ Ｓｔａｔｅ Ｆｏｒｅｓｔｒｙ Ａｄｍｉｎｉｓｔｒａｔｉｏｎ （国家林业局）． Ｔｈｅ
Ｄｅｔｅｒｍｉｎａｔｉｏｎ ｏｆ Ｔｏｔａｌ Ｐｈｏｓｐｈｏｒｕｓ ｉｎ Ｓｏｉｌ Ｆｏｒｅｓｔ． Ｆｏｒｅｓ⁃
ｔｒｙ Ｉｎｄｕｓｔｒｙ Ｓｔａｎｄａｒｄ ｏｆ ｔｈｅ Ｐｅｏｐｌｅ’ｓ Ｒｅｐｕｂｌｉｃ ｏｆ Ｃｈｉｎａ
（ＬＹ ／ Ｔ １２３２－１９９９）： Ｆｏｒｅｓｔ Ｓｏｉｌ Ａｎａｌｙｓｉｓ Ｍｅｔｈｏｄ． Ｂｅｉ⁃
ｊｉｎｇ： Ｃｈｉｎａ Ｓｔａｎｄａｒｄｓ Ｐｒｅｓｓ， ２０００： ８７ － ９０ （ ｉｎ Ｃｈｉ⁃
ｎｅｓｅ）
［２６］ 　 Ｚｈａｎｇ Ｃ⁃Ｓ （张朝生）， Ｚｈａｎｇ Ｓ （章 　 申）， Ｈｅ Ｊ⁃Ｂ
（何建邦）． Ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｈｅａｖｙ
ｍｅｔａｌｓ ｉｎ ｔｈｅ ｓｅｄｉｍｅｎｔｓ ｏｆ Ｃｈａｎｇｊｉａｎｇ Ｒｉｖｅｒ ｓｙｓｔｅｍ． Ａｃｔａ
Ｇｅｏｇｒａｐｈｉｃａ Ｓｉｎｉｃａ （地理学报）， １９９７， ５２（２）： １８４－
１９２ （ｉｎ Ｃｈｉｎｅｓｅ）
［２７］ 　 Ｗａｎｇ Ｚ⁃Ｑ （王政权）． Ｓｔａｔｉｓｔｉｃｓ ａｎｄ Ｉｔｓ Ａｐｐｌｉｃａｔｉｏｎ ｉｎ
Ｅｃｏｌｏｇｙ． Ｂｅｉｊｉｎｇ： Ｓｃｉｅｎｃｅ Ｐｒｅｓｓ， １９９９ （ｉｎ Ｃｈｉｎｅｓｅ）
［２８］　 Ｆａｎ Ｆ⁃Ｊ （范夫静）， Ｓｏｎｇ Ｔ⁃Ｑ （宋同清）， Ｈｕａｎｇ Ｇ⁃Ｑ
（黄国勤）， ｅｔ ａｌ． Ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｓｐａｔｉａｌ ｖａｒｉａｔｉｏｎ ｏｆ
ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｉｎ ｓｌｏｐｉｎｇ ｆｉｅｌｄ ｉｎ ａ ｇｏｒｇｅ Ｋａｒｓｔ ｒｅｇｉｏｎ，
ｓｏｕｔｈｗｅｓｔ Ｃｈｉｎａ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应
用生态学报）， ２０１４， ２５（１）： ９２－９８ （ｉｎ Ｃｈｉｎｅｓｅ）
［２９］　 Ｔａｎｇ Ｔ （唐　 涛）， Ｃａｉ Ｑ⁃Ｈ （蔡庆华）， Ｐａｎ Ｗ⁃Ｂ （潘
文斌）． Ａｐｐｌｉｃａｔｉｏｎ ｏｆ ｇｅｏｓｔａｔｉｓｔｉｃｓ ｉｎ ｆｒｅｓｈｗａｔｅｒ ｅｃｏｌｏｇｙ．
Ｊｏｕｒｎａｌ ｏｆ Ｌａｋｅ Ｓｃｉｅｎｃｅｓ （湖泊科学）， ２０００， １２（３）：
２８０－２８８ （ｉｎ Ｃｈｉｎｅｓｅ）
［３０］　 Ｒｅｎａｔｏ Ｖ， Ｒｏｂｉｎ ＢＦ， Ｇｏｒｋｙ Ｖ， ｅｔ ａｌ． Ｔｒｅｅ ｓｐｅｃｉｅｓ ｄｉｓ⁃
ｔｒｉｂｕｔｉｏｎｓ ａｎｄ ｌｏｃａｌ ｈａｂｉｔａｔ ｖａｒｉａｔｉｏｎ ｉｎ ｔｈｅ Ａｍａｚｏｎ：
Ｌａｒｇｅ ｆｏｒｅｓｔ ｐｌｏｔ ｉｎ ｅａｓｔｅｒｎ Ｅｃｕａｄｏｒ． Ｊｏｕｒｎａｌ ｏｆ Ｅｃｏｌｏｇｙ，
２００４， ９２： ２１４－２２９
［３１］　 Ｄｉｎｇ Ｊ （丁 　 佳）， Ｗｕ Ｑ （吴 　 茜）， Ｙａｎ Ｈ （闫 　
慧）， ｅｔ ａｌ． Ｅｆｆｅｃｔｓ ｏｆ ｔｏｐｏｇｒａｐｈｉｃ ｖａｒｉａｔｉｏｎｓ ａｎｄ ｓｏｉｌ
ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｎ ｐｌａｎｔ ｆｕｎｃｔｉｏｎａｌ ｔｒａｉｔｓ ｉｎ ａ ｓｕｂｔｒｏｐｉｃａｌ
ｅｖｅｒｇｒｅｅｎ ｂｒｏａｄ⁃ｌｅａｖｅｄ ｆｏｒｅｓｔ． Ｂｉｏｄｉｖｅｒｓｉｔｙ Ｓｃｉｅｎｃｅ （生物
多样性）， ２０１１， １９（２）： １５８－１６７ （ｉｎ Ｃｈｉｎｅｓｅ）
［３２］　 Ｚｈａｎｇ Ｚ⁃Ｈ （张忠华）， Ｈｕ Ｇ （胡　 刚）， Ｚｈｕ Ｊ⁃Ｄ （祝
介东）， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｈｅｔｅｒｏｇｅｎｅｉｔｙ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ａｎｄ
ｉｔｓ ｉｍｐａｃｔ ｏｎ ｔｒｅｅ ｓｐｅｃｉｅｓ ｄｉｓｔｒｉｂｕｔｉｏｎ ｉｎ ａ Ｋａｒｓｔ ｆｏｒｅｓｔ ｏｆ
Ｓｏｕｔｈｗｅｓｔ Ｃｈｉｎａ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ｐｌａｎｔ Ｅｃｏｌｏｇｙ （植
物生态学报）， ２０１１， ３５（１０）： １０３８－ １０４９ （ ｉｎ Ｃｈｉ⁃
ｎｅｓｅ）
［３３］　 Ｃａｏ Ｊ （曹 　 娟）， Ｙａｎ Ｗ⁃Ｄ （闫文德）， Ｘｉａｎｇ Ｗ⁃Ｈ
４１７ 应　 用　 生　 态　 学　 报　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 ２７卷
（项文化）， ｅｔ ａｌ． Ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｓｏｉｌ ｐｈｏｓｐｈｏｒｕｓ ｉｎ
ｄｉｆｆｅｒｅｎｔ ａｇｅｄ ｓｔａｎｄｓ ｏｆ Ｃｈｉｎｅｓｅ ｆｉｒ ｐｌａｎｔａｔｉｏｎｓ ｉｎ
Ｈｕｉｔｏｎｇ， Ｈｕｎａｎ Ｐｒｏｖｉｎｃｅ． Ａｃｔａ Ｅｃｏｌｏｇｉｃａ Ｓｉｎｉｃａ （生态
学报）， ２０１４， ３４（２２）： ６５１９－６５２７ （ｉｎ Ｃｈｉｎｅｓｅ）
［３４］　 Ｙａｎｇ Ｘ⁃Ｙ （杨小燕）， Ｆａｎ Ｒ⁃Ｙ （范瑞英）， Ｗａｎｇ Ｅ⁃Ｈ
（王恩姮）， ｅｔ ａｌ． Ｔｏｐｓｏｉｌ ｐｈｏｓｐｈｏｒｕｓ ｆｏｒｍｓ ａｎｄ ａｖａｉｌａ⁃
ｂｉｌｉｔｙ ｏｆ ｄｉｆｆｅｒｅｎｔ ｓｏｉｌ ａｎｄ ｗａｔｅｒ ｃｏｎｓｅｒｖａｔｉｏｎ ｐｌａｎｔａｔｉｏｎｓ
ｉｎ ｔｙｐｉｃａｌ ｂｌａｃｋ ｓｏｉｌ ｒｅｇｉｏｎ， Ｎｏｒｔｈｅａｓｔ Ｃｈｉｎａ． Ｃｈｉｎｅｓｅ
Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应用生态学报）， ２０１４， ２５
（６）： １５５５－１５６０ （ｉｎ Ｃｈｉｎｅｓｅ）
［３５］　 Ｌｉｕ Ｗ⁃Ｊ （刘文杰）， Ｃｈｅｎ Ｓ⁃Ｙ （陈生云）， Ｈｕ Ｆ⁃Ｚ （胡
凤祖）， ｅｔ ａｌ． Ｄｉｓｔｒｉｂｕｔｉｏｎｓ ｐａｔｔｅｒｎ ｏｆ ｐｈｏｓｐｈｏｒｕｓ， ｐｏ⁃
ｔａｓｓｉｕｍ ａｎｄ ｉｎｆｌｕｅｎｃｉｎｇ ｆａｃｔｏｒｓ ｉｎ ｔｈｅ ｕｐｓｔｒｅａｍ ｏｆ Ｓｈｕｌｅ
ｒｉｖｅｒ ｂａｓｉｎ． Ａｃｔａ Ｅｃｏｌｏｇｉｃａ Ｓｉｎｉｃａ （生态学报）， ２０１２，
３２（１７）： ５４２９－５４３７ （ｉｎ Ｃｈｉｎｅｓｅ）
［３６］　 Ｎａｔｉｏｎａｌ Ｓｏｉｌ Ｓｕｒｖｅｙ Ｏｆｆｉｃｅ （全国土壤普查办公室）．
Ｓｏｉｌ Ｓｕｒｖｅｙ ａｎｄ Ｔｅｃｈｎｏｌｏｇｙ ｏｆ Ｃｈｉｎａ． Ｂｅｉｊｉｎｇ： Ｃｈｉｎａ Ａｇ⁃
ｒｉｃｕｌｔｕｒｅ Ｐｒｅｓｓ， １９９２ （ｉｎ Ｃｈｉｎｅｓｅ）
［３７］　 Ｙａｎｇ Ｘ⁃Ｑ （杨秀清）， Ｈａｎ Ｙ⁃Ｚ （韩有志）， Ｌｉ Ｌ （李
乐）， ｅｔ ａｌ． Ｔｈｅ ｅｆｆｅｃｔ ｏｆ ｈｅｔｅｒｏｇｅｎｅｏｕｓ ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕ⁃
ｔｉｏｎ ｏｆ ｓｏｉｌ ｎｉｔｒｏｇｅｎ ｏｎ ｒｅｇｅｎｅｒａｔｉｏｎ ｏｆ Ｌａｒｉｘ ｐｒｉｎｃｉｐｉｓ⁃
ｒｕｐｐｒｅｃｈｔｉｉ ｓｅｅｄｌｉｎｇｓ ｉｎ ｔｙｐｉｃａｌ ｎａｔｕｒａｌｌｙ⁃ｒｅｇｅｎｅｒａｔｅｄ
ｍｏｎｔａｎｅ ｆｏｒｅｓｔｓ ｏｆ Ｎｏｒｔｈｅｒｎ Ｃｈｉｎａ． Ａｃｔａ Ｅｃｏｌｏｇｉｃａ Ｓｉｎｉｃａ
（生态学报）， ２００９， ２９（９）： ４６５６－４６６４ （ｉｎ Ｃｈｉｎｅｓｅ）
［３８］　 Ｌｉ Ｙ （李　 渊）， Ｇｏｎｇ Ｙ⁃Ｂ （宫渊波）， Ｓｕ Ｈ⁃Ｗ （苏宏
伟）， ｅｔ ａｌ． Ｓｏｉｌ ｐｈｏｓｐｈｏｒｕｓ ｖａｒｉａｔｉｏｎ ｏｆ Ｍａｓｓｏｎ ｐｉｎｅ
ｐｌａｎｔａｔｉｏｎ ｉｎ Ｓｏｕｔｈｅｒｎ Ｓｉｃｈｕａｎ． Ｊｏｕｒｎａｌ ｏｆ Ｎｏｒｔｈｅａｓｔ Ｆｏｒｅ⁃
ｓｔｒｙ Ｕｎｉｖｅｒｓｉｔｙ （东北林业大学学报）， ２０１４， ４２（６）：
６３－６７ （ｉｎ Ｃｈｉｎｅｓｅ）
［３９］　 Ｓｏｎｇ Ｘ （宋　 轩）， Ｌｉ Ｌ⁃Ｄ （李立东）， Ｋｏｕ Ｃ⁃Ｌ （寇长
林）， ｅｔ ａｌ． Ｓｏｉｌ ｎｕｔｒｉｅｎｔ ｄｉｓｔｒｉｂｕｔｉｏｎ ａｎｄ ｉｔｓ ｒｅｌａｔｉｏｎｓ
ｗｉｔｈ ｔｏｐｏｇｒａｐｈｙ ｉｎ Ｈｕａｎｇｓｈｕｉ Ｒｉｖｅｒ ｄｒａｉｎａｇｅ ｂａｓｉｎ． Ｃｈｉ⁃
ｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应用生态学报），
２０１１， ２２（１２）： ３１６３－３１６８ （ｉｎ Ｃｈｉｎｅｓｅ）
［４０］　 Ｇｕｏ Ｘ⁃Ｄ （郭旭东）， Ｆｕ Ｂ⁃Ｊ （傅伯杰）， Ｍａ Ｋ⁃Ｍ （马
克明）， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｂａｓｅｄ
ｏｎ ｇｅｏｓｔａｔｉｓｔｉｃｓ ｃｏｍｂｉｎｅｄ ｗｉｔｈ ＧＩＳ： Ａ ｃａｓｅ ｓｔｕｄｙ ｉｎ
Ｚｕｎｈｕａ Ｃｉｔｙ ｏｆ Ｈｅｂｅｉ Ｐｒｏｖｉｎｃｅ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐ⁃
ｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应用生态学报）， ２０００， １１（４）： ５５７－
５６３ （ｉｎ Ｃｈｉｎｅｓｅ）
［４１］　 Ｚｈｏｕ Ｈ⁃Ｚ （周慧珍）， Ｇｏｎｇ Ｚ⁃Ｔ （龚子同）， Ｌａｍｐ Ｊ．
Ｓｔｕｄｙ ｏｎ ｓｏｉｌ ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ． Ａｃｔａ Ｐｅｄｏｌｏｇｉｃａ Ｓｉｎｉｃａ
（土壤学报）， １９９６， ３３（２）： ２３２－２４１ （ｉｎ Ｃｈｉｎｅｓｅ）
［４２］　 Ｔａｔｅｎｏ Ｒ， Ｔａｋｅｄａ Ｈ． Ｆｏｒｅｓｔ ｓｔｒｕｃｔｕｒｅ ａｎｄ ｔｒｅｅ ｓｐｅｃｉｅｓ
ｄｉｓｔｒｉｂｕｔｉｏｎ ｉｎ ｒｅｌａｔｉｏｎ ｔｏ ｔｏｐｏｇｒａｐｈｙ⁃ｍｅｄｉａｔｅｄ ｈｅｔｅｒｏｇｅ⁃
ｎｅｉｔｙ ｏｆ ｓｏｉｌ ｎｉｔｒｏｇｅｎ ａｎｄ ｌｉｇｈｔ ａｔ ｔｈｅ ｆｏｒｅｓｔ ｆｌｏｏｒ． Ｅｃｏｌｏｇｉ⁃
ｃａｌ Ｒｅｓｅａｒｃｈ， ２００３， １８： ５５９－５７１
［４３］　 Ｌｉｕ Ｌ （刘　 璐）， Ｚｅｎｇ Ｆ⁃Ｐ （曾馥平）， Ｓｏｎｇ Ｔ⁃Ｑ （宋
同清）， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｈｅｔｅｒｏｇｅｎｅｉｔｙ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｉｎ
Ｋａｒｓｔ ａｒｅａ’ ｓ Ｍｕｌｕｎ Ｎａｔｉｏｎａｌ Ｎａｔｕｒｅ Ｒｅｓｅｒｖｅ． Ｃｈｉｎｅｓｅ
Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应用生态学报）， ２０１０， ２１
（７）： １６６７－１６７３ （ｉｎ Ｃｈｉｎｅｓｅ）
［４４］　 Ｃｈａｎｇ Ｃ （常 　 超）， Ｘｉｅ Ｚ⁃Ｑ （谢宗强）， Ｘｉｏｎｇ Ｇ⁃Ｍ
（熊高明）， ｅｔ ａｌ． Ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｏｆ ｄｉｆ⁃
ｆｅｒｅｎｔ ｖｅｇｅｔａｔｉｏｎ ｔｙｐｅｓ ｉｎ ｔｈｅ Ｔｈｒｅｅ Ｇｏｒｇｅｓ ｒｅｓｅｒｖｏｉｒ ａｒｅａ．
Ａｃｔａ Ｅｃｏｌｏｇｉｃａ Ｓｉｎｉｃａ （生态学报）， ２００９， ２９ （１１）：
５９７８－５９８５ （ｉｎ Ｃｈｉｎｅｓｅ）
［４５］　 Ｌｉｎ Ｂ （林　 波）， Ｌｉｕ Ｑ （刘　 庆）， Ｗｕ Ｙ （吴　 彦），
ｅｔ ａｌ． Ｅｆｆｅｃｔ ｏｆ ｆｏｒｅｓｔ ｌｉｔｔｅｒｓ ｏｎ ｓｏｉｌ ｐｈｙｓｉｃａｌ ａｎｄ ｃｈｅｍｉｃａｌ
ｐｒｏｐｅｒｔｉｅｓ ｉｎ ｓｕｂａｌｐｉｎｅ ｃｏｎｉｆｅｒｏｕｓ ｆｏｒｅｓｔｓ ｏｆ Ｗｅｓｔｅｒｎ Ｓｉ⁃
ｃｈｕａｎ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ ＆ Ｅｎｖｉｒｏｎｍｅｎｔａｌ Ｂｉｏ⁃
ｌｏｇｙ （应用与环境生物学报）， ２００３， ９（４）： ３４６－３５１
（ｉｎ Ｃｈｉｎｅｓｅ）
［４６］ 　 Ｘｕ Ｗ⁃Ｍ （徐武美）， Ｓｏｎｇ Ｃ⁃Ｙ （宋彩云）， Ｌｉ Ｑ⁃Ｍ
（李巧明）． Ｒｅｌａｔｉｏｎｓｈｉｐ ｂｅｔｗｅｅｎ ｓｏｉｌ ｒｅｓｏｕｒｃｅ ｈｅｔｅｒｏｇｅ⁃
ｎｅｉｔｙ ａｎｄ ｔｒｅｅ ｄｉｖｅｒｓｉｔｙ ｉｎ Ｘｉｓｈｕａｎｇｂａｎｎａ ｔｒｏｐｉｃａｌ ｓｅａ⁃
ｓｏｎａｌ ｒａｉｎｆｏｒｅｓｔ， Ｓｏｕｔｈｗｅｓｔ Ｃｈｉｎａ． Ａｃｔａ Ｅｃｏｌｏｇｉｃａ Ｓｉｎｉｃａ
（生态学报）， ２０１５， ３５（２３）： ７７５６ － ７７６２ （ ｉｎ Ｃｈｉ⁃
ｎｅｓｅ）
［４７］　 Ｌｕｎｄｈｏｌｍ ＪＴ． Ｐｌａｎｔ ｓｐｅｃｉｅｓ ｄｉｖｅｒｓｉｔｙ ａｎｄ ｅｎｖｉｒｏｎｍｅｎｔａｌ
ｈｅｔｅｒｏｇｅｎｅｉｔｙ： Ｓｐａｔｉａｌ ｓｃａｌｅ ａｎｄ ｃｏｍｐｅｔｉｎｇ ｈｙｐｏｔｈｅｓｅｓ．
Ｊｏｕｒｎａｌ ｏｆ Ｖｅｇｅｔａｔｉｏｎ Ｓｃｉｅｎｃｅ， ２００９， ２０： ３７７－３９１
［４８］　 Ｋｏｎｇ Ｑ⁃Ｂ （孔庆波）， Ｂａｉ Ｙ⁃Ｌ （白由路）， Ｙａｎｇ Ｌ⁃Ｐ
（杨俐苹）， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ ｃｈａｒａｃｔｅｒｉｓｔｉｃ ａｎｄ
ｉｔｓ ｉｎｆｌｕｅｎｔｉａｌ ｆａｃｔｏｒｓ ｏｆ ｓｏｉｌ ｐｈｏｓｐｈｏｒｕｓ ｉｎ ｒｅｇｉｏｎ ｏｆ ｔｈｅ
Ｈｕａｎｇ⁃Ｈｕａｉ⁃Ｈａｉ ｐｌａｉｎ． Ｓｏｉｌ ａｎｄ Ｆｅｒｔｉｌｉｚｅｒ Ｓｃｉｅｎｃｅｓ ｉｎ
Ｃｈｉｎａ （中国土壤与肥料）， ２００９（５）： １０－１４ （ｉｎ Ｃｈｉ⁃
ｎｅｓｅ）
作者简介　 胡瑞彬，男，１９９０年生，硕士研究生． 主要从事环
境生态学研究． Ｅ⁃ｍａｉｌ： ２４００６０９６１＠ ｑｑ．ｃｏｍ
责任编辑　 孙　 菊
胡瑞彬， 方晰， 项文化， 等． 中亚热带南酸枣落叶阔叶林土壤磷素空间变异及其影响因素． 应用生态学报， ２０１６， ２７（３）：
７０５－７１５
Ｈｕ Ｒ⁃Ｂ， Ｆａｎｇ Ｘ， Ｘｉａｎｇ Ｗ⁃Ｈ， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｈｅｔｅｒｏｇｅｎｅｉｔｙ ａｎｄ ｉｎｆｌｕｅｎｃｉｎｇ ｆａｃｔｏｒｓ ｏｆ ｓｏｉｌ ｐｈｏｓｐｈｏｒｕｓ ｃｏｎｃｅｎｔｒａｔｉｏｎ ｉｎ ａ ｍｉｄ⁃ｓｕｂｔｒｏｐｉｃａｌ
Ｃｈｏｅｒｏｓｐｏｎｄｉａｓ ａｘｉｌｌａｒｉｓ ｄｅｃｉｄｕｏｕｓ ｂｒｏａｄ⁃ｌｅａｖｅｄ ｆｏｒｅｓｔ， Ｃｈｉｎａ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ， ２０１６， ２７（３）： ７０５－７１５ （ｉｎ Ｃｈｉｎｅｓｅ）
５１７３期　 　 　 　 　 　 　 　 　 　 胡瑞彬等： 中亚热带南酸枣落叶阔叶林土壤磷素空间变异及其影响因素