1、第4 1卷第5期2023年5月水 电 能 源 科 学W a t e r R e s o u r c e s a n d P o w e rV o l.4 1 N o.5M a y 2 0 2 3D O I:1 0.2 0 0 4 0/j.c n k i.1 0 0 0-7 7 0 9.2 0 2 3.2 0 2 2 0 7 6 8干旱区内陆河流域植被恢复对水文过程的影响翟婧雅a,b,金 鑫a,b,c,金彦香a,b,c,傅 笛a,b(青海师范大学 a.地理科学学院;b.青海省自然地理与环境过程重点实验室;c.高原科学与可持续发展研究院,青海 西宁 8 1 0 0 1 6)摘要:为更好地理解植被恢
2、复对水文过程的影响,以地表水-地下水交互频繁且植被恢复明显的柴达木盆地东北部巴音河流域中下游为例,将S WA T-MO D F L OW模型与动态HR U相结合,开发了L U-S WA T-MO D F-L OW模型,克服了原始HR U在植被动态变化模拟中的缺陷,准确模拟了地表覆被变化及地表水地下水交互过程。结果表明,L U-SWA T-MO D F L OW模型模拟的L A I、E T、地下水位效果较好;研究区植被恢复后,不同子流域月蒸散发增加了01.5 mm,年蒸散发增加了06 mm;地下水补给量平均每月增加约1.2 7 mm,平均每年增加1 4.0 2 mm;研究区地表水-地下水转换关系
3、及水量在植被恢复类型为低覆盖度草地和裸地向林地转化的区域有较明显变化。关键词:植被恢复;L A I;E T;地下水中图分类号:T V 1 1;P 3 4 3.9 文献标志码:A 文章编号:1 0 0 0-7 7 0 9(2 0 2 3)0 5-0 0 2 6-0 4收稿日期:2 0 2 2-0 4-1 7,修回日期:2 0 2 2-0 8-0 9基金项目:国家自然科学基金项目(4 2 1 6 1 0 2 0);青海省科技厅自然科学基金项目(2 0 2 1-Z J-7 0 5)作者简介:翟婧雅(1 9 9 7-),女,硕士研究生,研究方向为水文过程模拟,E-m a i l:j i n g y a
4、 s e a 1 6 3.c o m通讯作者:金鑫(1 9 7 0-),女,博士、教授、博导,研究方向为生态水文,E-m a i l:j i n x 1 3 l z u.e d u.c n1 概况巴音河流域地处青藏高原柴达木盆地东北边缘,是柴达木盆地第四大河流。该流域属于高原荒漠半荒漠干旱性气候,风大沙多、干旱少雨、土地贫瘠,是我国生态环境最为严酷和脆弱的地区之一。巴音河出山后流入人口及工农业生产集中的中下游地区,该区域地表水地下水转换频繁。近2 0年来,随着“退耕还林”等一系列生态恢复措施的实施,巴音河流域植被恢复明显。但由于气候干旱、水资源时空分布不均匀,灌溉是该区域人工恢复植被的水量保证
5、。在此背景下,干旱区恢复植被的需水量、耗水特征等问题备受关注1。干旱区植被恢复的可持续性需更好地理解植被恢复及相应灌溉措施对水文过程的影响。在水文过程复杂区域,多模型耦合更能有效刻画水循环,如将地表水和地下水作为一个整体系统进行模拟的SWAT-MO D F L OW耦合模型更能反映区域水循环特点2。SWAT-MO D F L OW模型的最小计算单元HRU s为虚拟单元,且各HRU自身为集总式单元,同一HRU上的地表特征相同,这使得SWAT模型不能有效体现同一HRU上的土地覆被向不同类型转变及同一HRU上土地覆被部分变化等情景。克服静态HRU的缺点就成为模拟地表覆被变化对水文过程影响的关键。为此
6、,本文以地表水地下水交互频繁且植被恢复明显的巴 音 河 流 域 中 下 游 为 例,将SWAT-MO D F-L OW模型与能够反映流域真实土地覆被变化的动态HRU s相结合,模拟了植被恢复背景下流域水循环过程,研究结果可为干旱区植被恢复的可持续性、水资源合理配置等提供支撑。2 研究方法2.1 基于动态 HRUs 的 SWAT-MODFLOW 耦合模型 参照B A I L E Y R T等2的研究成果,将原始SWAT模型的基本计算单元HRU改进为动态HRU。基于动态HRU s的SWAT-MO D F L OW耦合模 型 运 算 流 程 见 图1。由 于 基 于HRU s的SWAT-MO D F
7、 L OW耦合模型能够反映土地利用/覆被 的 动 态 变 化,将 其 命 名 为L U-SWAT-第4 1卷第5期翟婧雅等:干旱区内陆河流域植被恢复对水文过程的影响开始HRU地表水分量与网络对应MODFLOW地下水模拟网络地下水分量与对应HRU年循环上一年为最后一年?SWAT地表水模拟日循环上一天为本年最后一天?SWAT后续模拟计算并输出结果程序结束是否是当年地表覆被高时空分辨率LAI动态HRUs图1 基于动态H R U s的S WA T-MO D F L OW耦合模型运算流程图F i g.1 T h e c a l c u l a t i o n f l o w c h a r t o f
8、t h e S WA T-MO D F L OW w i t h d y n a m i c H R U sMO D F L OW。2.2 模型建立建立SWAT模型所需的气象数据为德令哈气象站逐日降水、温度、相对湿度、风速及太阳辐射数据。所需农、林、草灌溉数据来自德令哈市水务局(表1)。所需D EM数 据为3 0 m分 辨 率S R TM数据(h t t p:/g d e x.c r.u s g s.g o v/g d e x/)。土壤数据(图2(a)裁剪自全国11 0 0万土壤类型图。图2(b)显示了子流域划分结果。表1 研究区灌溉水量T a b.1 T h e i r r i g a t
9、i o n w a t e r o f t h e s t u灌溉类型年灌溉水量/(m3h m-2)灌溉次数/次灌溉周期/月农业灌溉5 8 0 0631 0林地灌溉5 4 0 0641 1草地灌溉3 6 0 0541 1地下水观测井子流域划分5 km5 km内陆盐土寒漠土栗钙土灰棕漠土绿洲土(a)土壤类型(b)子流域N图2 S WA T模型输入数据F i g.2 I n p u t s d a t a o f t h e S WA T m o d e l 模拟时段内(2 0 0 02 0 1 8年)的土地利用类型数据为基于L a n d s a t 3 0 m遥感影像,其分类精度在9 0%以上
10、。图3为2 0 0 02 0 1 8年研究区植被覆盖情况。由于SWAT模型自带的土地利用/植被数据库中并无枸杞相关参数,本文研究将SWAT相关数据库中灌木的相关参数定义为枸05 km枸杞春小麦林地草地水域城镇裸地05 km春小麦林地草地水域城镇裸地N(a)2000年土地利用/覆被图(b)2018年土地利用/覆被图图3 研究区植被覆盖情况F i g.3 T h e l a n d u s e/c o v e r m a p o f t h e s t u d y a r e a杞的默认参数,利用L A I、E T等数据对相关参数进行校准3,基于SWAT模型的作物生长模块模拟其生长过程。其他土地利
11、用/覆被类型的相关参数采用模型自带数据库中的默认值或经校准得到。对于MO D F L OW模型,将SWAT模型划分的流域边界作为其隔水边界,将流域西部及东部出口作为定流量边界。将SWAT模型提取的河网作 为MO D F L OW的 河 流 边 界。将s t e a d y-s t a t e状态下模拟的水头作为模拟非稳定流的初始水头。此外,将研究区浅层含水层概化为非均质各向异性,将地下水运动概化为二维非稳定流。2.3 模型校准资料及方法将2 0 0 02 0 0 1年设置为模型预热期,2 0 0 22 0 1 8年为模型率定期。采用T P D C(h t t p:/d a-t a.t p d
12、c.a c.c n/)提供的2 0 0 22 0 1 8年3 0 m3 0 m分辨率逐月L A I数据集在HRU s尺度上校准SWAT模型中的植被生长参数。采用T P D C提供的0.1 0.1 分辨率地表蒸散发数据集在子流域尺度对L U-SWAT-MO D F L OW模型模拟的实际蒸散发进行校准。采用流域内逐月观测的地下水位数据对MO D F L OW模型相关参数进行率定,地下水位观测井的位置见图2(b)。文献4调整SWAT模型E T相关参数。选取纳什系数(EN S E)、误差百分数(eP B I A S)、均方根误差与决定系数(R2)评价L U-SWAT-MO D F L OW模型的E
13、T模拟效果。选取决定系数(R2)及绝对误差评价L U-SWAT-MO D F L OW模型的地下水位模拟效果2。3 结果及讨论3.1 研究区植被恢复情况表2为研究区2 0 0 22 0 1 8年主要的植被恢复情况。由表2可知,主要的植被恢复类型为低覆盖度草地、耕地和裸地转化为林地,转化面积分别为2 8 7 7、5 1 8、3 2 1 h m2。相应地,各植被恢复地块的灌溉水量也有较大变化。此外,耕地向林72地的转化主要发生在流域东南部的尕海灌区,低覆盖度草地和裸地向林地的转化主要发生在流域西北部的德令哈灌区。SWAT模型的人工灌溉模块,通过灌溉效率系数、地表径流率等定义进入土壤剖面的水量,能够
14、真实模拟土壤水剖面的漫灌方式3。表2 研究区主要植被恢复情况T a b.2 T h e m a i n v e g e t a t i o n r e s t o r a t i o n o f t h e s t u d y a r e a主要植被恢复类型恢复面积/h m2灌溉水量变化/m3(h m2a)-1低覆盖度草地林地2 8 7 705 4 0 0耕地林地5 1 85 8 0 05 4 0 0裸地林地3 2 105 4 0 03.2 LU-SWAT-MODFLOW 模型适用性评价3.2.1 L A I模拟效果图4为L U-SWAT-MO D F L OW模型模拟的2 0 0 22 0
15、1 8年多年7月(植物生长旺季)平均L A I与遥感反演的多年7月平均L A I的对比情况。由图4可知,L U-SWAT-MO D F L OW模型基于理想叶面积估算模型,在HRU上计算L A I。各个HRU作为同质单元,基于其计算的L A I与遥感反演的基于网格的L A I空间尺度不一致,但二者的空间分布仍有较好的一致性。100LAI010 km100LAI010 kmNN(a)遥感反演LAI(b)LU SWAT MODFLOWLAI-模拟的图4 遥感反演L A I及 L U-S W A T-M O D F L OW模型模拟L A IF i g.4 T h e r e m o t e l y
16、 s e n s e d a n d s i m u l a t e d L A I3.2.2 E T模拟效果图5为L U-SWAT-MO D F L OW模型校准后模拟的E T与遥感反演E T的对比情况。由图5研究区ET/mm302 329188 20.0 10 kmET/mm427.9302.329188.20010 kmN(a)遥感反演的多年平均ET(b)LU SWAT MODFLOWET-模拟的多年平均图5 L U-S WA T-MO D F L OW模型E T模拟结果F i g.5 T h e s i m u l a t e d E T可知,遥感反演的多年平均E T显示出东北部山区高
17、 而 西 南 部 平 原 底 的 空 间 分 布 特 点,而L U-SWAT-MO D F L OW模型模拟的E T并未显示出该特征。这可能是由于采用的基于遥感的实测E T数据在计算过程中并未考虑农、林、草灌溉5。而L U-SWAT-MO D F L OW模型较好考虑了植 被 恢 复 及 相 应 灌 溉 措 施。经 计 算,L U-SWAT-MO D F L OW模型模拟的各子流域月平均E T,其NN S E0.7 2、eP B I A S在-2 0%2 0%之间,R20.7 5。3.2.3 地下水位模拟效果对比研究区地下水位资料有限,观测井1仅有2 0 0 92 0 1 1年 观 测 数 据
18、,观 测 井24仅 有2 0 1 32 0 1 5年观测数据,观测井5仅有2 0 1 42 0 1 5年观测数据。基于有限观测数据,图6为L U-SWAT-MO D F L OW月地下水位模拟效果。由图6可知,模型能较好地刻画研究区地下水位变化,R2均在0.9 5以上,经计算,绝对误差均在0.5 m以内。2 8652 8622 8592 8562 8532 8502 847SWAT-MODFLOWm模拟地下水位/2 8652 8622 8592 8562 8532 8502 847观测地下水位/mR2=0.963 8R2=0.952 9R2=0.967 8R2=0.957 4R2=0.951
19、2图6 L U-S WA T-MO D F L OW月地下水位模拟结果F i g.6 T h e m o n t h l y g r o u n d w a t e r l e v e l s i m u l a t e d b y t h e L U-S WA T-MO D F L OW3.3 植被恢复对水文过程的影响3.3.1 蒸散发量为准确分析植被恢复对水文过程的影响,在L U-SWAT-MO D F L OW模型中设定2个情景:研究区未进行植被恢复,其他类土地利用/覆被按实际情况变化;研究区按实际情况进行植被恢复,其他类土地利用/覆被按实际情况变化。图7为植被恢复对流域E T的影响情况
20、。由图7可知,植被恢复及相应灌溉措施并未改变流域上逐月E T的变化趋势,多数月份蒸发量较无植被恢复情景高01.5 mm。在年际尺度上,植被恢复及相应灌溉措施亦未改变E T的变化趋势,自2 0 0 4年以后,蒸发量明显高于无植被恢复情景。4、1 2、1 3、1 4、2 6、3 3号子流域上的E T值大于1 0 mm,即在这些子流域上E T增大最为明显。研究时段内与研究区土地利用/覆被类型图综合比较后发现这些子流域植被恢复相对明显。3.3.2 地下水补给量图8为植被恢复对流域地下水补给量的影响82水 电 能 源 科 学 2 0 2 3年 第4 1卷第5期翟婧雅等:干旱区内陆河流域植被恢复对水文过程
21、的影响(a)ET无植被恢复及植被恢复情景下月际变化(b)无植被恢复及植被恢复情景下年际变化(c)ET两种情景下的年平均之差010203040506070801 5 9 1 5 9 1 5 9 1 5 9 1 5 9 1 5 9 1 5 9 1 5 9 1 5 9 1 5 9 1 5 9 1 5 9 1 5 9 1 5 9 1 5 9 1 5 9 1 5 9年-月无植被恢复ET植被恢复后ET150200250300350年份无植被恢复ET植被恢复后ETET/mmET/mm200220032004200520062007200820092010201120122013201420152016201
22、72018521183124622510152228301814331373232317269162741211292010010 kmN2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 20152017 20182016图7 植被恢复对流域E T的影响情况F i g.7 T h e i m p a c t s o f t h e v e g e t a t i o n r e s t o r a t i o n o n E T01020304020022003200420052006200720082009201
23、020112012201320142015201620172018年份05010015020022003200420052006200720082009201020112012201320142015201620172018地下水补给量/mm地下水补给量/mm年份无植被恢复植被恢复无植被恢复植被恢复补给量/m.d)3-1(00.060.060.420.422.502.5014.5114.5183.87010 km(a)无植被恢复及植被恢复情景下地下水补给量的月际变化(b)无植被恢复及植被恢复情景下地下水补给量的年际变化(c)植被恢复后的平均地下水补给分布图图8 植被恢复对流域地下水补给量的影响
24、情况F i g.8 T h e i m p a c t s o f t h e v e g e t a t i o n r e s t o r a t i o n o n g r o u n d w a t e r r e c h a r g e情况。由图8可知,植被恢复后,地下水补给量平均每月增加约1.2 7 mm,平均每年增加1 4.0 2 mm。树冠截留和根系吸水将导致土壤水分、地表径流和林地地下水补给减少6。然而,在巴音河中下游,虽然有相当一部分低覆盖草地、农田和裸地被改为林地,但植被恢复后的地下水补给量明显高于无植被恢复时的地下水补给量。可见农业灌区及被灌溉的植被恢复区地下水补给量最
25、高(1 4.5 1 m3/d)。3.3.3 地表水地下水交换量研究区地表水地下水转换频繁,在研究区水文过程中占主导。为探究植被恢复对其影响,将植被恢复和未恢复两个情景下的多年日平均地表水地下水转换水量相减,得到图9。由图94 8262 7656 677.69592.193 778.750.07河流高程/m地下水补给地下水排泄地下水排泄地下水补给地下水补给增量天(m/)3地下水补给减量天(m/)3地下水排泄增量天(m/)3地下水排泄减量天(m/)35.1401.00011.401.0001020 km图9 L U-S WA T-MO D F L OW模拟的多年日平均地表水地下水交换量F i g.
26、9 T h e L U-S WA T-MO D F L OW s i m u l a t e d s u r f a c e w a t e r-g r o u n d w a t e r e x c h a r g e可知,河段处于研究区上游,地表水大量补给地下水。河段处于下游,该河段地下水大量补给地表水。河道既有地表水对地下水的补给,亦有地下水对地表水的补给,该河段处于德令哈灌区内,受农业灌溉、林草灌溉影响较大,地表水地下水转换关系相对复杂。河段所在的区域亦是研究区主要的植被恢复区。河段处有6个网格的地表水地下水转换关系发生了改变。这与该区域植被恢复后,相应地块的灌溉水量发生变化有关。尕海
27、灌区与德令哈灌区均进行了植被恢复,但德令哈灌区地表水地下水转换关系及水量受植被恢复的影响更大,而尕海灌区几乎不受影响。这与植被恢复类型及相应灌溉水量差异有关。4 结论a.L U-SWAT-MO D F L OW模 型 模 拟 的L A I、E T、地下水位效果较好。b.研究区植被恢复后,不同子流域流域月蒸散发增加了01.5 mm,年蒸散发增加了06 mm;地下水补给量平均每月增加约1.2 7 mm,平均每年增加1 4.0 2 mm。c.研究区地表水地下水转换关系及水量在植被恢复类型为低覆盖度草地和裸地向林地转化的区域有较明显变化。参考文献:1 张宝庆,邵蕊,赵西宁,等.大规模植被恢复对黄土高原
28、生态水文过程的影响J.应用基础与工程科学学报,2 0 2 0,2 8(3):5 9 4-6 0 6.(下转第1 0 0页)92开启度的减压阀,建立了减压阀优化布置模型并采用遗传算求解模型。b.算例验证结果表明,本文模型及方法可行、有效,且优化效果显著。参考文献:1 张鹏飞,何芳,麦伟成,等.水务企业编制DMA规划工作的实践与探讨J.给水排水,2 0 1 4,4 0(1):1 5 6-1 6 0.2 GU L L O T T A AUR O R A,B UT L E R D AV I D,C AMP I S ANO A L B E R T O,e t a l.O p t i m a l l o
29、c a t i o n o f v a l v e s t o i m p r o v e e q u i t y i n i n t e r m i t t e n t w a t e r d i s t r i b u t i o n s y s t e m sJ.J o u r n a l o f w a t e r r e s o u r c e s p l a n n i n g a n d m a n a g e m e n t,2 0 2 1,1 4 7(5):0 4 0 2 1 0 1 6.3 S AN T ONA S T A S O G F,D I NA R D O A,C
30、R E A C O E.D u a l t o p o l o g y f o r p a r t i t i o n i n g o f w a t e r d i s t r i b u-t i o n n e t w o r k s c o n s i d e r i n g a c t u a l v a l v e l o c a t i o n sJ.U r b a n w a t e r j o u r n a l,2 0 1 9,1 6(7):4 6 9-4 7 9.4 A R AU J O L S,R AMO S H,C O E LHO S T.P r e s-s u r e
31、 c o n t r o l f o r l e a k a g e m i n i m i s a t i o n i n w a t e r d i s-t r i b u t i o n s y s t e m s m a n a g e m e n tJ.W a t e r r e s o u r c e s m a n a g e m e n t,2 0 0 6,2 0(1):1 3 3-1 4 9.5 唐鹏翔,许仕荣,盛炟.基于粒子群算法的阀门协同泵站控制供水管网压力J.湖南师范大学自然科学学报,2 0 2 0,4 3(3):7 0-7 5.6 俞焱曦.基于泵阀联合调度的供水管网漏
32、损控制及能耗优化D.杭州:浙江大学.7 赫彦文,刘紫阳,李建义,等.基于遗传算法的测试用例自动生成方法综述J.廊坊师范学院学报(自然科学版),2 0 2 0,2 0(4):1 6-2 0.R e s e a r c h o n O p t i m i z i n g P o s i t i o n o f P r e s s u r e R e l i e f V a l v e i n W a t e r S u p p l y N e t w o r k B a s e d o n I n t e l l i g e n t O p t i m i z a t i o n A l g o
33、r i t h mWANG X i n g,Z HONG Y u-f e i,L V M o u(S c h o o l o f E n v i r o n m e n t a l a n d M u n i c i p a l E n g i n e e r i n g,Q i n g d a o U n i v e r s i t y o f T e c h n o l o g y,Q i n g d a o 2 6 6 4 0 0,C h i n a)A b s t r a c t:T o a d d r e s s t h e p r o b l e m o f e x c e s s
34、 p r e s s u r e i n t h e w a t e r s u p p l y n e t w o r k s y s t e m,t h i s p a p e r p r o p o s e d t o a d d p r e s s u r e r e d u c i n g v a l v e s w i t h c e r t a i n o p e n i n g i n t h e s y s t e m.W i t h t h e o b j e c t i v e o f m i n i m i z i n g t h e s q u a r e s u m
35、 o f t h e d i f f e r e n c e b e t w e e n t h e o p e r a t i n g w a t e r p r e s s u r e a t e a c h n o d e a n d t h e m i n i m u m s e r v i c e h e a d a t t h e n o d e,t h e p r e s s u r e r e-d u c t i o n v a l v e o p t i m i z e d p r e s s u r e l e a k a g e c o n t r o l m o d e
36、 l w a s c o n s t r u c t e d b y c o m b i n i n g t h e h y d r a u l i c b a l a n c e c o n s t r a i n t o f t h e p i p e n e t w o r k a n d t h e n o d e p r e s s u r e c o n s t r a i n t.T h e m o d e l w a s s o l v e d b y u s i n g g e n e t i c a l g o r i t h m.T a k i n g t h e m a
37、 i n w a t e r s u p p l y n e t w o r k o f s o u t h e a s t c o a s t a l t o w n s a s a n e x a m p l e,t h e r e s u l t s s h o w t h a t t h e p r o p o s e d m o d e l a n d m e t h o d a r e e f f e c t i v e a n d f e a s i b l e,a n d i t h a s b e t t e r o p t i m a l e f f e c t i v e
38、 n e s s.K e y w o r d s:p i p e n e t w o r k s y s t e m;o v e r p r e s s u r e;c o m p r e s s i o n r e l e a s e v a l v e;m u l t i-o b j e c t i v e;g e n e t i c a l g o r i t h m(上接第2 9页)2 B A I L E Y R T,W I B L E T C,A R A B I M,e t a l.A s-s e s s i n g r e g i o n a l-s c a l e s p a t
39、 i o-t e m p o r a l p a t t e r n s o f g r o u n d w a t e r-s u r f a c e w a t e r i n t e r a c t i o n s u s i n g a c o u-p l e d S WA T-MO D F L OW m o d e lJ.H y d r o l o g i c a l p r o c e s s e s,2 0 1 6,3 0(2 3):4 4 2 0-4 4 3 3.3 傅笛,金鑫,金彦香,等.巴音河中下游农业灌溉对地下水补给量与排泄量的影响J.水电能源科学,2 0 2 1,3 9
40、(1 0):6 3-6 7.4 WH I T E K,CHAU B E Y I.S e n s i t i v i t y a n a l y s i s,c a l-i b r a t i o n,a n d v a l i d a t i o n s f o r a m u l t i s i t e a n d m u l t i v a-r i a b l e SWA T m o d e lJ.J a w r a j o u r n a l o f t h e Am e r-i c a n w a t e r r e s o u r c e s a s s o c i a t i o
41、n,2 0 1 0,4 1(5):1 0 7 7-1 0 8 9.5 MA N,S Z I L AG Y I J,Z HAN G Y,e t a l.C o m p l e-m e n t a r y-r e l a t i o n s h i p-b a s e d m o d e l i n g o f t e r r e s t r i a l e v a p o t r a n s p i r a t i o n a c r o s s C h i n a d u r i n g 1 9 8 22 0 1 2:V a l i d a t i o n s a n d s p a t i
42、o t e m p o r a l a n a l y s e sJ.J o u r n a l o f g e o p h y s i c a l r e s e a r c h a t m o s p h e r e s,2 0 1 9,1 2 4(8):4 3 2 6-4 3 5 1.6 QU B A J A R,AME R M,T A T A R I NOV F,e t a l.P a r t i t i o n i n g e v a p o t r a n s p i r a t i o n a n d i t s l o n g-t e r m e-v o l u t i o n
43、 i n a d r y p i n e f o r e s t u s i n g m e a s u r e m e n t-b a s e d e s t i m a t e s o f s o i l e v a p o r a t i o nJ.A g r i c u l t u r a l a n d f o r e s t m e t e o r o l o g y.2 0 1 9,2 8 1:1 0 7 8 3 1.T h e I m p a c t s o f V e g e t a t i o n R e s t o r a t i o n o n t h e H y d
44、r o l o g i c a l P r o c e s s e s i n A r i d I n l a n d R i v e r B a s i nZ HA I J i n g-y aa,b,J I N X i na,b,c,J I N Y a n-x i a n ga,b,c,F U D ia,b(a.S c h o o l o f t h e G e o g r a p h i c a l S c i e n c e;b.K e y L a b o r a t o r y o f P h y s i c a l G e o g r a p h y a n d E n v i r
45、o n m e n t a l P r o c e s s e s;c.A c a d e m y o f P l a t e a u S c i e n c e a n d S u s t a i n a b i l i t y,Q i n g h a i N o r m a l U n i v e r s i t y,X i n i n g 8 1 0 0 1 6,C h i n a)A b s t r a c t:I n o r d e r t o b e t t e r u n d e r s t a n d t h e i m p a c t o f v e g e t a t i
46、o n r e s t o r a t i o n o n h y d r o l o g i c a l p r o c e s s e s,t h e SWA T-MO D F L OW m o d e l a n d d y n a m i c HR U w a s c o u p l e d,a n d L U-S WA T-MO D F L OW m o d e l w a s d e v e l o p e d t o a c c u r a t e l y s i m-u l a t e t h e d y n a m i c c h a n g e o f v e g e t
47、a t i o n a n d g r o u n d w a t e r i n t h e l o w e r r e a c h e s o f t h e B a y i n R i v e r B a s i n.T h e r e s u l t s r e-v e a l e d t h a t t h e L A I,E T a n d g r o u n d w a t e r l e v e l s i m u l a t e d b y L U-S WA T-MO D F L OW m o d e l a r e b e t t e r;A f t e r v e g e
48、 t a t i o n r e s t o r a t i o n i n t h e s t u d y a r e a,t h e m o n t h l y e v a p o t r a n s p i r a t i o n o f d i f f e r e n t s u b-b a s i n s i n c r e a s e d b y 0-1.5 mm,a n d t h e a n-n u a l e v a p o t r a n s p i r a t i o n i n c r e a s e d b y 0-6 mm;T h e a v e r a g e m
49、 o n t h l y g r o u n d w a t e r r e c h a r g e i n c r e a s e d a b o u t 1.2 7 mm,a n d t h e a v e r a g e a n n u a l g r o u n d w a t e r r e c h a r g e i n c r e a s e d 1 4.0 2 mm.T h e e x c h a n g e o f t h e s u r f a c e w a t e r a n d g r o u n d w a t e r i n t h e s t u d y a r
50、 e a h a v e o b v i o u s c h a n g e s i n t h e a r e a s w h e r e t h e v e g e t a t i o n r e s t o r a t i o n t y p e i s l o w-c o v e r a g e g r a s s l a n d a n d b a r e l a n d c o n v e r t e d t o f o r e s t l a n d.K e y w o r d s:v e g e t a t i o n r e s t o r a t i o n;L A I;E