利用1961—2017年青藏高原109站降水量资料、NCEP全球逐月再分析资料,讨论了雨季期间高原的水汽输送特征。结果表明:高原雨季降水呈显著的年际变化特征,高原雨季降水主模态为南北反向型和全区一致型。气候态高原雨季的水汽输送路径为来自阿拉伯海的偏西风水汽输送,在孟加拉湾附近分为三支水汽输送气流:一支向北输送,自高原南缘输入;一支在南海附近转为偏南风水汽输送,自东南侧输入高原;一支受高原大地形的阻挡作用,转为偏西风水汽输送。在全区一致型降水偏多年,高原主体呈现北部弱辐合、南部辐散的分布形态;在全区一致型降水偏少年,来自阿拉伯海偏南水汽输送在25° N附近转为偏东水汽输送,在高原南缘形成较强的水汽辐合中心。雨季期间高原各边界的水汽收支情况表现为西、南、东边界以水汽输入为主,北边界为水汽输出,南边界水汽输入量最大,西边界次之。
土壤水文过程(水分运移和传热)及其对气候变化的响应是寒区水文学的前沿问题。然而,冻土的存在使得寒区土壤水文过程变得极其复杂。此外,寒区自然环境恶劣,较难获取长时间序列和高分辨率的野外观测资料。近年来,充分利用已有的观测数据,构建寒区土壤水热耦合模型,并开展相应的数值模拟研究,已成为理解寒区土壤水文物理过程,揭示其动力学机制的重要途径。基于寒区土壤水文物理过程和计算流体力学方法,构建了高分辨率、适用于完全饱和状态下的寒区土壤水热耦合模型,且自主研发了相应的数值求解器和软件包。随后,通过一系列完全饱和状态下的验证算例,如经典的一维传热方程解析解、被广泛应用的二维基准测试算例和室内土柱冻结实验等,对已构建的模型进行了系统的检验。模型模拟结果与解析解、基准算例的结果以及实验数据相比,均有较好的一致性,表明该模型较为准确且高效地模拟了寒区土壤在完全饱和状态下的水分运移和传热过程,尤其能够精细刻画冻土水-冰相态变化等关键过程,有望成为研究寒区土壤水文过程的有力工具。
作为大地测量的一种新兴空间技术,合成孔径雷达干涉(synthetic aperture radar interferometry, InSAR)具有全天时、高精度、大范围和速度快的优点,逐渐被应用于多年冻土区地表形变监测中。通过综述多年冻土形变原理及InSAR监测多年冻土形变的应用实例,研究表明:在气候变暖的背景下,多年冻土区地表年际形变以下沉为主,多年冻土上限附近地下冰含量的大小是影响年际形变量的主要因素;活动层内土壤含水量影响着地表季节形变量的大小,不同类型多年冻土区的地表年际形变量和季节形变量存在着较大的差异。研究还表明,不同波长的SAR产品在不同类型多年冻土区的适用性不同,下垫面特征对利用InSAR获取地表形变量有较大影响,L波段的SAR数据在植被覆盖度较好的区域有更好的效果。由于InSAR的失相干问题,加之目前还缺少长时间、多类型、高频率的实测形变结果作为验证和标校数据,获取准确且连续的大范围形变数据较为困难。针对目前寒区研究需求,布设野外长期观测站点,建立适用于不同多年冻土区的地表形变反演算法,构建具有较高精度和较高时空分辨率的地表形变数据集具有重要的实践和科学意义。
调查冰川资源的分布与变化,对区域乃至全球的自然环境与经济社会发展都具有十分重要的意义。基于315景Landsat 8 OLI遥感影像,结合中国第二次冰川编目数据与Google Earth软件,通过人工目视解译等方法调查了2018年中国冰川的分布与变化。结果表明:中国现存冰川53 238条,总面积为(47 174.21±19.93) km2,72%的冰川面积<0.5 km2,规模在1~32 km2的冰川的面积占中国冰川总面积的60%。2008—2018年,中国冰川总面积减少1 393.97 km2,面积变化率为-0.43%?a-1。冰川面积变化率表现出明显的空间差异,面积退缩最快的是冈底斯山,达-1.07%?a-1;最慢的是羌塘高原,为-0.05%?a-1。坡度上,各山系之间的冰川面积变化率差异较为明显。超过70%的山系位于正东和东南方向的冰川面积退缩快,2008—2018年退缩率为-5.0%;正北方向的冰川面积退缩相对缓慢,同时期退缩率为-3.8%。气温和降水变化率差异以及海拔、坡度、坡向等地形差异,共同影响中国冰川的变化。
Landsat-8 OLI因其空间分辨率较高、重复周期适中、高辐射分辨率、高图像获取率(图像质量)的特点,在北极地区大范围冰川流速监测研究中有较大优势。利用2017/2018年格陵兰岛、斯瓦尔巴群岛、北地群岛、法兰士约瑟夫地群岛、德文岛5处北极区域的Landsat-8全色波段数据,采用特征追踪方法提取入海冰川消融期流速。结合MEaSUREs冰川流速数据,分析了198条北极地区入海冰川流速的空间分布特征及其影响因素,同时探究了格陵兰岛Kangerlussuaq冰川流速随时间变化特征。结果表明:与北极其他区域相比,格陵兰岛前缘流速在5~10 m·d-1及10~20 m·d-1的入海冰川在数量上最多,最大流速达到了31.62 m·d-1。而格陵兰岛内部的冰川流速存在差异,北海岸入海冰川平均流速最慢(1.99 m·d-1),东海岸平均流速(6.13 m·d-1)大于西海岸(4.14 m·d-1)。这种流速空间分布差异可能由冰川规模、冰床地势、海流作用、冰盖消融情况等多种因素共同导致。2018年3—10月期间,Kangerlussuaq冰川前缘流速为21.02~22.87 m·d-1,整体流速为10.02~11.39 m·d-1。冰川流速在6—7月和9—10月出现峰值,在8—9月出现低谷,主要缘于冰川融水导致的运动加速和冰川物质平衡变化。
额尔齐斯河发源于中国新疆阿尔泰山南坡,是中国唯一流入北冰洋的河流。随着中俄“冰上丝绸之路”提出,额尔齐斯河通航具有重要意义,监测其河冰物候变化显得尤为必要。本研究基于2001—2018年MOD09GQ第2波段数据,分析了额尔齐斯河布尔津至斋桑泊段每年河冰物候信息,包括开始封河时间、结束开河时间、冰/雪最大覆盖时间、冰期持续时间、封河速率和开河速率。研究结果表明:(1)MOD09GQ的反射率随河冰开河而逐渐降低,经与南湾水文站实际观测的冰期持续时间相比较,其所得冰期持续时间与实测冰期持续时间变化一致,均略有延长趋势,这表明MOD09GQ能有效监测河冰物候信息;(2)开始封河时间与结束开河时间均有提前趋势,最早开始封河比最晚开始封河提前45天,最早结束开河比最晚结束开河时间提前31天,开始封河比结束开河提前时间更多,导致整个冰期持续时间有延长趋势;(3)封河速率有变慢趋势,开河速率有逐渐加快趋势,这意味着若未来该河段通航,每年通航期将延长,将带来更大的经济效益。
为了解北疆积雪中砷(As)和汞(Hg)的空间分布特征及来源,基于2018年1月在北疆地区58个采样点采集的积雪样品,采用原子荧光光谱法测定积雪中溶解性砷和汞的含量,用反距离加权插值法分析空间分布特征,并利用后向轨迹模型探讨了其来源。结果表明:积雪中溶解性砷和汞的浓度分别在0.21~2.69 μg?L-1和5.32~64.09 ng?L-1,均低于地表水I类水质标准限值,说明积雪中砷和汞元素污染较轻。反距离加权插值法分析表明,积雪中砷和汞污染物的空间分布存在差异性,其中准噶尔盆地和天山北坡地区的浓度分布较高。结合后向轨迹聚类分析结果,推断积雪中砷和汞污染物主要来自局地活动,周边国家气团远距离传输所产生的影响较小。研究结果可为评估北疆生态环境质量及区域污染防治提供科学支持。
基于2012年RTK-GPS、2015年三维激光扫描和2018年无人机航测数据,以天山乌鲁木齐河源1号冰川为研究区,分别从物质平衡、面积、末端等方面分析近期冰川变化。结果表明:乌鲁木齐河源1号冰川近年来呈快速消融趋势。2012—2018年冰川面积减少0.07 km2,年平均面积变化率为-0.01 km2·a-1;同期,冰川末端退缩率为6.28 m·a-1,且2015—2018年退缩速率大于2012—2015年;2012—2018年间表面高程下降,物质平衡为-1.13±0.18 m w.e.·a-1,物质损失主要发生在消融区。将2012—2018年间大地测量法冰川物质平衡与传统的花杆/雪坑法观测结果相比较,发现二者较为一致。而2012—2018年间物质平衡减小速率(-0.64 m w.e.·a-1)大于1980—2012年间(-0.47 m w.e.·a-1),表明近期乌鲁木齐河源1号冰川继续呈快速消融趋势。
利用东北地区121个气象站逐日冻土深度、积雪深度、平均气温、地表平均气温及降水量数据,分析了1964—2017年冬半年冻土的变化特征及气象要素对冻土的影响。结果表明:东北地区积雪深度、平均气温、地表平均气温与冻土深度相关系数较高,降水量相关性不大。20世纪60年代平均气温、地表平均气温及负积温最低,最大冻土深度为历年代最深;随着气候变暖,最大冻土深度以6.15 cm?(10a)-1的速率显著减小。冬半年平均最大冻土深度为123 cm,呈显著纬向分布,自辽东半岛向大兴安岭北部递增;随纬度和海拔高度的增加,平均气温和地表平均气温降低,负积温增加,且由北向南地气温差增大。最大冻土深度全区有90%以上的站点减少,减少速率以0.1~10 cm?(10a)-1为主。冻土持续时间随纬度升高而增加,月最大冻土深度和积雪深度最大值分别出现在3月和1月,最大冻土深度的增加要滞后于积雪深度的增加。由于积雪对地温的保温作用,积雪深度较浅时,冻土深度增加较明显,随着积雪深度的增加,冻土深度变化较小,积雪对冻土起到了保温的作用。对于高纬度地区站点,30 cm左右为积雪的保温界限值;对于沿海站点,积雪保温的界限值在5 cm左右;在相同地形下,冻土深度较浅区域积雪的保温值因海拔高度、气候特点而异。最大冻土深度对地表平均气温升温的响应更为显著,地表平均气温和平均气温每升高1 ℃,最大冻土深度将减小8.4 cm和10.6 cm,负积温每减少100 ℃?d,最大冻土深度减少4.9 cm。
为了研究秦岭高山积雪事件变化特点,利用秦岭陕西境内32个国家气象站1980—2016年度共37年高山积雪观测记录,统计分析5个或者5个以上气象站连续积雪日数≥3 d的区域性高山积雪事件、5个或者5个以上气象站连续积雪日数≥20 d的区域性长时间高山积雪事件、5个或者5个以上气象站连续积雪日数≥60 d的区域性稳定高山积雪事件,结果表明:秦岭1980—2016年度共出现区域性高山积雪事件114次,其中包含区域性长时间高山积雪事件29次,区域性稳定高山积雪事件6次。区域性高山积雪事件均发生在冷季(11月—次年4月),其中60%的区域性高山积雪事件发生在冬季(12月—次年2月)。6次区域性稳定高山积雪事件都发生在12月—次年3月。区域性高山积雪事件1980—1989年度最多(44次),其次是1990—1999年度和2000—2009年度(各29次),2010—2016年度最少(12次)。1980—2016年度,区域性高山积雪事件呈现明显减少趋势(通过α=0.01的信度检验),减幅为-0.86次?(10a)-1,区域性高山稳定积雪事件次数与冷季平均气温呈显著的负相关关系(通过α=0.01的信度检验)。
积雪是水文过程的重要环节,基于1979—2017年中国雪深长时间序列数据集、中国区域地面气象要素数据集提供的降水和气温数据,结合DEM数字高程模型等,运用Mann-Kendall检验法、Sen氏坡度法、Pearson相关分析法,分析了雅鲁藏布江流域雪深时空变化分布特征,并对雪深与气象因子(气温、降水)和地形因子(高程、坡度、坡向)的相关性进行了分析。结果表明:1979—2017年雅江流域多年平均雪深为1.95 cm,且以0.02 cm·a-1的速率呈现显著减少趋势;雪深空间分布特征差异性明显,呈现“二高二低”相间分布的特征,高值区为流域西部边缘和东部的山地区域,低值区为中游河谷、流域出口低谷区;气象因子对雪深的变化起决定性作用,其中年平均气温与雪深的相关系数数值为-0.63,二者相关性显著;雪深呈现出随着高程的增加而增加的变化趋势,但最大雪深并非出现在最高海拔处;雪深随坡度的变化呈现“减少—增加—减少”三段式分布规律,且东坡和南坡的雪深厚度高于西坡和北坡的雪深厚度。
利用1951—2016年努尔苏丹和乌鲁木齐的逐日最低气温资料, 以当年7月到翌年6月为统计年, 对比分析了中亚地区的努尔苏丹与乌鲁木齐的寒潮、 强寒潮和超强寒潮过程频数气候变化特征, 并初步分析了冬季北极涛动(AO)、 北大西洋涛动(NAO)与中亚两个城市同期寒潮活动关系。结果显示: 努尔苏丹的年平均寒潮过程频数(15.7次·a-1)是乌鲁木齐(4.1次·a-1)的3.88倍, 强寒潮和超强寒潮过程频数分别是乌鲁木齐的5.91倍和7.55倍; 努尔苏丹的超强寒潮过程频数月分布呈单峰型, 1月最多, 乌鲁木齐的呈双峰型, 峰值分别出现在11月和4月。努尔苏丹的寒潮过程持续日数普遍比乌鲁木齐短。65 a来, 努尔苏丹与乌鲁木齐的春季寒潮过程频数之间的相关关系显著。努尔苏丹和乌鲁木齐的年寒潮过程频数的线性递减率分别为-0.111次·(10a)-1和-0.445次·(10a)-1; 两城市秋、 冬、 春季寒潮、 强寒潮和超强寒潮过程频数线性趋势大多以递减为主, 但是努尔苏丹冬季和乌鲁木齐春季的超强寒潮过程频数均呈递增趋势。努尔苏丹的年超强寒潮过程频数在近5 a最多。努尔苏丹的年寒潮过程频数序列有31 a、 20 a左右的显著周期, 乌鲁木齐的年寒潮过程频数序列有39 a、 8 a左右显著周期。北大西洋涛动(NAO)与努尔苏丹和乌鲁木齐寒潮过程的关联程度比北极涛动(AO)与它们之间的关联程度更为密切; 努尔苏丹的寒潮过程频数受NAO以及AO的影响程度更深。
青藏高原多年冻土区冻融循环过程对地表能量及其分配的影响研究相对较少,青藏高原唐古拉站多年冻土的实测资料,依据10 cm土壤温度划分浅层土壤冻融循环的各个阶段并结合能量闭合率、地表能量各通量等数据探讨浅层土壤冻融循环过程与地气间水热交换过程之间的影响。结果表明:浅层土壤冻融循环过程各阶段均受气候变化的影响,其融化过程起始时间提前同时冻结过程起始时间推后,完全融化阶段持续时间增加,且逐渐接近完全冻结阶段持续时间;在浅层土壤不同冻融状态下,能量闭合率差值较大,其中完全融化阶段能量闭合状况普遍好于完全冻结阶段;净辐射值在完全融化阶段高于完全冻结阶段,净辐射在完全冻结阶段主要转化为感热通量,在完全融化阶段主要转化为潜热通量,地表土壤热通量在完全融化阶段为正值,在完全冻结阶段为负值。
研究发现,通电作用下冻土中的未冻水会发生迁移,这种持续的迁移是一个复杂的物理化学过程,并最终伴随着冻胀的过程。为了探究这一电场作用对冻土的影响,选取冻结冻胀敏感性较高的粉质兰州黄土作为研究对象,分析其在3、4和5 V?cm-1电势梯度作用下的阴阳极变形量、通电前后水分分布规律和电流及电能损耗。结果表明:随着电势梯度的增大,土体阴阳极变形量的差值逐渐增大,且阴极都发生膨胀,阳极发生沉降;随着电势梯度的增大,水分由阳极至阴极的迁移量增大,电流降低的幅度增大,两者变化特征与含水率差值变化特征类似;在5 V?cm-1的电势梯度下,电能总能耗最大,单位含水率能耗最小。
青藏高原地表辐射变化影响多年冻土地表能量收支平衡、热工计算温度场边界条件以及冻土热稳定性。降雨变化是影响地表辐射的重要因素,以青藏高原北麓河地区的地表辐射资料和浅层地表水热监测数据为基础,分析了降雨对多年冻土区地表辐射特征的影响。研究表明:受降雨影响,辐射分量具有明显的日变化和年际变化特征,短波辐射被严重削弱,地面长波辐射削弱程度较低,而太阳长波辐射增强。小雨作用后,太阳短波、地面短波、地面长波辐射日积分量分别减少24.6%、37.9%、4.2%,太阳长波增加了4.3%;中雨作用后,太阳短波、地面短波、地面长波辐射日积分量分别减少32.2%、43.4%、1.7%,太阳长波增加了11.6%;大雨作用后,太阳短波、地面短波、地面长波辐射日积分量分别减少56.3%、65.5%、4.4%,太阳长波增加了10.7%;降雨作用后地表反照率减小,且呈不对称“U”形变化,地表净辐射增加,长期时间尺度上降雨对净辐射影响不明显;降雨对辐射特征的影响程度与降雨强度呈正相关且地表辐射对降雨变化的响应有滞后效应。在暖湿化背景下,冻土热工计算中应考虑降雨对多年冻土区地表能量平衡和水热边界的影响。
为研究水泥土冻胀融沉特性,对南京地区典型粉质黏土进行了不同条件下的冻胀融沉试验。试验控制因素包括水泥掺入比、含水率、龄期、冷端温度和荷载等,各控制因素在包含基准数值的区间内选取代表性条件进行单一变量试验。试验采用自制冻胀融沉仪进行,冻胀融沉仪下部为制冷块,上部设置补水装置,模拟真实土体冻结时水分迁移过程。结果表明:粉质黏土水泥土冻胀率随水泥掺入比增大、龄期增加呈指数规律减小,随含水率增大、冷端温度升高和荷载减小呈线性增大;融沉系数随水泥掺入比变化规律与冻胀率类似,同一条件下融沉系数均大于冻胀率;水泥土完成融沉的时间随解冻温度升高而减少,但最终融沉系数相同;冻胀融沉存在合理水泥掺入比,粉质黏土的合理水泥掺入比为10%,此时冻胀率为1.56%,融沉系数为2.15%。据此提出,冻结法施工时,预先在土体中按合理水泥掺入比掺入水泥,既可有效抑制土体冻胀融沉,又可提高施工效率。
在多年冻土地区,路基工程的修筑将对下伏多年冻土的热力稳定性产生显著的影响。为保护多年冻土,保证路基的长期稳定性,通风管作为一种对流换热类主动冷却降温措施被广泛应用。通风管内的对流换热强度与管内风速密切相关,针对这一问题,采用现场监测与数值模拟相结合的方式开展了通风管内风速特征及影响因素的研究。结果表明:随通风管管径的增加,管内风速呈抛物线型增加,当管径达到0.6 m后,管内风速增加不再明显;通风管外伸长度对管内风速的影响较小,但随环境风速的增加这一影响逐渐加强;管内风速随通风管的埋设高度的增加呈线性增加,当通风管的埋设高度达到2 m后,管内风速受路基高度的影响较小;此外,在两幅路基并行条件下,受上风向路基的遮挡作用,下风向路基管内风速明显降低,以两幅路基管内风速差值不超过0.4 m·s-1为标准,路基高度为3 m时两幅路基最小间距为50 m。
盐渍土相变温度是判断土体中水分冻结与融化、盐分结晶与溶解的重要参数。不同盐分含量相变温度的差异,给盐渍土在降温过程中的水盐迁移过程及变形规律的模拟带来极大的不确定性。通过降温试验,研究了降温过程中氯盐和硫酸盐综合作用盐渍土中水盐相变温度的变化情况。结果表明:全盐量相同时,盐结晶温度随NaCl和Na2SO4比例的不同而不同。随NaCl的加入,在Na+同离子效应的影响下,Na2SO4更容易结晶,但土体的冰和芒硝共晶点温度下降,使得冰含量显著减少,从而降低了孔隙溶液中固相的产生比例,起到抑制Na2SO4盐渍土盐冻胀变形的作用。当土中只含Na2SO4盐时,随Na2SO4浓度的增加,冰和芒硝共晶点的温度先上升而后缓慢下降,二次相变前冰盐的累积量是导致冰和芒硝共晶点产生这种变化的主要原因。盐渍土三相共晶点温度随NaCl含量的增加呈现上升趋势,这是因为随着NaCl的加入,在发生三相共晶前,孔隙溶液发生相变的固相含量减少,从而使孔隙结构对三相共晶点的影响减小。此外,含有NaCl与Na2SO4双组分的盐渍土,水分和盐分可能以单固相、双固相以及三固相状态析出。研究结果可为深入认识盐渍土的相变规律及物理性质提供理论支撑。
水文模拟不确定性长期以来是制约寒区水文发展的瓶颈问题。水稳定同位素示踪为认识冰川流域径流过程提供了重要“指纹”信息,但仍缺乏有效模型将该信息与冰川水文模型耦合,而且同位素信息对冰川流域水文模型不确定性的约束效果也有待检验。将水稳定同位素信息(δ18O)与冰川流域水文模型FLEXG相耦合,实现对冰川流域水稳定同位素和径流过程的耦合建模(FLEXG-iso),并在乌鲁木齐河源1号冰川流域进行模拟检验和径流分割。结果表明:模型不仅对2013—2016年径流过程有良好的模拟效果,还可以重现水稳定同位素、冰川物质平衡等重要过程。利用水稳定同位素这一辅助数据,提高了模型参数的识别能力,减少了模拟过程中各水源的相互妥协效应和不确定性范围。2013—2016年1号冰川断面径流32%~34%来自融雪,48%~51%来自融冰,0%~7%来自地下水,12%~15%来自降雨径流。水稳定同位素对雪和冰川相关中间过程有明显的约束能力,原有模型对融冰贡献的模拟偏高约7%。FLEXG-iso模型的建立有助于推动寒区水文学理论和方法的发展,以及寒区水资源、生态环境保护等相关决策制定。
冻土覆盖率高的小流域的径流形成受温度因素控制明显,普通水文模型不适用,而常规冻土水文模型因需要较多的气象观测要素而难以应用。考虑冻土流域产流机制,利用青藏高原腹地风火山小流域2017—2018年逐日降水、气温、径流观测数据,以降水、气温为输入,径流为输出,基于长短期记忆神经网络(LSTM)建立了适用于小流域尺度的冻土水文模型,并利用2019年观测数据进行验证。模型得益于LSTM特殊的细胞状态和门结构能够学习、反映活动层冻融过程和土壤含水量变化,具有一定的冻土水文学意义,能很好地模拟冻土区径流过程。模型训练期R2、NSE均为0.93,RMSE为0.63 m3·s-1,验证期R2、NSE分别为0.81、0.77,RMSE为0.69 m3·s-1。同时,为了验证模型可靠性,将模型应用于邻近的沱沱河流域,模型训练期(1990—2009年)R2、NSE均为0.73,验证期(2010—2019年)R2、NSE分别为0.66、0.64,模拟结果较好。考虑到未来气候变化,通过模型对风火山流域径流进行了预测:降水每增加10%,年径流增加约12%;气温每升高0.5 ℃,年径流增加约1%;春季融化期、秋季冻结期径流增幅明显,而由于蒸发加剧、活动层加深,径流在8月出现了减少。模型经训练后依靠降水、气温作为输入能较好地模拟、预测青藏高原冻土区小流域径流,为缺少土壤温度、水分等观测数据的冻土小流域径流研究提供了一种简单有效并具有一定物理意义的方法。
利用新疆自动气象站资料、NCEP/NCAR再分析资料、GDAS数据,基于拉格朗日方法的气流轨迹模式(HYSPLITv4.9),对2018年5月16—22日南疆西部特大暴雨过程的水汽源地和输送路径进行了分析。此次HYSPLIT后向轨迹模拟结果表明,暴雨水汽源于巴伦支海、喀拉海、挪威海和地中海,两地的水汽先分别沿西北气流和偏西气流向下游地区输送。水汽轨迹在哈萨克丘陵汇聚后进入北疆,再绕过天山东侧到达罗布泊地区后随低层的偏东急流抵达暴雨区上空。虽然从巴伦支海、喀拉海、挪威海出发的水汽轨迹略多于地中海,但两地的水汽贡献率分别占62%和38%。因此,巴伦支海、喀拉海、挪威海的水汽对于此次暴雨起到明显的增强作用。另外,后向轨迹模拟结果中并未出现从孟加拉湾出发的低层东方路径水汽,说明该区域的水汽并不是构成南疆西部暴雨水汽的必要条件。北极挪威海、巴伦支海、喀拉海地区海表温度异常偏高,有利于水汽源地蒸发增强和上升气流背景场的形成,北极海域的温度异常对于区域气候变化有很好的指示意义。
基于1961—2017年北疆47站夏季逐日降水资料、NOAA海温资料和NCEP/NCAR再分析资料,利用90%分位确定北疆夏季单站降水阈值,得出极端降水量,讨论了春季热带海温与北疆夏季极端降水的关系。结果表明:北疆夏季极端降水和春季热带印度洋(20° S~15° N,50°~110° E)及赤道东太平洋(15° S~15° N,90°~180° W)海温呈正相关,两个关键区春季热带海温异常偏暖时,北疆夏季极端降水偏多,仅春季热带印度洋关键区海温异常偏暖时,北疆西北地区夏季极端降水偏多。当春季热带印度洋和赤道东太平洋关键区海温同时异常偏暖时,200 hPa西风急流轴明显偏南,500 hPa西西伯利亚和中亚地区低值系统活跃,南方路径输送的水汽增加,有利于北疆夏季极端降水的发生;仅春季热带印度洋关键区海温异常偏暖时,200 hPa西风急流强度增强,500 hPa西西伯利亚地区低值系统活跃,配合偏东路径的水汽输送,北疆西北地区夏季极端降水偏多。
利用1960—2017年日降水量资料,采用线性倾向趋势分析、滑动分析和泰森多边形法等,对河西地区多年降水时空变化特征及不同量级降水日数及降水强度的变化趋势进行了研究。结果表明:河西地区年均降水量为99.0 mm,呈现明显的逐年上升趋势,平均倾向率为8.72 mm?(10a)-1,月降水量为单峰分布,5—10月夏秋汛期降水量占年降水量的89.2%,各季节降水量均呈现显著上升趋势;年均降水日数为36.7天,呈现明显的上升趋势,增幅为3.18 d?(10a)-1,降水日数主要分布在夏季,约占总降水日数的54.6%;平均降水强度为2.70 mm?d-1,呈现减弱趋势,变化速率为-0.04 mm?d-1?(10a)-1;零星小雨和小雨降水日数均呈现增加趋势,而二者平均降水强度均为下降趋势,小到中雨降水日数和降水强度呈现增加趋势,中雨及以上的降水变化趋势不明显。
我国西北内陆干旱区水资源匮乏,生态环境脆弱,在全球气候变化和人类活动干扰背景下,采用同位素方法进行精细尺度地表水-地下水交互作用研究是探求当地水循环变化和水资源管理的基本要求。通过测量格尔木河流域河水、地下水样品2019年5月和8月的222Rn浓度和典型断面流量,结果发现:山区河段河水222Rn浓度最高,平均值为948.72 Bq·m-3,指示基岩裂隙水是山区河段重要补给来源;山前冲洪积扇河水222Rn浓度最低,平均值为76.71 Bq·m-3,地下水补给较少;溢出带地区河水222Rn浓度上升至平均676 Bq·m-3,地下水溢出补给河水,向下至细土平原,河水222Rn浓度呈下降趋势。时间变化上,8月与5月相比,河水222Rn浓度下降,表明地下水补给减少。溢出带S1~S2河段河水与地下水交互关系以双向转化为主,基于质量守恒原理计算河水与地下水交互通量,5月和8月累积河水渗漏通量分别为3.87 m3?s-1和0.9 m3?s-1,地下水补给通量分别为0.51 m3?s-1和0.47 m3?s-1,河水渗漏强度大于地下水补给,二者交互通量存在时空差异。
全球气候变暖将加剧水循环,增大洪水风险。阿克苏河流域位于天山南坡,是北半球中纬度典型的高山流域。本流域不仅有暴雨洪水、冰川和积雪融水造成的洪水,而且还有冰川溃决突发洪水。以阿克苏河的两条支流库玛拉克河和托什干河为研究对象,利用块最大值抽样方法(block maximum)和超定量阈值(peak over threshold,POT)抽样方法提取出1958—2011年的洪水事件,其中基于POT方法在两条支流分别提取106次和112次洪水事件,主要集中在6月、7月和8月份。研究表明:阿克苏河的两条支流库玛拉克河和托什干河的年最大洪水强度分别以8.48 m3?s-1和3.40 m3?s-1的速率增加;在洪水发生时间上,以冰川融水补给为主的库玛拉克河,洪水发生时间有提前的趋势,而以降雨和融雪水补给为主的托什干河,洪水发生时间变得更加分散,表现为春季最大洪水提前、秋季最大洪水推后。
玉龙雪山是我国季风海洋型冰川发育最为典型的代表性地区,也是北半球最靠近赤道的现代冰川区。丽江市的部分水源来自于玉龙雪山冰川融水,同时玉龙雪山景区每年吸引着上千万的游客来欣赏现代冰川。旅游活动与水环境息息相关且相互作用,水环境的准确刻画是指导与优化区域旅游活动的基础。以世界著名旅游城市丽江市为例,基于实地采样与数据分析,对2018年丽江市的旅游活动与水环境进行耦合分析。结果表明:丽江市旅游活动与水体污染耦合度极高,旅游活动与水体富营养化和水体毒理的耦合度较低。旅游活动与水环境存在较强的关联,特别是与水环境指标中的高锰酸盐指数,存在强烈的耦合关系。总体上,旅游活动与水环境存在相关关系,旅游旺季水环境质量较旅游淡季更差。
乡镇水权分配是乡镇水资源承载力研究的前提,是乡镇落实最严格的水资源管理和水权交易制度的基础,为新时期乡村振兴与美丽乡村建设提供物质保障。以最严格水资源管理制度确定的甘州区用水总量控制目标为可分配水权总量,在优先保障生活和生态基本用水基础上,选取2017年为规划水平年,采用人口、面积、产值单指标和混合分配模式,以及综合指标分配模式5种方法,开展了乡镇尺度水权分配对比研究。结果表明:领域专家对水权分配的公平性的重视程度略大于效率性,区域水资源的管理者更注重水权分配的公平性,而研究学者更倾向于效率性。5种模式对同一乡镇的水权分配结果范围较为均衡,但不同模式因侧重点各异对具体乡镇的分配比例存在差异。各乡镇不同模式下最大水权分配比例处于10%以上,而最小比例不超过0.3%。基于AHP的综合指标模式的分配结果更为合理,更贴合各乡镇水资源利用现状与经济社会未来发展潜力,其分配结果为各乡镇水资源承载力研究奠定了基础。
地貌学和沉积学证据表明,山东蒙山地质灾害遗迹分布广泛,通常以砾石堆积堤及河谷堆积阶地的形式保留在现代河床之上。其中,峨峪口砾石堤在蒙山具有典型性和代表性。探明峨峪口砾石堤的成因类型和发生历史,对中国东部(105° E以东)第四纪环境变化、气候地貌过程具有重要的科学研究和科普教育价值。2020年10月和2021年2月,笔者团队两次赴蒙山考察,野外工作主要是砾石组构测量,寻找砾石埋压的断木等有机14C测年材料,以及历史洪水调查、地方文献资料收集等。研究结果表明:(1)峨峪口砾石堤不连续、不对称分布于凹岸(弯道),砾石组构、石线构造、叠置空洞构造和砾石群构造等成因专属性特征,以及砾石表面擦痕的明度,揭示其可能为数百年前所形成的一次性历史泥石流堆积垄岗,且被后期山洪改造过;(2)年代学和遥感资料显示,1991年夏季峨峪口等地曾暴发大型山洪泥石流,这一结果与蒙山泥石流的历史记载具有一致性,从证据链上进一步增强了峨峪口砾石堤泥石流成因的确然性;(3)峨峪口砾石堤可能是低频、短历时快速地貌灾害过程,属于暴雨直接类特大型历史泥石流堆积,其物源主要是峨峪口上游谷地的重力崩积物;(4)峨峪口砾石堤缺少冰缘现象的地貌组合证据,不存在全新世早期冰缘地貌的年代学依据,当地不具备发育全新世早期“石河”和“冻土环境”的气候条件;(5)蒙山地区分布有丰富的碎屑物质和比较陡的沟谷纵剖面,暴雨频繁、水分充足,具有泥石流暴发的形成条件,周边居民点、公路和登山步道应加强山洪泥石流地质灾害防治工作。
利用唐古拉站2004—2012年气象观测资料,基于KNN算法,结合机器学习思想,建立了一个气象回归模型,模拟了2005年唐古拉地区表层土壤水热变化趋势,结合实测数据,将模拟值与观测值进行对比,并对模型模拟效果进行了评估。结果表明:KNN模型能够较好地模拟活动层表层土壤水热状况,各层土壤温度的模拟值与观测值的相关系数均在0.99以上,均方根误差在1.25 ℃以内;不同深度土壤水分的模拟值与观测值的相关系数均在0.95以上,均方根误差在0.02 m3?m-3以内。总体上,KNN模型能够对青藏高原多年冻土区唐古拉地区表层土壤水热状况进行较为精确地模拟,该模型对于青藏高原其他地区的适用性有待进一步研究验证。
气候变暖导致的滑雪场营业时间变短和造雪成本升高、大型体育赛事用雪缺乏保障等问题日益突出,储雪是解决这些问题的有效措施之一。然而目前基于高效隔热结构和新材料的场地储雪技术研究稀缺,无法满足日益增长的滑雪赛事与娱乐需求。为筛选储雪的最优技术方案,基于不同的隔热材料和隔热结构,在河北省张家口市崇礼区万龙滑雪场设置了四个建造成本可比的储雪堆,开展储雪试验。通过2019年2月21日—3月1日不同储雪堆隔热层的温度、雪面热流通量等数据的分析,对雪堆的储雪效果及性能进行先期评估。试验结果表明:从主要隔热结构外侧、内部、内侧到积雪表面,温度波动逐渐减小,不同隔热结构的隔热性能差异逐渐凸显;在四个储雪堆中,真空夹层储雪堆隔热效果最好,厚层聚氨酯储雪堆次之,多层中空储雪堆再次之,真空绝热板储雪堆隔热效果最差。根据试验结果,建议储雪使用真空夹层隔热结构以及较厚的隔热材料。