钱 争，冯绍元，庄旭东，2，袁成福，3
（1.扬州大学水利科学与工程学院，江苏扬州 225009；
2.宿迁市新沂河调度工程管理处，江苏宿迁 223800；
3.江西农业大学国土资源与环境学院，江西南昌 330045）
河套灌区是我国重要的商品粮、油生产基地，也是北方干旱地区的重要生态屏障［1］。目前河套灌区正面临着土壤盐碱化和水资源短缺等制约灌区农业可持续发展的关键问题。河套灌区盐碱地分布广泛，盐分组成复杂，可通过渗透胁迫、离子毒害和矿质营养缺乏等方式造成盐分胁迫［2］，影响作物生长发育。灌区常采用水利、农业、生物和化学等措施改良盐碱化土地［3］。其中，水利措施中的暗管排水技术是根据“盐随水来，盐随水去”的水盐运移特点，布设适宜埋深和间距的排水暗管，通过控制地下水位和抑制潜水蒸发，改善土壤水盐环境。
近年来，国内外有关土壤盐分分布和组成等方面的研究逐渐增多［4-6］。刘延峰［7］等通过对不同深度土壤盐离子进行相关性分析和主成分分析，评价了焉耆盆地耕地和荒地土壤盐渍化状况。徐友信等［8］在河北近滨海高水位盐碱区治理的研究得出暗管排水可加强盐分离子的淋洗，减缓土壤碱化趋势。郑复乐等［9］研究的表明盐碱地微咸水滴灌较常规黄河水灌溉会促使土壤盐分钠质化发展，生物炭和脱硫石膏可以抑制钠质化问题。窦旭等［10］在乌拉特灌域针对春秋季土壤盐分和离子分布的研究得出Na+、Cl-、盐分含量、Ca2+和Mg2+代表了土壤盐碱化状况。侯晨丽等［11］通过紫花苜蓿和玉米根系对水盐分布的研究得出盐分离子分布及浓度随作物种类而改变。综上可知，已有部分研究从灌溉方式、改良材料、季节变化和作物品种等不同角度分析土壤盐碱化与盐离子组成，但针对河套灌区田间尺度不同土层深度作物生育期内土壤盐碱化及其盐离子特征的研究相对薄弱，有待进一步的研究。
本文基于2020年田间试验观测资料，通过对葵花生育期内土壤样品进行全盐量与盐离子分析，以期得到盐离子与全盐量之间的相关关系以及主导盐碱化的特征因子。旨在为当地暗管排水区域合理布设和改良土壤盐碱化提供理论依据。
1.1 试验区概况
本研究于2020年在内蒙古巴彦淖尔市中国农业大学永济试验基地（107°16′E，40°44′N）进行。该试验基地所在的河套灌区紧邻黄河，平均海拔为1 043.4 m，为中温带半干旱大陆性气候。该地区拥有丰富的光热资源，年日照时间约3 200 h，年蒸发量2 032～3 179 mm，而降雨匮乏，年降雨量90～300 mm，故常年采用黄河水灌溉。2020年葵花生育期内总降雨量121.2 mm，地下水位埋深在0.7～2.0 m范围内变化。
1.2 试验设计
2020年试验区种植作物为食用向日葵RH1号，采用宽窄行种植，种植密度为22 725株/hm2，种植日期为6月2日，收获日期为9月20日。暗管排水工程采用1.5 m埋深、45 m间距、13.3 cm直径的PE 打孔波纹管进行布置。灌溉方式为畦灌，灌溉用水为引黄河水，灌水量用水表控制，灌水定额为1 200 m3/（hm2·次），来水时间由当地管理局统一调配，2020年灌水日期为7月10日，其离子组成见表1。其余田间管理措施与当地一致，具体参考《内蒙古河套灌区盐碱地食用葵花抗盐高产栽培技术规范》［12］。
表1 灌溉及降雨的离子组成 mg/LTab.1 Salt ion composition of irrigation and rainfall
1.3 测定项目与方法
土壤数据主要通过采集并分析各试验小区土壤样品，分析项目主要包括土壤全盐量、颗粒级配等。其中，采用马尔文帕纳科激光粒度仪（Mastersizer 3 000）测定颗粒级配；
采用环刀法测定干体积质量；
采用电导率仪（DDSJ-308A）测定纯水与土壤样品配制的土水比1∶5 土壤浸提液的电导率值EC1∶5，并运用内蒙古河套灌区经验公式（TS=2.882EC1∶5+0.183）换算成相应的土壤全盐量TS［13］。试验区土壤机械组成及水力参数见表2。主要盐离子测定方法见表3。
表2 试验区土壤机械组成Tab.2 Mechanical composition of soil in experimental area
表3 主要盐离子测定方法Tab.3 Method for determining salt-based ions
1.4 数据处理和分析方法
采用Excel 2016 进行数据处理、初步分析和绘图，使用SPSS 23.0对数据进行显著性分析、相关性分析和主成分分析。
2.1 土壤盐离子剖面分布特征
将研究区的0～1 m 剖面土壤划分为表层（0～20 cm）、中层（20～40 cm）和深层（40～100 cm），葵花生育期内各层土壤全盐量TS、pH及盐离子含量如图1所示。
图1 各土层土壤全盐量、pH值及盐离子含量Fig.1 Total salt content，pH value and salt ion content of each soil layer
生育期内各层土壤pH 值均呈碱性，其变化趋势一致，均有所增大，平均值从生育初期的7.79 增大到生育末期的8.22，均为碱性土。土壤全盐量TS呈“Γ”型分布，在生育期内土壤盐分积聚性加强，且各土层TS均呈现波动增大。
在生育初期，各层土壤阴离子以HCO3-为主。就整个0～1 m剖面而言，HCO3-、Cl-和SO42-含量占阴离子百分比分别为45.02%、28.34%和26.63%；
在葵花生长进程中，各层土壤同一阴离子变化趋势相同，HCO3-含量减少且占比降低，SO42-含量增加且占比升高，而Cl-含量增加但占比基本不变；
在生育末期，各层土壤阴离子以SO42-为主，0～1 m剖面SO42-、Cl-和HCO3-占比分别为55.62%、27.59%和16.79%，表明研究区土壤盐分组成类型经历了从生育初期的苏打型到生育末期的硫酸盐的转变。
各层土壤阳离子始终以K++Na+为主，生育初期0～1 m 剖面土壤K++Na+、Ca2+和Mg2+的阳离子占比分别为65.26%、23.52%和11.22%，生育末期分别为74.76%、12.52%和12.72%；
葵花生长进程中，各层K++Na+和Mg2+总体表现为含量增加且占比升高，Ca2+含量基本不变而占比降低。考虑到K+含量远低于Na+含量，因此在葵花生育期内研究区的土壤盐分以钠盐为主。
2.2 土壤盐离子相关性分析
通过对土壤盐离子进行相关性分析，可以揭示盐分在土体中的分布和存在状态，在一定程度上反映盐离子的同源性、差异性和组合情况。研究区的各土层TS 和不同离子的相关性分析结果如表4所示。
表4 各土层土壤盐离子相关系数表Tab.4 Correlation coefficients of soil ions in different soil layers
在各土层中，土壤全盐量均与阴离子中的SO42-均呈极显著正相关（P＜0.01），相关系数最大且相关系数随深度增加而减少，说明硫酸盐在表层上积聚强烈，盐分以向上运动为主［14］；
各土层土壤全盐量与Cl-均呈极显著正相关，说明Cl-离子是促进1 m 土体中土壤盐分变化的重要原因之一；
各层土壤全盐量与Cl-的相关系数均小于SO42-，说明土壤盐分中硫酸盐较氯化物有更高的含量，这与研究区蒸发强烈、有效灌溉与降雨的频率集中以及地质构造等因素有关，多次的积盐和脱盐过程使得更易溶于水的氯化物并未在土壤表层有效积聚；
全盐量与HCO3-的相关性较差，且在20～40 cm 土层深度中呈显著负相关（P＜0.05），说明HCO3-离子并不是促进土壤盐分变化的主导因素。在各土层中，土壤全盐量均与K++Na+离子的相关性最高，均呈极显著正相关，说明K++Na+离子是促进各层土壤盐分变化的主导因素之一；
全盐量与Mg2+的相关性较好，且在0～40 cm 土壤深度中呈极显著正相关；
全盐量与Ca2+的相关性较弱，仅在20～40 cm土层深度呈显著相关。
阴离子中，各土层土壤中的HCO3-与其他离子的相关性较差，其原因是HCO3-在水中有较高的占比，在降雨和灌溉后，其他离子受淋洗作用含量降低，而土壤碳酸钙溶解和生物厌氧呼吸导致HCO3-含量上升，这也导致HCO3-和Ca2+相关性大于其他离子；
在0～20 cm 土层，Cl-与SO42-、K++Na+和Mg2+均具有较好的相关性，而随着土壤深度增加，Cl-仅与K++Na+呈显著相关。阳离子中，仅K++Na+和Mg2+在0～20 cm 土层深度呈极显著相关，而在其余土层深度阳离子之间的相关性较弱。
2.3 土壤盐离子主成分分析
主成分分析法是将众多的指标化为少数几个相互独立的综合指标（主成分）的一种多元统计方法［15］。本文采用累计方差贡献率≥85%的标准确定各土层的主成分数目，找出具有代表性的主导因子，使得选取的主成分能够较为全面地表达原有数据信息，以正确评价研究区土壤盐碱化状况，结果如表5所示。
表5 各土层土壤盐离子主成分分析表Tab.5 Principal component analysis of soil ions in different soil layers
在表层土壤中，Cl-、SO42-、K++Na+、Mg2+在第一主成分上的荷载较大，在0.878 以上，特征值为4.956，方差贡献率为70.806%，第一主成分主要包含土壤盐渍化（中性盐）的信息，可作为土壤盐渍化的特征因子。在第一主成分中，与其相关程度最高的阴离子为SO42-、阳离子为K++Na+，说明该土层盐分特征在一定程度上受Na2SO4的影响最大［16］。而第二主成分的特征值为1.215，方差贡献率为17.354%，主成分为与第一主成分离子相关性较差的HCO3-，HCO3-主要表征土壤碱性离子，说明第二主成分可作为土壤碱性的特征因子。同理可知，中层和深层土壤的第一主成分可作为土壤盐渍化的特征因子，第二主成分可作为土壤碱性的特征因子，在一定程度上受Na2SO4的影响最大。
综合各土层各项指标对各主成分的系数大小，各主成分方差贡献率，以及研究区的实际情况，可以认为全盐量、SO42-、Cl-、HCO3-、K++Na+作为研究区盐碱化状况的特征因子。
R-分析可用于性状或指标的分类［17］。利用R-分析，对土壤盐碱化状况指标在第一、第二主成分构成的空间作散点图，直观反映各项指标间的相互关系，指标之间越靠近表示关系越密切，结果如图2所示。
从图2 中可根据各项指标的离散情况进行归类。总体来看，HCO3-可以单独归类，代表土壤碱化特征；
全盐量TS、SO42-、Cl-、K++Na+可以归为一类，代表土壤盐分组成及土壤盐化趋势。
图2 各土层R-分析图Fig.2 R-analysis diagram of each soil layer
各土层主成分荷载向量除以相应主成分特征值的算术平方根得到不同土层土壤盐离子系数矩阵，结果如表6 所示。第一主成分得分F1、第二主成分得分F2和第三主成分得分F3可用于表征研究区2020年葵花生育期内各土层土壤可溶性盐离子特征，按公式（1）计算各主成分得分，并按公式（2）计算综合得分Y。
表6 各土层土壤盐离子系数矩阵Tab.6 Matrix of ion coefficient of soil salt base in each soil layer
式中：F为各主成分得分；
i为主成分的指标个数；
xi为第i个主成分指标系数；
Xi为第i个主成分指标的含量；
Y为综合得分；
j为主成分个数；
Fj为第j个主成分的得分；
fj为第j个主成分的方差贡献率在总累计方差贡献率中的占比。
以0～20 cm 表层土壤为例，各主成分条件下的得分与综合得分表达式如下所示：
计算得到2020年研究区葵花生育期内不同取样时间下的土壤得分，结果如图3所示。
图3 生育期内各土层土壤主成分得分图Fig.3 Soil PCA score of each soil layer during growth period
在表层中，由第一主成分得分F1可知，生育期内土壤受盐化影响程度总体增加；
由第二主成分得分F2可知，土壤受碱化影响程度总体降低，土壤受碱性影响最小的时间为生育末期；
由综合得分Y可知，土壤受盐化和碱化的综合影响总体呈波动上升。同理可知，生育期内中层土壤受盐碱化综合影响程度总体呈波动上升。深层土壤受盐化影响程度总体呈波动上升，受碱化影响程度总体下降，受盐碱化综合影响程度呈先增加后平稳的状态。上述结果表明，通过各主成分得分和综合得分可以反映不同取样时间的各土层土壤受盐化和碱化的影响程度，为治理该地区土壤盐碱化提供理论依据。
（1）研究区葵花生育期内，0～1 m 剖面的土壤盐分组成类型以钠盐为主，经历了从生育初期的苏打型到生育末期的硫酸盐的转变；
生育期始末各层土壤中，阴离子从HCO3-为主转变为以SO42-为主，阳离子始终以K++Na+为主。
（2）随土壤剖面深度增加，土壤全盐量与可溶性盐离子的相关性发生改变，相关性最强的阴离子为SO42-，阳离子为K++Na+。
（3）在各层土壤中，HCO3-可代表土壤碱化特征，SO42-、Cl-、K++Na+可代表土壤盐分组成及土壤盐化趋势。
猜你喜欢 盐分生育期葵花 基于单片机原理的土壤盐分实时观测系统安全与环境工程(2021年5期)2021-10-08葵花田阅读与作文（小学高年级版）(2020年8期)2020-09-12葵花加油站小天使·二年级语数英综合(2019年10期)2019-11-08海水这么咸的原因祝您健康·文摘版(2019年3期)2019-06-11葵花点穴手小演奏家(2018年9期)2018-12-062014年峄城区花生生育期气象条件分析新农村(2016年12期)2017-01-12辣椒不同生育期光合作用比较长江蔬菜·学术版(2016年12期)2017-01-12夜晚的葵花扬子江(2016年3期)2016-06-17摄影欣赏读者·校园版(2016年6期)2016-03-07盐碱地变良田/盐碱地为什么不立苗？中国青年(1966年2期)1966-08-17