南方典型丘陵山区不同高程耕地土壤养分变化特征分析

林建平 邓爱珍 赵小敏 江叶枫 韩 逸 谢 雨

(1.江西农业大学林学院， 南昌 330045； 2.江西农业大学江西省鄱阳湖流域农业资源与生态重点实验室， 南昌 330045；3.江西应用技术职业学院测绘地理信息学院， 赣州 341000)

0 引言

耕地土壤养分是农作物生长的必要条件，是表征耕地质量与肥力的重要指标，其含量直接影响作物的整个生长过程[1]。高程被认为是影响土壤养分流失与累积的关键因素之一[2]。高程变化导致水热资源在空间分布上存在差异，这些差异会影响土壤养分的矿化、降解、迁移和累积，进而导致不同高程下土壤养分含量存在较大空间异质性[3-4]。具体地，在低海拔地区，水热条件优越，土壤物质循环加速，养分的利用率高而沉积较少；在高海拔地区，气温下降，微生物分解速度减慢，造成有机质以及氮、钾等养分含量累积较多[5]。因此，深入了解土壤养分空间变异与高程的耦合关系，对于实现土壤养分的精准管理和耕地地力的提升具有重要意义。

近年来，国内外学者针对不同高程土壤属性空间变异特征做了大量研究，研究尺度主要集中在县域、小流域和一些特定山区，土地利用类型包括森林[6]、草地[7]和农田土壤[8]，研究内容包括不同高程土壤理化性质[9-10]、土壤侵蚀[11]和土壤微生物及酶活性变化[12]等。如TANG[13]在对峨眉山的研究发现，不同高程下植被和土壤性质均存在一定差异。丛静等[14]以神农架自然保护区为研究对象，研究发现，高程上升，土壤微生物和微生物生物量碳、氮增加。黄婷等[15]对云南省曲靖县植烟区不同高程下土壤属性分析表明，土壤pH值及土壤养分与高程显著相关。吴崇书等[16]分析了杭州市滨海平原、水网平原、河谷平原和丘陵山地4个地貌区耕地地力，结果显示，不同高程地貌区耕地土壤肥力差异明显。上述结果表明，高程对土壤属性空间变异会产生重要影响。

南方红壤丘陵区是我国重要的农业生产基地，该区土地总面积约占全国土地面积的1/5，人口占全国的40%，耕地占全国的30%，负担了我国近一半人口的粮食[17]，耕地分布具有明显的垂直地带性特征，揭示该区不同高程下耕地土壤养分的空间变异特征及其影响因素，对保障我国粮食安全、防止水土和养分流失、实现土壤养分的精细管理具有重要参考意义。为此，本文以我国南方典型红壤丘陵区——江西省赣州市为研究区，基于2012—2013年测土配方项目不同高程采集的3 098个耕层(0～20 cm)土壤数据，通过测定土壤有机质(Soil organic matter，SOM)、全氮(Total nitrogen，TN)、全磷(Total phosphorus，TP)、有效磷(Available phosphorus，AP)、全钾(Total potassium，TK)、有效钾(Available potassium，AK)养分含量，运用相关分析、单因素方差分析和协同克里格(Co-Kriging)插值法，探讨不同高程耕地土壤养分含量空间变异特征，以期为南方丘陵地区耕地表层土壤养分调控和耕地的可持续利用提供理论依据。

1 材料与方法

1.1 研究区概况

赣州市位于江西省南部，又称“赣南”，地理坐标介于24°29′～27°09′N，113°54′～116°38′E之间，土地总面积39 362.96 km2，辖18个县(市、区)。赣州市属中亚热带湿润季风型气候，气候温和，日照充足，雨量充沛，年平均气温18.9℃，年平均日照时数1 748.7 h，年降雨量1 639.8 mm。研究区土地利用类型主要为林地、耕地和园地(图1a)。地形地貌以山地、丘陵为主，总体格局为周高中低、南高于北，四周群山环绕，中部丘陵与平原交错分布。全市耕地空间分布存在明显的垂直地带性，主要分布在100～650 m高程范围内(图1b)。研究区成土母质主要为酸性结晶岩类风化物、石英岩类风化物和泥质岩类风化物(图1c)；土壤类型主要为水稻土、红壤、草甸土、石灰(岩)土、紫色土和山地黄壤6个土类(图1d)。作为典型的红壤区，该区水土流失较为严重，土壤侵蚀以水力侵蚀为主(图1e)。

图1 研究区土地利用类型、样点分布、成土母质、土壤类型和土壤侵蚀程度图Fig.1 Distribution maps of land use type, soil sampling site and digital elevation model (DEM) data, parent material, soil types and soil erosion

1.2 样品采集及分析

土壤样品采集按照NY/T 1634—2008《全国耕地地力调查与质量评价技术规程》要求，在考虑研究区地形地貌、土壤类型和成土母质等自然条件的前提下确定采样点位(图1b)。2012—2013年，通过实地采集0～20 cm耕层土样，每个采样点采集4个点，充分混匀后按照四分法留取1 kg样品，初始采集表层土壤样点3 098个，每个采样点均以GPS记录其经纬度坐标，并详细记录该点的高程、成土母质、土地利用类型和土壤类型等信息。土壤样品经室内自然风干、剔除杂物和研磨，统一过0.25 mm尼龙筛。分别采用高温重铬酸钾氧化容量法测定SOM含量；全自动凯氏定氮仪(Kjeltec8400型，Foss，瑞典)测定TN含量；钼锑抗比色法测定TP、AP含量；火焰光度法测定TK、AK含量。

1.3 数据处理

受采样及实验室分析误差的影响，土壤养分含量测量结果可能存在离群值，本文采用域值法(3倍标准差)[18]对离群值进行剔除，剔除后总样点为2 928个。采用SPSS 22软件，将高程与土壤养分含量进行Pearson相关性分析[19]；采用单因素方差分析(One-way ANOVA)进行差异性检验；采用GS+7.0进行半变异函数计算和模型拟合；运用ArcGIS 10.2进行协同克里格插值生成土壤养分含量空间分布图。高程数据来源于地理空间数据云(http:∥www.gscloud.cn)下载的赣州市30 m的数字高程模型(Digital elevation model，DEM)数据。土壤侵蚀图由水利部调查的江西省土壤侵蚀图处理提取。

2 结果与分析

2.1 土壤养分含量描述性统计

对照第二次全国土壤普查土壤养分分级标准，从表1可以看出，总体上，研究区AP含量处于2级(0.02～0.04 g/kg)丰富水平，养分含量均值0.02 g/kg，值域范围为0.003～0.08 g/kg；SOM、TN含量处于3级(20～30 g/kg、1～1.5 g/kg)中等水平，均值分别为28.91、1.47 g/kg，值域范围分别为12.10～61.10 g/kg、0.33～3.52 g/kg；TP、TK、AK含量处于4级(0.4～0.6 g/kg、10～15 g/kg、0.05～0.1 g/kg)缺乏水平，均值分别为0.51、11.16、0.07 g/kg，值域范围分别为0.22～1.04 g/kg、2.70～22.80 g/kg、0.02～0.26 g/kg。

从变异系数看，土壤养分呈中等变异[20]，养分变异系数由大到小依次为AP(66.47%)、AK(56.86%)、TK(39.61%)、TN(31.97%)、SOM(31.68%)、TP(31.37%)。总体上，TP相对稳定，变异系数较小，受成土母质的影响较大；而AP和AK变异系数最大，这可能与AP和AK为速效态养分有关。郭宏等[21]研究表明，速效态养分由于迁移性强和稳定性较差的特性，导致其变异程度较大。从分布类型看，TN养分分布的K-S检验系数为0.24，符合正态分布；SOM、TP、AP、TK、AK对数转换后K-S检验系数分别为0.19、0.17、0.09、0.11、0.21，也符合正态分布，满足地统计学研究要求。

表1 赣州市土壤养分含量(质量比)描述性统计特征Tab.1 Descriptive statistics of soil nutrient content in Ganzhou City

2.2 土壤养分含量与高程的相关性分析

将高程与土壤养分含量进行Pearson相关分析(表2)，结果表明，高程与SOM、TN、TK呈极显著正相关，相关系数分别为0.159(p<0.01)、0.156(p<0.01)、0.048(p<0.01)，可解释为海拔越高，SOM、TN、TK含量越高，这主要是海拔越高、温度越低，低温阻碍了土壤矿化作用，为养分的累积提供了条件。江叶枫等[22]对南方丘陵区土壤养分的研究也得出了类似结果。高程与TP、AP呈极显著负相关，相关系数分别为-0.090(p<0.01)和-0.222(p<0.01)，表明地势越低，TP、AP含量越高，这主要是由于磷含量随高程受成土母质和施肥水平共同影响，易在低洼处积累[23]。高程与AK相关性不强，AK作为易淋溶养分，受人为活动干扰较大。

表2 土壤养分与高程Pearson相关性分析Tab.2 Correlation analysis between elevation and soil nutrients

注：** 表示极显著相关(p<0.01)。

2.3 不同高程分组土壤养分含量的方差分析

高程与土壤养分含量相关性分析虽然能在总体上反映土壤养分随高程的变化关系，但局部信息描述不够。为详细地分析高程与耕地土壤养分含量的关系，利用DEM数据，按照自然断点法，将高程由低至高分成6个组别：0～150 m、150～200 m、200～250 m、250～300 m、300～350 m和大于350 m，分别对不同高程组别下的土壤养分均值进行单因素方差分析(表3)。从表3可以看出：

(1)SOM、TN养分含量均值随高程组别上升而增加，与表2相关性分析结果一致。SOM、TN养分含量在小于150 m组别最小，均值分别为26.88、1.36 g/kg，随着高程的增加，SOM、TN均值在大于350 m组别时达到最大，均值分别为32.26、1.62 g/kg，分别比最低高程组别(0～150 m)均值高5.38、0.26 g/kg。从不同组别养分含量变异系数看，SOM、TN变异系数分别为28.11%～33.56%、28.65%～33.24%，属中等变异。方差分析结果表明，当高程小于250 m时，SOM和TN含量随高程上升而呈上升趋势，但大于250 m时，差异性并不显著。

(2)TP和AP养分含量均值随高程上升呈先下降后上升趋势，与表2结果较为一致。其中，TP在小于300 m各组别养分含量均值随高程上升而下降，在250～300 m高程组别达最低值，为0.48 g/kg；当高程达到300～350 m和大于350 m组别时，养分均值又分别上升至0.50、0.52 g/kg；AP与TP类似，小于350 m各组别养分均值随高程上升而下降，在300～350 m高程组别达最低值，为17.80 mg/kg，当大于350 m时，养分均值又上升至18.90 mg/kg。从变异系数看，TP、AP变异系数分别为26.65%～32.27%、50.85%～72.77%，属中等变异。方差分析结果表明，当高程小于200 m时，TP和AP含量随高程上升而呈下降趋势，但大于200 m时，差异性并不显著。

(3)TK和AK养分含量均值随高程组别上升有所波动。TK养分总体随高程组别上升而上升，养分均值由高至低依次为大于350 m(12.65 g/kg)、300～350 m(11.90 g/kg)、0～150 m(11.42 g/kg)、200～250 m(11.20 g/kg)、150～200 m(10.79 g/kg)、250～300 m(10.64 g/kg)；AK养分均值随高程组别上升波动较大。从变异系数看，TK、AK变异系数分别为32.07%～42.91%、49.22%～57.91%，属中等变异。方差分析结果表明，TK在200～350 m高程范围内组间有显著差异，而在0～150 m、150～200 m和大于300 m组间无显著差异；AK在高程大于350 m与其他高程组间均有显著差异。

2.4 不同高程下土壤养分的空间分布特征分析

为直观反映土壤养分含量在不同高程的空间分布，首先进行半方差函数分析，确定土壤养分变异的最优理论空间插值模型；其次，根据最优理论空间插值模型，以高程为协变量进行协同克里格插值，得到土壤养分空间分布预测图；最后，将预测图与高程图进行叠加分析，得到不同高程土壤养分分布图。

表3 不同高程组土壤养分含量均值统计Tab.3 Statistics of mean value of soil nutrient contents in different elevation groups

注：同一列平均值后的小写字母相同，表示数据间未达到p=0.01水平上的显著性差异，下同。

2.4.1土壤养分含量空间分布拟合模型

半方差函数能够很好地反映引起土壤养分空间变异的随机性和结构性特征。结构性因素主要是指高程、成土母质、土壤类型等自然环境因素，随机性因素则是指土地利用方式等人为活动因素。通常用块金效应值(C0/(C0+C))说明样本间空间的相关性及由随机性因素引起的空间变异程度[24]。从图2可看出，TP、TK的最优模型分别为高斯模型和球状模型，其他养分最优模型均为指数模型。从模型的参数来看(表4)，研究区SOM、TN块金效应小于25%，说明其具有较强的空间相关性，空间变异主要由高程等结构因素引起；TK、AK、TP、AP块金效应在25%～75%之间，具有中等强度的空间自相关性，变异由随机和结构因素共同引起，结构因素大于随机因素。模型决定系数均在0.72以上，模型总体拟合较好。SOM、TN、TP、AP、TK、AK的变程分别为15.00、16.20、22.52、39.60、41.80、34.50 km，变程相对较大，这与赣州山地丘陵复杂的地形地貌空间结构相符合。

表4 土壤养分半方差函数模型参数Tab.4 Parameters of semi-variogram models for soil nutrients

2.4.2不同高程下土壤养分的空间分布

将土壤养分进行协同克里格插值，交叉验证结果表明(表5)：SOM、TN、TP、AP、TK、AK预测值与实测值的相关系数R分别为0.713、0.676、0.899、0.864、0.962、0.853，相关系数越大，表明预测精度越高[25]。均方离比差(MSDR)分别为1.395、1.814、1.072、0.950、0.839、0.963，MSDR系数越接近于1，说明精度越高[26]。总体上，协同克里格插值结果较为可信。土壤养分含量插值后，按照第二次全国土壤普查土壤养分分级标准进行重分类并与高程进行叠加，得到不同高程土壤养分含量分布图(图3)。从图3可以看出：①不同土壤养分各等级分布面积不同。SOM、TN、AP集中分布在2～3级(丰富～中等)，分别占研究区总面积的95.69%、95.46%、93.46%；TP养分含量主要集中在3～4级(中等～缺乏)，占研究区总面积的80.95%。TK、AK集中分布在 4～5级(缺乏～很缺乏)，分别占研究区总面积的80.70%、86.30%，与表1分析结果一致。②不同高程土壤养分含量有不同的空间分布格局。SOM、TN、TK养分含量高值区主要分布在高海拔地区，养分含量分布格局由外向内随高程降低呈周边高、中间低、逐渐下降的分布趋势。而TP、AP正好相反，养分含量总体随高程上升而下降，分布格局由外向内随高程降低呈周边低、中间高、逐渐上升的分布趋势。AK养分含量随高程分布规律不明显，这与表3分析结果比较一致。其中，SOM和TN的空间分布较为相似(图3a、3b)，高低分布趋势相当，这主要是因为两者有较强的空间耦合关系[27]。同时，大部分土壤养分含量的空间分布为突变而非渐变，能在一定程度上体现养分含量分布的空间异质性。

表5 土壤养分含量空间插值交叉验证结果Tab.5 Cross-validation results of spatial estimation of soil nutrients

图3 不同高程土壤养分含量分布图Fig.3 Soil nutrient distribution maps at different elevations

3 讨论

3.1 丘陵山区土壤养分的空间变异

受地形、母质、气候、生物和人类活动等因素的共同影响，丘陵山区土壤养分存在空间变异性。从研究区来看， SOM、TN均处3级中等水平，这与南方长期高温多雨的亚热带季风气候环境一致，这种气候条件有利于植物生长并促进植物凋落物转化为SOM，潮湿的空气也易产生嫌气条件，使微生物的活动降低，促进SOM、TN的累积；TP、TK、AK普遍处于4级缺乏水平，TP、TK剔除人为因素的影响，主要来源于土壤中含磷、钾矿物，南方砖红壤和红壤粘粒中含磷、钾矿物以高岭石为主[28]，磷、钾含量和释放量较低，加上南方地区高温高湿多雨易造成钾素的大量流失；AP处于2级含量丰富水平，一个重要原因可能是水田较多，占研究区耕地的80%。在淹水条件下，磷的释放率提高且不易淋失，与杨文等[29]对脱甲河小流域土壤碳氮磷研究结论一致；同时，研究区磷肥和有机肥(如猪粪等)施用量增加也是AP含量累积的重要原因。从半方差函数分析可以看出，土壤养分的块金值与基台值之比由小到大分别为SOM(12.130%)、TN(12.775%)、TK(25.189%)、AK(29.143%)、 TP(40.741%)、AP(42.744%)，其中，土壤SOM、TN块金效应均小于25%，说明两者含量具有强烈的空间相关性，其空间变异主要受结构性因子影响；而TK、AK、TP、AP块金效应在25%～75%之间，说明其具有中等空间相关性，其空间变异受结构性因子和随机性因子共同影响，但受高程等结构因素影响更大。从变异系数看，土壤养分变异系数由大到小依次为AP(66.47%)、AK(56.86%)、TK(39.61%)、TN(31.97%)、SOM(31.68%)、TP(31.37%)，其中以AP和AK变异系数最大，这可能与AP和AK均为速效态养分迁移性强和稳定性较差的特性有关。

3.2 高程对土壤养分空间变异的影响

不同高程土壤养分空间变异规律和程度存在差异，影响因素也不同，可能的原因如下：

(1)高程对SOM、TN含量的影响。南方丘陵山区耕层土壤SOM、TN含量随海拔升高而增加，养分含量在大于350 m高程组平均值达到最大值。究其原因，主要是SOM、TN含量主要受植物残体、动物废弃物归还量和矿化速率的影响：①不同高程下秸秆还田量不同。秸秆作为植物残体是农田有机物来源的重要组成部分，还田后可积累土壤有机质。但不同高程耕地秸秆还田量不一，海拔较低处，耕地集中连片程度高，耕作、交通等条件较海拔较高处要好，因而，大量秸秆收集用于工业造纸或喂养牲畜，而海拔较高处田块分散且细碎，加上交通不便，根茬残留量和秸秆就地还田量要比海拔较低处要多，从而使更多秸秆转变为土壤有机质。②不同高程植物残体矿化的速率不同。植物残体矿化的速率受温度、湿度和微生物等的影响较大，而高程通过影响光、温、水、热资源的分配，导致局部小气候发生变化，从而影响土壤中SOM、TN的分解和积累。一般而言，海拔越高，土壤温度越低，土壤湿度越大，土壤微生物活动会受到抑制，分解速度减慢，从而矿化作用减弱，因而，海拔较高处有利于有机质的积累。李丹维等[30]对太白山以及陈桂秋等[31]对湖南省、广西壮族自治区不同温度长期定位实验也得出相似结论。③少耕、免耕技术的推广，也是农田土壤SOM、TN增加的重要原因。吴建富等[32]对江西省南部稻区的研究表明，免耕1年后，较翻耕处理土壤有机质含量增加1.22 g/kg。研究区海拔较高处，很多耕地由于气温及人类活动的影响，只种植一季水稻，减少了耕作次数，因而有利于有机质的累积。

(2)高程对TP、AP含量的影响。研究区TP、AP含量随高程上升而下降， TP、AP养分均值均在小于150 m达到最大值。出现上述规律的可能原因是：①不同高程成土母质不同，高程与成土母质存在一定的耦合关系。从图1c可以看出，赣南土壤成土母质主要为酸性结晶岩类、泥质岩类风化物，紫色岩类、红砂岩类、石英类风化物以及第四纪红色粘土。研究表明，不同成土母质含磷量变化较大，并随岩土酸性程度的增加而下降[33]。将高程(图1b)与成土母质图(图1c)进行叠加分析得出，不同高程成土母质分布存在明显差异，研究区海拔较高的地区，广泛分布着花岗岩及石灰性岩类母质，由于含云母、长石矿物较多，且土壤质地较粗，因此，磷的含量较低。研究区海拔较低的中部地区，由河湖沉积物发育形成的土壤磷含量显著高于第四纪红土和泥质岩类风化残积—坡积物发育形成的土壤，赣州市土壤志表明该类成土母质含磷量较多，由此形成的土壤熟化程度也较高，且质地为壤土，对磷的富集起着一定作用。总体上，不同高程成土母质发育的土壤磷含量均表现出明显的从上往下依次递减规律，并具有较强的表聚性。②不同高程下土地利用方式不同，对土壤磷的分布产生重要影响[34]。研究区耕地的土地利用方式主要为水旱轮作、一季旱地、两季旱地、一季水田和两季水田，耕地的土地利用方式与高程密切相关，随着高程的上升，研究区土地利用方式逐渐由两季水田、两季旱地或水旱轮作过渡到一季旱地或一季水田。相比两季种植，一季种植施肥量少，田间管理措施、投入管理水平没有两季精细，且一季种植撂荒时易生长杂草，通过杂草生长也带走部分养分，因而，总体上，研究区磷的含量表现出随高程增加而逐渐降低的趋势。

(3)高程对TK、AK含量的影响。研究区TK含量随高程上升而增加，可能原因为：不同高程发育不同土壤类型、成土母质和土壤类型对钾含量的影响较大，从图1d可以看出，研究区土壤类型主要有水稻土、红壤、草甸土、石灰(岩)土、紫色土和山地黄壤，各种土壤的供钾能力大小不一[35]。从高程与土壤类型叠加分析结果可以看出，海拔较低处，分布着一定数量的花岗岩红壤和第四纪红土红壤，花岗岩土壤砂性强、保水保肥能力较弱、钾离子相对容易淋失，第四纪红土红壤虽然粘性强、土壤吸附性能好，但由于土壤胶体几乎全部被H+或AL3+饱和，从而使盐基离子很难在土壤胶体中保持下来，土壤对钾素养分的保持能力亦相对较弱。在海拔较高处，分布着一定数量的微碱性紫色土，该土壤阳离子交换量较大、交换性盐基离子含量较高，对钾素的保持和供应能力亦较强，故土壤中出现了钾的累积。AK含量随高程变异较大，规律性不明显，主要原因是AK为易淋溶养分，易受到施肥、水流和土壤侵蚀的影响。南方丘陵区是典型的红壤区，降雨充沛且汛期集中，是我国水土流失最严重的区域之一[36]，土壤侵蚀带走大量泥土的同时，也带走了土壤养分。将DEM高程图与土壤侵蚀图(图1e)进行空间叠加分析可以看出，随着高程的增加，高程与土壤侵蚀大致呈倒“U”形关系，即随着高程的上升，土壤侵蚀加剧，上升到一定高度后又下降，由于南方丘陵山区复杂的地形地貌，以及速效钾的易淋溶特征，使AK含量随高程变化规律异常复杂，与赵明松等[37]对安徽省速效钾养分流失评估、邓欧平等[38]对四川省紫色土速效钾的研究结论基本一致。

3.3 高程对土壤养分含量空间变异影响的启示

高程是土壤养分含量空间变异的重要影响因素之一，高程通过改变生态环境的光、热、水、气等条件，影响土壤的物理、化学和微生物性质，从而使不同高程土壤养分含量存在明显的差异并呈现一定的规律变化。研究高程与耕地土壤养分的相关性及其垂直分布特性，分析海拔高度对土壤养分的影响，可为不同地形地貌耕地地力评价、农用地分等定级、土地利用潜力评价、土壤改良措施制定、水土流失评价等提供依据。从不同高程土壤养分空间变异特征分析，可以得出如下启示和建议：①根据山地丘陵区不同海拔高程耕地土壤养分分异结果，可对土壤养分实行分高程梯度精准管理，不断提升耕地质量，实现农作物高产稳产。②土壤侵蚀在研究区大于350 m高程带以中度及以上侵蚀为主，易淋溶养分流失最严重，是土壤侵蚀治理的重点区域，建议对该区域耕地土地利用方式进行转型，如将坡耕地转变为果园、林地，或在耕作的基础上进一步采取坡改梯等工程措施和生物等水保措施，防止水土和养分流失。③为进一步细化高程对土壤养分的影响，掌握不同高程耕地土壤养分的运移及累积规律，建议在以后的研究中，在不同高程建立土壤养分含量长期定位调查，并与土地利用方式、土壤类型等相结合， 针对土壤养分变化与作物生长状况，制定适宜的耕地利用方式，维持山地丘陵区耕地土壤养分和生态平衡。

4 结论

(1)统计分析表明，南方丘陵山区AP养分含量丰富，均值为0.02 g/kg；SOM、TN 含量处中等水平，均值分别为28.91、1.47 g/kg；TP、TK、AK含量缺乏，均值分别为0.51、11.16、0.07 g/kg。

(2)相关分析结果表明，高程与SOM、TN、TK养分含量呈极显著正相关，相关系数分别为0.159(p<0.01)、0.156(p<0.01)、0.048(p<0.01)；与TP、AP呈极显著负相关，相关系数分别为-0.090(p<0.01)和-0.222(p<0.01)；与AK显著性不明显。

(3)方差分析结果表明，不同高程分组下土壤养分含量均值存在显著差异。SOM、TN、TK养分含量随高程组上升而增加；TP、AP养分含量均值随高程上升而降低；AK随高程分组无明显规律性。

(4)空间插值与高程叠加分析表明，SOM、TN和TK养分含量高值区分布在研究区周边高海拔区形成的连片区域；TP、AP养分含量高值主要分布在研究区中部低海拔区域；AK养分含量分布规律不明显。