基于Copula的洞庭湖-流域-长江系统水文干旱概率分析*

李 珍，李相虎，张 丹，蔺亚玲

(1：中国科学院南京地理与湖泊研究所，中国科学院流域地理学重点实验室，南京 210008)(2：中国科学院大学，北京 100049)

干旱是全球发生最普遍且损失最严重的自然灾害之一，具有发生频率高、影响范围广、持续时间长等特点[1-3].在全球气候变暖背景下，干旱事件的发生频率和强度都在呈显著增加趋势，而且根据预测，未来将会有更多陆地区域出现极端干旱[4-5].我国地处东亚季风区，气候条件复杂，降水的时空变率大，特殊的气候变化特征与地理环境导致大范围干旱灾害频繁发生，对当地粮食安全和人民生产生活造成严重危害[6-7]，已成为制约经济以及生态可持续发展的重要因素之一[8].研究发现，近年来我国北方干旱频发的同时，南方湿润区干旱事件也在不断增加[9-10]，尤其是季节性干旱[11-12]和骤发性干旱[13]突出.例如西南地区2009年秋季至2010年春季遭遇百年一遇的特大旱灾，致使2500万人生活用水困难，大量农作物歉收或绝收[14]；2011年长江中下游地区春旱造成江河及湖泊水位异常降低，湖区生态环境恶化，威胁水生动植物的生存和繁衍，造成了严重的经济损失[15].因此，对南方湿润地区季节性干旱的研究也已引起国内外学者的广泛关注，这对深入认识气候变化和人类活动影响下的干旱变化规律、影响因素以及区域防旱抗旱减灾等都具有重要的科学意义.

根据研究对象或应用领域的不同，干旱可分为气象干旱、水文干旱、农业干旱以及社会经济干旱4种类型[16].干旱指数是识别区域干旱最常用的指标[17]，它不仅可以量化干旱的强度、持续时间等特征，而且在干旱评价以及监测预警方面也具有重要作用.常用的干旱指标包括Palmer干旱指数(PDSI)、标准化降水指数(SPI)、标准化降水蒸散指数(SPEI)以及标准化径流指数(SRI)等[18-20].国内外诸多学者也从不同角度基于不同的干旱指标进行研究，例如莫兴国等[21]采用PDSI评估了未来气候变化情景下我国干旱事件发生频次、强度和持续时间的可能变化；李明等[22]计算了不同时间尺度的SPI，结合游程理论识别中国东部季风区气象干旱事件.由于变化环境下气象水文时间序列受到气候变化和人类活动的影响呈现非平稳特征[23]，而传统干旱指数的计算很大程度上依赖平稳性假设[24]，因此一些学者尝试构建非平稳干旱指数来识别干旱事件.例如温庆志等[25]在SPEI和非平稳性理论基础上构建了非平稳性标准化降水蒸散指数，用来评估不同排放情景下中国气象干旱的时空格局；李敏等[26]通过GAMLSS模型构建时变SRI以评价流域变化环境下的水文干旱.此外，干旱具有多个特征变量，基于单一变量的分析无法完全解释干旱各特征变量之间的关系，因此，越来越多的研究更注重基于多变量的联合分析.其中，Copula函数能够通过边缘分布和相关性结构来构造多维联合分布[27]，既能进行特征的联合，又可以进行区域联合，特别适用于多变量水文分析，已在极端降水、洪水频率分析和干旱特征分析等方面广泛应用[28-30].如陈永勤等[31]基于游程理论识别了鄱阳湖流域水文干旱特征变量，采用Copula函数分析了主要支流的干旱历时和干旱烈度的联合概率特征；Xu等[32]在考虑干旱事件时空变化基础上，构建了干旱持续时间、影响面积以及严重程度的三变量联合分布.Zhang等[33]利用河流径流和湖泊水位构建水文干旱指数对鄱阳湖、鄱阳湖流域和长江的联合水文干旱进行了研究.Copula函数已被证明是多元水文分析和模拟的有效工具.

洞庭湖是我国第二大淡水湖，也是三峡水库坝下第一个大型通江湖泊[34]，素享“鱼米之乡”的美誉.湖区接纳流域“四水”(湘江、资水、沅江、澧水)和长江“三口”(松滋口、太平口、藕池口)的径流，经湖泊调蓄后于城陵矶汇入长江，形成了复杂的江-湖-河水系格局[35].受气候变化与人类活动的影响，洞庭湖在过去几十年间多次发生不同程度的旱灾[36].尤其是2000年以来在流域径流变化及长江上游大型水利工程的共同影响下，江湖关系发生变化，水量平衡关系改变，湖泊干旱事件更为频繁，洞庭湖湖容明显减小，枯水期提前[37]，季节性水文干旱不断加剧[38].目前有关洞庭湖干旱的研究多侧重于湖泊或流域干旱特征的描述以及三峡水库运行后江湖关系改变对洞庭湖水情变化的影响等方面，未将长江、洞庭湖及其流域作为一个系统考虑，也未揭示系统内洞庭湖干旱与流域、长江干流水文干旱的内在关联与相互作用关系.因此，本文基于传统和时变标准化水位指数(SWI)和标准化径流指数(SRI)，通过Copula函数分析洞庭湖-流域-长江系统1964-2016年水文干旱事件的联合概率分布特征，明确具有不同属性的水文干旱事件在洞庭湖-流域-长江系统的相关性，探讨洞庭湖水文干旱的成因及其对江湖关系变化的响应.本研究对进一步认识洞庭湖干旱的发生机理具有重要意义，同时也为洞庭湖防旱抗旱减灾以及科学实施干旱预警提供重要的科学依据.

1 研究区域与方法

1.1 研究区域与数据资料

洞庭湖(28°03′～30°20′N，111°40′～113°30′E)地处长江中游荆江河段以南，湖南省东北部，由西洞庭湖、南洞庭湖和东洞庭湖3个湖区组成，是国际重要的湿地保护区.洞庭湖流域水系发达，流域面积约26万km2，占长江流域总面积的12%[39].洞庭湖是我国典型的大型调蓄湖泊，长江径流由松滋、太平、藕池“三口”分流入湖，与流域“四水”在湖区汇合调蓄后又由唯一的通道城陵矶注入长江[40]，流域与长江汇水以及长江的顶托作用共同决定了湖泊的水位.洞庭湖流域与长江中游整体上属于亚热带湿润季风气候，四季分明，降水季节集中且年际变化大.因此长江与流域入湖径流量的年内分配也不均匀，流域与长江的水情变化直接影响洞庭湖的水量变化，季节性降水导致洞庭湖具有明显的季节性水文特征.通常流域汛期为4-6月，洞庭湖受流域补给在4月水位开始上升；长江汛期为7-9月，湖水位继续抬升至年内最高[41].流域与长江之间的汛期转换使洞庭湖成为一个复杂的江-湖-河系统.

本研究选取洞庭湖南咀、小河咀、城陵矶3个水文站1964-2016年水位数据，以3站平均水位反映洞庭湖的整体水位变化情况.以洞庭湖流域桃江、桃源、湘潭水文站径流量代表流域入湖水量变化情况(图1).选取长江干流宜昌站径流量代表长江来水变化情况.并且以2003年三峡大坝运行为时间断点，在总研究时段中重点关注2003年之后的水文干旱.上述数据均来源于长江水利委员会水文局.流域内同期降水数据来源于国家气象中心，主要用于水文干旱事件的成因分析.

图1 洞庭湖-流域-长江系统水系和水文站、气象站分布

1.2 研究方法

1.2.1 标准化水位/径流指数(SWI/SRI) 标准化水位指数SWI和标准化径流指数SRI是参照标准化降水指数SPI这一概念而提出，用于识别区域水文干旱和量化干旱等级[42]，两者的计算过程与SPI类似.考虑到长时间水位和径流序列可能存在非平稳特征，因此在干旱指数计算前首先进行MK趋势检验[43]，趋势变化显著的季节则在GAMLSS框架下计算时变标准化水位指数和径流指数[26,44-45].水文干旱根据SWI和SRI大小共划分为5个干旱等级[46](表1)，当SWI和SRI值低于-1时，即认为发生水文干旱，并且其数值越小，表明其干旱强度越大，本研究关注中度及以上程度的水文干旱(即SWI≤-1、SRI≤-1).同时，为了反映洞庭湖-流域-长江系统水文干旱变化的总体特征，分别从年和季节时间尺度分析SWI与SRI的变化情况，其中选用12个月尺度下的SWI与SRI值代表年际变化，选用3个月尺度下5、8、11以及次年2月的值分别代表春、夏、秋、冬季变化.

表1 水文干旱等级划分

1.2.2 边缘分布函数 建立联合分布函数之前首先需要确定每个水文干旱指标的边缘分布函数[30]，水文科学领域存在多种随机变量分布函数，本文选择在干旱分析中广泛应用的Weibull、Logis和Gamma分布函数来研究SWI和SRI概率分布特征，并采用最大似然函数和线性矩法进行参数估计.

1.2.3 Copula函数 Copula理论最初由Sklar[47]提出，目前已广泛应用于二元和多元干旱频率分析.Copula函数根据随机变量之间的依赖结构，连接一维边缘分布以形成概率区间在[0，1]上的多元联合分布[48].在研究中常用的Copula函数分为4种类型：阿基米德型(Frank、Clayton、Gumbel)，椭圆型(t、Gaussian)，极值型(Husler-Reiss、t-EV)和混合型(Plackett)[32].本文采用Gumbel、Clayton、Frank、Gaussion、t和Plackett 6种Copula函数进行拟合，同时根据相关性指标法，建立起Kendall秩相关系数τ与Copula函数的参数θ之间的关系来进行参数估计[49].

根据Sklar定理，假设两个相关的水文事件X和Y分别具有FX(x)、FY(y)分布，则它们的联合分布F(x,y)如下：

F(x,y)=P(X≤x,Y≤y)=Cθ(FX(x),FY(y))

(1)

式中，C表示变量X和Y的二元依赖结构的Copula函数，θ为待定参数.在本研究中，X和Y指SWI和SRI的时间序列.

Kendall秩相关系数[50]计算公式如下：

(2)

式中，τ为Kendall秩相关系数，(xi,xj)为对应时间段SWI和SRI的值，n为系列长度.

某一水文干旱事件的联合重现期[51]可由下式计算：

(3)

式中，s是序列中出现的时间间隔.当SWI和SRI时间尺度为12个月时，s=1年；当时间尺度为3个月时，s=0.25年.

在选择最优Copula函数之前，利用R语言软件包Copula对长江、洞庭湖及其流域水文干旱指数的可交换性进行检验[52]，以获得统计量的近似P值，用来评估未知极值CopulaC是否满足C(u，v)=C(v，u)，u、v∈[0，1].如果不拒绝可交换性检验，就可以排除不对称的Copula函数，从而减少Copula函数在拟合优选中的数量.

1.2.4 边缘分布函数和Copula函数的优选 边缘分布函数通过Kolmogorov-Smirnov(K-S)检验法进行假设检验，达到95%置信水平则采用AIC信息准则[53]评价分布函数的拟合优度.统计量D和AIC的值越小，表明函数拟合效果越好.Copula函数同样通过AIC信息准则进行拟合优度评价，最终从6种Copula函数中选择拟合最优的联合分布函数.各检验方法的表达式如下：

D=maxij|F(X(i))-G(Y(j))|

(4)

(5)

AIC=2k+nln(MSE)

(6)

式中，F(X)为经验分布函数，G(Y)为理论分布函数，MSE为均方误差；AIC为赤池信息准则，k为参数个数，n是样本个数.

2 结果分析

2.1 边缘分布函数和Copula函数的选择

2.1.1 边缘分布函数的选择 基于Weibull、Logis和Gamma 3种边缘分布函数，分别对洞庭湖SWI、流域SRI以及长江SRI的分布进行拟合，K-S检验结果如表2所示.3种分布函数对不同时间尺度上SWI和SRI的拟合均通过了K-S检验.根据统计量D与AIC最小值原则，选择拟合优度最高的分布函数作为不同区域各个季节水文干旱的边缘分布函数(表3).洞庭湖SWI在年尺度和夏季、秋季Logis函数对应的AIC值最小，遂选Logis函数为该时段的最优边缘分布函数.在春季和冬季，Weibull与Gamma函数对应的AIC值最小，则Weibull与Gamma函数分别为春季和冬季的最优边缘分布函数；而对于流域SRI，在冬季最优边缘分布函数为Gamma，其他3个季节和年尺度均选择Weibull函数为最优边缘分布函数；长江SRI在年尺度和夏季最优边缘分布函数为Logis，在春季、秋季为Weibull函数，而在冬季为Gamma函数.

表2 基于K-S检验的SWI和SRI边缘分布拟合优度值

表3 基于AIC信息准则的SWI和SRI边缘分布函数拟合优选*

2.1.2 Copula联合分布函数的选择 可交换性检验是基于经验Copula与潜在二元Copula可互换性的评估测试[54].表4显示各联合水文干旱指标在95%置信水平下均通过交换性检验，意味着以上6种Copula函数可用于洞庭湖-流域-长江系统水文干旱联合概率的研究.基于AIC最小值原则，从Gumbel、Clayton、Frank、Gaussion、t和Plackett Copula中选择拟合最优的Copula函数.如表5所示，洞庭湖-流域系统在年尺度和4个季节水文干旱联合概率拟合效果最优的Copula函数分别为Clayton、Frank、t、Gumbel和Gumbel；洞庭湖-长江系统则分别为t、t、Gaussion、Clayton和Clayton Copula.

表4 SWI和SRI可交换性检验

表5 基于AIC选择Copula函数*

2.2 各区域水文干旱概率分布特征

基于选定的不同区域最优边缘分布函数，分别分析洞庭湖、流域及长江的水文干旱概率分布特征(图2).在年尺度上，洞庭湖水文干旱的发生概率(SWI<-1)为14.01%，重现期为7.14 a.1964-2016年间，洞庭湖共发生9次水文干旱事件，其中4次出现在2003-2016年(表6)，最严重的水文干旱发生在2011年，SWI为-2.47(极端干旱).流域水文干旱的发生概率为16.20%，重现期为6.17 a，研究期内共识别出10次水文干旱事件，有2次发生在2003-2016年间，其中SRI最小值-2.55(极端干旱)也出现在2011年.对于长江干流，1964-2016年水文干旱的发生概率为14.35%，重现期6.97 a，过去53年间共发生7次水文干旱事件，其中2003年之后有3次，SRI最小值-3.39(极端干旱)发生在2006年.从3个区域对比来看，虽然整体上流域发生水文干旱的概率最高，但2003-2016年间发生水文干旱的频次较少，而洞庭湖及长江在2003年以后发生水文干旱的频次明显增多，表明在2003年以后洞庭湖及长江水文干旱呈加剧态势，且二者的加剧态势具有同步性.

图2 年尺度洞庭湖SWI概率分布(a)；流域SRI概率分布(b)；长江SRI概率分布(c)

表6 1964-2016年洞庭湖-流域-长江系统水文干旱概率对比

季节尺度上(图3)，洞庭湖春季水文干旱的发生概率为15.58%，研究时段内共发生11次水文干旱事件，其中发生在2003年以后的有2次.春季流域水文干旱的发生概率均大于洞庭湖和长江，研究时段内共识别出10次水文干旱事件，其中有4次发生在2003-2016年间(表6)，SRI最小值-2.60(极端干旱)出现在2011年.春季长江在2003年以后未识别出水文干旱事件，其SRI最小值-1.90(重度干旱)出现在1979年.在夏、秋季，洞庭湖水文干旱概率均大于流域和长江，尤其秋季洞庭湖水文干旱的发生概率高达 32.46%，重现期0.77 a，在研究时段内共识别出15次水文干旱事件，其中有7次发生在2003年以后，SWI最小值-3.30(极端干旱)发生在2006年.秋季长江中游2003年以后识别出5次水文干旱事件，占总次数的一半以上.冬季洞庭湖水文干旱概率最小，为4.60%，长江冬季水文干旱的发生概率为16.27%，2003年以后长江与洞庭湖冬季均未识别出水文干旱事件，流域仅识别出1次(表6).

图3 季节尺度(a～d)洞庭湖SWI概率分布(a～d)；流域SRI概率分布(e～h)；长江SRI概率分布(i～l)

从不同季节对比来看(表6)，总体上洞庭湖秋季水文干旱频率最高、频次最多，表明洞庭湖秋季水文干旱最为严重且2003年以后呈加剧的态势.其次是春季水文干旱次数较多，但2003年之后所占比重相对较小，表明2003年以后洞庭湖春季水文干旱有所减轻；流域春季水文干旱发生概率高于洞庭湖和长江，且2003年以后干旱频次所占比重较大，表明流域春季水文干旱与其它2个区域相比更为严重，并仍有持续加剧的趋势；对于长江干流，2003年以后春冬两季均未识别出中度以上水文干旱事件，但秋季水文干旱的发生次数占总次数的一半以上，表明长江秋季水文干旱进一步加剧.

2.3 洞庭湖-流域-长江系统水文干旱联合概率分布特征

根据优选的Copula函数，进一步研究洞庭湖-流域-长江系统水文干旱的联合概率分布特征.图4为年尺度上基于Clayton Copula的洞庭湖-流域系统和基于t Copula的洞庭湖-长江系统水文干旱的联合概率分布.在1964-2016年间，洞庭湖-流域系统的水文干旱联合概率为9.65%，重现期10.36 a，其中发生在2003年后的联合水文干旱仅2次，其对应的洞庭湖SWI和流域SRI值分别为-2.47和-2.55，均达到极端干旱程度.而洞庭湖-长江系统的水文干旱联合概率为8.58%，重现期为11.66 a，年尺度上共识别出4次联合水文干旱事件，其中有3次发生在2003-2016年间，水文干旱最严重的年份为2006年，其对应的洞庭湖SWI为-2.15，长江SRI为-3.39，均达到极端干旱程度.对比发现，总体上洞庭湖-流域系统水文干旱的联合概率大于洞庭湖-长江系统，表明年尺度上流域补给减少对洞庭湖水文干旱的贡献更大.但2003年后洞庭湖-长江系统水文干旱同时发生的频次明显多于洞庭湖-流域系统，且所占比重较大，表明2003年之后长江来水量的减少对洞庭湖水文干旱的影响增强.

图4 年尺度1964-2016年洞庭湖-流域系统水文干旱联合概率(a)；洞庭湖-长江系统水文干旱联合概率(b)

从季节分布来看，春、夏、秋和冬季洞庭湖-流域系统水文干旱的联合概率分别为9.52%、7.64%、7.98%和2.63%(图5和表7)，春季洞庭湖与流域同时发生水文干旱的概率均高于其它季节.此外，1964-2016年间，春季洞庭湖-流域系统同时发生水文干旱事件7次，其中有2次出现在2003-2016年.虽然洞庭湖与长江也同时发生水文干旱事件达7次，但全部出现在2003年以前，表明洞庭湖春季水文干旱与流域及长江径流变化均存在响应关系，但2003年后，流域径流对洞庭湖水文干旱的贡献增强.对于洞庭湖-长江系统，其秋季水文干旱联合概率为15.38%(图6)，是秋季洞庭湖-流域系统水文干旱联合概率的近2倍，远远高于其它季节.研究时段内共识别出7次水文干旱事件，其中有5次发生在2003-2016年间，占总频次的一半以上，表明2003年以后秋季洞庭湖-长江系统同时发生水文干旱事件的频次明显增多.总体而言，春季洞庭湖-流域系统水文干旱联合概率大于洞庭湖-长江系统，表明春季流域补给对洞庭湖水位具有主导作用.夏、秋和冬季洞庭湖-长江系统水文干旱联合概率高于洞庭湖-流域系统，尤其在2003年之后秋季洞庭湖-长江系统水文干旱呈加剧态势.

表7 1964-2016年洞庭湖-流域-长江系统水文干旱联合概率

图5 季节尺度1964-2016年洞庭湖-流域系统水文干旱联合概率(a～d分别表示春、夏、秋、冬季)

图6 季节尺度上洞庭湖-长江系统水文干旱联合概率(a～d分别表示春、夏、秋、冬季)

3 讨论

洞庭湖水文干旱受流域和长江来水的共同影响，但其主导因素随时间发生变化[55].本文研究发现，在年和春季时间尺度上，洞庭湖-流域系统的水文干旱联合概率高于洞庭湖-长江系统，尤其是春季洞庭湖-流域系统的水文干旱联合概率达9.52%，这表明春季洞庭湖水文干旱与流域水文干旱具有较好的同步性，流域入湖补给减少是春季洞庭湖水文干旱的主导因素.孙占东等[34]通过对比不同季节长江与流域来水变化与典型干旱事件的对应关系，也认为湖泊春季干旱基本上是由洞庭湖流域来水偏少造成的；贺秋华等[56]研究发现流域“四水”与洞庭湖年均入湖径流量变化曲线呈同步增减趋势，且二者的相关系数为0.91，而长江“三口”年均分流入湖与洞庭湖入湖水量无明显相关性，说明洞庭湖年均入湖水量主要受流域来水影响，这与本文的研究结果较为一致.洞庭湖流域来水变化的主要影响因素之一是降水的改变[57].由图7所示，洞庭湖流域年际和季节性降水分配不均匀，季节性水文干旱与降水减少相对应，其中春季的每一次水文干旱事件均与降水量骤减相吻合.比如，2011年春季流域降水量达到历史最低值，此时SRI也达到最低值-2.6(极端干旱).吉红霞等[15]利用遥感影像和水文气象观测资料对2011年洞庭湖极端干旱事件的成因进行了分析，发现2011年洞庭湖流域4-5月降水量比多年平均偏少50%以上，降水减少是洞庭湖春季干旱发生的主导因素.

图7 洞庭湖流域SRI与降水量季节分布(a～d分别表示春、夏、秋、冬季)

夏季、秋季和冬季的洞庭湖水文干旱主要受到长江来水变化的影响，尤其是秋季洞庭湖-长江系统同时发生水文干旱的概率达15.38%.特别是1964-2016年间秋季洞庭湖-长江系统同时发生的7次水文干旱事件有5次发生在2003年以后，这表明2003年后秋季洞庭湖-长江系统同时发生水文干旱事件的可能性增加，这主要由长江干流水利工程调度引起的长江来水变化导致的，以三峡水库蓄水最为典型[58-59].根据三峡水库优化调度方案，水库在9月15日后开始蓄水，此时长江中下游干流水位快速降低，对湖泊的顶托作用减弱[60].Gao等[61]研究了三峡水库运行前后长江中下游的流量变化，发现三峡水库秋季蓄水是长江中下游河道流量显著减少的主要原因.赖锡军等[62]研究发现，三峡水库秋季蓄水促使湖水加速流入长江，造成洞庭湖水位快速下降、低枯水位提前出现.值得注意的是，秋季三峡水库开始蓄水，同时该时期洞庭湖流域降水也相对较少，流域来水对湖泊的补给减弱，两者的叠加效应极易引发湖泊极端水文干旱事件[63-64]，如Ou等[65]研究表明，2006年秋季洞庭湖的极端干旱事件是由长江“三口”和流域“四水”来水锐减共同引起的，其中长江“三口”和流域“四水”来水量较同期多年平均值分别减少了74.7%和30.4%.

除降水变化等气候因素和三峡蓄水对河湖水文要素的季节变化产生了显著影响外，其它人类活动的影响也不容忽视.据湖南省水利发展统计公报显示[66]，到2016年末，全省已建成各类水库14098座，总库容514亿m3，其中包括大型水库45座和中型水库359座，各类蓄水工程年末蓄水量242.5亿m3，水库的调度运行直接改变了径流的季节性分配[67].同时，洞庭湖流域也是我国植被变化最为显著的地区之一，该区域的森林覆盖率从1989年的36.7%上升到2016年的59.6%，在一定程度上对流域径流过程产生影响[68-69].此外，洞庭湖平原作为我国重要的商品粮基地，农业耗水量在2016年为101.22亿m3，占各类用水消耗量的70%以上，耗水率也高达50%以上[70]，尤其是农田用水量的峰值出现在秋季，与水库秋季蓄水、流域降水减少相叠加，更容易引发极端水文干旱事件[71-72].总的来说，人类活动对水资源的影响过程和机制是复杂的，在江-湖-河关系发生变化的背景下，一方面需要加强流域生态环境建设、优化水利工程调度方案，满足生产生活用水和生态用水需求[73]，另一方面也应发展节水农业，提高农业生产用水效率.而对定量评估不同人类活动对洞庭湖干旱的影响程度仍是目前研究的重点和难点，还需进一步通过定量模拟等手段开展深入研究.

4 结论

本文基于标准化水位指数SWI和标准化径流指数SRI，通过Copula函数对1964-2016年洞庭湖-流域-长江系统水文干旱联合概率分布特征进行研究.主要结论如下：

1)年尺度上，洞庭湖流域发生水文干旱的概率最高，但洞庭湖及长江在2003年以后发生水文干旱的频次明显较多，所占比重较大.季节尺度上，洞庭湖和长江2003年之后秋季水文干旱加剧，洞庭湖流域春季水文干旱概率远高于其它季节.

2)年尺度上, 洞庭湖-流域系统水文干旱联合概率大于洞庭湖-长江系统，但2003年后洞庭湖-长江系统联合水文干旱频次明显多于洞庭湖-流域系统.季节尺度上，洞庭湖-流域系统春季水文干旱联合概率最高，且两者同时发生水文干旱事件的次数最多；而洞庭湖-长江系统，其秋季水文干旱联合概率最大，并且一半以上的水文干旱事件发生在2003年以后.

3)洞庭湖春季水文干旱与流域的水文干旱具有较好的同步性，流域入湖补给减少对洞庭湖春季水文干旱的影响更大.而洞庭湖-长江系统秋季水文干旱自2003年以后更加极端和频发，这一方面受秋季降水减少和流域内人类活动的影响，造成流域径流对湖泊的补给减弱，另一方面三峡水库秋季蓄水使长江中下游干流水位降低，长江对湖泊顶托作用减弱也是重要原因之一.