掌握土壤养分的空间变异性,是实现土壤养分分区管理、作物优化施肥,指导基本农田建设,进而提高作物产量和品质的重要基础。在土壤养分空间变异研究方面,既有对区域大尺度的研究[1-2],如通过土壤有效磷时空变异研究进行烟草施肥分区[3];也有小尺度的研究[4-6],如探讨果园土壤肥力空间分布与作物产量的关系[7]。这种结合地统计学和GIS对土壤养分的空间变异特征进行的研究已成为土壤科学研究的热点。关于土壤养分空间变异的养分指标以氮、磷、有机质等居多[8-9],钾被誉为植物的品质元素与抗逆元素,且土壤速效钾(AK)作为植物可以直接吸收利用的元素形态,探究其在土壤中的空间变异对区域种植规划、提高作物品质等具有重要意义。目前专门对AK空间分布的研究主要有川西南山区[10]、新疆艾比湖湿地区[11]、鄱阳湖平原区[12]及秦巴中部山区[13]等区域,各研究表明耕作层速效钾空间变异多受结构性因素与随机性因素共同影响,且其空间插值斑块大小、形状及分布等具有显著的区域差异。然而,上述研究仅选取0~20 cm耕作层土壤进行水平空间异质性的研究,对垂直空间异质性研究不足。
重庆潼南国家农业科技园区地处以紫色土为主的典型渝西浅丘带坝区,其基本农田建设规模大,农业综合开发强度高,是以市场为导向、以科技为支撑、以企业为主导的现代农业建设新模式。然而,基本农田建设中表土剥离与回填、深松深翻等措施将改变原有土壤养分的空间分布格局。本文选取该园内小岭村2017年集中建设区,结合普通克里格法和聚类分析法对其不同土层的AK空间分布及影响因素进行研究,以期为园区种植规划、基本农田建设和农业综合开发等提供参考。
1 材料与方法 1.1 研究区域概况潼南国家农业科技园区位于渝西丘陵区,核心区约10 km2,涉及罐坝、鱼溅、郭坡、小岭四个行政村。小岭村2017年集中建设区内兼有平坝(北部及中部)和丘陵(南部),面积97.55 hm2,海拔250~272 m。土壤以侏罗系遂宁组泥岩发育而成的紫色土和以冲洪积砾石及砂土母质或遂宁组泥岩母质发育的水稻土为主。潼南区属亚热带湿润季风气候,年均温17.9℃,年降水1 006 mm。区内以水田、旱地为主,旱地多种玉米、花生。
1.2 土壤样品采集与测定采用S形取样法于2017年春通过GPS定位,对41个采样点分0~20、20~40和40~60 cm采集3层土样,共获得土样104个(图 1)。由于部分采样点土层较薄,三层土样分别为41、34和29个。野外采回的土样经风干、过筛、混匀后,采用乙酸铵浸提—火焰光度法[14]进行AK含量的测定。
经正态检验(K-S检验法)剔除异常值(3δ法、Grubbs检验法)后,比较不同土壤类型、地形、土地利用类型及土层间土壤AK分布(单因素ANOVA分析、皮尔森(Pearson)相关性分析)及Q型聚类分析均在SPSS 25.0中完成。常规统计分析均在Excel 2016中完成。
空间自相关分析中的全局莫兰指数(Global Moran’s I)是在定量基础上对研究目标空间自相关程度进行衡量的一种方法,同时也可判别区域化变量在指定区域内是否存在集聚区与孤立区[15]。双变量空间自相关分析是利用定量指标来描述统计关系间强弱的一种方法,其目的是分析两变量之间有无空间相关性,进而探究它们之间是否可以进行相互预测或者解释[16]。本文选取基于距离的空间权重,以最小距离阈值的空间自相关分析、地形因子与AK双变量空间自相关分析,分析过程均在GeoDa 1.10完成。
半变异函数(Semi variation function)计算在GS+ 9.0完成,对应最佳函数模型以决定系数(R2)越大、残差平方和(RSS)越小为原则确定[17]。结合GS+ 9.0半变异函数最优模型参数,采用ArcGIS 10.2普通克里格方法进行空间插值分析与制图。采用ArcGIS 10.2区域统计分析得到研究区各土地利用类型斑块在不同土层深度的速效钾均值,供后续进行聚类分析和统计分析。
2 结果 2.1 不同土层AK描述性统计特征如表 1所示,研究区0~20、20~40和40~60 cm土层AK平均值分别为111.6、96.1和90.2 mg·kg-1。0~20 cm土层AK含量显著高于其他两层(P < 0.05),AK随土层深度的增加而降低。就变异系数而言,各土层AK均为中等变异。
与全国第二次土壤普查土壤养分分级标准[18]进行对比(图 2),0~20 cm土层AK含量大部分集中于中等水平(100~150 mg·kg-1),20~40 cm和40~60 cm土层AK含量大部分处于缺乏水平(50~100 mg·kg-1)。
各土层AK的全局莫兰指数依次为-0.184、-0.213、0.318,表明研究区AK具有显著的空间自相关性(P < 0.05)。图 3表明,0~20和20~40 cm土层速效钾主要分布在第2、4象限且对应全局莫兰指数均为负值,属于中等偏低的“低-高”和“高-低”孤立类型,表明该两层AK较高或较低的地区交错随机分布;40~60 cm主要分布于第1、3象限且全局莫兰指数为正,属于“高-高”和“低-低”集聚类型,表明该层AK较高的地区之间相邻,含量较低地区之间相邻。
如图 4所示,0~20 cm土层AK的空间自相关随距离增大起伏较明显,全局莫兰指数在500~650 m范围内为正值,表明在该范围内AK在结构上呈正相关,存在空间集聚。20~40 cm土层AK在小于850 m范围内全局莫兰指数均为正值且随距离增大出现明显波动,具有一定空间集聚性;大于850 m范围内其指数均为负值且随距离增大逐渐减小,空间孤立性显现且逐渐明显。40~60 cm土层AK的全局莫兰指数变化表现为随距离增加而减小的趋势,小于750 m范围内莫兰指数均为正值且随距离增大逐渐减小,表明空间集聚性极弱;大于750 m范围内其指数转为负值且随距离增大而逐渐减小,表明空间孤立性显现并趋于平稳。该研究结果可对土壤野外采样提供指导[19],表明较合理采样点间距应设置在500~750 m之间。
经检验数据符合正态分布,符合半方差函数分析要求,表 2为最优模型的块金值、基台值等相关参数。块金值(C0)是由实验误差和小于实际取样尺度引起的变异,其大小直接限制了空间内插的精度,基台值(C0+C)表示系统内的总变异量,是区域化变量的最大变异。块基比(C0/(C0+C))用于判定区域化变量空间异质性的程度[20]。0~20和20~40 cm土层最优模型均为高斯模型,而40~60 cm土层为球型模型。其中,0~20和40~60 cm土层AK的块基比均高于75%,表明AK的空间变异以随机因素变异为主,且0~20 cm土层随机性更强。而20~40 cm土层AK块基比处于25%~75%之间,表明其空间变异是结构因素(如土壤母质、地形因子与水文等)和随机因素(如耕作措施、施肥管理和土壤改良等人类活动[13])共同造成的。变程是指变异函数在有限步长上达到基台值时对应的步长,可反映出变量空间自相关范围大小。从表 2可知,各土层AK的变程介于74~452 m之间,对应AK空间连续性由大到小依次为20~40、40~60、0~20 cm。0~20 cm变程仅为74 m,表明随机因素占其空间变异主体。
各土层AK水平空间分布主要呈现出条带状和斑块状格局(图 5)。0~20 cm土层属于中等水平(100~150 mg·kg-1)的面积占绝大部分,高值集中于研究区中部,低值分布较少且位置随机。20~40和40~60 cm的AK分布趋势大致相同且缺乏水平(< 100 mg·kg-1)面积均超过一半。20~40 cm土层AK高值集中于区域东南部,在区域西部边缘仅有少量分布;低值集中于东北部。40~60 cm土层高值区仅表现为个别岛屿分布于中部且面积较小,低值集中于东北部且范围更广。0~20 cm土层AK高低分布更为随机,这与空间自相关分析结果(图 3)一致。
整体上0~20 cm土层AK分布较其余两层差异较大,且后两层高低分布相对均匀(图 5)。Q型聚类分析(图 6)表明,由欧氏距离系数[21]可将研究区各斑块土壤速效钾剖面明显分为3类,即均匀型、低钾递减型和高钾递减型(图 7),三者占斑块总数的比例依次为58.7%、35.7%和5.6%,占总面积的比例依次为49.4%、47.9%和2.7%。
如图 7所示,均匀型剖面的主要特征是各土层AK含量差异较小,无明显递增、递减趋势。低钾递减型剖面的主要特征是整体偏低(70~110 mg·kg-1),且中下层较表层明显降低,平均减幅达35.5%。高钾递减型剖面的主要特征是AK含量整体较高(120~150 mg·kg-1),并随土层深度增加AK含量递减,平均减幅达16.1%。
3 讨论 3.1 土壤类型对AK的影响研究区主要土壤类型包括由冲洪积砾石及砂土(Q2)发育的黄潮田(U34),由遂宁组泥岩(J3s)发育的红棕紫泥土(N32)和紫泥田(U42)。研究区北部主要为紫泥田,中部主要为黄潮田,南部主要为红棕紫泥土。紫泥田三层AK含量较其他两种土壤类型均较低(表 3),该结果充分解释了AK中下层水平空间分布主体南高北低的原因,而表层因受施肥等随机因素影响,南高北低趋势不明显。
一般而言,黄潮田对应母质含钾量较低[22],而本文研究中黄潮田AK在三层土层中均为最高(表 3),高钾递减型剖面类型均对应黄潮田。这可能与黄潮田对应区域为平坝,土地利用类型变化较频繁有关。红棕紫泥土表层AK显著高于中下层,这可能受作物耕作补充表层钾素有关。均匀型剖面在三种土壤类型中均有分布,低钾递减型多对应紫泥田。
3.2 地形地貌对AK的影响研究区丘陵和平坝0~20 cm土层AK分别为117.0±38.3、108.8±31.9 mg·kg-1,该结果与杨柳等[23]对川中丘顶、丘坡及冲沟土壤养分研究结果相似,由于林地主要集中于丘陵,故其AK较高。有研究表明,在中小尺度区域内地形因子如海拔、坡度、平面与剖面曲率等对土壤养分的空间变异影响显著[19]。皮尔森(Pearson)相关性分析结果表明,研究区土壤AK与海拔、坡度、剖面曲率、平面曲率及坡向均无显著性相关。为进一步明确地形地貌因子在空间上与AK的相关性,本文采用双变量空间自相关分析法进一步分析(表 4),其中坡向由于是环形变量,故进行余弦转化为-1~1新变量。结果表明,五种地形因子与各土层AK双变量空间自相关系数绝对值范围在0.002~0.336之间,海拔和坡度与各土层AK含量存在极显著的双变量空间自相关性(P < 0.01),即二者在局部空间上仍与AK含量的变化趋势一致,坡向与底层土壤AK含量也存在上述现象(P < 0.05)。综合皮尔森相关性及双变量空间自相关分析,由于本研究区为浅丘带坝地貌,海拔变化较小(250~272 m),从而使得总体上各地形因子对AK空间变异影响均较弱,未表现出显著相关性,这解释了0~20和40~60 cm土层AK空间变异主要是由随机因素所导致;但在空间分布上,海拔、坡度和坡向地形因子与AK仍具有一定的显著双变量空间自相关性,该结果也解释了20~40 cm土层AK空间变异是由结构性变异和随机性因素共同导致。
本研究区内除建设用地和水域外主要为草地、旱地、水田和林地4种土地利用类型,总体上各土地利用类型土壤AK均以0~20 cm土层最高,各土层中均以林地最高(图 8),这与林地土壤养分是以充足草木枯萎微生物分解得以补充有关[24]。林地主要集中于研究区中部及北部少数区域(图 1),对应0~20 cm土层AK水平空间分布中高值斑块区域。0~20 cm土层旱地AK显著高于其他两层,这主要是由于施肥与秸秆还田[25]等随机因素造成的;水田由于长期处于淹水状态AK淋溶损失较大,故其含量低于林地,这与张倩等[26]对重庆铜梁的研究结果相似。同土层旱地与水田对应AK含量较为相近,且二者面积占研究区主体部分,使得各土层空间分布中出现较大的面状斑块。
土壤AK三种剖面类型中,均匀型剖面斑块总数最多,其中各土地利用类型斑块占该剖面类型斑块总数的比例由大到小依次为:旱地55.4%、水田16.4%、草地16.0%、林地12.2%。表明大部分旱地由于施肥补充土壤钾素与作物吸收消耗钾素相互作用导致表层土壤与中下层AK含量相近。低钾递减型速效钾剖面主要包括旱地、草地及小部分水田,表明该部分由母质风化发育而来的土壤含钾较低,故中下层AK含量较低,而表层AK由于施肥等随机因素影响而明显高于中下层。高钾递减型速效钾剖面仅含6个水田斑块与1个草地斑块(水田撂荒地),这可能由于水田施钾肥后钾素溶于水中随水下渗使得所在深层土壤含钾量也较高,同时也解释了20~40与40~60 cm空间分布中出现的少数高值岛屿。
4 结论渝西丘陵区潼南国家农业科技园核心区在0~60 cm土层中AK平均含量均处于中等-缺乏水平且均为中等变异,AK含量具有显著的空间自相关性,水平和垂直空间变异明显。不同土层表现为不同的水平空间分布特征,中上层具有明显的空间孤立性,尤以表层更为突出,而下层具有明显的空间集聚性;上述特征受土壤类型、地形地貌和土地利用类型的综合影响,其中,土地利用类型对中层影响显著,土壤类型对底层影响显著。垂直空间分布特征表现为均匀型、高钾递减型和低钾递减型三种剖面类型,以均匀型为主。土壤类型中红棕紫泥土、土地利用类型中旱地和林地的各剖面AK含量均为显著递减。在基本农田建设中,应根据AK含量空间分布特征采取相应的措施,研究区大部分区域需重视表土剥离和回填,东部和南部区域可适度进行深松、深翻改土,并在整治后尤其重视低钾递减型剖面区域的钾素施用。
[1] |
张智, 任意, 鲁剑巍, 等. 长江中游农田土壤微量养分空间分布特征. 土壤学报, 2016, 53(6): 1489-1496. Zhang Z, Ren Y, Lu J W, et al. Spatial distribution of micronutrients in farmland soils in the mid-reaches of the Yangtze River (In Chinese). Acta Pedologica Sinica, 2016, 53(6): 1489-1496. (0) |
[2] |
Paterson S, Minasny B, McBratney A. Spatial variability of Australian soil texture:A multiscale analysis . Geoderma, 2018, 309: 60-74. DOI:10.1016/j.geoderma.2017.09.005
(0) |
[3] |
李强, 戴美玲, 向德明, 等. 湘西喀斯特区植烟土壤有效磷时空变异及风险评估. 土壤, 2018, 50(1): 181-189. Li Q, Dai M L, Xiang D M, et al. Temporal and spatial variability risk assessment of soil phosphorous in Karst tobacco-planting region of Xiangxi, Hunan (In Chinese). Soils, 2018, 50(1): 181-189. (0) |
[4] |
Wang S, Zhuang Q L, Jia S H, et al. Spatial variations of soil organic carbon stocks in a coastal hilly area of China . Geoderma, 2018, 314: 8-19. DOI:10.1016/j.geoderma.2017.10.052
(0) |
[5] |
Li T, Zheng W W, Zhang S R, et al. Spatial variations in soil phosphorus along a gradient of central city-suburb-exurban satellite . Catena, 2018, 170: 150-158. DOI:10.1016/j.catena.2018.06.011
(0) |
[6] |
刘文全, 卢芳, 徐兴永, 等. 废弃盐田复垦利用后土壤盐分与有机质含量空间变异特征. 生态学报, 2018, 38(4): 1311-1319. Liu W Q, Lu F, Xu X Y, et al. Spatial variation of soil salinity and organic matter under reclamation in an abandoned salt pan (In Chinese). Acta Ecologica Sinica, 2018, 38(4): 1311-1319. (0) |
[7] |
Dos Santos R O, Franco L B, Silva S A, et al. Spatial variability of soil fertility and its relation with cocoa yield . Revista Brasileira de Engenharia Agrícola e Ambiental, 2017, 21(2): 88-93. DOI:10.1590/1807-1929/agriambi.v21n2p88-93
(0) |
[8] |
江叶枫, 孙凯, 郭熙, 等. 南方红壤区不同侵蚀程度下土壤有机质空间变异的影响因素研究. 自然资源学报, 2018, 33(1): 149-160. Jiang Y F, Sun K, Guo X, et al. Spatial variability of organic matter and its influencing factors in red soil with different erosion degrees in South China (In Chinese). Journal of Natural Resources, 2018, 33(1): 149-160. (0) |
[9] |
张泽宇, 李占斌, 李鹏, 等. 汉江小流域土壤氮素空间分布特征及影响因素. 水土保持研究, 2017, 24(4): 46-52, 58. Zhang Z Y, Li Z B, Li P, et al. Spatial variability of soil nitrogen and the impact factors of land-use types in Yujiehe Watershed, Han River (In Chinese). Research of Soil and Water Conservation, 2017, 24(4): 46-52, 58. (0) |
[10] |
王永豪, 王昌全, 李启权, 等. 川西南山区耕地土壤速效钾空间特征及影响因素分析. 西南农业学报, 2016, 29(3): 651-657. Wang Y H, Wang C Q, Li Q Q, et al. Geostatistics and GIS analysis of soil available potassium in mountainous southwestern Sichuan (In Chinese). Southwest China Journal of Agricultural Sciences, 2016, 29(3): 651-657. (0) |
[11] |
艾尤尔·亥热提, 王勇辉, 海米提·依米提. 艾比湖湿地土壤速效钾空间变异性分析. 土壤通报, 2015, 46(2): 375-381. Ghayur·G, Wang Y H, Hamit·Y. Spatial variability of soil available K concentrations in Ebinur lake wetland (In Chinese). Chinese Journal of Soil Science, 2015, 46(2): 375-381. (0) |
[12] |
江叶枫, 郭熙. 基于协同克里格的耕层土壤速效钾空间异质性研究. 浙江农业学报, 2019, 31(1): 139-148. Jiang Y F, Guo X. Spatial heterogeneity of available potassium in topsoil based on co-Kriging (In Chinese). Acta Agriculturae Zhejiangensis, 2019, 31(1): 139-148. DOI:10.3969/j.issn.1004-1524.2019.01.18 (0) |
[13] |
陈洋, 齐雁冰, 王茵茵, 等. 秦巴中部山区耕地土壤速效钾空间变异及其影响因素. 环境科学研究, 2017, 30(2): 257-266. Chen Y, Qi Y B, Wang Y Y, et al. Spatial variability and factors affecting soil available potassium in the central Qinling-Daba mountain area (In Chinese). Research of Environmental Sciences, 2017, 30(2): 257-266. (0) |
[14] |
杨剑虹, 王成林, 代亨林.土壤农化分析与环境监测.北京: 中国大地出版社, 2008. Yang J H, Wang C L, Dai H L. Analytical methods for soil, agro-chemistry and environmental monitoring. Beijing: China Land Press, 2008. (0) |
[15] |
杨奇勇, 蒋忠诚, 袁道先, 等. 广西典型岩溶区土壤水分含量空间自相关分析. 中国岩溶, 2015, 34(3): 260-265. Yang Q Y, Jiang Z C, Yuan D X, et al. Spatial autocorrelation analysis of soil water content in a karst region of Guangxi Province (In Chinese). Carsologica Sinica, 2015, 34(3): 260-265. (0) |
[16] |
姚小薇, 曾杰, 李旺君. 武汉城市圈城镇化与土地生态系统服务价值空间相关特征. 农业工程学报, 2015, 31(9): 249-256. Yao X W, Zeng J, Li W J. Spatial correlation characteristics of urbanization and land ecosystem service value in Wuhan urban agglomeration (In Chinese). Transactions of the Chinese Society of Agricultural Engineering, 2015, 31(9): 249-256. (0) |
[17] |
吴敏, 刘淑娟, 叶莹莹, 等. 喀斯特地区坡耕地与退耕地土壤有机碳空间异质性及其影响因素. 生态学报, 2016, 36(6): 1619-1627. Wu M, Liu S J, Ye Y Y, et al. Spatial variability of surface soil organic carbon and its influencing factors in cultivated slopes and abandoned lands in a Karst peak-cluster depression area (In Chinese). Acta Ecologica Sinica, 2016, 36(6): 1619-1627. (0) |
[18] |
四川省农牧厅, 四川省土壤普查办公室.四川省第二次土壤普查数据资料汇编.成都: 成都地图出版社, 1992. Department of Agriculture and Animal Husbandry of Sichuan Province, Office of Soil Survey of Sichuan Province. Compilation of data of the second soil census in Sichuan Province. Chengdu: Chengdu Cartographic Publishing House, 1992. (0) |
[19] |
宋轩, 李立东, 寇长林, 等. 黄水河小流域土壤养分分布及其与地形的关系. 应用生态学报, 2011, 22(12): 3163-3168. Song X, Li L D, Kou Z L, et al. Soil nutrient distribution and its relations with topography in Huangshui River drainage basin (In Chinese). Chinese Journal of Applied Ecology, 2011, 22(12): 3163-3168. (0) |
[20] |
张敏, 贺鹏飞, 陈伟强. 基于GIS和地统计学的土壤养分空间变异分析. 东北农业大学学报, 2010, 41(3): 53-58. Zhang M, He P F, Chen W Q. Spatio-temporal variability analysis of soil nutrients based on GIS and Geostatistics (In Chinese). Journal of Northeast Agricultural University, 2010, 41(3): 53-58. (0) |
[21] |
李二焕, 胡海波, 鲁小珍, 等. 苏北滨海盐土区土壤盐分剖面特征及其理化特性. 水土保持研究, 2016, 23(4): 116-119, 127. Li E H, Hu H B, Lu X Z, et al. Characteristics of soil salinity profiles and physical and chemical properties of coastal saline soil in the estuary of North Jiangsu (In Chinese). Research of Soil and Water Conservation, 2016, 23(4): 116-119, 127. (0) |
[22] |
四川省农牧厅, 四川省土壤普查办公室.四川土种志.成都: 四川科学技术出版社, 1994. Department of Agriculture and Animal Husbandry of Sichuan Province, Office of Soil Survey of Sichuan Province. Soil taxonomy of Sichuan. Chengdu: Sichuan Science and Technology Press, 1994. (0) |
[23] |
杨柳, 胡玉福, 许宗林, 等. 川中丘陵区土壤养分时间变异特征研究——以中江县玉兴镇为例. 土壤通报, 2009, 40(5): 1057-1062. Yang L, Hu Y F, Xu Z L, et al. Temporal variation characteristics of soil nutrient in hilly region of mid-Sichuan Basin:A case study in Yuxing Town of Zhongjiang County (In Chinese). Chinese Journal of Soil Science, 2009, 40(5): 1057-1062. (0) |
[24] |
刘钊, 魏天兴, 朱清科, 等. 黄土丘陵沟壑区典型林地土壤微生物、酶活性和养分特征. 土壤, 2016, 48(4): 705-713. Liu Z, Wei T X, Zhu Q K, et al. Microbes, enzyme activities and nutrient characteristics of rhizosphere and non-rhizosphere soils in forests of loess hilly region (In Chinese). Soils, 2016, 48(4): 705-713. (0) |
[25] |
隋鹏祥, 张文可, 梅楠, 等. 不同秸秆还田方式对春玉米产量、水分利用和根系生长的影响. 水土保持学报, 2018, 32(4): 255-261. Sui P X, Zhang W K, Mei N, et al. Effects of different straw returning methods on spring maize yield, water use and root growth (In Chinese). Journal of Soil and Water Conservation, 2018, 32(4): 255-261. (0) |
[26] |
张倩, 胡胜勇, 高明, 等. 紫色丘陵区土壤养分空间变异特征研究——以重庆市铜梁县为例. 中国农学通报, 2011, 27(20): 216-223. Zhang Q, Hu S Y, Gao M, et al. Study on spatial variability of soil nutrients in purple soil hilly region——Take Tongliang in Chongqing as an example (In Chinese). Chinese Agricultural Science Bulletin, 2011, 27(20): 216-223. (0) |