《Frontiers in Sustainable Food Systems》：Modeling cascading drought effects in a national food supply network using minimum cost flow
编辑推荐：
干旱日益威胁粮食安全，然而评估往往忽视了局部农业损失如何通过复杂供应链产生级联传播。研究人员开发了一个空间显式的网络优化模型，采用最小成本流（Minimum Cost Flow, MCF）算法，模拟在一个由开放地理空间数据构建的包含超过830,000个节点的全
干旱日益威胁粮食安全，然而评估往往忽视了局部农业损失如何通过复杂供应链产生级联传播。研究人员开发了一个空间显式的网络优化模型，采用最小成本流（Minimum Cost Flow, MCF）算法，模拟在一个由开放地理空间数据构建的包含超过830,000个节点的全国尺度瑞典食品供应网络中的食品分配。为测试系统对极端水文气候胁迫的响应，研究人员施加了一个与2018年经历条件类似的严重干旱情景。模型揭示了显著的空间错位：农业产量损失集中在南部农业区，而未满足的消费需求则主要出现在北部城市。此外，网络分析表明，干旱胁迫下的流量重新配置产生了新的关键瓶颈，实质上改变了系统的整体脆弱性。尽管这些具体地理结果仅用于说明该框架的诊断能力，而非经验验证的实际运行故障，但它们证明了干旱影响通过供应网络呈非线性传播。这强调了在粮食安全评估中考虑网络拓扑结构的关键重要性。最终，该方法为识别国家食品供应基础设施中的隐藏脆弱性并制定韧性策略提供了一个可扩展的框架。
干旱等极端水文气候事件对农业产量的影响可能通过高度互联的食品供应链传播，使远离初始冲击区域的消费者面临食品短缺。然而，现有干旱风险评估大多聚焦于直接农业损失，忽略损失在物流网络中的级联放大效应；而拓扑网络研究又常以随机节点删除模拟扰动，难以刻画产量部分减少导致的持续供应削减。为填补这一空白，研究人员开发了一种数据驱动的网络优化方法，以瑞典为例，构建了一个包含832,442个节点、约2280万条潜在边的全国尺度合成食品供应网络，并应用最小成本流（Minimum Cost Flow, MCF）算法模拟正常与干旱胁迫下的食品最优分配。研究采用2018年欧洲干旱情景，利用标准化降水蒸散指数（Standardized Precipitation Evapotranspiration Index, SPEI-12）转换为市政级产量减少系数，得到全国平均产量损失13.86%，并引入虚拟节点处理需求超过供给时的不可行问题。结果表明，尽管产量损失集中在南部和中部的生产区域，未满足的消费需求却主要出现在北部偏远城市，全国总未满足需求达14.59%，形成明显的“空间错位”。同时，流量重新配置使部分原本次要的节点（如H?ssleholm）中心性显著上升，成为新的关键瓶颈；而原有枢纽（如Arvika和Boden）则被绕过。基于产量损失与未满足需求的四象限分类进一步揭示了“传播型”（低产量损失、高未满足需求）和“缓冲型”（高产量损失、低未满足需求）两种典型区域模式，表明节点在网络中的位置而非灾害邻近性决定其脆弱性。该研究通过空间显式的大规模网络优化，展示了干旱影响非线性传播的机制，为识别国家食品供应基础设施中的隐藏脆弱性、制定针对性韧性策略提供了可扩展的分析框架。
研究人员主要使用了以下关键技术方法：一是基于开放地理空间数据构建网络，数据来源包括开放街道地图（OpenStreetMap, OSM）、谷歌地图位置服务（Google Places API）、瑞典统计局（Statistics Sweden）、瑞典农业委员会（Jordbruksverket）、联合国粮食及农业组织统计数据库（FAOSTAT）以及由SPEI-12推导的产量减少系数。网络被划分为生产、加工、分销、零售和消费五个层级，以欧氏距离作为边成本系数，并采用k近邻（k-NN）规则连接生产-加工（k=20）和零售-消费（k=15）节点。二是采用最小成本流（MCF）优化模型，基于流守恒约束最小化总运输成本，并使用Google OR-Tools求解。三是在干旱情景分析中，将SPEI-12数据通过线性回归模型转换为市政级产量减少系数，对生产节点的供应能力进行削减。四是通过虚拟节点法解决供给不足时的模型不可行问题，虚拟节点与所有需求节点相连并赋予极高惩罚成本，其流量被解释为该节点的未满足需求。五是采用网络中心性分析与统计检验，包括流加权介数中心性（Flow-Weighted Betweenness Centrality）、Wilcoxon符号秩检验、Kruskal-Wallis H检验和Dunn事后检验，以量化节点关键性的显著变化。
### 3.1 供应网络的拓扑结构
研究人员构建的基线瑞典食品供应网络呈现高度集中化的拓扑特征，具有重尾度分布和显著的空间集聚。生产与加工节点集中于南部农业区（如斯科讷和西约塔兰），消费节点则沿人口分布向南部落差递减。分布网络从南部纵向延伸至全国，少量大型分销中心承担了大部分流量，符合无标度网络特征。
### 3.2 压力下的优化：空间错位
在2018干旱情景下，MCF优化产生了14.59%的总未满足需求。产量损失集中在南部和中部的Uppsala、V?stmanland等地，而最严重的消费者短缺（>95%未满足需求）出现在北部的V?sterbotten和Norrbotten等市政区。这一空间错位表明，成本最小化物流优先满足人口密集的南部需求，导致短缺被推向北部偏远节点。
### 3.3 网络重构与关键性转移
流量差异分析显示，干旱影响区域内向外流量显著减少，同时出现新的关键物流走廊。介数中心性分析表明，H?ssleholm等分销节点中心性增加超过300%，成为压力网络中的新瓶颈；而Arvika和Boden等原枢纽中心性降至接近零，被新的优化路由绕过。统计检验确认分销和零售节点的中心性变化显著（p < 0.0001），说明网络重构并非随机波动。
### 3.4 传播类型学
基于本地产量损失与未满足需求的组合，市政区被划分为四种类型：韧性型（低损失、低短缺）、缓冲型（高损失、低短缺）、传播型（低损失、高短缺）和脆弱型（高损失、高短缺）。北部大量市政区属于传播型，尽管本地产量损失很小，却因网络位置偏远而遭受极端短缺；南部部分市政区尽管产量损失高达49%，仍因连接到未受影响产区而保持低短缺水平，表现出缓冲行为。
### 讨论与结论
讨论部分指出，成本最优化的基础设施设计存在效率-脆弱性权衡：MCF算法为满足南方人口密集地区需求而优化物流，却造成北部边缘地区严重短缺，说明仅考虑总产量的粮食安全评估具有局限性。流量重路由还产生了“适应不良”效应——旨在绕过干旱区域的优化路径将大量商品集中到单个节点，形成新的关键瓶颈，使网络面临二次级联失效的风险。该框架可作为压力测试工具，帮助决策者识别基线条件下不显著、但危机时关键的节点，并制定基础设施加固、缓冲库存布局和运输优先级等策略。研究的主要局限性在于网络基于开放地理空间数据而非实际物流数据，且采用欧氏距离而非真实道路距离，模型为静态均衡，并将食品视为单一同质商品；但敏感性分析表明空间错位是网络拓扑的稳健属性。结论部分指出：第一，干旱影响源（农业产量损失）与级联影响位置（消费短缺）之间可能存在显著空间错位；第二，供应冲击导致供应路线重构，次要节点可能成为意外瓶颈，增加次级失效风险；第三，大规模图优化为关键基础设施压力测试提供了数据驱动方法，可作为传统经济模型的空间替代方案。研究强调，国家粮食安全评估应同时考虑生产损失和供应网络在压力下的表现。
