松嫩平原苏打盐渍化旱田面积约80万hm²,主要以轻、中度盐渍化为主。土壤空间分布上表现为盐渍化程度存在较大的空间异质性^[1],土壤盐渍化的空间异质性导致了盐渍化旱田作物产量具有较大差异。尤其是旱田中存在盐渍化程度较高的盐碱斑块,增加了农田耕作和作物种植的难度,多年来一直是“种地不见苗”。盐渍化旱田土壤的这种不均一性,严重影响了当地的农业生产,极大地限制了盐渍化旱田的大面积利用和机械化管理,是限制该区旱田作物高产的主要障碍因子^[2]。因此分析苏打盐渍化旱田土壤盐渍化空间分布特征,实现盐渍化程度快速诊断和土壤盐渍化空间异质性评价,提出有针对性的改良和管理措施,对于提高当地农业粮食产量和保持农业可持续发展具有重要的理论和现实意义。

盐渍化旱田土壤改良目标是耕层土壤盐渍化指标达到作物正常生长范围。就一具体田块或区域而言,要实现大规模盐碱地农业高产高效,必须解决耕层土壤盐渍化空间异质性问题。诸多学者开展了土壤盐渍化性质的空间变异研究^[3-5],提出了针对土壤性质空间差异实施有针对性的管理技术^[6-8]。有研究发现,土壤表观电导率与盐渍化程度等理化性质密切相关,并用土壤表观电导率指示土壤性质空间变异,指导田间灌排等管理^[9-10]。因此,土壤表观电导率作为一种非扰动、快速、可靠、操作简单的土壤田间实时测试参数,已经逐渐被广泛应用,成为精准农业研究的重要工具^{[4, 11]}。针对不同盐碱地类型,建立土壤表观电导率与土壤盐碱性质相关方程,并将盐碱空间分布特征与田间土壤改良、水肥管理等农艺技术紧密结合,形成精准化管理技术是国际上研究的热点^{[7, 9]}。然而以往研究主要集中在氯化物盐渍土区^[12]。由于氯化物盐渍土的土壤特性与苏打盐渍土完全不同,其表观电导率所指示的盐碱指标也完全不同。苏打盐渍土旱田空间变异研究报道较少,仅利用经典统计方法研究旱田空间变异^[13],该方法对空间变异研究具有较大的随机性。如何通过方便、快捷的土壤表观电导率测定,对苏打盐渍土旱田盐碱空间异质性进行定量表征,是盐渍化旱田精准管理的关键问题。因此,本研究采用磁感应电导率仪(EM38)对典型苏打盐渍化旱田土壤表观电导率进行监测,利用地统计方法量化土壤盐碱的空间变异特征,并结合定位采样,建立苏打盐渍土壤表观电导率与土壤盐碱化指标方程,对松嫩平原苏打盐渍化旱田盐碱程度进行快速诊断,将空间变异理论与盐碱地改良实践相结合,实施局部清除或降低土壤盐碱危害,实现苏打盐渍土均质化改良和利用,对于扩大有效耕地面积具有重要现实意义。

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

研究区位于松嫩平原吉林省西部大安市乐胜乡。地理坐标为123°21′31″—123°40′58″E,45°26′28″—45°29′27″N。该区盐渍土广泛分布,土壤属内陆型苏打盐渍土,盐分组成中以苏打(NaHCO₃与Na₂CO₃)为主,pH值多在8.5以上,是松嫩平原轻度和中度盐渍化低产旱田的典型分布区。土壤类型为壤土和壤黏土。该区属于中温带大陆性季风区,多年平均降水量413 mm,年蒸发量较大,为降雨的2—3倍。年均气温4.3℃,无霜期约140 d^[14]。

1.2 方法与试验设计

选取盐碱程度不均一的苏打盐渍化旱田典型地块作为调查小区。小区长300 m,宽160 m,面积为4.8 hm²。利用GPS和EM38大地电导仪相结合,将调查小区的盐碱旱田画成网格,网格大小为20 m×20 m,在网格节点设定观测点,作为土壤表观电导率的观测点,共144个观测点。本研究采用EM38电导仪(加拿大Geonics公司)对调查区的表观电导率进行测定。 EM38利用电磁感应的原理,可以实现非接触地快速精确测量土壤的电导率。它有两种探测模式,垂直偶极模式的探测深度达1.5 m,水平偶极模式的探测深度达0.75 m。依此测得的土壤表观电导率也分为垂直方向电导率和水平方向电导率^[15]。

调查在2013年4月20日进行。在调查区内,充分考虑不同土壤表观电导率的空间分布,使样点在调查区内合理散布的原则下,利用GPS定位,在土壤表观电导率的监测点上共布置了46个测试点,采集土壤样品,测定土壤样品的盐碱指标。土壤取样点的取样深度为0—30 cm。对土壤样品进行室内分析,采用1∶5土水比例浸提液，测定土壤可溶性盐分离子(K⁺,Na⁺,Ca²⁺,Mg²⁺,CO₃^2-,HCO₃^-,SO₄^2-,Cl^-),土壤可溶性离子总量(SC),EC_1∶5和pH_1∶5^[16],计算土壤钠吸附比(SAR)^[17]。其中SAR计算公式如下：

钠吸附比(SAR)是溶液中Na⁺离子同Ca²⁺、Mg²⁺离子的平均浓度的平方根的比值,即

(1)

式中,Na⁺、Ca²⁺、Mg²⁺离子浓度单位为mmolc/L,SAR的单位为(mmolc/L)^1/2。

1.3 数据处理方法

样本的描述性分析采用SPSS 11.5软件进行,用Kolmogorov-Smirnov (K-S) 检验数据是否呈正态分布;半方差函数应用地统计学处理软件GS⁺for windows 5.3b进行计算。首先在GS⁺软件中对不符合正态分布的数据源进行对数转换;其次进行半方差函数的计算、模拟、分析和检验;利用Surfer软件进行Kriging插值,生成2维空间分布图。

2 结果与分析 2.1 经典统计分析旱田土壤盐渍化空间变异 2.1.1 土壤表观电导率统计特征

EC_h和EC_v分别为水平和垂直方向测定的土壤表观电导率。水平方式敏感程度在地表最高,并随深度降低。垂直方式在近地表敏感程度很低,随深度增加,敏感程度增高。用EC_h/EC_v反映两种方式在近地表敏感程度的区别。如果水平方式的读数大于垂直方式的读数,即(EC_h/EC_v)＞1,说明近地表土壤盐碱大于下层土壤盐碱。如果两种方式的读数相差很小,说明土壤剖面盐碱较一致。由表 1可知,EC_h最大值为74 mS/m,是最小值的18倍。而EC_v最大值为64 mS/m,是最小值的6倍。EC_h和EC_v均值分别为15.9 mS/m和23.5 mS/m,说明该调查区以轻度盐渍土为主,在轻度盐渍土上分布着中、重度盐渍土斑块。由表 1可知,调查区土壤表观电导率EC_v的平均值为23.5 mS/m,较EC_h高。从EC_h和EC_v的比值可知,其均值为0.63,标准偏差为0.16,说明调查区域土壤表层的电导率小于深层土壤电导率,说明该区深层盐碱高于表层盐碱。变异系数(CV) 反映相对变异,即随机变量的离散程度。已有研究将空间变异系数分为弱空间变异(CV≤ 0.1),中等空间变异(0.1<CV<1)和强空间变异(CV≥1)^[18]。调查区中EC_v和EC_h的变异系数范围为0.1 <CV< 1,为中等空间变异。EC_v的变异系数为0.31,较EC_h小0.24,说明表层土壤盐渍化空间变异大于深层土壤盐渍化空间变异。因此,在生产实践中,应主要评估表层土壤(作物根区土壤)盐碱空间分布,将不同盐碱化土壤进行定位分区,有针对性的降低根区土壤盐碱。同时,应注意盐渍化旱田的田间管理,防止深层土壤盐碱向上运移,造成表层土壤盐碱加重,影响作物产量。

2.2 地统计分析旱田土壤盐渍化空间变异 2.2.1 盐渍化旱田土壤表观电导率的半方差分析

地统计方法有效的前提是数据服从正态分布。因此,数据必须进行正态分布检验。由表 2可知,通过Kolmogorov-Smirnov法对土壤表观电导率分布进行正态检验(P＜0.05,2-tailed),发现水平方向和垂直方向土壤表观电导率均不符合正态分布,对EC_v和EC_h进行对数转换,转换后EC_h(P=0.338)、EC_v(P=0.119)均符合正态分布,满足地统计学分析的要求,可以用地统计学的方法来定量化研究土壤表观电导率的空间变异。变异函数计算均采用对数转化后的数据。

在半方差模型模拟过程中,最优模型要求残差平方和(RSS)最小,同时决定系数(R²)最大。由表 3可知,综合考虑半方差指标,球状模型为EC_h的半方差模拟最优模型。而EC_v半方差模拟的最优模型为指数模型。由于模型的块金值均较小,该模型受实验误差等随机性因素影响较小。对于EC_h有效变程为91.9 m,观测值之间大于91.9 m时,观测值是相互独立的,而观测值之间在91.9 m以内存在着空间自相关。EC_v的有效变程较EC_h大,有效变程为198.9 m,说明观测值之间大于198.9 m时,观测值是相互独立的,而观测值之间在198.9 m以内存在着空间自相关。由块金值与基台值之比可知,当比例>75%时,说明变量具有弱空间自相关;比例在25%—75%之间,变量具有中等空间自相关;当比例<25%时,变量具有强空间自相关^[19]。EC_h半方差模型的块金值与基台值的比值为6.5%(<25%),说明变量具有强空间自相关,即水平方向土壤表观电导率EC_h主要由地形地貌和水文状况等结构性因素引起的空间变异。而EC_v半方差模型的块金值与基台值的比值为47.8%(在25%—75%之间),因此变量具有中等的空间自相关,即垂直方向土壤表观电导率EC_v的空间变异是由结构性因素和随机性因素共同影响的。因此,对旱田土壤盐渍化程度评估时,应选用水平方向土壤表观电导率EC_h作为评估指标,来研究盐渍化空间变异特征,分析地形地貌和水文状况等结构性影响因素与土壤盐渍化的关系,进而实现对盐渍化旱田有针对性的改良。

2.2.2 基于Kriging插值的盐渍化旱田表观电导率空间分布

依据半方差函数模拟模型,应用克立格法(Kriging)进行插值,绘制经对数转换后的土壤表观电导率EC_v和EC_h的空间分布图,见图 1。经对数转换后的水平与垂直表观电导率EC_h和EC_v空间分布非常相似,均表现为不同表观电导率的土壤呈条带状或斑块状镶嵌分布。试验小区土壤水平和垂直方向均表现为西北和南部各有一较大条带状的高盐碱分布区,另在东北和北部中间各有一较小的高盐碱分布区。而在调查区的中部和东南部有斑块状的低盐碱区。盐渍化旱田地表呈现微地貌特征,盐碱严重的区域地势稍高,无作物生长,为光板地;地势稍低区域盐碱化程度稍高,作物产量较低;在调查区最低部位,为轻度盐渍化和非盐碱化区,作物长势良好。因此土壤盐渍化程度的不同直接影响作物长势,使作物产量表现为极大地不均匀性。

2.3 土壤表观电导率与相关因子分析 2.3.1 土壤表观电导率与土壤盐碱指标相关分析

将水平方向和垂直方向土壤表观电导率与土壤盐碱指标进行Pearson分析。土壤盐碱指标包括土壤可溶性盐离子(K⁺,Na⁺,Ca²⁺,Mg²⁺,CO₃^2-,HCO₃^-,SO₄^2-,Cl^-),EC_1∶5,pH_1∶5,SAR和SC,土壤盐碱指标与土壤表观电导率经检验后具有相关性的见表 4。EC_h和EC_v均与盐碱指标呈极显著相关(P＜0.01)。土壤表观电导率能够反映土壤的盐碱程度,可以用土壤表观电导率来指示土壤盐碱程度,实现土壤盐渍化的快速诊断。比较EC_v和EC_h,结果表明,水平方向土壤表观电导率(EC_h)与盐碱指标相关系数均高于垂直方向土壤表观电导率(EC_v),该结果与其他研究结果是一致的^[20]。由于EM38水平方式敏感程度在地表最高,土壤表观电导率EC_h与近地表层(0—30 cm)的土壤盐碱指标相关性好。因此,在松嫩平原利用EM38测定土壤表观电导率,只需采用水平方向土壤表观电导率对土壤进行监测,即可达到指示土壤盐碱程度的目的。

2.3.2 土壤表观电导率与相关因子的回归分析

为了更好的研究土壤表观电导率与土壤盐碱指标的关系,对土壤表观电导率与相关因子进行回归分析。由表 5可知,EC_h与EC_1∶5、pH_1∶5、SAR、SC和Na⁺呈线性分布关系,决定系数范围为0.6—0.82,具有较好的模拟效果(P=0.000)。因此,可以利用EC_h来计算SAR、SC、EC_1∶5、pH_1∶5和Na⁺的值。由表 5可知,EC_v与EC_1∶5、pH_1∶5、SAR、SC和Na⁺也呈线性分布关系,决定系数在0.5—0.7之间,具有较好的模拟效果(P=0.000)。比较EC_h和EC_v与盐碱指标回归模型的决定系数发现,EC_h模型的决定系数均大于EC_v。实际应用中,只需使用EC_h来计算表层土壤(0—30 cm)盐碱指标,指示土壤盐碱程度。

3 讨论与结论 3.1 讨论

EM38测得的水平方向土壤表观电导率能够指示苏打盐渍化旱田土壤盐碱程度。土壤表观电导率具有空间变异性,表现为不同表观电导率的土壤呈条带状或斑块状镶嵌分布,体现盐渍化旱田土壤盐碱的空间分异特征。松嫩平原盐渍化旱田多数由盐碱退化草地开垦而成的,地表呈现微地貌差异,进而影响到水盐运动和再分配过程,在相对高起的微地形上,同时存在着纵向和横向两个方面的湿度差。根据土壤毛管水由湿度大的土层向湿度小的土层移动规律,微高地上既有毛管上升水流的补给,也有毛管侧向水流的补给,当水分沿土壤毛细管由下而上,由缓坡低处向缓坡高处移动过程中,盐分也随之迁移,并通过蒸发而表聚。蒸发量愈大,水分的补给愈快,盐分积累也愈多^[21],成为盐碱化严重区域^[22],甚至为雨养旱田内的“光板地”;而在地势最低区域为轻度盐渍化和非盐渍化区。因此,调查小区在微高地处土壤盐渍化程度明显高于低地,形成所谓“一步三换土”^[23]。松嫩平原旱田土壤盐渍化程度空间异质性,给土壤改良、规模化利用和作物生产带来严重障碍。明确需要改良土壤的面积和比例,是经济、合理、高效利用盐渍化旱田的关键。

将土壤盐分空间变异理论研究与盐碱地改良实践相结合,利用水平方向土壤表观电导率对旱田盐渍化进行评估,有利于实现盐渍化旱田快速低成本改良。不同作物适应盐碱能力不同。对于盐碱化程度较低、空间变异较小的盐渍化旱田,可以利用土壤表观电导率指示不同作物的耐盐碱阈值,针对土壤盐碱程度种植相应作物,是合理利用盐渍化土壤的一种重要方式。对于盐渍化较重、空间变异显著,斑块盐碱程度较重的区域,种植耐盐碱作物无法实现调节,这部分区域应进行定位分区改良,定点清除或降低土壤盐碱,使土壤达到适合作物生长。局部盐碱斑块,由于大量交换性Na⁺的存在,导致土壤黏粒分散,土壤渗透性差,且随着交换性Na⁺的增加,土壤性质不断恶化,影响作物生长。针对盐碱斑块,最好的方法是加入含有二价阳离子的改良剂,置换土壤中的交换性Na⁺,达到改良目的。当然,种植的作物不同,耐盐碱能力不同,需要改良剂用量也不同。因此,改良利用此类盐渍化旱田应综合考虑土壤表观电导率、作物、土壤盐碱之间的关系,以期对不同盐渍化程度的区域进行改良,达到作物正常生长,实现精准化管理技术。

3.2 结论

经典统计分析土壤水平方向表观电导率EC_h与垂直方向表观电导率EC_v均为中度空间变异。地统计分析表明经对数转换后的土壤水平方向表观电导率EC_h具有强空间自相关,其变异特征主要是由地形地貌和水文状况等结构性因素引起的;而对数转换后的垂直方向土壤表观电导率ECv具有中等的空间相关性,即空间变异是由结构性因素和随机性因素共同影响的。EC_h、EC_v空间变异的半方差模拟最优模型为球状模型和指数模型。

相关分析和回归分析均表明土壤表观电导率与盐碱指标呈正相关关系(P＜0.05)和线性关系,且EC_h测量土壤表观电导率效果好于EC_v,在生产实践中用EC_h能够很好的反映土壤盐碱状况。因此,可以用水平方向土壤表观电导率EC_h来计算苏打盐渍化旱田土壤盐碱指标,对土壤盐渍化程度进行评估,指导旱田土壤盐渍化管理。