气象干旱触发水文（农业）干旱的概率及其触发阈值的动态变化-贝叶斯copula模型

---

前言

在干旱研究中，一个关键的科学问题是：在某一地区发生不同等级的气象干旱时，气象干旱会以何种概率引发不同等级的水文干旱、农业干旱和地下水干旱？换句话说，气象干旱的不同程度会分别引发其他类型干旱的哪种等级？尤其是在大区域尺度上，并且涉及多种干旱类型时，干旱之间的传播过程和其动态变化对于触发不同类型干旱的概率及阈值至关重要。

为了回答这一问题，我们首先评估了在不同等级气象干旱胁迫下，水文干旱的诱发风险；接着，分析了水文干旱胁迫下对农业干旱的诱发风险；最后，探讨了农业干旱胁迫下地下水干旱的诱发概率。通过这一系列的评估，我们逐级构建了触发不同等级水文、农业和地下水干旱的阈值模型，并分析了这些模型随干旱发生过程中的动态变化规律。

这些研究结果为建立基于气象、水文、农业和地下水干旱之间相互作用的预警系统提供了科学依据。该系统对于减轻人为用水与自然生态系统水资源压力，并帮助制定适应性管理策略，具有重要的实践意义。

本研究复现了韩知明博士在其博士论文《中国多类型干旱时刻演变特征及其传播过程研究》中的相关工作，旨在为干旱风险预测和管理提供理论支持。

一、评估条件概率

为了评估在不同等级气象干旱胁迫下诱发不同等级气象、水文、农业和地下水干旱的风险，我们采用了贝叶斯概率模型。这一模型通过结合条件依赖的概念，推导了不同随机变量之间的关系。特别是，通过分析水文干旱相对于不同程度气象干旱的条件概率，我们可以估计气象干旱对水文干旱（SRI变化）的影响。

在多种干旱情景下，基于干旱等级的划分标准，我们进一步推导了不同等级气象干旱诱发不同等级水文干旱的概率计算公式。总体而言，这些情景涵盖了16种干旱组合情形，本文仅展示其中的两种情景计算。

二、求解触发阈值

对于给定等级的水文干旱，随着气象干旱等级的加重，理论上诱发水文干旱的概率应趋近于1。根据干旱等级划分，当SRI（气象干旱指标）小于等于-0.5时，认为发生水文干旱。因此，我们通过迭代SPEI（气象干旱指数）从-0.5开始，每次间隔-0.1，直到-3，来估算每次迭代所对应的条件概率。具体做法是，当条件概率大于或等于0.5时，我们返回对应的SPEI区间，区间的右侧值即为触发水文干旱的阈值。

当触发阈值较低（即SPEI值较小）时，说明只有较高等级的气象干旱才能引起SRI小于等于-0.5的情况，意味着水文系统对气象干旱的抵抗能力较强。相反，当触发阈值较高（即SPEI值较大）时，说明较低等级的气象干旱就可能引发SRI小于等于-0.5，表明该地区水文系统对气象干旱的抵抗能力较弱，值得相关部门关注并采取相应的应对措施。公式如下

气象干旱到水文干旱的触发阈值模型

1.代码片段

%Matlab version>=2021
close; % 关闭所有图形窗口
clear; % 清除工作区的变量
tic; % 计时开始

% 该代码的目的是分析不同等级气象干旱引发其他类型干旱的概率及阈值。
% 输入数据为一个包含气象干旱和水文干旱指数（例如SPI和SSMI）的CSV文件。

%% 数据读取
data = importdata('copula-data.csv');  % 导入数据
data = data.data;  % 获取数据矩阵
SPI = data(:, 2);  % 获取SPI数据列
SSMI = data(:, 3); % 获取SSMI数据列

%% 计算边缘分布并拟合
% 计算气象干旱（SPI）和水文干旱（SSMI）的边缘分布
[D_U1, PD_U1] = marginfitdist(SSMI);  % SSMI的边缘分布拟合
[D_U2, PD_U2] = marginfitdist(SPI);   % SPI的边缘分布拟合

% 计算SSMI和SPI的累积分布函数值
EP1 = cdf(PD_U1{
   
   1}, SSMI);
EP2 = cdf(PD_U2{
   
   1}, SPI);

% 组合SSMI和SPI的累积分布函数，避免负值
EP = [EP1, EP2];
EP(EP <= 0) = 1e-4;

% 选择最优的Copula模型（通过AIC准则）
[Family1, thetahat1, loglik1] = doublecopulaselect(