考虑融雪量的多变量综合干旱指数的构建与评估.pdf

Water Resources and Hydropower EngineeringVol.54No.879CEditorial Department of Water Resources and Hydropower Engineering.This is an open access article under the CC BY-NC-ND license.JJ.Water Resources and Hydropower Engineering,2023,54(8):79-90.et al.Development and evaluation of a multivariate composite drought index accounting for snowHUYingying,ZHOUYuyan,LUFan90.胡莹莹，周毓彦，鲁帆，等考虑融雪量的多变量综合干旱指数的构建与评估J水利水电技术（中英文），2 0 2 3，54（8）：7 9-第54卷2023年第8 期水利水电技术（中英文）考虑融雪量的多变量综合干旱指数的构建与评估胡莹莹，周毓彦，鲁帆，肖伟华，侯保灯，刘建伟1，薛伟，张佳悦（1.中国水利水电科学研究院流域水循环模拟与调控国家重点实验室，北京100038;2.河海大学水文水资源学院，江苏南京210098)摘要：【目的】为解决大多数干旱评估中雪因素常被忽略的问题，需要开发适用于高寒融雪流域的干旱指数。【方法】根据潜在蒸散量、融雪量、降雨量计算水分增益/损失，然后利用Copula函数计算水分增加/损失与土壤水分的联合分布概率，从而构建了一个基于降雨、融雪量、潜在蒸散量和土壤水分新的标准化融雪综合干早指数（SSCDI），并比较和评估了3月和6 月时间尺度的SSCDI在黄河源区的干早监测性能。【结果】结果表明，Fischer-Hinzmann Copula是构建SSCDI的最佳连接函数；新构建的SSCDI与SPEI和MSDI相比，同SSI的拟合程度更好（R?高于0.7）；SSCDI比SPEI和MSDI更符合观测到的干旱变化。【结论】特别是在高寒流域，对雪因素的考虑提高了SSCDI的干旱监测能力，并且SSCDI能够同时监测出SPEI和MSDI所监测到的干旱情况，总体而言，SSCDI可用于识别和监测高寒区不同时间尺度的气象和农业综合干旱。关键词：干旱指数；融雪量；Copula函数；黄河源区与作者互动D0I:10.13928/ki.wrahe.2023.08.007开放科学（资源服务）标志码（OSID）：中图分类号：P429文献标志码：A文章编号：10 0 0-0 8 6 0(2 0 2 3)0 8-0 0 7 9-12Development and evaluation of a multivariate composite drought index accounting for snowHU Yingying,ZHOU Yuyan,LU Fan,XIAO Weihua,HOU Baodeng,LIU Jianwei-2,XUEWei,ZHANG Jiayue(1.State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin,China Institute of Water Resourcesand Hydropower Research,Beijing100038,China;2.College of Hydrology and Water Resources,Hohai University,Nanjing210098,Jiangsu,China)Abstract:ObjectiveJIn order to address the issue of neglecting snow factors in most drought assessments,it is necessary to de-velop a drought index that is suitable for alpine snowmelt watersheds.Methods Water gain/loss is calculated based on potentialevapotranspiration,snowmelt and rainfall.The Copula function is used to calculate the joint distribution probability of water gain/收稿日期：2 0 2 3-0 2-10；修回日期：2 0 2 3-0 3-11；录用日期：2 0 2 3-0 4-13；网络出版日期：2 0 2 3-0 6-2 0基金项目：2 0 2 2 年度科技智库青年人才计划（2 0 2 2 0 6 15ZZ07110156）；青海省中央引导地方科技发展资金项目（2 0 2 2 ZY020）；第二次青藏高原综合科学考察研究（2 0 19 QZKK0207-02）；国家自然青年科学基金项目（519 0 9 2 7 5）；中国水利水电科学研究院基本科研业务费专项项目（WR110145B0052021）；流域水循环模拟与调控国家重点实验室开放研究基金项目（IWHR-SKL-KF202204）作者简介：胡莹莹（19 9 0 一），女，博士研究生，主要从事生态水文研究。E-mail：h o o _y i n g 16 3.c o m通信作者：周毓彦（19 9 1一），男，高级工程师，博士，主要从事水文水资源研究。E-mail：z h y i w h r.c o m80水利水电技术（中英文）第54卷2023年第8 期胡莹莹，等/考虑融雪量的多变量综旱指数的构建与评估loss and soil water.A new standardized snow melt drought index(SSCDI)based on rainfall,snowmelt,potential evapotranspira-tion,and soil moisture is constructed.The drought monitoring performance of SSCDI at 3-month and 6-month time scales in thesource region of the Yellow River is compared and evaluated.Results J The results show that the Fischer-Hinzmann Copula isthe best connection function for constructing SSCDI.Compared with SPEI and MSDI,the newly constructed SSCDI fit SSI better(R?was higher than O.7).SSCDI is more consistent with observed changes in meteorological and agricultural drought than SPEIand MSDI.Conclusion JEspecially in alpine basins,consideration of snow improves the drought monitoring capability of SSC-DI.SSCDI can monitor the drought situation monitored by SPEI and MSDI at the same time,and can be used to identify and mo-nitor meteorological and agricultural drought at different time scales in alpine regions.Keywords:drought index;snowmelt;Copula function;the source regions of the Yellow River0引言干旱是常见的自然灾害之一，由于其持续时间长、影响范围广，往往造成巨大的经济损失1-3。随着全球平均气温升高，极端天气事件发生频率增加，导致气候系统的稳定性下降，全球大部分地区干旱化趋势日益显著，随之而来的干旱灾害增加，严重威胁到农业和生态安全4-5。2 0 2 1年政府间气候变化专门委员会（IPCC)第六次评估报告（AR6）指出，气候变化预测表明许多地区将经历更频繁和更严重的干旱6 。因此，有效的干旱监测、预测和早期预警对于决策和减轻负面影响的缓解措施至关重要。温度和大气蒸发需求是干旱变化的重要驱动因素，它们的增加预计会减少未来的降雪量7-9 ，降水将更有可能以雨而不是雪的形式落下，这将对以融雪为主的流域水文产生重大影响7 。由雪变化引起的长期地表水失衡和水分亏缺可能导致更频繁的水文和农业干旱，进一步威胁水安全10-1。特别是对于中高海拔地区，未来降雪量的减少可能会导致更严重的干旱12 。考虑到雪在干旱演变中的重要性和复杂性，有必要将雪因素纳人干旱监测因为干旱的定义复杂，许多干旱指数被逐渐发展出来表征干旱。干旱指数是研究干旱的基础和关键，只有合理的干旱指数才能准确监测和量化干旱，反映出干旱规律和趋势的变化13。然而，仍然很难建立一个能够监测和识别各种干旱的通用干旱指数14。AR6将干旱分为气象干旱、农业干旱、水文干旱和生态干旱6 。干旱类型的划分是相对的，不同类型之间的干旱在一定条件下可以相互转化。为了量化不同类型的干旱，过去已做出一些努力来开发各种干旱指标，依靠这些干旱指标，干旱评估已经有了长足的进步：标准化径流指数（SRI)15】、归一化植被指数（ND VI）【16 、植被状况指数(VCI)【17 、标准化土壤湿度指数（SSMI）【18 、标准化降水指数（SPI)【19 、帕尔默干旱指数（PDSI）2 0 、标准降水蒸散指数（SPEI)2 1 已经在世界范围内广泛使用2 2 。SPI和SPEI通常用于表征气象干旱，SRI用于表征水文干旱，NDVI和VCI用于表征生态干旱，SSMI用于量化农业干旱，PDSI用于评估基于降水、蒸散量和径流的综合干旱。然而，他们通常忽略了雪对干旱变异性的影响2 3。这些方法忽略了雪因素可能会降低它们监测和识别干旱的能力。因此，干旱指数需要考虑积雪和/或融雪6 。虽然STAUDINGER等人2 3 提出了标准化融雪和雨水指数，该指数考虑了影响干旱的雨水和融雪赤字，但它只评估了水文干旱，而且是利用温度和降水来模拟雪，并不是直接采用雪数据；ZHANG24最近建立了结合雪动力学的标准化水分异常指数（SZIsnow），但它具有计算复杂和困难等局限性。HAO和AGHAKOUCHAK等2 5 提出了一种多变量多指标干旱模拟方法一多变量标准化干旱指数（M SD I），该方法利用降水和土壤水分的联合分布函数来综合这两个变量的干旱信息。Javed等2 6 比较了SPI、标准土壤水分指数（SSI）和MSDI，以确定最适合中国的农业干旱指数，结果表明MSDI在监测中国大陆农业干旱方面表现更好。然而，由于MSDI没有考虑蒸散量和雪因素，因此其在干旱监测中存在一定的缺陷。VICENTESERRANO等人【2 1 在SPI的基础上开发了SPEI，由于它物理意义和简单的计算方法、多标量和多类别特征，SPEI已经广泛应用于表征气候变化下干旱监测2 7 。SPEI考虑了供水和需求，并利用降水和潜在蒸散量之间的差异来评估水分亏缺或盈余。然而，潜在蒸散量并不是影响需水量和干旱条件的唯一因素，因为土壤水分在估算需水量方面也起着重要作用2 0 黄河源区位于青藏高原和黄土高原的过渡地带，是典型的旱地高寒草地生态系统，对气候变化响应最为敏感，也对生态环境变化最为脆弱。近年来全球变81水利水电技术（中英文）第54卷2023年第8 期胡莹莹，等/考虑融雪量的多变量综合干旱指数的构建与评估暖的趋势对黄河源区生态环境产生了重要影响。因此，本研究以黄河源区为研究对象，综合考虑了MS-DI和SPEI的局限性，主要目标是：（1）将融雪量纳人干旱评估，利用Copula函数建立一种新的干旱指标SSCDI，它基于MSDI的联合概率方法，同时具有MSDI和SPEI的优点，并且可用于多时间尺度干旱评估；（2)评价新构建的干旱指数SSCDI在黄河源区干旱监测中的性能。1研究区概况黄河源区是指为以唐乃亥水文站为流域出口的汇水区域。如图1所示，黄河源区位于东经9 550 一10330，北纬32 10 36 15，面积为1.2 2 10 5km，海拔2 56 9 6 0 7 2 m。黄河源区流域年均径流量为2.0 510 m，占全流域的35%，是黄河流域重要的水源产流区2 8 。草地是黄河源区最重要的土地覆被，占整个源区的8 0%2 8 ，主要类型为高寒草甸、高寒沼泽草甸和高寒草原。黄河源区属典型的高原大陆性气候，兼具高寒气候和干旱气候的特点。黄河源区长期年平均降水量为50 2.7 2 mm，时空变化极不均匀2 9 ，其中西北部较为干燥，年降水量不足350mm，而东南部气候较为湿润，年降水量约7 50mm30)。气候变暖对黄河源区产生了重大的生态水文影响。特别是黄河源区的平均气温自19 8 0 年代以来以0.48/10 a的速度升高31。近40 年来，黄河源区降雪量以2.43mm/10a的速率增加，降雪率（降雪量占降水量的比例)多年平均约为2 6%31，黄河源区积雪期集中在11月到次年4月【32】，融雪期主要为每年的46 月和9 10 月332研究方法和数据2.1SPEI、SSI 和MSDI 的计算SPEI是一种广泛使用的气象干旱指数，用于说明大气湿度与平均状态的偏差34。SPEI计算的基本步骤为：（1)计算降水量与潜在蒸散量之差；（2）通过拟合一定时间尺度下差异的累积值来求出概率分布；（3)归一化累积概率密度函数以获得SPEI值，本研究采用Python计算gamma分布的SPEI。SSI 通过充分考虑土壤水分作为输人来描述干旱状况19 SSI是根据正常降水百分比和PalmerZ指数的概念以960000E980000E1000000E1020000E1040000EN.00.00.96N.00.00.96共和县贯南县泽库县山麻莱县同德县玛多县河南蒙古族自酒县称多县N.00.00.tN.00.00.t石渠生升德县玛出县图例达山果久治县若尔益县河流班书阿坝县湖泊松通县红原具黄河源区边界N.00.002cN.00.00.20海拔/m60722569km050100200960000E980000E1000000E1020000E1040000E图1研究区位置Fig.1Location of the study area82第54卷2023年第8 期水利水电技术（中英文）胡莹莹，等/考虑融雪量的多变量综合干旱指数的构建与评估及SPI的统计结构开发的35。MSDI是HAO等2 5 开发的基于SPI的指数的扩展，进一步将SPI扩展为基于降水和土壤水分的双变量模型。MSDI可以根据降水和土壤水分的联合概率计算得出2.2SSCDI 的构建2.2.1SSCDI的构建框架开发和评估SSCDI的框架如图2 所示。首先对NASA Global Land Data Assimilation System(GLDAS)数据进行预处理，以提取降水、降雨、潜在蒸散量、融雪和土壤水分。其次，根据潜在蒸散量、融雪量和降雨计算水分增减，然后用Copula函数计算水分增减与土壤水分之间的联合分布概率。在此基础上，本文选择了5个Copula函数来构建联合分布函数，通过赤池信息准则（AIC）、贝叶斯信息准则（BIC）、均方根误差（RMSE）和纳什-萨特克利夫效率（NSE）对每个候选Copula函数的性能进行评估，以选择最佳拟合的分布函数。最后，通过对联合分布概率的归一化得到 SSCDI。同时计算了SPEI、M SD I、SSI干旱指数，并将其用于SSCDI的可靠性评估。此外，根据从黄河源区实际干旱统计资料中选取的典型干旱年份，对SPEI、M SD I和SSCDI监测的干旱结果进行比较，验证SSCDI在干旱监测中的适用性。2.2.2SSCDI的构建方法本研究首先根据潜在蒸散量、融雪量、降雨量、土壤水分计算水分增益/损失，然后采用Copula函数计算水分增加/损失与土壤水分的联合分布概率，Copula函数可以导出两个或多个变量的联合分布函数而不管它们的原始边际分布如何。计算公式为P(D,d,S,s)=CF(D,),G(S,)=p(1)D,=R,+M,-PET,(2)式中，C是Copula函数；F(D)和G(S)分别是随机变量D，和S，的边际累积分布函数；p表示两个变量的联合分布概率；d和s表示需要满足的特定值；D,为水分增益/损失；S,是土壤水分；M，是融雪量；R;是降雨量；PET，是潜在蒸散量；i代表月。SSCDI=-i(p)(3)式中，为标准正态分布函数。通过将联合概率标准化，可以将其与不同的干旱指数进行交叉对比，与MSDI相似，它可以用于监测不同的时间尺度（1个月、2 个月、3个月、6 个月、9 个月等），不同的时间尺度可以反映不同的干旱类型，3月的尺度更能准确监测到黄河源区的季节性干旱事件，6 月的时间尺度能较好地反映出黄河源地区干季和雨季的变化。因此，在本研究中使用3月和6 月尺度的SSCDI表征黄土壤水分降水潜在蒸散量降雨融雪量MSDISPEI计算水分增加/损失111构建SSCDI111SSI1AIC1111111Copula函数优选BIC计算SPEI&SSI、M SD I&SSI、SSCD I&SSI的线性回归拟合结果1RMSESSCDI与SPEI和MSDI的干旱监测性能对比分析SSCDINSE气象灾害的统计数据各指数在典型干旱年份的干旱监测能力比较分析图2SSCDI的构建流程Fig.2Flowchart of the SSCDI construction83第54卷2023年第8 期水利水电技术（中英文）胡莹莹，等/考虑融雪量的多变量综干旱指数的构建与评估河源区的干旱状况。如表1所列，SSCDI负值代表干旱，而正值表示潮湿表1SSCDI值的等级划分Table 1Classification of SSCDI valuesSSCDI值的范围干早分类等级-0.5SSCDI0.5正常D1-1 SSCDI-0.5轻度干旱D2-1.5SSCDI-1中度干旱D3-2 0sign=(x;-x,)(y;-yj)=0(7)-1,(x;-x)(y:-y)0式中，n代表数据系列的总长度表2本研究中使用的5个不同Copula函数的数学描述Table2Classification of SSCDI values Mathematical descriptions of five different Copula functions used in this study名称Copula函数参数范围Claytonmax(u0+-1,0)-1/00E-1,8(exp(-Qu)-1)(exp(-u)-1)FrankQERexp(-0)-10,E 0,1,10,min(u,0)+(1-0,)um/0Fischer-Hinzmann02 ER0,E0,BB11/662=1,1,02E0,1,Tawnexp/In(u-(-0,ln(u),E1,84水利水电技术（中英文）第54卷2023年第8 期胡莹莹，等/考虑融雪量的多变量综合干旱指数的构建与评估(3)Spearmans 相关系数pZ(R,-R)(Q.-Q)ni=1P(8)2(R,-R)22(Q.-Q)n式中，R,和Q.分别是样本值；x，y:在序列中取值等级的顺序；R，Q 分别是是样本R,和Q,的平均等级2.4GLDAS陆地表面数据集GLDAS-2是通过使用全球气象强迫数据集作为输人，驱动Noah陆面过程模型（LSM）来模拟全球陆面过程而获得各地表数据4。GLDAS-2 能够较好地捕捉到土壤水分的时间变化，与观测土壤水分的相关系数在0.5以上45，同时，GLDAS-2中的融雪量数据也是根据Noah陆面模式的算法计算出来的，该数据为推进与雪相关的水文效应研究提供了独特的机会，因为该数据集提供了高质量、长期、全球覆盖和多变量的地表和地表条件估计46 ，它准确再现雪过程的能力是我们干旱指数发展的重要基石。本研究中采用19 48 年1月至2 0 14年12 月的CLDAS-2.0NoahLSM月度数据和2 0 0 0 年1月至2 0 2 1年12 月的GLDAS-2.1NoahLSM月度数据，它们的分辨率均为0.25度。GLDAS-2.0和CLDAS-2.1数据集都可以通过戈达德地球科学数据和信息服务中心获得(http:/disc.sci.gsfc.nasa.gov/hydrology/data-hold-ings）。分析中使用的各种参数包括降水量、降雨量、潜在蒸散量、融雪量和土壤水分，对其进行距离加权平均重映射进行插值，插值精度为0.0 5，该插值方法是对所有输人字段的四个最近的邻居值进行距离加权平均重映射。GLDAS提供四个土层的水分数据，分别为0 10 cm、10 40 c m、40 10 0 c m 和10 0 2 0 0cm，在本研究中我们采用顶层（即0 10 cm）的土壤水分数据2.5气象灾害统计数据青海省历年干旱灾害情况数据来源于19 45一1995中国灾情报告中国气象灾害统计年鉴气象干旱年鉴和中国气象灾害大典青海卷。这些数据用于反映现实世界干旱的发生情况，因此用于验证本研究中提出的SSCDI在黄河源区干旱监测中的性能。中国气象灾害大典-青海卷数据时间范围为1950一2 0 0 0 年，下载自国家青藏高原科学数据中心(https:/ 0 0 4一2 0 2 0 年，气象干旱年鉴时间范围为20112017年。3研究结果与分析3.1最佳拟合Copula优选结果验证变量之间的相依性结构是正确进行Copula分析的关键，如表3所列，3月和6 月尺度的D,和S；的Pearson，K e n d a l l 和Spearman的相关系数值高达0.79，0.6 3和0.8 3，这表明D,和S,具有较强的相关关系，同时可以看出，3月尺度的三个相关系数均大于6 月尺度的相关系数。表3D,和S,在3月和6 月时间尺度上的相关性结果Table 3Correlation between D,and S,at 3-month and 6-monthscales时间尺度Pearson相关系数Kendall 相关系数Spearman 相关系数3月0.790.630.836月0.720.530.745种Copula函数用于推导D,和S,，的联合概率分布的拟合优度检验结果如表4所列。从表中可以看出，5种候选Copula函数的3月和6 月尺度RMSE，NSE，A IC，BIC检验结果非常接近，各时间尺度的平均值分别为0.2 2，0.9 9，-6 7 13，-6 7 0 4，表明所选Copula函数对变量的拟合效果都较好，具体来看，Fischer-HinzmannCopula在各时间尺度上的RMSE值、BIC值和AIC值最小，且NSE值最大，因此Fischer-HinzmannCopula被选为最佳拟合分布，故而本研究构建SSCDI时使用Fischer-HinzmannCopula作为连接函数。表45个Copula函数的拟合度测试评估结果Table 4 The goodness-of-fit test assessment of fiveCopula functions时间排Copula家族RMSENSEBICAIC尺度名Clayton0.3060.990-6403.973-6408.7624Frank0.3160.990-6381.024-6385.81353月Fischer-Hinzmann0.0830.997-7 307.620-7 317.1981BB10.1870.993-6683.1876692.7653Tawn0.1600.994-6 746.010-6760.3772Clayton0.1630.974-5 346.714-5351.5035Frank0.3760.993-6 472.601-6 477.39046月Fischer-Hinzmann0.0740.999-7760.670-7770.2481BB10.2820.996-6.976.806-6986.3833Tawn0.2820.996-6 970.663-6985.03023.2SSCDI与其他干旱指数的对比验证为了解SSCDI在量化干旱事件方面的表现，我们将其与文献中其他已经建立的干旱指数进行比较85水利水电技术（中英文）第54卷2023年第8 期胡莹莹，等/考虑融雪量的多变量综合干旱指数的构建与评估采用了SSI来初步评估SSCDI，过去虽然已经开发了多种能够衡量农业干旱的指数，例如PALMER开发的作物湿度指数（CMI），CA M M A L L ERI等人47 开发的土壤水分干旱严重程度指数（DSI）等等，然而这些指数只能采用固定时间尺度进行分析，而SSI可以在多时间尺度上衡量农业干旱，因此可以在不同尺度上与SSCDI进行对比。基于Copula函数的优选结果，以3月和6 月尺度为例，计算了黄河源区SPEI、SSI、MSDI和SSCDI的指标值，并记为SPEI-3、SSI-3、MSDI-3和 SSCDI-3，以及SPEI-6、SSI-6、M SD I-6 和SSCDI-6。计算一定时间尺度下干旱指数之间的线性回归拟合程度的目的是评价在相同条件下，哪种干旱指数在干旱识别方面表现最好，不同时间尺度的SSI与SPEI、M SD I和SSCDI的散点关系图如图3所示，其中3-月SSI-SSCDI的R?值为0.7 11，高于SSI-SPEI（0.52 5）和SSI-MSDI(0.630），6 月尺度SSI-SSCDI的R?值为0.7 32，同样高于SSI-SPEI(0.694）和SSI-MS-DI（0.57 1），并且两种时间尺度上SSI-SSCDI的RMSE值都最小。因此，SSCDI在各时间尺度上和空间位置上的误差是最低的，表明SSCDI具有更好的干旱监测性能贵南县和玛多县SPEI、M SD I 和SSCDI在3月尺度的时间序列如图4和图5所示。其中，玛多县平均海拔430 0 m，多年平均融雪量15.8 3mm，而贵南县平均海拔310 0 m，多年平均融雪量3.16 mm，远远低于玛多县，因此本文称玛多县为融雪区域，贵南县为非融雪区域，将两者进行对比分析。如图4所示，非融雪区域的各干旱指数走势和变化规律具有一致性，SSCDI-3折线与大多数MSDI-3折线重叠，SPEI-3、M SD I-3和SSCDI-3的干旱开始和结束时间相似，这意味着SPEI-3、M SD I-3和SSCDI-3的干旱监测性能相似，原因在于非融雪区域雪对干旱表征的影响较小。融雪区域中三种干旱指数显示的干湿状况并不完全一致（见图5），根据气象干旱年鉴中2 0 16 年中国各省干旱事件统计表可知，2 0 16 年8 月下旬至9月上旬青海省发生了特级干旱，如图5(b）所示，SS-CDI指数值小于-2，监测出了该特干旱事件，而SPEI和MSDI仅监测出了中等干旱。说明融雪量的加人对干旱的评估具有一定的改善作用。3.3典型干旱年份的比较分析为了进一步测试SSCDI的干旱监测性能，选择了典型的干旱年份作为案例研究。据中国气象灾害大典-青海卷和中国气象灾害年鉴记载，青海省19601961年全省发生持续干旱早。2 0 13年34月，青海省遭遇严重干旱，全省重旱、特旱面积居全国之首。因此，本研究中，以19 6 0 年和2 0 13年为案例，评估SPEI、M SD I 和SSCDI在3月和6 月时间尺度上的干旱监测性能，因为研究区域在这些年份具有典型的干旱状态48 1960年7 月和2 0 13年4月SPEI、M SD I 和SSCDI444R2=0.525R2=0.630R2=0.711333RMSE=0.611RMSE=0.539RMSE=0.478222yuoul-s111IdSIASW0001-1-2-220-3-3-3-3-2-1012344-3-2-1012344-3-2-10134SSISSISSI1444R?=0.571R2=0.694R2=0.732333RMSE=0.580RMSE=0.488RMSE=0.458-2222yuow-g111IdSIASWOSS0003-1-2-22-3-3-31-4-3-2-101234-4-3-2-1012344-3-2-10134SSISSISSI注：颜色是指点的密度，从最高（红色）到最低（蓝色），黑色实线代表1：1线。图3SPEI、M SD I、SSCD I与SSI在3月和6 月时间尺度上的散点图Fig.3Scatter plots of SPEI,MSDI,and SSCDI versus SSI on 3-month and 6-month time scales86水利水电技术（中英文）第54卷2023年第8 期胡莹莹，等/考虑融雪量的多变量综合干旱指数的构建与评估3220O2231960-031960-071960-101961-031961-06 1961-092016-012016-042016-072016-10时间时间SPEIMSDISSCDI43201231954-101965-091976-081987-081998-072009-072020-06时间图4贵南县19 48 2 0 2 1年SPEI-3、M SD I-3和SSCDI-3时间序列干旱指数比较Fig.4Comparison of three drought indices of SPEI-3,MSDI-3,and SSCDI-3 time series for Guinan County from 1948 to 2021在3月和6 月时间尺度上的干旱监测结果如图6 和图7 所示。总体而言，与气象干旱相比，农业干旱影响的区域更大。SSCDI成功监测出了气象干旱和农业干旱，此外，SSCDI的干旱空间格局与MSDI相似。图中蓝色的干旱等级为D1，表示正常或偏湿润，颜色越深代表干旱等级越高。从图6 中紫色圆圈标记的区域可以看出，19 6 0 年对于MSDI-3未能监测出干旱而SPEI-3监测出D3干旱的区域，SSC-DI-3同样监测出了D3干旱，2 0 13年对于MSDI-3未能监测到干旱而SPEI-3监测出D2干旱的区域，SSCDI-3同样监测出了D2干旱。从图7 中紫色圆圈标记的区域可以看出，19 6 0 年对于SPEI-6未能监测出干旱，而MSDI-6监测出干旱的区域，SSCDI-6同样监测出了干旱情况，2 0 13年对于MSDI-6未能监测出干旱，而SPEI-6监测出D2干旱的区域，SS-CDI-6监测出了D3干旱因此可以得出，SSCDI能够同时监测出SPEI或MSDI监测到的干旱，且SSCDI量化的干旱等级大于或等于MSDI，可能原因在于MSDI没有考虑蒸散发因素。由此我们可以得出结论，SSCDI在黄河源区干旱监测中表现良好4结论基于Copula理论，本文将融雪量纳人干旱评估体系，提出了一种新的综合干旱指数SSCDI。虽然本研究是基于GLDAS数据分析的，但SSCDI同样适用于站点数据，它是多变量的，可以在多个时间尺度上进行计算，并且计算过程简单清晰。该指数综合考虑了降雨、融雪量、潜在蒸散量和土壤水分，为干旱量化（特别是高纬度地区的干旱量化）带来了更广泛的视角。（1）SSCD I是在SPEI和MSDI的基础上制定的，并且对5种Copula函数进行拟合度优选，得出了构建SSCDI时使用Fischer-HinzmannCopula作为连接函数。（2）通过与SSI进行线性回归及误差分析，结果显示在3月和6 月时间尺度上SSCDI与SSI拟合度更高（R分别达到了0.7 11和0.7 32）。SSCDI在融雪87水利水电技术（中英文）第54卷2023年第8 期胡莹莹，等/考虑融雪量的多变量综合干旱指数的构建与评估3210-1-2-23-341960-041960-101961-041961-102016-012016-042016-072016-10时间时间SPEIMSDI-SSCDI432E-IASSO1954-101965-091976-081987-081998-072009-072020-06时间图5王玛多县19 48 2 0 2 1年SPEI-3、M SD I-3和SSCD/-3时间序列干旱指数的比较Fig.5 Comparison of three drought indices of the SPEI-3,MSDI-3,and SSCDI-3 time series for Maduo County from 1948 to 2021SPEIMSDISSCDILO-0961Jkm0170340D1D2D3D4D5图61960年7 月和2 0 13年4月SPEI、M SD/和SSCD/在3个月时间尺度上的空间格局Fig.6Spatial patterns of SPEI,MSDI,and SSCDI on a 3-month time scale in July 1960 and April 201388第54卷胡莹莹，等/考虑融雪量的多变量综合干旱指数的构建与评估SPEIMSDISSCDIL0-0961Jkm0170340D1D2D3D4D5图71960年7 月和2 0 13年4月SPEI、M SD/和SSCD/在6 个月时间尺度上的空间格局Fig.7Spatial patterns of SPEI,MSDI,and SSCDI on a 6-month time scale in July 1960 and April 2013较多的区域表现更佳，而在非融雪区域与SPEI和MSDI具有较高的一致性，因此，一个流域受雪的影响越多，就越值得用SSCDI来评价其干旱情况（3)在黄河源区将这三者的结果与实际统计资料的干旱情况进行了对比，验证了SSCDI在干旱监测中的适用性。以19 6 0 年7 月和2 0 13年4月为例，分析了SSCDI、SPEI、M SD I 在3月和6 月尺度上的干旱监测性能，结果表明，SSCDI能够同时监测出SPEI和MSDI监测到的干旱，纳人融雪量使SSCDI在高寒区能够产生更可靠的干旱监测结果综上所述，构建的干旱指数SSCDI在高寒地区的干旱量化方面具有良好的监测效果，适用于气象和农业干旱的综合评估。然而，尽管构建的干旱指数很新颖，但仍存在一些局限性，本文仅使用GLDAS数据对SSCDI进行了验证和分析，未来还需要利用多种遥感数据和实测数据进行研究，进一步验证SSCDI的可靠性。参考文献(References)：1MUKHERJEE S,MISHRA A,TRENBERTHK E.Climate changeand drought:a perspective on drought indices J.Current ClimateChange Reports,2018,4(2):145-163.2FUNK C,SHUKLA S.Drought early warningdefinitions,challen-ges,and opportunitiesM.Elsevier,2020:23-42.3MUKHERJEE S,MISHRA A,TRENBERTH K E.Climate changeand drought:a perspective on drought indices J.Current