《Agronomy Journal》：GIS–TOPSIS assessment of soil and water constraints on sunflower in arid Algeria
编辑推荐：
在日益干旱的背景下，识别能够支撑油料作物生产的土地，是增强北非旱地农业韧性的核心任务。研究人员构建了一个经统计验证、空间显式的土地适宜性评估框架，将地理信息系统（GIS）与技术逼近理想解排序法（TOPSIS）相结合，覆盖阿尔及利亚西南部Ain Sefra盆地M
在日益干旱的背景下，识别能够支撑油料作物生产的土地，是增强北非旱地农业韧性的核心任务。研究人员构建了一个经统计验证、空间显式的土地适宜性评估框架，将地理信息系统（GIS）与技术逼近理想解排序法（TOPSIS）相结合，覆盖阿尔及利亚西南部Ain Sefra盆地Maader灌区（15.8 km2）内34个地理定位样点。八项土壤–水文指标——包括阳离子交换量（CEC）、有机质（OM）、交换性钠百分比（ESP）、质地、石膏含量、坡度、pH值及灌溉水质指数（IWQI）——经层次分析法（AHP）赋权并标准化（一致性比率 = 0.07）。TOPSIS贴近系数（Ci；相对逼近理想解的程度）介于0.37至0.76之间，据此将研究区划分为1个高度适宜样点（占2.9%）、29个中度适宜样点（占85.3%）及4个边缘适宜样点（占11.8%）；无任何样点被判定为不适宜。统计验证——包括Kruskal–Wallis检验（CEC：p = 0.00056；OM：p = 0.00313；ESP：p = 0.0272）、Pearson相关分析（CEC：r = +0.65；OM：r = +0.62；ESP：r = ?0.70）、主成分分析（PCA）及方差膨胀因子（VIF < 2.5）——表明CEC、OM和ESP是主导且相互独立的适宜性驱动因子。在空间分布上，中部钙镁质地带具有高CEC、天然石膏含量及平缓坡度，呈现出向日葵（Helianthus annuus L.）集约化种植的“即扩即用”潜力（实测产量约2.8 t ha?1），而邻近wadi河床的钠化洼地则需要定向石膏改良（≥6 t ha?1）及地下排水措施。敏感性测试在因子权重±20%扰动下证实了模型稳健性，适宜性等级变化低于15%。总体而言，该GIS–TOPSIS方法不仅具有诊断功能，更提供了一个透明、可扩展且可直接投资的平台，用于指导气候胁迫下的油料作物扩张。通过将定量严谨性与可操作的空间输出相结合，该方法支撑了阿尔及利亚农业自给自足议程，并为该区域旱地适宜性规划提供了可移植的模板。
论文解读
一、研究背景与问题提出
全球粮食安全正面临日益严峻的压力。人口增长、气候变化、土地退化与快速城镇化不断加剧有限农业用地和水资源之间的竞争。联合国粮农组织（FAO）的预测表明，到2050年全球农产品需求可能增长60%–70%，部分估计认为作物产量需近乎翻倍才能养活90–100亿人口。在此背景下，摆脱高投入单一作物种植模式、转向可持续的气候智慧型农业系统成为迫切需求，而空间显式的土地资源规划则是实现这一转型的关键支撑。
阿尔及利亚西南部的Ain Sefra盆地是典型的干旱–半干旱地区，面临水资源匮乏、土壤贫瘠和气候变率显著等多重农业约束。该地区地下水和土壤资源以往多被独立评估，缺乏将两者综合纳入农业适宜性评价的系统框架。尤其值得关注的是，灌溉水质与土壤属性之间的交互作用——特别是盐渍化动态——在干旱农业生态系统中至关重要，但迄今未得到充分研究。同时，多准则决策分析（MCDA）方法虽在土地适宜性研究中日益普及，但其与灌溉专项指标（如IWQI）在GIS框架内的整合在阿尔及利亚尚属空白，针对向日葵等特定作物的应用更为罕见。
向日葵（Helianthus annuus L.）凭借其强耐旱性和在食品、生物燃料及工业领域的广泛经济用途，成为干旱环境中的优选作物。阿尔及利亚农业部已将油料作物种植列为国家农业议程的优先事项，设定到2028年实现30万公顷油料作物的目标，以减少进口依赖、增强粮食安全。在此国家战略背景下，对Maader灌区向日葵种植土地适宜性开展严格评估，具有重要的现实意义和政策价值。
二、研究设计与核心发现
本研究在阿尔及利亚西南部Ain Sefra盆地Maader灌区（15.8 km2，位于Tiout绿洲西北2 km处，海拔1130–1210 m）开展。研究区覆盖钙镁质土壤（包括黑色石灰土、石灰性褐土及含石膏结壳的石灰性褐土），下伏古老及中更新世 glacis 阶地，并伴有砾质坡积物和沿wadi河床分布的近代砂质冲积物。长期（1990–2022）气候记录显示年均降水量仅165.6 mm，1月均温5.9°C，7月均温30.4°C，夏季FAO-56参考蒸散量（ETo）约达256 mm/月，三种独立干旱指数均将该地区归入干旱等级。
研究人员于2024年3月至4月间，采用分层随机样线网络布设34个地理定位采样点，覆盖从高钙镁质glacis、坡积裙到wadi沿岸远端砂质冲积物的全部主要地貌单元。每样点采集5个0–20 cm土芯混合样，并同步采集34口井及渠道取水口的水样。土壤样本经风干、过筛后测定颗粒组成、饱和浸提液电导率（ECe）、pH、CEC、ESP、有机质（OM）及石膏含量；水样测定电导率（EC）、钠吸附比（SAR）、残余碳酸钠（RSC）及主要阴阳离子，并按Meireles等（2010）方法整合为灌溉水质指数（IWQI）。
三、主要研究方法
本研究采用的核心技术方法包括：其一，层次分析法（AHP）——经两轮德尔菲专家咨询，由9位区域土壤与水资源专家在Saaty 1–9标度上对八项指标进行两两比较赋权，一致性比率为0.07；其二，技术逼近理想解排序法（TOPSIS）——通过构建决策矩阵、向量归一化、加权归一化、确定正负理想解、计算欧氏距离及贴近系数（Ci）等九个步骤，将多维土壤–水文数据融合为单一适宜性评分；其三，地理信息系统（GIS）空间插值——采用反距离权重（IDW）法将离散样点数据生成连续专题面，坡度图层由数字高程模型（DEM）经ArcGIS坡度函数提取；其四，统计验证框架——包括Kruskal–Wallis非参数检验、Pearson相关分析、主成分分析（PCA）及方差膨胀因子（VIF）诊断，用以量化各准则对Ci的独立及联合影响。
四、研究结果
标准化决策矩阵。 经向量归一化后，34个样点的八项指标呈现出清晰的分异格局。CEC归一化值从0.054延伸至0.279，IWQI从0.146至0.212，而ESP展现出数量级的跨度（0.036–0.481），成为总方差的主要来源。CEC最高的五个样点（m33、m2、m18、m27、m3）空间集中于中部钙镁质glacis，黏粒含量35%–42%，以蒙脱石–伊利石间层矿物为主；CEC最低的三个样点（m13、m12、m19）位于wadi河床边缘的砂质–石膏质透镜体上，以石英和海泡石为主导矿物。IWQI最高值出现在深部碳酸盐含水层（>90 m）供水的m18和m29，EC < 1.2 dS m?1，SAR ≈ 2.4；而浅层潜水系统（<30 m）供水的m12和m13虽IWQI > 0.18，却伴随极端ESP。仅m11、m29和m30三个样点OM归一化值超过0.30，对应十年以上>10 Mg ha?1 yr?1的堆肥羊粪施用历史。ESP峰值（m13为0.481、m12为0.371、m29为0.307）均出现在地下水位1.2 m以内的点位，证实毛细上升是钠的主要运移途径。质地呈双峰分布：坡积裙上的粗砂质壤土（归一化0.068，坡度5%–8%）与远端冲积地上的细壤土至黏壤土（≥0.204）形成鲜明对照。所有土壤均为碱性（pH 7.8–9.0），石膏浓度在m1–m4处达到峰值（归一化0.357）。
加权归一化矩阵。 应用AHP权重（CEC和IWQI各0.145；OM、ESP和pH各0.127；坡度、质地和石膏各0.109）后，高CEC与优质灌溉水组合的样点（m33、m18、m27、m3、m2）从这两项准则获得超过40%的总权重；而m13的极端ESP经加权后将其综合值压至0.01以下，尽管其IWQI尚可。加权引起的排名变动分为三组：m18和m33因CEC与IWQI协同上升4位；m12、m13和m29因高ESP下降3–5位。弹性分析表明，加权ESP每增加0.01对综合得分的压低效应约为加权OM同等增幅的3倍。
贴近系数与适宜性等级。 Ci值介于0.37至0.76之间，将34个样点划分为1个高度适宜（2.94%）、29个中度适宜（85.29%）和4个边缘适宜（11.76%），无样点落入不适宜等级。样点m18以Ci = 0.76成为唯一高度适宜点位，其优势源于高CEC（0.268）、深井灌溉水（IWQI 0.212）和低ESP（0.075）的三重协同，2024年实测向日葵籽粒产量2.8 t ha?1，较盆地均值高出40%。中度适宜区呈现两个亚群：上段（m33、m11、m27、m30、m7、m32）以高CEC配合良好水质或石膏缓冲为特征；下段（m5、m6、m9、m15、m25）ESP尚可（<0.15）但因OM偏低或质地偏粗而受限。边缘适宜样点m12、m13、m22和m20的ESP均超过0.25，其中两个超过0.37，按Schoonover石膏需求量公式计算需施用6 t ha?1石膏方可将ESP降至15%基准以下。±20%权重扰动敏感性测试显示仅不到15%的样点跨等级边界。
Kruskal–Wallis检验。 非参数检验结果显示三项指标对Ci具有统计显著的判别力：CEC（p = 0.00056）、OM（p = 0.00313）和ESP（p = 0.0272）；而IWQI（p = 0.51）、坡度（p = 0.4135）、质地（p = 0.47）、石膏（p = 1）和pH（p = 0.2388）均未达到显著水平。箱线图显示：CEC从“极低”到“中等”等级，Ci中位数从约0.45升至约0.60；OM“低”等级中位数约0.58，“极低”等级降至0.50以下；ESP从“低钠化”（<5%）的约0.65逐级下降至“重度钠化”（>15%）的0.45以下，且高ESP等级区间宽度增大，反映局部排水和石膏缓冲的差异。
综合统计验证。 Pearson相关分析确认CEC（r = +0.65）、OM（r = +0.62）与Ci呈强正相关，ESP（r = ?0.70）呈强负相关；其余变量相关系数绝对值均低于0.35。VIF诊断显示所有变量均远低于5的阈值（最高为ESP的2.21和CEC的1.95），确认各准则贡献独立信息。PCA双标图中，第一主成分（Dim1，占方差32.6%）清晰对置了高CEC/OM（正载荷）与高ESP（负载荷），构成主导的“肥力–钠化”权衡轴；第二主成分（Dim2，占16.9%）反映地形（坡度）、矿物学（石膏）和化学活性（pH）的次要效应。样点在双标图中形成两个明确聚类：高适宜性组（Group 1）与高CEC/OM、低ESP相关联，边缘组（Group 2）以钠化主导、肥力偏低为特征。
空间分布与分区管理。 适宜性图揭示三个空间格局：第一，中部钙镁质坡积裙呈高至中上适宜性带状分布，得益于深厚坡积覆盖层（≥1 m）、高CEC和缓坡（3%），样点m18锚定该带；第二，占研究区85%的中度适宜区呈现功能均衡状态，CEC介于0.15–0.24 cmol(+) kg?1，IWQI处于低限制范围，ESP一般低于0.25阈值；第三，四个边缘适宜样点集中于wadi轴线附近的洼地，表现出极端钠化（ESP最高达0.34）、潜育层和3–5 m厚黏土–粉砂隔水层，电阻率层析成像证实其阻碍垂直淋溶。
五、讨论与结论
讨论部分指出，土壤相关限制因素主导了研究区的适宜性格局，而非水质因素。CEC、OM和ESP被确认为统计上主导且农艺上意义明确的Ci驱动因子；IWQI、坡度、质地、石膏和pH虽在局部或操作层面仍有参考价值（如在高ESP区域定向施用石膏、在坡度受限点位调整灌溉量），但并未主导整体适宜性景观。研究同时承认若干局限：输入数据的空间分辨率和可得性可能影响适宜性输出的精度；AHP赋权过程固有的主观性；缺乏长期时间序列数据限制了对土壤和地下水季节性及年际变化的捕捉；有限的田间验证限制了模型输出的直接核验。
研究结论可归纳为：在Maader灌区，85%以上的土地至少达到中度适宜等级，为向日葵扩张提供了坚实的农艺基础，前提是配套定向石膏改良和排水升级。钠害管理——而非内在肥力——是决定该区域适宜性的关键约束。中部钙镁质地带具备“即扩即用”的集约化潜力，而wadi邻近的钠化洼地需以≥6 t ha?1石膏配合地下排水作为扩张前提。该GIS–TOPSIS框架不仅提供了诊断工具，更构建了一个透明、可扩展、可直接投资的决策平台，服务于阿尔及利亚油料自给自足议程，并为北非及更广泛旱地区域的土地适宜性规划提供了可移植的方法论模板。研究同时支撑了联合国可持续发展目标（SDGs），通过促进明智的作物选择和可持续土地管理，助力提高农业生产力（SDG 2）、优化灌溉用水（SDG 6）及缓解土壤退化和盐渍化（SDG 15）。
