淮河小流域多源降水数据产品径流模拟对比分析.pdf
《淮河小流域多源降水数据产品径流模拟对比分析.pdf》由会员分享,可在线阅读,更多相关《淮河小流域多源降水数据产品径流模拟对比分析.pdf(10页珍藏版)》请在咨信网上搜索。
1、第54卷第8 期2023年8 月文章编号:10 0 1-417 9(2 0 2 3)0 8-0 10 9-0 9引用本文:渠畅,魏玲娜,董建志,等.淮河小流域多源降水数据产品径流模拟对比分析J.人民长江,2 0 2 3,54(8):10 9-117,130.人民长江YangtzeRiverVol.54,No.8Aug.,2 0 2 3淮河小流域多源降水数据产品径流模拟对比分析渠畅12,魏玲娜1.2,董建志3.4,徐士惠,玮6,7(1.南京信息工程大学水利部水文气象灾害机理与预警重点实验室,江苏南京2 10 0 44;2.南京信息工程大学水文与水资源工程学院,江苏南京2 10 0 44;3.天津
2、大学地球系统科学学院,天津30 0 0 7 2;4.天津大学表层地球系统科学研究院,天津30 0 0 7 2;5.中国气象科学研究院浙江分院,浙江杭州31130 0;6.抚河水文水资源监测中心,江西抚州3440 0 0;7.鄱阳湖水文生态监测研究重点实验室,江西南昌330 0 0 2)摘要:为研究多源降水数据产品的地区适用性,对比不同产品的水文模拟效果,选取淮河上游紫罗山、中游黄泥庄两个流域,对比分析CHIRPS、C M O RPH、ERA 5、I M ERG 4种产品模拟流域平均降水的误差特征,利用集总式FLEX模型对产品模拟结果进行比较,探讨了产品在准河流域不同位置对径流模拟效果的影响。结
3、果表明:两个流域CHIRPS降水误差最大,CMORPH次之ERA5和IMERG的精度较高,上游紫罗山流域多源数据集整体高估;不同降水数据产品径流模拟性能不一,模拟效果与产品精度评估结果基本一致,总体ERA5和IMERG模拟径流与实测吻合较好;总体而言,4种数据集有效反映了流域的空间水文特性,上游紫罗山流域径流模拟效果优于中游黄泥庄流域。关键词:径流模拟;FLEX模型;卫星降水产品;淮河流域中图法分类号:P3390引言降水是一种时空变异特性显著的水文气象要素,陆地径流的产生、土壤水分的变化和蒸散发的形成等流域水文循环过程无一不与降水分配密切相关。作为水文模型的关键输人变量2-3,降水数据的准确性
4、对水文模拟的不确定性影响极大4。因此,高精度大范围的连续降水数据对流域水文过程模拟及预报具有重要意义。目前获取降水数据的途径主要有地面雨量计观测、雷达探测、数值天气模式模拟以及卫星遥感产品5-6 等。传统的降水信息利用地面站点雨量计实测获得,例如自动雨量站等,尽管时间分辨率较高,但仍不能满足空间覆盖率的需求,尤其对于时空变化剧烈收稿日期:2 0 2 2-10-2 2基金项目:国家自然科学基金项目(5190 912 1,52 17 90 2 1)作者简介:渠畅,男,硕士研究生,研究方向为水文气象。Em a i l:q c s t r u g g l e 16 3.c o m通信作者:魏玲娜,女,
5、副教授,博士,主要从事流域水循环模拟与水文气象研究。E-mail:文献标志码:AD0I:10.16232/ki.1001-4179.2023.08.015的热带地区7 ,地面站点观测更是无能为力。卫星遥感技术能够提供具有一定时空精度的连续降水数据,其空间覆盖范围广、获取便捷,可在一定程度上弥补其他来源降水数据的缺陷 ,为水文气象研究与实际应用提供有力支持,尤其给缺资料地区提供了新的数据来源9-10 1。卫星遥感技术快速发展促进了卫星降水数据产品的日趋成熟和完善,同时大大推动了卫星产品在水文、气象、农业、环境等诸多领域的应用。现有的卫星降水数据产品已显示了全球化的监测与预测潜力,例如IMERG的
6、时空分辨率已能达到0.1和30min。同时,随着观测手段的多样化和数据同化技术的成熟,对地面、卫星、雷达、探空、飞机、船舶等多种观测资料进行质量控制,利用全球气候系统模式对实际的大气状况进行模拟,可获得高分辨率、长时间序列的全110球降水数据,例如欧洲中期天气预报中心的ERA-In-terim和ERA5,以及美国国家环境预报中心的CFSR等等国内外大量研究是针对降水数据产品的精度和地区特征进行的12-14。例如,Xin 等15 在粤港澳大湾地区比较了ERA5-Land和ERA5-HRES两种产品,发现后者在探测降水强度和降水事件方面更优,同时发现ERA5在城市化水平较高的地区降水探测表现不佳。
7、Huang等16 以40 0 多个台湾当地雨量计观测值为基准,对比发现IMERG产品在数量和质量上都优于TRMM,可以更准确地描述台湾夏季降水在多个时间尺度的变化。Tang等17 对IMERG在中国日尺度和小时尺度的产品质量综合分析后发现IMERG产品优于除GSMaP以外的其他数据集(TRMM3B42、CM O RPH、PERSIAN-CDR、G SM a P、C H I RPS、SM 2 RA I N、ERA 5、ERA-Interim和MERRA2),且产品质量随时间推移不断提高。在此基础上,学者们也在不断尝试将各种数据产品应用到水文模拟中,对其适用性能加以探讨。例如,邹磊等18 曾利用分
8、布式时变增益模型对元江流域多套降水产品的可替代性进行研究;许心怡等19 利用滦河流域IMERG和CMFD(中国区域地面气象要素驱动数据集)开展定量多尺度时空精度评估,借助SWAT模型对比分析了两种产品的水文效用。尽管降水数据产品驱动水文模拟领域的研究已积累了一定的成果,但对于具体流域不同降水产品的误差如何,其误差又是如何影响水文模拟结果的仍有待进一步研究。淮河流域位于中国南北气候过渡带,气候条件复杂,旱涝灾害严重,区域内降水集中,导致在地势起伏较大的山丘区,小流域山洪源近流急、峰高量大,破坏力强,但雨量站网信息不足,降水量分布估计存在较大误差影响洪水预报的精度。因此,有必要针对淮河小流域的多源
9、降水产品的可替代性和水文模拟的适应性加以研究。本文选取淮河上游紫罗山和中游黄泥庄两个典型流域,对比分析CHIRPS、CM O RPH、ERA 5、I M E-RG这4种常用降水数据产品面雨量的误差特征,分析其对径流模拟效果的影响,以此评估这4种数据集在紫罗山与黄泥庄两流域水文过程模拟中的精度。1石研究流域与数据1.1沃流域概况选择淮河上游的紫罗山和中游的黄泥庄作为典型流域进行研究,研究区域如图1所示。其中紫罗山流域为紫罗山站(112 31 E,3410N)以上集水区,面积1800km位于淮河西部沙颖河水系北汝河上游,属人民长江半干旱大陆性气候,多年平降水量7 58 mm,降水主要集中在汛期,水
10、面蒸发量9 54.3mm,流域地形起伏较大,海拔2 90 2 12 0 m。黄泥庄流域为黄泥庄站(11537E,3128N)以上集水区,面积8 0 5km,位于大别山北麓,淮河流域史河源头,受东亚季风活动控制,属于湿润地区,年均气温12 2 1,多年平均降水量约1380mm,年内分配不均,主要集中在6 9月。流域呈不规则扇形,地形以丘陵和山地为主,海拔130 1500m,河谷狭长,水系发育良好。两流域植被覆盖良好,受人类活动影响较小,适合开展天然状态下的水文过程模拟研究。1120E11210E11220E11230E11520E11530E11540EN,O.tEN,OS.5EN,Ot.c5(
11、b)紫罗山流域(b)ZiluoshanBasin125250km(a)淮河流域(a)HuaiheBasin图1研究流域位置、DEM与站点分布Fig.1 Location of the study watersheds,their DEM and gauges1.2数据收集与处理20102018年的逐日流量数据和地面雨量数据来源于中华人民共和国水文年鉴中的淮河流域水文资料。黄泥庄流域内有禅堂站、关庙站、吴店站、斑竹园站、西河站、银沙站、徐坳站、马鬃岭站8 个雨量站以及黄泥庄水文站共9个站点;紫罗山流域内有孙店站、龙王庙站、两河口站、禅螳站、木植街站、黄庄站、排路站、沙坪站、付店站、娄子沟站、十八
12、盘站、秦亭站、王坪站、三屯站14个雨量站以及紫罗山水文站共15个站点。收集整理紫罗山流域内汝阳站与黄泥庄流域临近气象站点金寨站2 0 10 2 0 18 年逐日的日平均气温、日最高气温、日最低气温,分别用于模型蒸发量计算。选用CHIRPS、C M O RPH、I M ERG 卫星降水产品和ERA5再分析降水数据,时间为2 0 10 年1月1日至2018年12 月31日,采用最邻近法【19 将空间分辨率统一至0.10.1以便对比不同数据产品的流域平均面雨量精度,同时计算流域内所有数据集格点的算术平均值作为流域平均降水量。各降水数据产品的信息如表1 所列。2023年N.O.16N,OI.1S(c)
13、黄泥庄流域N(c)Huangnizhuang Basin高程:m水文站22200。雨量站130口降水数据产品格点流域边界水素第8 期Tab.1Information of the four precipitation data products产品名称研发背景CHIRPS由美国地质调查局数据时间序列长(198 1https:/www.chc.ucsb.edu/(U SG S)和加利福尼亚年至今),是多时间尺data/chirps大学气候危害小组开发度、多空间分辨率降水数据集,全球覆盖范围为50 S50N,空间分辨率较高,广泛应用于中低纬度区域内的降水、干早研究CMORPH基于美国NOAA气候预
14、集合了多平台卫星制作https:/www.cpc.ncep.noaa.测中心(CPC)创立发展而成,目前在国内多应gov/products/janowiak/cmor-的一种利用时间空间联用于数据反演、区域大ph_description.html合内插的全球高时空分尺度气候预报分析,降辨率降水产品的技术水产品评估等领域19IMERC基于GPM卫星群获取较高的空间分辨率以及https:/disc.gsfc.nasa.gov/的大量遥感信息、微波、覆盖率、多时间分辨率,红外和雷达数据反演所可精确探测微量降水,产生的降水数据集,来产品可靠性高【2 0 】自美国NASA主导的TRMM的后续全球卫星降水
15、观测计划ERA5欧洲中期天气预报中心采用改进的数据同化系htps:/w w w.e c mw f.i n t/e n/(ECM W F)发布的第五统(IFSCycle41r2)、四维forecasts/datasets/reanalysis代再分析数据集2研究方法2.1FLEX模型与构建基于 FLEX(Flu x Ex c h a n g e h y d r o lo g ic a l m o d e l)模型构建刻画不同水文过程响应单元、适应不同流域所需具体应用目标的水文模型2 0-2 。该模型在搭建过程中可以根据研究区域的实际气候、地形等下垫面特征对水箱进行增删或对参数加以修改,其结构灵活
16、度高,区域适应性强。考虑到研究流域的面积较小,水文要素的空间差异较小,集总式概念模型比分布式模型有更好的参数化灵活性和更强的误差过滤能力2 3,本文利用FLEX模型框架构建了包含:植被冠层截留水箱(S.)、非饱和土壤层水箱(S.)、快速退水水箱(S.)以及慢速退水水箱(S.)的四水箱FLEX水文模型,分别用以模拟植被截留、非饱和土壤水产汇流等环节,反映降水经过植被的截留、土壤与植被的蒸散发(见图2)。模型包含了11个主要参数,参数的含义及取值范围如表2 所列2 4 针对气象站点高程(eleu_CM)和流域平均高程(e le v _P4)的差别对流域进行温度修正,用式(1)计算渠畅,等:准河小流
17、域多源降水数据产品径流模拟对比分析表1四种降水数据产品信息特点变分同化分析(4D-d a t a s e t s/e r a 5Var)以及辐射变分偏差校正技术2 1,相较于前代再分析数据集降水,数据集的质量得到系统性提升111修正因子(Tcalt)。潜在蒸散发(E,)计算采用Har-greaves 公式 2 5 来源参数S iMax植被截留能力,用于表明最大截留量SuMax最大储水量,用于表明包气带的土壤层根区存储水能力C包气带水分中消耗于实际蒸腾的比例系数包气带蓄水容量曲线抛物线型指数,用于表征景观空间异质性,进行径流系数的计算D分径流成分比例参数,将径流分成直接径流和深层地【0.1,0.
18、9下径流分别流入S,和S两个水箱Kf浅层地下径流(壤中流)退水参数,为水箱S,中的壤中流流人河网的时间K深层地下径流(壤中流)退水参数,为水箱S。中的基流汇人河网的时间TiagF水流在S水箱中的滞时水流在S水箱中的滞时SmMax直接径流(包括坡面流和壤中流)阈值参数,可理解为S水箱的最大储水量Kf地表径流退水参数,为超过S,储水能力的地表径流流人河网的时间降水经过S,水箱首先进行植被冠层截留,当截留量超过该水箱的最大储水量SiMax时,多余的水到达非饱和土壤层S.水箱,此处采用新安江模型的不饱和蓄水层的水量计算方法,流域平均蓄水容量W为W.=1+非饱和土壤层S.中的水量部分消耗于实际蒸散发,这
19、里的实际蒸发蒸腾量E。与潜在蒸散发E。与土壤相对湿度S./SuMax成正比,计算公式为E.=min(E,E,3S.uMa当流入S.水箱的水量大于最大储水量SuMax时,多elev_P4-elev_CMTcalt=100 0.6E.个PSCPuSaMaxS.图2 FLEX模型结构Fig.2Structure of FLEX model表2FLEX模型使用的参数Tab.2Parameters of FLEX model参数含义SuMax(1)TuneFRDRsTuesRsIS.RaLKQnS,K,QrK.Q取值范围0.1,5 20,8000.2,0.80.1,31,20 10,4000.8,32,
20、50 10,2001,10(2)(3)112余的水量由比例系数D分为两部分,一部分在流向河网过程中发生人渗,经滞时Tlags后流向慢速退水水箱S,最终成为深层地下径流汇人河网;另一部分经滞时Tlagr后流向快速退水水箱Sr,当快速退水水箱的储水量超过其储水能力SrMax时,多余的水量经过Kr时间后转化为地表径流或壤中流流人河网。最终流域出口断面出流就是地表径流、壤中流和地下径流之合。模型的输人数据为流域平均逐日的降水数据、气象站逐日温度(日平均气温、日最高气温、日最低气温)数据和流域出口断面逐日流量数据。选取2 0 10 2014年为模型率定期,2 0 15 2 0 18 年为验证期,模拟两个
21、流域的日流量变化过程。由MOSCEMU A(Multi-objective Shuffled Complex Evolution Metropo-lis)算法2 6 ,选取KCE(K l i n g-G u p t a e f f i c i e n c y)作为目标函数(KGE最优值为1),同时以开方形式与绝对值形式的KGE作为辅助判据,优化确定模型参数。2.2附降水误差评估采用相对偏差BIAS、均方根误差RMSE和相关系数r三项指标来定量评价ERA5,C H I RPS,I M ERG,CMORPH4种产品与地面雨量站观测算得的流域平均降水量在日尺度上的一致性。以地面雨量站观测为基准对比不
22、同数据集的降水量误差,BIAS值反映了数据产品的系统偏差程度;RMSE值用来衡量产品偏离地面观测的程度,它对误差极值敏感,可反映产品的精密度;r值表征数据集与地面观测之间的线性相关程度。BIAS值和RMSE值反映了偏差误差的大小。BIAS=RMSE=Cov(S,C)r=VVar S VarC式中:C,表示地面雨量站观测日降水量;S,表示产品的日降水量数据,i=1,2,3,n 表示时间序列;Cou(S,C)为S与C的协方差,VarS为S的方差,VarC为 C 的方差。采用探测率Pd、空报率RrA和Heidke技巧评分指数SH三项分类评价指标来反映产品对日降水事件的探测能力。P.越高,表明数据集对
23、日降水事件的漏报程度越小2 7 ;Rr越小,空报程度越小;S综合反映了卫星对降水事件是否发生的估计能力,表示实际预报的准确率比随机预报的准确率好多少,其值范围在-1,1之间,预报完全正确时,S=1。人民长江式中:n.表示降水数据产品和地面观测基准数据同时有雨的频数;n1o表示降水数据产品有雨地面观测无雨的频数;noi表示降水数据产品无雨地面观测有雨的频数;noo表示降水数据产品和地面观测同时无雨的频数。本次研究重点在于理解模型利用多源降水数据导致的输人误差对径流模拟结果影响,而非关注流域内部降水产品的空间精度,本文并未考虑多源降水产品在空间分布上的精度,而将主要篇幅集中于评估多源降水产品对流域
24、面平均降水量的刻画能力。2.3径流模拟效果评价由于水文模型参数的优化常会出现“异参同效”问题,这种参数的不确定性会对评判参数优化标准、提高水文模型模拟精度造成阻碍。Beven认为如果采用单一的目标函数进行模型模拟结果的评判,容易弱化部分模拟结果的特性,具有局限性2 8 。为了尽可能避免模拟结果受某组模型参数影响,同时兼顾径流模拟效果评价指标的合理性,本文采用比较流域出口断面的径流变化过程线图以及纳什效率系数NSE和Kling-Gupta效率系数KGE两种方式评价模拟效果。Z(S:-C.)NSE是常用于评价验证水文模型模拟结果的目标n(4)(5)(6)2023年nlPa=nm+noln10RrA
25、=n+nio2(nnoo-nonol)函数。KCE充分考虑模拟值和实测值的均值和误差,是对NSE的一种补充描述,研究表明,使用KGE更能充分反映流域水文过程2 9NSE=1-nKGE=1-/(r-1)+(-1)+(-1)sim二uobs0sim=Qobs式中:Qsim和Qobs.i分别是模拟和实测径流量,mm/d;r为模拟和实测径流间的线性相关系数;sim和obs分别是模拟和实测径流量的平均值;sim和obs分别是模拟和实测径流量的均方差。(7)(8)(9)(10)(11)(12)(13)第8 期3结果和讨论3.17不同降水产品的面雨量精度评价图3给出了紫罗山流域4种产品与地面站点流域平均日降
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 淮河 小流域 降水 数据 产品 径流 模拟 对比 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。