土壤是一种形态和演化过程都十分复杂的自然综合体, 受气候、生物、母质、地形、成土时间等成土因素的影响, 具有复杂性和时空变异性^[1]。即使在土壤类型和质地相同的区域内, 同一时刻土壤特性在空间上也有明显的差异^[2]。土壤养分是土壤提供的植物生长发育所必需的营养元素, 是土壤肥力的物质基础, 也是土地评价和管理的重要指标之一。由于受到自然因素和人为因素的共同作用, 土壤养分也具有高度的空间变异性^[3-6]。土壤P、K不仅是植物生长发育所必需的营养元素, 也是影响区域水体生态环境的重要属性, 因此, 研究其时空变异对喀斯特生态系统的恢复和重建等具有重要的理论与现实指导意义, 以实现喀斯特地区的生态、经济、社会的和谐和可持续发展^[7]。

一些学者^[8-9]利用土壤学家Fisher所创立的传统统计方法对土壤性状的空间分布进行了定量描述, 但由于此方法要求样本的取值相互独立, 给土壤空间变异的研究带来了很大的局限性。20世纪60年代Matheron创立了地统计学理论^[10], 为土壤空间变异性研究提供了有效的方法, 掀起了国外学者对土壤空间变异性研究的热潮^[11-13]。目前, 土壤空间变异性研究已成为土壤科学的重要研究内容和热点^[14-16]。如：Corss和Schlesinger^[17]认为小尺度上的土壤异质性分布与群落优势种分布一致；Farlely等^[18]在研究阔叶林地时发现, 在20 cm小尺度上, 土壤的养分含量差异非常大, 同时也有研究表明^[19-21]土层深度对森林土壤养分含量的影响差异显著, 且随着土层的加深土壤有机质和速效钾的含量降低；宋同清等^[22]认为中小尺度上微生境土壤理化特性等因素可能是推动植物群落结构变化的主力军。有关土壤养分空间分布的研究颇多, 但国内外的研究多数集中在小尺度、单因素^[23-25], 而在较大尺度上的喀斯特土壤养分空间分布研究较少^[26-38]。喀斯特地区人口密度高、增长速度快、可利用土地少、生产方式落后等不利因素使得该区域生态系统承受压力大, 水土流失严重, 基岩裸露率高, 最终导致愈来愈严重的石漠化现象。本文以广西木论国家级自然保护区典型森林(常绿落叶阔叶混交林)土壤为研究对象, 应用经典统计和地统计法学方法分析了该森林土壤磷钾养分TP、TK、AP、AK的空间格局、相关性与变异性等特征；同时, 通过Pearson相关性分析该区域内坡度、坡向、土层深度等环境因子对土壤养分分布的影响, 旨在为喀斯特地区的土壤养分库的评价、生态保护和恢复研究提供理论支持。

1 研究区概况与研究方法 1.1 研究区概况

研究区位于广西环江毛南族自治县西北部木论国家级自然保护区25°07′01″—25°12′22″N, 107°54′01″—108°05′51″E, 南缘云贵高原, 北连贵州茂兰国家级自然保护区, 总面积8969 hm², 海拔400—1000 m。该区域为中亚热带季风气候, 年平均气温15—18.7℃, 最冷月(1月)平均气温3.4—8.7℃, 最热月(7月)平均气温23—26.7℃, ≥10℃, 年积温4700—6300℃；年降雨量1530—1820 mm, 集中于4—8月, 年平均相对湿度为80%—90%, 无霜期235—290 d, 气候温暖, 雨量充沛。土壤主要为石灰土和零星分布的硅质土, 河流及水时空分布不均, 地域性来水差异大。研究区植被为中亚热带石灰岩常绿落叶阔叶混交林, 属隐域性喀斯特森林植被顶极群落类型, 是目前世界上同纬度地区残存下来的仅有的、原生性强、相对稳定的喀斯特森林生态系统, 也是喀斯特区原生性森林分布面积最大的区域。植被成层现象比较明显。

1.2 样品采集与分析

通过勘测, 于2014年3月在研究区内设置一块500 m×500 m的大样地, 作为长期监测样地(图 1)。参照CTFS标准采用规则网格法(20 m×20 m)进行样点布置, 用全站仪和GPS基站相结合的测量方法设置和划分样地, 参考木论自然保护区的森林状况和植被、地形分布状况, 共设置625个小样方, 样方4个角用水泥桩作永久标记。在每个20 m×20 m样方内, 按照其中心点半径2 m范围内五点取样法, 用直径5 cm的土钻取表层(0—10 cm)土壤5个混合成一个样品, 代表该样点样本。同时记录各样方坡向、坡位、海拔、土层厚度等相关立地因子。将所有混合土样带回实验室风干、过筛, 以备分析测定全磷(TP)、全钾(TK)、速效磷(AP)、速效钾(AK)。坡向视为坡度方向或山体所面对的罗盘方向, 以度为单位按逆时针方向进行测量。岩石裸露率以样点周围2米范围内岩石出露的面积估算。凹凸度为样方的海拔减去该样方相邻的8个样方海拔的平均值。全磷(TP)的测定采用氢氧化钠熔融—钼锑抗显色—紫外分光光度法；全钾(TK)的测定采用氢氧化钠熔融—火焰原子吸收分光光度法；速效磷(AP)的测定采用0.5 mol/L碳酸氢钠浸提—钼锑抗显色—紫外分光光度法；速效钾(AK)的测定采用乙酸铵浸提—火焰原子吸收分光光度法^[29]。

1.3 数据处理

使用SPSS 18.0(统计分析与检验、Pearson相关性)、Minitab 16(Box-Cox法)、GS+(地统计分析)、ArcGIS 9.3(克里格插值图)等软件, 采用经典统计学和地统计学方法对样本数据进行分析处理。

1.3.1 半方差函数分析

半方差函数是地统计学研究中特有的基本工具, 用来描述属性变量在空间上的差异程度, 它是距离的函数, 也是应用最广泛的空间格局描述工具, 其公式如下:

(1)

式中, γ(h)为半方差函数值；h为两样本点的空间距离；N(h)为间距等于h的点对数；Z(x_i)为区域化变量Z在位置x_i的实测值；Z(x_i+h)为与x_i距离为h处样点的值。一般认为半方差函数只有在最大间隔的l/2内才有意义。块金值(C₀)、基台值(C₀+C)和变程(A)为半方差分析的重要参数, 块金值与基台值的比值用来表示空间变异程度, 一般认为块基比C₀/(C₀+C) < 25%时, 空间变量为强烈的空间自相关, 在25%—75%之间时, 为中等空间自相关, >75%为弱空间自相关。本研究中所有半方差函数的有效滞后距均设为最大滞后距的1/2, 采用指数模型进行拟合。

1.3.2 克里格插值

克里格插值是建立在半方差函数理论的基础上, 在一定区域内对变量的取值进行无偏最优估计的一种方法, 也是最简单、应用最广泛的一种插值方法。在普通克里格插值中, 对于任意待估点(块段) x₀的实际值Z(x₀), 其估计值Z^#(x0)通过该待估点(块段)周围的n个有效样点Z(x_i)的线性组合得到。

(2)

式中, λ_i为权重系数。在满足无偏条件下, 要使预测值的估计方差最小, 即:

(3)

2 结果与分析 2.1 土壤磷钾养分经典统计描述

采用样本均值±3倍标准差来识别特异值, 在此区间外的数据均定为特异值, 分别用最大和最小值代替, 后续计算均采用处理后的原始数据。由表 1可知, 土壤中TP、TK、AP、AK的含量分别为(1.60±0.76) g/kg、(5.42±2.74) g/kg、(5.74±3.63) mg/kg、(5.20±2.96) mg/kg。从变异系数看, 土壤磷钾养分的变异系数范围为47.66%—63.22%, 均属中等变异水平。其中变异系数最大的指标为AP, 为63.22%。通过对偏度、峰度的观测以及采用K-S法对4种土壤磷钾养分数据进行非参数检验, 在5%的检验水平下, 结果均不符合正态分布。用Box-Cox法对不符合的样本进行转换, 转换后数据再进行单样本Kolmogorov-Smirnov检验, 结果表明均符合正态分布。

2.2 土壤磷钾养分空间自相关分析

Moran′s I的数值大小反应了土壤养分指标的空间自相关性大小。由图 2可以看出, 研究区土壤磷钾养分呈现出一定的结构性, TP与AP之间具有相似的空间结构。AP的Moran′s I系数在-0.124—0.497之间, I系数达0.497, 表明相对于其他养分指标, AP的对空间的依赖程度较大；TP、TK、AK的Moran′s I在-0.043—0.313之间；各养分空间自相关的排序大小为：AP>TP>AK>TK。所有养分指标的Moran′s I系数均随着滞后距的增大而减小, 在220 m左右出现拐点, Moran′s I下降为负相关并趋于平缓, 说明在该研究区220 m范围内森林原生性较好。

2.3 土壤磷钾养分半方差函数

土壤特性空间异质性的经典统计分析中将土壤特性看作随机的、彼此相互独立的。然而实际中, 在一定的空间范围内, 土壤性状存在一定的空间依赖性。由表 2和图 3可看出, 研究区土壤磷钾养分TP、TK、AP、AK变异函数值的最佳拟合模型均为指数模型, 决定系数R²均很高, 在0.671—0.995之间, 能很好地反映土壤磷钾养分的空间结构特征。研究区各土壤磷钾养分指标的块金值/基台值C₀/(C₀+C)依次为TK>AK>AP>TP, 其中TP、AP块金值/基台值C₀/(C₀+C)为50.4%、64.4%, 空间自相关中等, 主要受结构因子影响；TK、AK块金值/基台值分别为88.4%、86.6%, 为弱空间自相关。变程反映土壤性状的有效空间自相关距离, 研究区TP、AP的变程分别为336.00、373.50 m, 表明空间连续性较好；TK、AK的变程在33.30—64.50 m之间, 空间依赖性较强。

2.4 土壤磷钾养分空间分布格局

为了能更直观地反映研究区土壤参数的空间分布特征, 根据半方差函数模型Arc GIS 9.2软件进行Kriging插值分析, 分别绘制全磷(TP)、全钾(TK)、速效磷(AP)、速效钾(AK)4种土壤养分的空间分布图(图 4)。从图中可看出, 研究区土壤养分TP表现为坡下(含洼地)含量高, 坡上含量较低, 在偏西北方向出现最大值；AK表现为中坡含量高于洼地含量；AP呈斑块状分布, 规律性不强, 样地正北方向出现最大值；TK斑块破碎化程度较高。

2.5 土壤磷钾养分空间变异的影响因素

经过Pearson相关性分析计算4个土壤磷钾养分指标的相关性矩阵(表 3), 可以看出研究区4种养分指标TP、TK、AP、AK均与海拔呈负相关, 其中TP、TK、AP与海拔呈极显著负相关；各指标与坡向和地面凹凸度均为负相关, 其中坡向与TP呈极显著负相关关系, 地面凹凸度与TP、TK呈极显著负相关；土层深度与AK呈极显著负相关；岩石裸露率与AP呈显著负相关, 与TP、TK、AK呈正相关关系；坡度与TP、AP呈极显著负相关。说明土壤磷钾养分主要受海拔、坡度、地面凹凸度的影响。海拔较低、坡度较缓的点, 土壤磷钾养分的含量也较高。

3 讨论

研究区位于木论自然保护区缓冲带内, 为一典型的峰丛洼地单元, 森林保存完好、人为干扰少、植被多样性高、地形地貌复杂多变、微生境丰富^[7]。经典描述性统计分析发现, 研究区土壤各磷钾养分指标差异较大, 这可能与喀斯特地区土壤的高异质性及营养元素的淋溶特征密切相关。喀斯特常绿落叶阔叶林土壤磷钾养分变异系数在47.66%—63.22%之间, 均为中等变异；TP、TK、AP、AK养分含量的变化范围与刘璐等^[30]的研究结果相符, 即原生性较好的喀斯特森林的养分水平皆呈现出较高的水平, 同时也与自然保护区内的生态结构功能和保护区的顶极群落优势密切相关。研究区常绿落叶阔叶混交林土壤因子的基本特征与其他亚热带森林有所不同。与浙江常绿阔叶林土壤相比^[31], 本研究区土壤TP含量均值高于浙江常绿阔叶林水平；与鄂西南亚热带常绿落叶阔叶混交林相比^[32], 研究区土壤TP、AP较高, AK较低。

基于前期对喀斯特土壤主要养分和水分的空间异质性研究^[33-36], 进一步对喀斯特常绿落叶阔叶混交林土壤磷钾养分含量进行分析发现, 研究区土壤磷钾养分呈现出一定的结构性, 所有养分指标的Moran′s I系数均随着滞后距的增大而减小。研究区TP、TK、AP、AK变异函数值的最佳拟合模型均为指数模型, 对土壤空间变异具有极高的拟合度, 表明研究区土壤具有良好的空间结构。土壤性状的空间异质性是结构因子和随机因子共同作用的结果。影响土壤空间变异的结构因素为母质、气候、地形等自然地理要素, 随机因素则为取样设计、测定误差与干扰等人为原因。喀斯特常绿落叶阔叶林土壤TP、AP呈空间自相关中等, 说明其变异是由结构性变异(环境因子等)和随机因素共同引起的；TK、AK为弱空间自相关, 说明随机因素对这两项土壤养分指标空间分布的的影响较大, 可能与干扰、试验误差有关。研究区TP、AP的变程较高, 空间连续性较好；TK、AK空间依赖性较强。与湖北^[32]及浙江常绿落叶阔叶混交林^[37]的同类研究相比具有相同点与不同点：相同点在于均发现土壤AP的变异系数相对其他因子较高, 这表明相比其他养分指标, 土壤AP易受其他因素影响；不同点在于本文发现土壤AP的半方差变程最大, 而湖北与浙江同类森林中土壤AP较其他因子的半方差拟合的变程相对较小, 这可能与喀斯特生态系统具有独特的二元水文结构及复杂的微生境有关, 生态系统的结构与功能对土壤空间分布格局产生重要影响。

空间异质性导致空间格局的存在, Kriging等值线图能更深刻、全面和直观地反映土壤磷钾的空间分布。研究区土壤TP表现为坡下(含洼地)含量高, 坡上含量较低, 这是由于土壤中的P主要来源于基岩的风化, 并通过植物的表聚作用在表层土壤积累, 因而凋落物P是表层土壤P的直接来源^[38-39]；潘复静等^[40]研究表明, 在峰丛洼地, 凋落物P沿坡面具有明显的空间分异规律, 随坡位的升高而降低, 这也是本研究土壤TP坡上含量较低的主要原因。AP、TK呈斑块破碎化, 这可能与土壤的淋溶特征有关。由于土壤养分受多种因子的综合作用, 不同养分的破碎化程度不一样, 还需进一步分析地形因子在土壤养分变化中的作用, 以便对土壤养分空间格局机理进行合理的解释。

环境因子对喀斯特森林土壤磷钾养分空间分布特征和变异规律起着非常重要的作用^[41]。研究区TP、TK、AP、AK均与海拔、坡向和地面凹凸度呈负相关关系；AK与土层深度极显著负相关；TP、AK与坡度呈极显著负相关关系；岩石裸露率与AP负相关。海拔、坡度和地面凹凸度对土壤磷钾养分的影响最大, 说明海拔越高、坡度越陡的点, 土壤磷钾养分的含量较低。这可能是由于淋溶作用强烈、植物的表聚作用以及凋落物相关。

以贵州为中心连带成片的中国西南喀斯特地区(世界三大岩溶地区之一)受地球内动力、强烈的地质运动、高温多雨且分布不均、碳酸盐岩溶蚀性强、水文二维结构明显以及其适生植物具有嗜钙性、耐旱性和石生性等限制特点的综合影响, 生境和植被具有高度异质性, 生态系统的稳定性和抗干扰性差, 退化容易恢复难^[42]。土壤特性空间变异研究可为土壤质量的恢复与改善, 喀斯特脆弱生态系统植被的迅速恢复和生态重建提供理论依据。

4 结论

(1) 喀斯特常绿落叶阔叶混交林土壤各磷钾养分含量顺序依次为AP>TK>AK>TP, 变异程度均为中等变异。

(2) 各磷钾养分变异函数值的最佳拟合模型均为指数模型, 决定系数均很高；TP、AP空间自相关中等, TK、AK为弱空间自相关。

(3) 土壤TP表现为坡下(含洼地)含量高, 坡上含量较低；AK表现为中坡含量高于洼地含量；AP、TK呈斑块破碎化分布。

(4) 海拔、坡度和地面凹凸度是影响各土壤磷钾养分分布的主要因子。