2. 环境保护部南京环境科学研究所,南京 210042;
3. 江苏省耕地质量与农业环境保护站,南京 210036
钾与氮、磷是植物生长必需的三大营养元素。速效钾作为植物能够直接吸收的钾素形态,通常被视作衡量土壤供肥能力的重要指标之一。土壤速效钾主要来源于土壤矿物的释放、钾肥投入以及植物秸秆中钾的返还,植物吸收和钾素径流是土壤速效钾减少的主要途径[1-2]。长期以来,学术界通过实验或统计分析研究了成土母质、气候、土壤质地、土壤pH、施肥、种植制度等影响因子与土壤速效钾转化循环的关系[3-6],为提升土壤钾素的有效性,增加土壤速效钾含量提供了科学依据。但这些研究大多是针对土壤钾素有效性或土壤速效钾含量变化的单一影响因子研究,对不同影响因子之间交互作用,特别是自然因素与施肥等人为因素之间的交互效应探讨较少。
研究表明,第二次土壤普查后江苏省农田土壤速效钾含量一直处于下降趋势,直到20世纪90年代中期推广“补钾工程”,这一趋势才得到遏制[7-9]。为了对土壤钾素资源进行科学管理,不少学者围绕江苏农田土壤速效钾的时空变异及其影响因子开展了大量研究。徐茂[10]研究发现,江苏环太湖地区21世纪初的土壤速效钾含量与第二次土壤普查相比总体持平,自然因素和人为因素对该地区土壤速效钾的影响力相当;杨艳丽等[11]和沈德福等[12]分别对江苏北部和沿江地区的土壤养分空间变异进行研究,认为土壤速效钾受成土母质影响显著,同时受生物、气候等条件影响深刻。孙瑞娟等[13]研究认为,钾肥投入不足是导致苏南地区土壤速效钾含量下降的重要原因。这些研究大多集中于较小的行政区域或特定流域,研究方式一般为近年来的肥力数据与第二次土壤普查的结果进行对比,缺乏对近年来全省农田土壤速效钾含量变化的长期定位监测研究,对土壤速效钾的各影响因子交互效应也较少涉及。
本文基于江苏基本农田质量长期监测点土壤速效钾的实测数据,分析了2008—2015年全省农田土壤速效钾的时空变异特征,并运用多因素方差探明了成土母质、气候、种植制度、化肥对土壤速效钾含量变化的影响及不同影响因子间的交互效应,为提升江苏农田土壤钾素资源的功能潜力和利用效率提供科学依据。
1 材料与方法 1.1 研究区域概况江苏省位于我国东部,地处长江、淮河下游,介于116°22′ ~ 121°55′E、30°46′ ~ 35°07′N,东滨黄海,西邻安徽省,北接山东省,南与上海市和浙江省毗邻,主要由苏南平原、江淮平原、黄淮平原和东部滨海平原组成。全省成土母质按其成因可分为9类,分别为残积母质、坡积母质、河流冲积母质、黄土和黄土状堆积母质、第四纪红土母质、火山喷发堆积母质、泻湖相沉积母质、三角洲相沉积母质和海相沉积母质[14-15]。根据热量和降水量的地区分布差异,全省农业气候区划可划分为湿润亚热带、湿润暖温带和半湿润暖温带3个气候带类型[16]。江苏经济发达,耕地资源约束矛盾突出,农业集约化程度高,化肥单位面积投入是全国平均水平的1.5倍[17];随着现代农业发展,全省农业种植结构发生了显著变化,粮食种植面积比重有所下降,蔬菜、瓜果种植面积比重上升较快,2000年以来蔬菜类成为江苏除水稻和小麦以外的第三个种植面积最多的作物[18]。
1.2 数据来源数据来自于300个江苏省基本农田质量长期监测点2008年和2015年的监测数据(图 1)。监测点设在永久性耕地的典型地块上,中心点用GPS定位。每年最后一季作物收获后,在监测地块采集土壤样品进行测试分析,每个样品由20个以上取样点采土混合。土壤速效钾含量采用醋酸铵浸提-火焰光度计法进行测定。
![]() |
图 1 江苏省基本农田质量监测点分布 Fig. 1 Geographic location of soil monitoring sites |
每个监测点的采集内容还包括成土母质、种植制度、气候类型和化肥投入等肥力影响因子信息(表 1、表 2)。全省300个监测点的主要成土母质可分为泻湖相沉积母质、三角洲相沉积母质、黄土和黄土状堆积物、河流冲击母质和海相沉积母质5类(残积母质、坡积母质、第四纪红土母质、火山堆积母质发育的土壤监测点均少于3个,予以舍去),种植制度分为稻麦轮作、棉玉旱作、蔬菜、茶果4类;气候类型依据监测点的GPS经纬度,参照买苗等[16]对江苏省气候带划分及其地理分布的研究,将300个监测点归类为湿润亚热带、湿润暖温带和半湿润暖温带。化肥投入的等级以全省稻麦轮作化肥投入量为参考,运用统计软件SPSS“可视分箱”对监测点的化肥投入量进行等级自动划分。
![]() |
表 1 研究区成土母质、气候类型、种植制度 Table 1 Types of soil parent material, climate and cropping system in the study area |
![]() |
表 2 研究区化肥投入种类及等级划分 Table 2 Fertilizer inputs and grading in the study area |
为保证分析结果的可靠性,用“平均值加二倍标准差法”剔除异常值,最终用于分析的样点数为298个。
利用GS+9.0软件对速效钾含量进行半方差分析,运用ArcGIS10.1软件进行基本克里金插值。采用多因素方差分析对土壤速效钾时空变异的影响因子进行研究。以实测点2008—2015年速效钾含量的变化值为因变量,以成土母质、气候类型、种植制度、化肥为自变量,运用SPSS 23.0中一般线性模型TypeⅢ平方和分析模型,在全变量模型的基础上,依据F测验,采用显著水平P = 0.05,逐步淘汰不显著的影响因子,直至模型中所含因子均显著为止,得到最终的优化模型。
2 结果与分析 2.1 土壤速效钾含量变化特征统计结果显示,2008年、2015年江苏省农田土壤速效钾含量均呈偏态分布,经对数变换后符合正态分布。2008年、2015年土壤速效钾含量分别为106.4 mg/kg和107.9 mg/kg,二者差异不显著,两个年份变异系数分别为42.67%和36.43%(表 3)。
![]() |
表 3 江苏省农田土壤速效钾含量统计特征(2008年和2015年) Table 3 Descriptive statistics of soil available potassium in farmlands in Jiangsu Province(2008 and 2015) |
采用对数转换数据对土壤速效钾含量进行半方差分析,2008年土壤速效钾的半方差拟合模型为球状模型,2015年为指数模型,两个模型的拟合决定系数均达到显著水平(表 4)。与2008年相比,2015年土壤速效钾的块金值(C0)和基台值(C0+C)均有所下降,说明速效钾的随机变异和总变异在强度上均有较大程度的减小。就块金值占基台值比(C0 /(C0+C))而言,速效钾在两个时期均具有中等的空间自相关性,但2015年的空间自相关性有所降低,变程也有所减小。
![]() |
表 4 土壤速效钾含量对数值半方差函数理论模型及相关参数(2008年和2015年) Table 4 Semivariance theoretic models and parameters of log available potassium (2008 and 2015) |
图 2A、图 2B分别是2015年和2008年江苏农田土壤速效钾含量的空间分布图,两个年份的空间分布规律基本一致,总体而言,东部高于西部、苏北高于苏南。速效钾含量高于130 mg/kg的土壤集中分布在苏北沿海地区;苏中泰州、如皋等地的土壤速效钾含量最低,低于80 mg/kg。与2008年相比,2015年江苏农田土壤速效钾含量下降的地区主要集中在沿海地区,其中盐城部分地区下降最多,达10.0 ~ 30.5 mg/kg;土壤速效钾含量增加较多的区域主要集中在徐淮部分地区,达10.1 ~ 20.5 mg/kg(图 2C)。
![]() |
(A. 2015年土壤速效钾分布;B. 2008年土壤速效钾分布;C.土壤速效钾含量变化(2008—2015年)) 图 2 2008年和2015年江苏省农田土壤速效钾含量空间分布及变化示意图 Fig. 2 Spatial distribution and change of soil available potassium(AP) in farmlands in Jiangsu Province(2008 and 2015) |
表 5是2008—2015年江苏农田土壤速效钾含量变化影响因子的多因素优化模型。结果显示,成土母质、气候、钾肥、钾肥×成土母质、钾肥×气候对速效钾的含量变化存在显著影响,氮肥、磷肥、种植制度以及作为协变量的土壤有机质对速效钾的含量变化影响不显著。从模型中各因子的Ⅲ类平方和来看,各因子对土壤速效钾含量变化的影响力依次为:钾肥 > 钾肥×成土母质 > 成土母质 > 钾肥×气候 > 气候,钾肥是土壤速效钾含量变化的主要影响因素。
![]() |
表 5 2008—2015年土壤速效钾含量变化影响因子多因素优化模型 Table 5 Multifactor optimal model of soil available potassium content (2008—2015) |
钾肥分别与成土母质、气候存在交互效应,说明在不同的成土母质和气候类型下,钾肥使用量的不同对土壤速效钾含量的变化存在显著差异。总体而言,在钾肥超低投入和低投入下土壤速效钾含量呈现下降,在钾肥中等投入和高投入下土壤速效钾含量呈现上升,说明江苏农田钾肥投入只有达到中等投入以上,才能实现土壤速效钾含量的增加。
图 3显示,钾肥超低投入时,河流冲击物母质土壤的速效钾含量几乎没有变化,海相沉积母质土壤有所下降,其他3种成土母质发育的土壤速效钾含量下降较多;在钾肥低投入、中等投入和高投入下,不同成土母质条件下土壤速效钾含量的变化差异不显著。图 4显示,钾肥超低投入时,湿润亚热带土壤速效钾的减少含量显著低于湿润暖温带和半湿润暖温带;在钾肥低投入、中等投入和高投入下,不同气候类型下土壤速效钾含量的变化差异不显著。
![]() |
(估算边际平均值是指因变量受多个自变量交互影响的情况下,控制其他自变量,单纯在一个自变量作用下的因变量均值[19]) 图 3 钾肥和成土母质对2008—2015年土壤速效钾含量变化的影响 Fig. 3 Effects of potassic fertilizer and parent material on changes in soil available potassium (2008—2015) |
![]() |
图 4 钾肥和气候对2008—2015年土壤速效钾含量变化的影响 Fig. 4 Effects of potassic fertilizer and climate on changes in soil available potassium (2008—2015) |
研究结果显示,2008—2015年江苏农田土壤速效钾含量无显著性变化,这与郁洁等人[20]的研究结果较为一致。2008年和2015年江苏农田土壤速效钾均具有中等的空间自相关性,说明成土母质、气候等自然因素和施肥等人为因素对速效钾空间变异的影响大致相当,这与该地区土壤速效钾含量变化影响因子的多因素方差分析结果较为吻合。与2008年相比,2015年江苏农田土壤速效钾的空间自相关性有所降低,表明人为因素对土壤速效钾含量的影响在增加,自然因素的影响在下降,这与孙永健等[21]的研究结果一致。
钾肥与成土母质存在交互效应,是影响江苏省农田速效钾变化的主要因素。钾肥超低投入时,成土母质是影响土壤钾含量的主要因子[22-23]。土壤矿物中钾长石、微斜长石、水云母富含钾素[24],含有这几类矿物的成土母质发育的土壤钾含量均较高。如苏北地区河流冲击母质以富含钾素的水云母为主;苏南地区主要为三角洲相、泻湖相、黄土和黄土状堆积母质,水云母含量较低,这是苏南、苏北土壤速效钾含量差异的主要原因之一[14]。随着钾肥投入量的逐步增加,不同成土母质和气候带土壤速效钾含量的变化差异不明显,且所有成土母质和气候带土壤的速效钾含量均随钾肥的投入量增加而上升,而钾肥投入低的地区土壤速效钾含量低,因此钾肥成为影响土壤速效钾含量变化的主要因素。本研究中,2008—2015年间江苏盐城地区土壤速效钾含量下降较多,这与钾肥投入不足相关。如盐城建湖县域范围,同期钾肥年均投入折纯仅为0.36 kg/hm2,远远低于全省0.51 kg/hm2的平均水平。
气候也是影响土壤速效钾含量的重要因子,土壤钾的影响主要体现在气温和降水对土壤交换性钾释放。金继运等[25]研究表明,在5、25和40℃3种温度下,随着温度的升高,土壤供钾能力增加、缓冲性能下降。有研究认为,在淹水条件下土壤溶液中有大量的可溶性Fe2+能从黏土复合体中取代出交换性钾,使交换性钾保持在较高的水平[4-5];吴利杰等[26]定量分析了近万年以来气候变化与土壤钾素演变规律,发现年均降水量越大或年均气温越高,土壤中速效钾、缓效钾含量也越高。本研究中,钾肥超低投入下气候对土壤速效钾含量变化的影响显著,其中湿润亚热带土壤速效钾减少量显著低于湿润暖温带和半湿润暖温带,这可能是因为湿润亚热带高温高湿环境更有利于土壤缓效钾的释放,当土壤供钾能力无法满足作物需求时,土壤中难溶性矿物钾及缓效钾会逐渐转化为速效钾,这一现象与金继运等[25]和吴利杰等[26]的研究结果比较吻合。此外,江苏湿润亚热带地区以稻麦轮作制度为主,在干湿交替条件下,由于无定型铝硅酸盐或铁铝氧化物结晶度的提高和水化度的降低,土壤对速效钾的固钾量会降低[1]。徐国华等[3]报道,经过干湿交换处理能促进低钾水平的土壤中交换性钾的释放。这可能也是江苏湿润亚热带地区在钾肥超低投入下土壤速效钾含量下降较少的原因之一。
4 结论与2008年相比,2015年江苏省域农田土壤速效钾平均含量无明显变化;就不同地区而言,同期江苏沿海地区农田土壤速效钾含量下降明显,其中盐城部分地区下降最多。
成土母质、气候、钾肥对2008—2015年江苏农田土壤速效钾含量的变化存在显著影响,其中成土母质与气候、钾肥分别存在交互效应。总体而言,土壤速效钾含量的变化值随着钾肥投入量的增加而增加。在钾肥超低投入下速效钾含量的变化主要受成土母质和气候影响,在钾肥低投入及以上水平,速效钾含量的变化主要受钾肥投入影响,说明钾肥是土壤速效钾含量变化的主要影响因素。
[1] |
徐晓燕, 马毅杰, 张瑞平. 土壤中钾的转化及其与外源钾的相互关系的研究进展[J]. 土壤通报, 2003, 34(5): 489-492 DOI:10.3321/j.issn:0564-3945.2003.05.025 ( ![]() |
[2] |
Feng X, Simpson M J. Temperature and substrate controls on microbial phospholipid fatty acid composition during incubation of grassland soils contrasting in organic matter quality[J]. Soil Biology and Biochemistry, 2009, 41(4): 804-812 DOI:10.1016/j.soilbio.2009.01.020 ( ![]() |
[3] |
徐国华, 鲍士旦, 史瑞和. 生物耗竭土壤的层间钾自然释放及固钾特性[J]. 土壤, 1995(4): 182-185 DOI:10.3321/j.issn:0253-9829.1995.04.004 ( ![]() |
[4] |
Dobremann A, Cassman K G, Stacruzp C, et al. Fertilizer inputs, nutrient balance, and soil nutrient supplying power in intensive, irrigated rice systems. Ⅱ. Effective soil K-supplying capacity[J]. Nutrient Cycling in Agroecosystems, 1996, 46: 11-21 DOI:10.1007/BF00210220 ( ![]() |
[5] |
Dobremann A, Stacruzp C, Cassman K G. Fertilizer input, nutrient balance, and soil nutrient supplying power in intensive, irrigated rice system. Ⅰ. Potassium uptake and K balance[J]. Nutrient Cycling in Agroecosystems, 1996, 46: 1-10 DOI:10.1007/BF00210219 ( ![]() |
[6] |
Srinivasarao C, Kundu S, Ramachandrappa B K, et al. Potassium release characteristics, potassium balance, and fingermillet (Eleusine coracana G.) yield sustainability in a 27-year long experiment on an Alfisol in the semiarid tropical India[J]. Plant and Soil, 2014, 374(1): 315-330 ( ![]() |
[7] |
王绪奎. 江苏耕地质量建设的实践与思考[J]. 江苏农村经济, 2013(12): 29-30 ( ![]() |
[8] |
钱剑峰, 王红. 二十年来无锡市耕地土壤养分含量的变化及培肥措施[J]. 上海农业科技, 2003(4): 13-14 DOI:10.3969/j.issn.1001-0106.2003.04.011 ( ![]() |
[9] |
巩万合, 郁伟, 王志强, 等. 南通市耕地土壤肥力演变趋势及培肥对策[J]. 现代农业科技, 2014(8): 189-191 DOI:10.3969/j.issn.1007-5739.2014.08.119 ( ![]() |
[10] |
徐茂.基于地统计学的江苏省环太湖地区土壤肥力质量演变特征研究[D].南京: 南京农业大学, 2006 http://cdmd.cnki.com.cn/Article/CDMD-10307-2007125653.htm
( ![]() |
[11] |
杨艳丽, 史学正, 于东升, 等. 区域尺度土壤养分空间变异及其影响因素研究[J]. 地理科学, 2008, 28(6): 788-792 DOI:10.3969/j.issn.1000-0690.2008.06.013 ( ![]() |
[12] |
沈德福, 史学正, 吕成文, 等. 江苏沿江地区土壤肥力空间分布及其区域对比研究[J]. 长江流域资源与环境, 2005, 14(3): 316-321 DOI:10.3969/j.issn.1004-8227.2005.03.010 ( ![]() |
[13] |
孙瑞娟, 王德建, 林静慧. 太湖流域土壤肥力演变及原因分析[J]. 土壤, 2006, 38(1): 106-109 DOI:10.3321/j.issn:0253-9829.2006.01.019 ( ![]() |
[14] |
周传槐, 彭补拙. 江苏省的成土母质及其与土属类型关系的研究[J]. 南京大学学报(地理学), 1983(s1): 53-65 ( ![]() |
[15] |
江苏省土壤普查办公室. 江苏土壤[M]. 北京: 中国农业出版社, 1995
( ![]() |
[16] |
买苗, 査书评, 朱宝, 等. 江苏省农业气候资源区划[J]. 安徽农业科学, 2011, 39(8): 4692-4695 DOI:10.3969/j.issn.0517-6611.2011.08.111 ( ![]() |
[17] |
闫湘.我国化肥利用现状与养分资源高效利用研究[D].北京: 中国农业科学院, 2008 http://cdmd.cnki.com.cn/Article/CDMD-82101-2008130416.htm
( ![]() |
[18] |
唐惠燕.基于GIS江苏种植结构演变研究(1949— 2011)[D].南京: 南京农业大学, 2014 http://cdmd.cnki.com.cn/Article/CDMD-10307-1016037632.htm
( ![]() |
[19] |
杜子芳. 多元统计分析[M]. 北京: 清华大学出版社, 2016
( ![]() |
[20] |
郁洁, 王绪奎. 江苏省耕地土壤肥力概述[J]. 江苏农村经济, 2017(4): 45-47 ( ![]() |
[21] |
孙永健, 周蓉蓉, 王长松, 等. 稻麦两熟农田土壤速效钾时空变异及原因分析——以江苏省仪征市为例[J]. 中国生态农业学报, 2008, 16(3): 543-549 ( ![]() |
[22] |
Luebs R E, Stanford G, Scott A D. Relation of available potassium to soil moisture[J]. Soil Science Society of America Journal, 1956, 20(1): 45-50 DOI:10.2136/sssaj1956.03615995002000010011x ( ![]() |
[23] |
章明清, 颜明娟, 林琼, 等. 土壤钾素解吸模型及其特征值的作物效应[J]. 土壤学报, 2008, 45(1): 120-129 DOI:10.3321/j.issn:0564-3929.2008.01.016 ( ![]() |
[24] |
黄昌勇. 土壤学[M]. 北京: 中国农业出版社, 2000
( ![]() |
[25] |
金继运, 高广领, 王泽良, 等. 温度对土壤钾素容量和强度(Q/I)关系的影响[J]. 土壤学报, 1992, 29(2): 137-141 ( ![]() |
[26] |
吴利杰, 于娟, 王成敏, 等. 长尺度气候变化下土壤钾素演变规律的研究——以泾河中游将军村黄土剖面为例[J]. 上海交通大学学报(农业科学版), 2006, 34(1): 69-73 ( ![]() |
2. Nanjing Institute of Environmental Sciences, Ministry of Environmental Protection, Nanjing 210042, China;
3. Jiangsu Station of Farmland Quality and Agro-Environmental Protection, Nanjing 210036, China