索马里2001–2025年水文气候变率与谷物生产对干旱的阈值响应

《Cogent Food & Agriculture》：Hydroclimatic variability and threshold responses of cereal production to drought in Somalia, 2001–2025

　　

摘要专业翻译 索马里的干旱通常被解读为侵蚀粮食生产的干燥化与变暖趋势。本研究探讨的决定性变化是否在于方差而非均值，以及谷物生产如何对此作出响应。将73个县的月度综合干旱指数（Combined Drought Index, CDI）与2001–2025年45个作物种植县的玉米和高粱（两种主粮谷物）的季节性生产数据相连接，随后采用趋势与方差检验、面板固定效应及阈值回归进行分析。平均干旱严重程度并未恶化，指数中的地表温度（land surface temperature）分量亦无趋势。变率发生了决定性变化。该指数在记录前后两半期之间的离散度大约翻倍，且变化源于降水而非温度。谷物生产呈非线性响应。在生长季指数低于约0.88时，损失累积速度是高于该值时的六倍；从有利条件转入干旱条件可使县级产出减少54%–61%。近半数生长季低于该阈值。这一灾害并非干燥化的均值，而是急剧的阈值敏感性遭遇日益加宽的降水摆动。早期预警应继续以降水为中心，但应从关注平均值转向关注变率与阈值跨越。

论文解读文章

一、研究背景与问题提出

索马里地处非洲之角干旱东翼，依赖两个短暂且不可靠的雨季，数十年的冲突削弱了民众在雨季失收时可依赖的制度保障。2011年饥荒、2017年近乎饥荒以及2020至2023年四十年最严重干旱，均被解读为降雨失败并造成严重人道代价。谷物承载了大部分此类风险：索马里农民种植玉米、高粱、水稻和小麦，但仅玉米和高粱构成国内谷物生产的实质份额，全国产出实际为两者之和。高粱为雨养作物，集中于中部干旱区，尤以拜州为主；玉米则在河岸区灌溉种植，谢贝利河是最重要的产区。谷物合计供应索马里膳食能量的三分之一至二分之一，玉米和高粱产出的变动直接关联粮食可及性。该产出长期未见改善：自1960年代以来无上升趋势，且年际波动极端剧烈，人均谷物产量从1972年的91公斤降至2012年的30公斤。本研究所考察的记录呈现同样模式：45个作物种植县的总产出从2001年的263,670吨降至2025年的145,901吨，序列以季节性摆动为主而非趋势主导，介于2002年峰值390,111吨与2023年谷值82,410吨之间。

针对这一记录的科学与政策回应长期聚焦于降雨是否呈下降趋势。这一追问因东非气候悖论而更加尖锐——观测到的长雨季下降与模型预测的更湿润未来并存。近期另一种框架兴起：大气变暖通过蒸发需求上升成为干旱的独立驱动因素。变暖提高饱和水汽压差和潜在蒸散量，使同等降雨总量维持的土壤水分少于以往；热量亦直接作用于谷物，因为一旦温度在开花和灌浆期超过作物特定极限，产量即急剧下降。这对索马里的实际含义是：若热力强迫现已构成干旱条件的实质份额，监测和早期预警需在降雨之外同时追踪温度。该热力信号是否存在于索马里使用的业务化干旱信息中，因此是需要检验而非假定的事项。

两种框架均关注均值条件的变化，均未涉及均值周围的分布本身是否在变化。变暖气候不仅预期改变平均条件，还预期使其失稳，加宽并加速干湿极端之间的摆动，此行为被描述为水文气候波动性或降水鞭打效应，预计在北部和东部非洲增强。索马里近期的序列与此行为一致：从2019和2020年洪水到2020至2023年干旱，再回到2023年破坏性洪水。这是否反映可测量的变率上升而非少数极端年份留下的印象，尚未在县级尺度上得到证实。趋势框架同样使影响侧基本未受检验。干旱之所以重要在于其对粮食的作用，而该联系不大可能是线性的。作物科学已表明产量对水分和热量的响应存在急剧非线性，损害集中在临界点之外而非平滑累积。若索马里谷物生产以此方式运行，则给定干旱变化的后果取决于生长季在该曲线上所处位置。使季节跨越临界点的变率，将比均值的缓慢移动更为重要。

由此产生三个研究空白：第一，索马里干旱研究以均值趋势检验为主，水文气候方差尚未在县级尺度上、以足以比较前后两半期的记录长度直接检验；第二，降水与温度对任何此类变化的相对贡献尚未针对索马里干旱监测所依赖的业务化指数加以量化；第三，将气候与谷物产出相联系的索马里证据局限于全国年度总量和线性设定。因此，不存在对索马里谷物生产开始失效的干旱水平的估计。本研究将两侧结合起来：使用FAO-SWALIM综合干旱指数及其分量覆盖73个县，以及FSNAU季节性谷物生产数据覆盖45个作物种植县，均涵盖2001至2025年。研究回答三个问题：平均干旱严重程度是否如趋势框架所假设的那样加剧，还是主导变化是水文气候变率的上升；县级谷物生产对干旱的响应是否非线性，是否包含产出急剧下降的阈值；上升的变率和阈值暴露如何在县与地区间分布，该分布对干旱监测和预见性行动有何含义。

二、主要关键技术方法

研究人员构建了2001至2025年索马里73个县的月度综合干旱指数（CDI）及其三个分量（降水干旱指数PDI、温度干旱指数TDI、植被干旱指数VDI）的面板数据，以及45个作物种植县的FSNAU季节性玉米和高粱生产数据。关键方法包括：采用Mann–Kendall检验和Sen斜率评估均值趋势，并应用Hamed和Rao方差修正及Benjamini–Hochberg错误发现率控制；采用Brown–Forsythe中位数中心形式的Levene检验比较前后两半期（2001–2013与2014–2025）的方差变化，辅以24个月滚动标准差和极端态占比分析；利用CDI为固定加权和的性质对方差变化进行精确分解；采用双向固定效应面板模型（含二次项）估计生产对干旱的非线性响应；采用Hansen回归扭结框架的阈值回归估计未知阈值，以集中最小二乘法搜索、Davies检验和通配聚类自助法（999次重复）进行推断；采用对子聚类自助法构建置信区间；以Driscoll–Kraay标准误和县—季节双向聚类标准误进行稳健性检验。数据来源为FAO-SWALIM和FSNAU的公开业务数据。

三、研究结果

干旱严重程度的趋势：2001至2025年间，平均干旱严重程度无加剧证据。没有任何县在用于构建指数的地表温度上呈现显著趋势，全国序列亦无趋势（p=0.69）。指数的温度分量无全国趋势（p=0.50）。降雨无全国趋势（p=0.20），尽管20个县显示显著增加。仅有的全国信号是降水分量（+0.23/十年，p=0.042）和CDI（+0.18/十年，p=0.027）朝向更湿润的弱上升趋势；排除2018和2023两个极端湿润年份后，CDI趋势不再显著（p=0.10）。这些结果仅涉及均值条件。

水文气候变率的变化：均值条件未加剧的同时，均值周围条件的离散度显著上升。全国CDI异常的标准差从前半期的0.315升至后半期的0.588，方差增加约3.5倍（Brown–Forsythe检验，p<0.001）。73个县中72个方差上升，经多重检验校正后54个显著。24个月滚动标准差呈显著上升趋势（p=0.009），12和36个月窗口结果相同。全国年降雨变异系数从0.21升至0.40。变宽在两端之间不均衡：严重干旱的县月占比从18.0%降至12.5%，而非常湿润占比从9.0%升至22.6%，变率上升主要是湿极端的扩张。该上升对分割年份的选择稳健（2008至2017年任意分割年均显著，方差比2.9–4.3），删除任意单年仍显著，但集中于记录末端的异常湿润年份——排除2019和2020年将全国方差比降至1.6，排除2018至2020及2023年则消除全国上升。此变率上升是2001至2025年记录的明确特征，但主要由2018年以来的一簇湿润极端驱动，此长度记录无法确定其是否持续。

指数变率变化的分解：由于CDI是其降水与温度分量的固定加权和，方差变化可精确分解。全国CDI异常方差增加中，95.1%源于降水分量方差增加，1.2%源于温度分量，3.8%源于二者协方差。降水项在方差上升的全部72个县中均为最大，中位份额98%。分量直接比较显示：降水分量异常标准差上升1.99倍，温度分量1.52倍，未进入指数的植被分量上升2.06倍。降水与温度异常间的相关性从0.66降至0.34。CDI变率上升主要与降水分量变率上升相关；温度变率亦上升，但其较低权重和较小方差使其对指数变化的贡献甚微。

研究期内的谷物生产：45个作物种植县的玉米高粱合计产量从2001年的263,670吨降至2025年的145,901吨。2001至2013年平均年产量为259,477吨，2014至2025年为166,747吨，下降36%。序列介于2002年峰值390,111吨与2023年谷值82,410吨之间。高粱占期间产量的53.5%，玉米占46.5%。此下降反映干旱以外的诸多因素，分析并未将其归因于干旱，而是估计在去除县和年效应后产量如何随干旱条件变化。

谷物生产对干旱的响应：对数谷物产量对生长季CDI呈凹形响应。线性项为正（+0.96），二次项为负（?0.36），均在p<0.001显著。在CDI为0.5时，指数增加0.1与产量提高14.3%相关；CDI为1.0时相关10.3%；1.5时相关6.4%。反向解读，产量随条件从正常移向干旱而加速下降。拟合曲线在CDI为2.36处达峰，但仅2.5%的生长季超出该值。玉米和高粱分别、以及两种生计系统中均呈现相同凹形模式。河岸系统响应较弱（+0.48和?0.19），因灌溉降低了对降雨的依赖；农牧系统响应较强（+1.03和?0.34），因其依赖降雨。季节最大CDI与县级谷物产量之间无统计可检测关联（p=0.34），季节内CDI标准差与产量之间亦无（p=0.87），但此检验功效有限，因季节最大值与季节均值高度相关（r=0.98）。

干旱—生产阈值：阈值回归将生产响应陡增点定位于生长季CDI为0.88（95%置信区间0.82–1.23）。阈值以下，对数产量对CDI的斜率为+1.94（1.28–2.38）；阈值以上为+0.33（?0.03至0.48）。阈值以下，CDI下降0.1与产量降低17.7%相关（12.0–21.2）；阈值以上，相同下降与3.2%降低相关（?0.3至4.6）。阈值以下斜率约为阈值以上的六倍，二者差异为1.61（1.16–2.07）。阈值模型显著优于线性模型（F=48.5，通配聚类自助法p=0.001），无第二阈值证据（p=0.37）。从正常条件（CDI=1.0）移至严重干旱边界（CDI=0.5）与谷物产量减少54%相关（47–60）；从有利范围湿端（CDI=1.5）移至同一边界与61%减少相关（56–68）。阈值接近观测到的生长季条件分布中心：2001至2013年44.3%的生长季低于阈值，2014至2025年为47.7%。估计对标准误计算方法稳健，5%或15%修剪下不变，排除2025年Deyr季阈值移至0.83（在置信区间内）。

变率的区域格局：变率上升在县间不均衡。中位县后半期方差为前半期的2.4倍。上升最大在西北部：Awdal四个县的方差上升7.9至9.4倍，为全国最大增幅，地区层面Woqooyi Galbeed为4.7倍、Togdheer为4.0倍。中部地区中，Galgaduud（3.9倍）和Middle Shabelle（3.4倍）亦远高于中位数。增幅最小在东北部的Bari（1.5倍）和Nugaal（1.7倍），以及南部产粮区Bay和Lower Juba（各1.8倍）。这些为描述性区域格局，未检验高变率县是否形成统计显著的空间聚类。

四、讨论总结与研究结论

讨论部分指出，这些结果重新定位了索马里干旱问题。过去四分之一世纪中，两种主流框架在县级数据中鲜获支持：平均干旱严重程度未恶化，进入指数的地表温度无县呈显著趋势，温度在指数中的权重未变。变化在于均值周围条件的离散度，且此变化几乎完全与指数的降水分量相关而非热力分量。变化亦为近期且集中：主要依托2018年以来一连串异常湿润年份，25年记录无法确定其为分布的持久转变还是极端簇。使此变化具有后果的是生产响应的形状：谷物产出在生长季指数约0.88以下陡峭响应、以上则平缓得多，且该点已接近观测分布中心——前后两半期均有近半数生长季低于此值。指数的给定变动因此取决于季节所处位置而意义迥异，更宽的摆动使更多季节跨越损失加速点。这是不同于干燥化趋势的问题，因其关乎方差和非线性而非水平。研究发现均值条件未恶化与记录严重近期干旱的文献并存而非对立：结果说明跨73县25年业务化干旱指数的均值无下降趋势，而非个别干旱不严重或未因变暖而加剧。变率上升与北非和东非水文气候波动性和降水鞭打效应增强的预测一致。与近期区域研究（如肯尼亚1981至2021年降雨变率变化、萨赫勒变率上升预测）方向一致。生产侧证据表明变率比均值更重要：亚季节性干湿格局变化可解释撒哈拉以南非洲玉米产量年际变异的50%–60%，约为均值降水与温度解释份额的两倍，非洲之角约达75%。本研究结果为索马里谷物的县级对应物。与东非面板阈值研究相比，本研究的增量在于提供了以业务化指数单位表达的阈值，可直接与监测触发值比较。对实践的三个含义：监测应继续以降雨为中心；阈值是正确的工具，且现已有可与业务值比较的经验估计——预见性行动计划在CDI为0.8或以下时启动，而损失加速点估计为0.88（置信区间0.82–1.23），业务触发值恰位于区间之下，提示以0.8为关键的行动可能在生产已进入陡峭响应段后才开始；由于灾害是阈值跨越而非均值水平，以典型条件校准的监测产品和应急融资与之匹配不佳，一个季节的相关统计量不是其相对于正常的干燥程度，而是其是否低于损失加速点。谷物生产是粮食安全的若干决定因素之一，结果涉及可及性和稳定性而非可获取性或利用。阈值接近生长季条件分布中心意味着县级产出在季节间于非常不同的水平间移动，更宽的指数摆动使这些移动更不规则。变率上升最大的县并非种植大部分谷物的县，但这不削弱论点，因作物区方差仍约翻倍且近半数生长季已低于阈值。湿极端方面，季节最大指数与县级谷物生产之间无统计可检测关联，不应解读为洪水不损害作物——季节最大值与季节均值几乎完全相关，湿极端的独立效应无法在这些数据中识别。最后，所估计的关系是关联而非因果证明：面板为观测性，县级产出响应降雨以外的诸多因素，县固定效应吸收县间永久差异、年固定效应吸收同年共同冲击，但不消除县内随时间变化且与干旱条件相关的因素。

研究结论部分总结：第一，平均干旱严重程度在73个县中无加剧证据，指数地表温度分量无显著趋势；变率则显著上升，全国指数异常方差在两半期之间上升约3.5倍，73个县中72个上升，该上升几乎完全与降水分量相关，主要由2018年以来异常湿润年份驱动。第二，县级谷物生产对干旱的响应是非线性的：生长季指数0.88以下，指数下降0.1与产量降低17.7%相关，以上为3.2%；从正常条件移至严重干旱边界与约54%产出损失相关；近半数生长季已低于阈值。第三，这些发现指向围绕降雨变率和阈值跨越而非平均条件组织的监测与预见性行动。两条超越证据的建议：检验业务触发值0.8与估计阈值0.88之间的关系，以确认行动是否被关键到损失已加速的点；在记录产量之外同时记录收获面积，以分离响应的面积和单产分量。