雅鲁藏布江大拐弯地区河流水位日变化时空分异特征
2024-05-13余国安岳蓬胜侯伟鹏
余国安 岳蓬胜 侯伟鹏
摘要：径流（或水位）日变化是反映河流水文情势的重要方面，解析径流（水位）日变化特征有助于明晰河流水文动态过程和规律，揭示径流来源和产汇流机制。以藏东南雅鲁藏布江大拐弯地区的易贡藏布、帕隆藏布、曲宗藏布、拉月曲、金珠曲和白马西路河为对象，基于典型河段2022—2023年小时分辨率实测水位数据，采用数理统计和表征参数分析河流水位日变化时空分异特征，并结合降雨、冰川/积雪数据探究河流水位日变化影响因素。结果表明：研究区河流水位日变幅在汛期多高于非汛期，日水位数据分布在汛期多呈正偏（均值大于中值），汛后多呈负偏（均值小于中值）;除白马西路河外，各河段水位日涨落过程在非汛期相对于汛期有所延迟，且汛后延迟趋势更明显;汛期各河段日水位上涨历时均小于回落历时，非汛期则多相反（拉月曲和金珠曲除外）;白马西路河水位日变化主要受降雨过程扰动，而其他河流水位日变化主要受冰雪消融过程影响。
关键词：水位日变化;实测水位;时空分异;影响因素;雅鲁藏布江
中图分类号：TV11
文献标志码：A
文章编号：1001-6791（2024）02-0274-15
收稿日期：2023-09-25;网络出版日期：2024-01-03
网络出版地址：https：//link.cnki.net/urlid/32.1309.P.20240103.1102.002
基金项目：第二次青藏高原综合科学考察研究项目（2019QZKK0903）;国家自然科学基金资助项目（42371015）
作者简介：余国安（1978—），男，安徽怀宁人，副研究员，博士，主要从事泥沙运动、河流地貌及灾害研究。
E-mail：yuga@igsnrr.ac.cn
通信作者：岳蓬胜，E-mail：yuepengsheng3913@igsnrr.ac.cn
径流（或水位）日变化是反映河流水文情势的重要方面，解析径流（水位）日变化特征有助于明晰河流水文动态过程和规律，揭示径流来源和产汇流机制，进而促进水文模型修正及参数率定［1-2］，推动气候变化影响下河流水文情势演变研究［3］。
狭义上的径流日变化也称径流昼夜波动，由蒸散发、植物蒸腾以及冰雪消融等多种过程驱动，本质上受控于太阳日升日落而产生的日尺度上的辐射和气温变化。自20世纪70年代末以来，相关成果在多个水文科学领域的主流学术期刊报道，内容涉及径流昼夜波动特征［4-6］、影响因素（山坡植被、融水过程、地下水补给过程等［7-9］）及对流域冰雪消融过程的反演［1-2，10］等。径流日变化通常表现为2种模式：① 白天河流流量减小（水位回落）;② 白天河流流量增加（水位上涨）。其中，模式①主要与蒸散发及植物蒸腾作用有关，Graham等［3］根据既往研究成果总结了3个理论假设：饱和楔形假说［11］，认为蒸散发改变土壤非饱和基质电位，导致白天流向河流的水力梯度降低，从而流量减小;滨河截流假说［7］，认为河流的侧向流动由滨河植被控制，白天由于植被蒸腾作用消耗潜水层蓄水，导致进入河道的流量减小;流路迁移假说［12］，认为岸坡植被蒸腾作用导致土壤湿度降低，上游来水的一部分将补给土壤水，导致河流补给径流减小。也有研究认为［9，13］，在径流补给地下水的河流中，由于下午水温高，水的黏度降低、导水率增加，流量在下午更小，其与蒸散作用共同影响径流日变化过程。与模式①不同的是，模式②多出现在冰雪消融集水区，这些区域融冰/融雪对径流的影响显著超过蒸散作用，土壤水分始终处于饱和状态，因此，白天流量得以增加，夜晚融冰/融雪过程减弱或停止，流量则相应减小［6］。
雅鲁藏布江（简称雅江）大拐弯地区位于藏东南，是青藏高原的降水和产流高值区，也是高原气候变化的敏感区。这一区域水能和水量资源十分丰富，生态和景观资源得天独厚，且区位优势显著，国道318线和新219线贯穿其中，是连接西藏和川、滇两省的战略交通廊道，也是未来10～20 a中国多项重大战略工程施工建设和运行区［14］。同时，这一地区也是第二次青藏高原综合科学考察研究和“西南河流源区径流变化和适应性利用”国家自然科学基金重大研究计划关注的典型区域之一［15-16］。因此，深入认识雅江大拐弯地区径流（水位）日变化特征及其影响因素，有利于理解和揭示藏东南地区产汇流机制和演变规律，也有利于区域水能水资源合理开发利用和工程长期安全运行，有重要的科学价值和实践意义。不过，由于区域地形复杂多变，野外监测困难，气象水文数据资料匮乏［17］，尚无有关径流（水位）日变化的研究报道。本文基于区域主要河流典型河段高时间分辨率水位监测数据，解析水位日变化特征和时空分异，并尝试探讨其主控影响因素。
1 研究区概况
藏東南雅江大拐弯地区地处喜马拉雅山东端构造结（图1），北与念青唐古拉山、东与横断山脉交接［18］，地理位置大致为27°N—31°N、92°E—97°E。该区域隶属印度大陆与欧亚大陆碰撞最前缘，地质环境脆弱，山高谷深坡陡，气候差异大，冰雪活动强烈，海洋性冰川广泛分布，是青藏高原海洋性冰川发育中心之一［19］，冰崩、滑坡、泥石流、山洪等自然灾害多发。区域降水充沛、水系密布，除雅江干流外，还发育帕隆藏布、易贡藏布、拉月曲、金珠曲、白马西路河等众多支流，河流补给主要以降雨、冰雪融水和基流为主。
2 数据来源与研究方法
2.1 数据来源
选取雅江大拐弯地区易贡藏布等6条主要河流的典型河段为研究对象（图1，表1），开展高时间分辨率水位和降雨监测（表2）。其中，水位监测采用HOBO U20L-02型自记式压力水位计（量程为0～30.6 m，数据分辨率为0.41 cm，标定精度为±3.0 cm），其外形为两头呈锥形的圆柱体（外径3 cm，长15 cm）。为尽可能减小水流紊动对水位监测结果的影响，水位计预先布置固定于钢管（管内径6 cm，长30 cm）内形成相对静水环境。在枯水季节将水位计（连同外置保护钢管）安装于河床底部基岩或漂石背水面（或石缝）等受水流扰动和直接冲击较小的位置，监测获取小时/半小时分辨率水位数据（本文以水位计在河床的固定布设位置为高程基准，这里的“水位”即河流自由水面相对于水位计布设位置的高程，即相对水深）。雨量监测基于HOBO RG3型雨量计（最大雨强为127 mm/h，数据分辨率为0.2 mm，标定精度为±1%）。各监测点具体位置、监测时间序列及监测频率（时间分辨率）列于表2。研究区冰川数据［20-21］来源于国家青藏高原科学数据中心（http：∥data.tpdc.ac.cn/zhhans）的中国第二次冰川编目数据集（2006—2011年）;归一化植被指数（NDVI）数据来源于资源环境科学数据注册与出版系统（http：∥www.resdc.cn/DOI），数据时间为2018年，空间分辨率为1 km。
2.2 分析方法
根据藏东南地区气候和水文特点，结合2022年3月至2023年2月各河段水位季节变化特征（图2），本文将区域径流过程分为3个阶段，即：汛前（3—5月）、汛期（6—10月）和汛后（11月至次年2月）。
为定量分析区域径流日变化特征，以实测小时分辨率水位数据分析水位日变化过程。图3为易贡藏布监测段2022年汛期和非汛期典型日的水位日变化过程。
根据研究需要采用以下指标，以定量反映水位日变化特征。
（1） 统计参数。日水位标准差（均方差，δ）表示日小时监测水位（Zh）偏离日均水位值（每日24 h水位的均值，Zh）的离散程度。由于δ难以体现均值不等数据系列的离散程度，采用变差系数（Cv）比较不同河段日水位动态特征。偏态系数（Cs）反映每日小时水位数据在日均值两侧偏离对称的程度，Cs>0，数据系列呈右偏分布，Zh多位于Zh以下;Cs<0则反之，且Cs绝对值越大日小时水位偏斜程度越明显。
（2） 水位日变幅（Zamp）。1个水位涨落周期（一般分布在1日及其相邻日）最高水位（Zmax）与最低水位（Zmin）之差的一半，单位为m，即
水位变幅受各河段河道自身断面特征（如河宽和坡降等）影响，不同河段不具直接可比性。为比较各河段水位变幅情况，采用相对水位日变幅（Zramp，即各河段水位日变幅与其年内最大水位日变幅（Zyamp）之比，为量纲一变量）：
（3） 日最高、最低水位出现时间。由于水位日变化（即日尺度上的水位涨落周期）很少正好完整地分布在1日，而多是横跨相邻日（如图3）;且研究区水位日涨落过程总体特征为：最高水位出现时间多在夜间至次日凌晨（甚至上午）;而最低水位出现时间多在中午至傍晚（甚至次日凌晨）。因此，在统计分析日最高/最低水位时间分布时，为保证数据的连续性，作如下处理：① 最高水位出现时间在相邻日分界点（0时）之前（即18—24时），记为-6—0时;② 最低水位出现时间在相邻日分界点（0时）之后（即0—6时），记为24—30时;③ 最高/最低水位出现在其余时段按采集时间值正常表示，不作调整。
（4） 日水位涨落历时。trise反映1个水位上涨过程的持续时间，即从日最低水位上升至日（或次日）最高水位的历时，h;tdecline反映1个水位回落过程的持续时间，即从日最高水位回落至日（或次日）最低水位的历时，h，由实测水位数据获得。
（5） 水位小时变化值。对各河段实测水位数据进行滑动平均处理（窗口为24 h），日水位小时变化值（ΔZi）为平滑后水位之差（后1 h水位与前1 h水位之差）：
式中：ΔZi为1 d中第i时的日水位小时变化值，m;Zi+1和Zi分别为1 d中第i+1和第i时的平滑后水位，m。
为了对比不同河段间的小时水位变化值，采用相对小时水位变化值（ΔZi∧）：
式中：ΔZmax为各河段年最大小时水位变化值的绝对值，m。
3 结果及分析
3.1 日水位分布特征
基于实测数据得到各河段汛前、汛期、汛后日水位标准差与变幅季节变化统计特征（表3），可以看出，各河段两参数具有相似的季节差异特征，汛期相较于非汛期（汛前和汛后）日水位标准差和变幅更大。
各河段相对水位日变幅季节变化（图4）与表3所体现的季节差异吻合。汛期Zramp相对非汛期呈增加趋势，易贡藏布（勒曲藏布）、曲宗藏布、金珠曲3個河段趋势更加明显，帕隆藏布、拉月曲2个河段趋势稍弱，而白马西路河则无明显季节差异。
各河段日水位Cv与Cs季节变化如图5和图6所示。由图5可见，各河段日水位Cv值大致呈“W”型分布，非汛期Cv值高于汛期。尽管各河段汛期也存在波峰，但除白马西路河外，其他河段汛期Cv波峰值明显小于非汛期。不过，这一结果的主要原因在于非汛期和汛期水位整体特征（日水位均值）的差异，非汛期水位低，导致水位的较小日变化（波动）也会产生较大的变差系数。
图6显示，易贡藏布（勒曲藏布）、帕隆藏布及金珠曲3个河段汛期日水位Cs大多为正值，说明日水位数据分布多右偏，Zh大多较当天Zh小;曲宗藏布和拉月曲2个河段汛期Cs正负值出现概率相当但正值绝对值较大，说明2个河段汛期日水位数据分布左偏和右偏天数相当，但右偏数据系列中Zh总体上小于所在日Zh的趋势更显著。汛后各河段Cs值为负值的天数更多，说明汛后各河段Zh多大于当日Zh。白马西路河日水位4—10月无明显季节差异，11月后日水位数据多左偏（Zh多高于所在日Zh）。
3.2 日水位涨落趋势
图7为各河段日小时水位变化，显示日水位涨落过程和季节变化趋势。虽然各河段日水位涨落过程受降雨扰动，但每日水位涨落过程和季节差异仍较为明显（白马西路河除外）。就易贡藏布（勒曲藏布）河段而言（图7（a）），汛前日水位回落时段在凌晨至午后（约4—16时），水位上涨时段在午后至次日凌晨（约16时至次日4时）;汛期日水位涨落过程有提前趋势，即午夜至正午（约0—12时）为水位回落时段，而正午至午夜（12—24时）为水位上涨时段;汛后河段日水位涨落趋势与汛前类似，但日水位上涨/回落开始时间有所延迟。
帕隆藏布、曲宗藏布、拉月曲、金珠曲4个河段日水位涨落过程与易贡藏布类似（图7（b）—图7（e）），非汛期相对于汛期日水位涨落过程有所延迟，且汛后的延迟趋势更明显;但季节变化的总体特征与易贡藏布有所不同：汛期8—9月时，正午至午夜为日水位上涨时段，而汛期其他阶段（6—7月，10月）水位上涨过程多出现在每日午后至次日凌晨。白马西路河日水位漲落过程无明显季节差异，水位回落时段多出现在凌晨至午后。
比较6个河段不同季节日水位涨落过程，9月水位涨落趋势总体上较为杂乱，结合在此区域的雨量监测数据（2022年主汛期（7—8月）降雨较少，而9月降雨较多），初步推测为降雨扰动，这与2022年西藏气候影响评价［22］中“区域夏季降水量为1981年以来历史同期最少;秋季降水量为1981年以来历史同期最多”的结论相吻合。同时，白马西路河与金珠曲的日水位涨落过程更为杂乱，推测出现此现象的原因为两河地处喜马拉雅山南侧，气候更加暖湿，降雨更加充沛（尤其白马西路河），频繁的降雨扰动日水位涨落过程，使其变化趋势无一定规律。
3.3 日最高和最低水位出现时间
图8为各河段日最高、最低水位出现时间散点图。与图7结果相呼应，各河段汛期日最高/最低水位出现时间相对于非汛期有所提前（图8）。就易贡藏布（勒曲藏布）河段（图8（a））来看，其汛期日最高和最低水位出现时间相对集中，分别大致为22时至次日2时和11—15时2个时段;非汛期时间分布则相对分散，日最高水位除分布在23时至次日6时外，也有部分点据分布在10—12时，日最低水位出现时间则大致在14—18时。不同河段日最高/最低水位时间分布特征存在空间差异：易贡藏布、帕隆藏布河段在汛期相对集中，非汛期较为分散;曲宗藏布、拉月曲河段在8月较为集中，其余时段较为分散;金珠曲、白马西路河则在全年都较为分散，无明显季节差异。
图9为各河段年内日最高、最低水位出现时间箱型图。在年尺度上，日高水位出现时间多在凌晨（中值约0—2时）;而最低水位时间多在下午（中值约15—17时），但存在明显空间差异。
3.4 水位涨落历时
表4、图10显示各河段日水位涨落历时季节变化。总体上，各河段日水位上涨历时由汛期至非汛期有增加趋势，回落历时相应呈减小趋势。易贡藏布（勒曲藏布）、帕隆藏布、曲宗藏布汛期日水位上涨历时小于回落历时，而非汛期大致相反;拉月曲、金珠曲全年日水位上涨历时均小于回落历时，且汛期更显著；白马西路河4—10月日水位上涨历时明显小于回落历时，11月则大于回落历时。
4 水位日变化的主要影响因素
河流水位日变化的影响因素主要表现为植被蒸散发［3］、冻融（冰川/积雪融水）过程［5］和降雨，蒸散发与冻融过程本质上受制于太阳辐射引起的光照或气温变化，具有日变化特征，但降雨或其他偶发事件（滑坡、冰崩、泥石流堵江等）则多存在偶然性，白天夜间均有可能发生，一般不具有日尺度变化特征。实际上降雨或其他偶发事件会在一定程度上扰动狭义上（由植物蒸腾/蒸散发及冻融过程引起）的水位日变化。
4.1 冰雪融水过程对区域水位日变化的影响
由前文分析可知，除白马西路河外，其他5条河流水位日变化特征大致相似（图7），非汛期（汛前与汛后）最高水位出现时间多在0—4时，最低水位多在14—18时，即凌晨至下午为水位回落阶段;汛期日水位涨落过程相对非汛期有所提前，最高水位出现时间多在22时至次日2时，最低水位多在11—15时，即午夜至次日正午为水位回落阶段。根据现有资料与国内外河流水位日变化特征的相关研究成果，初步判断研究区典型河段（除白马西路河外）河流日水位变化特征主要影响因素为冻融过程。
冻融过程主导型河流水位日变化主要受温度影响，表现为白天水位上涨而夜晚水位回落。本研究中，除白马西路河外其余各河流上游均存在大面积冰川和积雪（表1）。由于缺乏冰川/积雪融水以及太阳辐射等具体量化数据，本文将1 d中融水补给波动起始时间按当地昼夜变化（即日出时间）考虑。研究区非汛期日出时间大致在8时（7：30—8：30），汛期大致在7时（6：30—7：30），即汛期日出时间较非汛期提前约1 h。由于水位监测位置与融水补给源距离较远，水位涨落过程存在延迟，延迟时间与监测点上游河道长度正相关，与水流流速负相关。以易贡藏布（勒曲藏布监测河段Ⅰ）为例：上文分析已知，其非汛期水位上涨过程大致自14—18时开始，相对于非汛期融水补给变化出现时间延后约6～10 h;而汛期水位上涨过程大多自11—15时开始，相对于汛期融水补给变化的出现时间延后约4～8 h，即非汛期和汛期延后时长（监测点Ⅰ水位上涨起始时间相对于上游融水补给变化起始时间）不同，其主要原因在于非汛期与汛期河流流速的差异。本文采用德卡托地面测速雷达（SVR）电波流速仪及浮漂法实测了易贡藏布河段（监测点Ⅰ）非汛期（2月中旬）和汛期（7月中旬）的河道中央表面流速（表5），并基于实测流速数据对此差异作了初步估算。由于山区河流坡降较大，河床及岸坡由漂石和卵砾石等粗颗粒构成（表2），水流紊动强烈，其断面平均流速与表面流速的比值随流量大小不同而存在差异，中低流量时比值一般约0.7［23］（波动范围为0.5～1.0）［24］，而高流量时比值常小于0.7，洪峰过程可能低至约0.3［23-26］。结合易贡藏布河段情况，折中选取0.7和0.5分别作为其非汛期与汛期表面流速换算为断面平均流速的折算系数;另外，考虑到沿程汇流作用使得断面流速逐渐有所增大，则监测点Ⅰ上游河段沿程平均流速应小于监测点I的断面平均流速，折算系数初略取0.7（表5）。可以看出，非汛期易贡藏布河段融水径流延迟时间约8 h，汛期约6 h，与前文分析结果大致对应。
蒸散发主导型河流日水位变化多表现为白天水位回落而夜晚水位上涨。本研究中，虽然河流非汛期日水位变化特征与蒸散发特征有部分类似，但细节方面仍存在差异。研究区非汛期日出时间大致在8时，汛期在7时，若河流水位日变化由蒸散发控制，则其回落时段应在7—8时后即逐渐开始（蒸散过程主要受河道两侧植被影响，响应较快，不似冰雪融水过程其径流从冰雪融水补给源流动至监测断面存在明显的时间延迟），但实际上各河段日水位回落过程非汛期大致在凌晨至下午（4—16时），而汛期大致发生午夜至次日正午（0—12时），这与蒸散过程控制的日水位回落时间都不吻合。由于研究区数据资料匮乏，各河段蒸散发与太阳辐射未能定量，但从图7显示的河流日水位变化特征仍能推断：研究区蒸散发对各河段河流日水位变化的影响相对于冰雪融水作用居于次要地位。
4.2 降雨对水位日变化的影响
为分析降雨对水位日变化的影响，以易贡藏布（勒曲藏布）河段为例，结合实测降雨数据对水位数据进行筛选，剔除受降雨影响的水位数据：将日降水量>1 mm的日期作为降雨日并剔除;同时，将上述降雨日的后一日数据一并剔除，以消除降雨对后续水位回落过程的影响。图11为易贡藏布未剔除与剔除降雨影响的水位日变化过程对比，可见去除降雨影响后的水位日变化过程更加稳定，季节分异更明显，但整体趋势未改变。图12为易贡藏布河段去除降雨影响前后日最高、最低水位出现时间对比，可以看出，剔除降雨影响后，数据分布更加集中。
从表6也可以看出，剔除受降雨直接影响的水位数据后，易贡藏布河段水位涨落历时差异增大，汛期日水位上涨历时进一步缩短，回落历时相应增加。总体来说，降雨对易贡藏布河段水位日变化过程有一定扰动，但未改变其季节差异的整体趋势。
综上，虽然降雨对研究区易贡藏布（勒曲藏布）、帕隆藏布、曲宗藏布、拉月曲、金珠曲典型河段水位日变化有一定扰动，但各河段水位日变化与冻融过程主导型水位日变化特征更贴合。因此，雅江大拐弯地区墨脱上游五河段水位日变化主要受冻融过程主导，而墨脱下游的白马西路河水位日变化则主要受降雨过程扰动。应该指出，目前对区域河流水位日变化影响因素的探讨还较为定性和初步，更详细深入的定量分析有赖于高时空分辨率降雨、气温、太阳辐射、蒸散发、冰雪动态等气象水文要素的系统监测。
5 结论
本研究基于2022年3月至2023年2月雅江大拐弯地区6条主要河流小时分辨率水位监测数据，解析了区域水位日变化时空分异特征，并初步探讨了其主要影响因素，得到主要结论如下：
（1） 各河段汛期水位日变幅多高于非汛期，日水位数据分布汛期多呈正偏（日小时水位多低于当日水位均值），汛后多呈负偏（日小时水位多高于当日水位均值）。白马西路河水位波动强烈，无明显季节差异。
（2） 非汛期水位日涨落过程（涨/落起始时间）相对于汛期有所延迟，并且汛后延迟趋势更明显。非汛期水位上涨过程多发生在午后至次日凌晨，回落过程多发生在凌晨至午后;而汛期水位上涨过程多发生在正午至午夜，回落过程多发生在午夜至次日正午。非汛期日最高/最低水位出现时间相对于汛期同样有所延后，且汛期日最高/最低水位分布时段更集中。白马西路河无明显季节差异。
（3） 易贡藏布、帕隆藏布、曲宗藏布汛期日水位上涨历时小于回落历时，而非汛期上涨历时大于回落历时;拉月曲、金珠曲全年日水位上涨历时均小于回落历时，且汛期趋势更明显;白马西路河4—10月日水位上涨历时明显小于回落历时，而11月水位上涨历时明显大于回落历时。
（4） 易贡藏布、帕隆藏布、曲宗藏布、拉月曲、金珠曲五河段水位日变化主要受冻融过程主导，而白马西路河水位日变化主要受降雨过程扰动。
參考文献：
［1］KROGH S A，SCAFF L，KIRCHNER J W，et al.Diel streamflow cycles suggest more sensitive snowmelt-driven streamflow to climate change than land surface modeling does［J］.Hydrology and Earth System Sciences，2022，26（13）：3393-3417.
［2］WEBB R W，FASSNACHT S R，GOOSEFF M N.Defining the diurnal pattern of snowmelt using a beta distribution function［J］.Journal of the American Water Resources Association，2017，53（3）：684-696.
［3］GRAHAM C B，BARNARD H R，KAVANAGH K L，et al.Catchment scale controls the temporal connection of transpiration and diel fluctuations in streamflow［J］.Hydrological Processes，2013，27（18）：2541-2556.
［4］GRIBOVSZKI Z，SZIL?GYI J，KALICZ P.Diurnal fluctuations in shallow groundwater levels and streamflow rates and their interpretation：a review［J］.Journal of Hydrology，2010，385（1）：371-383.
［5］MUTZNER R，WEIJS S V，TAROLLI P，et al.Controls on the diurnal streamflow cycles in two subbasins of an alpine headwater catchment［J］.Water Resources Research，2015，51（5）：3403-3418.
［6］LUNDQUIST J D，CAYAN D R.Seasonal and spatial patterns in diurnal cycles in streamow in the Western United States［J］.Journal of Hydrometeorology，2002，3（5）：591-603.
［7］BREN L J.Effects of slope vegetation removal on the diurnal variations of a small mountain stream［J］.Water Resources Research，1997，33（2）：321-331.
［8］CAINE N.Modulation of the diurnal streamflow response by the seasonal snowcover of an alpine basin［J］.Journal of Hydrology，1992，137（1/2/3/4）：245-260.
［9］SCHWAB M，KLAUS J，PFISTER L，et al.Diel discharge cycles explained through viscosity fluctuations in riparian inflow［J］.Water Resources Research，2016，52（11）：8744-8755.
［10］CAUVY-FRAUNI?S，CONDOM T，RABATEL A，et al.Technical note：glacial influence in tropical mountain hydrosystems evidenced by the diurnal cycle in water levels［J］.Hydrology and Earth System Sciences，2013，17（12）：4803-4816.
［11］BURT T P.Diurnal variations in stream discharge and throughflow during a period of low flow［J］.Journal of Hydrology，1979，41（3/4）：291-301.
［12］BOND B J，JONES J A，MOORE G，et al.The zone of vegetation influence on baseflow revealed by diel patterns of streamflow and vegetation water use in a headwater basin［J］.Hydrological Processes，2002，16（8）：1671-1677.
［13］CONSTANTZ J，THOMAS C L，ZELLWEGER G.Influence of diurnal variations in stream temperature on streamflow loss and groundwater recharge［J］.Water Resources Research，1994，30（12）：3253-3264.
［14］余国安，岳蓬胜，张晨笛，等.藏东南地区的河流水文研究：进展与挑战［J］.科学通报，2023，69（3）：394-413.（YU G A，YUE P S，ZHANG C D，et al.River hydrology studies in Southeast Tibet：progress and challenges［J］.Chinese Science Bulletin，2023，69（3）：394-413.（in Chinese））
［15］胡春宏，鄭春苗，王光谦，等.“西南河流源区径流变化和适应性利用”重大研究计划进展综述［J］.水科学进展，2022，33（3）：337-359.（HU C H，ZHENG C M，WANG G Q，et al.Reviews of the major research plan “runoff change and its adaptive management in the source region of major rivers in Southwestern China”［J］.Advances in Water Science，2022，33（3）：337-359.（in Chinese））
［16］徐宗学，班春广，张瑞.雅鲁藏布江流域径流演变规律与归因分析［J］.水科学进展，2022，33（4）：519-530.（XU Z X，BAN C G，ZHANG R.Evolution laws and attribution analysis in the Yarlung Zangbo River basin［J］.Advances in Water Science，2022，33（4）：519-530.（in Chinese））
［17］徐宗学，周祖昊，姜瑶，等.西南河流源区径流量变化规律及其未来演变趋势［J］.水科学进展，2022，33（3）：360-374.（XU Z X，ZHOU Z H，JIANG Y，et al.Variation laws and future evolution trends of runoff in the headwaters region of Southwestern rivers［J］.Advances in Water Science，2022，33（3）：360-374.（in Chinese））
［18］楊逸畴.雅鲁藏布江大拐弯峡谷的地貌特征和成因［J］.地理研究，1982，1（1）：40-48，111-114.（YANG Y C.The topographic features and the origin of the Great Bend Valley at the Yarlung Zangbo River in Tibet［J］.Geographical Research，1982，1（1）：40-48，111-114.（in Chinese））
［19］施雅风，郑本兴，姚檀栋.青藏高原末次冰期最盛时的冰川与环境［J］.冰川冻土，1997，19（2）：97-113.（SHI Y F，ZHENG B X，YAO T D.Glaciers and environments during the last glacial maximum （LGM） on the Tibetan Plateau［J］.Journal of Glaciology and Geocryology，1997，19（2）：97-113.（in Chinese））
［20］刘时银，姚晓军，郭万钦，等.基于第二次冰川编目的中国冰川现状［J］.地理学报，2015，70（1）：3-16.（LIU S Y，YAO X J，GUO W Q，et al.The contemporary glaciers in China based on the Second Chinese Glacier Inventory［J］.Acta Geographica Sinica，2015，70（1）：3-16.（in Chinese））
［21］GUO W Q，LIU S Y，XU J L，et al.The second Chinese glacier inventory：data，methods and results［J］.Journal of Glaciology，2015，61（226）：357-372.
［22］西藏自治区气候中心.西藏气候影响评价（2022年）［R/OL］.拉萨：西藏自治区气象局.（2023-02-24）［2023-08-25］.http：∥xz.cma.gov.cn/zfxxgk_85277/zwgk/qxbg/202302/t20230224_5326504.html.（Tibet Autonomous Region Climate Center.Climate impact assessment of Tibet （2022）［R/OL］.Lhasa：Meteorological Bureau of Tibet Autonomous Region.（2023-02-24）［2023-08-25］.http：∥xz.cma.gov.cn/zfxxgk_85277/zwgk/qxbg/202302/t20230224_5326504.html.（in Chinese））
［23］MARCHAND J P，JARRETT R D，JONES L L.Velocity profile，water-surface slope，and bed-material size for selected streams in Colorado［M］.Lakewood：US Department of the Interior，Geological Survey，1984.
［24］ASANO Y，UCHIDA T.Detailed documentation of dynamic changes in flow depth and surface velocity during a large flood in a steep mountain stream［J］.Journal of Hydrology，2016，541：127-135.
［25］WIBERG P L，SMITH J D.Velocity distribution and bed roughness in high-gradient streams［J］.Water Resources Research，1991，27（5）：825-838.
［26］THOME C R，ZEVENBERGEN L W.Estimating mean velocity in mountain rivers［J］.Journal of Hydraulic Engineering，1985，111（4）：612-624.
Spatiotemporal variation of diurnal cycle of water level in the grand bend
area of the Yarlung Zangbo River
The study is financially supported by the Second Tibetan Plateau Scientific Expedition and Research of China （No.2019QZKK0903） and the National Natural Science Foundation of China （No.42371015）.
YU Guoan1，YUE Pengsheng1，2，HOU Weipeng1，2
（1. Key Laboratory of Water Cycle and Related Land Surface Processes，Institute of Geographic Sciences and Natural Resources Research，
Chinese Academy of Sciences，Beijing 100101，China;
2. University of Chinese Academy of Sciences，Beijing 100049，China）
Abstract：Diurnal streamflow （or water level） cycle is an important indicator of hydrologic regime.Analyzing the diurnal variation of streamflow （water level） can reveal the source of runoff and its formation mechanism，and clarify the hydrological processes/dynamics of the river.The reaches of six important tributaries （namely，Yigong Zangbo，Parlung Zangbo，Quzong Zangbo，Layuequ，Jinzhuqu，and Baima Xilu rivers） in the lower Yarlung Zangbo River located in Southeast Xizang，were collected to investigate the spatiotemporal characteristics of the daily water level changes.Based on the monitored hourly-resolution water level data of the six river reaches during 2022—2023，the diurnal water level characteristics in the study area were studied using mathematical statistics and parameter analyses.Factors （including rainfall，glacier，and snow melting） influencing the diurnal variation of water level were also investigated.The diurnal variation of water level during flood season is mostly higher than that during non-flood season;the distribution of daily water level is mostly positively skewed during the flood season （mean value higher than median value） and is negatively skewed during post-flood season （mean value lower than median value）.Except for the Baima Xilu River，where the daily water level changes have no marked seasonal difference，the start time of daily water level rise/fall during non-flood season is normally later than that during flood season in each analyzed river reach;this delay is more pronounced in the post-flood period.The occurrence time of daily maximum and minimum water level in each reach is more centrally distributed during the flood season，and usually earlier than that during the non-flood season.The duration of daily water level rising （i.e.，the time span for the water level to go from the daily lowest to highest） is shorter than falling （from the daily highest to lowest） during flood season in each reach;however，the trend is almost opposite in non-flood season except for the Layuequ and Jinzhuqu rivers.The diurnal variation of water level of the Baima Xilu River is mainly disturbed by rainfall，while that of other rivers is primarily influenced by snow/ice melting.
Key words：diurnal water level cycles;measured water level;spatiotemporal variation;influencing factor;Yarlung Zangbo River
