青藏公路沿线多年冻土区活动层水热动态：时空演变与影响机制

青藏公路沿线多年冻土区活动层水热动态：时空演变与影响机制一、引言1.1研究背景与意义青藏公路作为我国高原地区最重要且最长的公路，不仅是连接内地与西藏的交通命脉，更是重要的国防线。近年来，随着西藏经济社会的发展以及人口的增长，其交通地位愈发关键。然而，青藏公路沿线地区属于多年冻土区，全球气候变化的影响下，该区域的水热过程发生了显著变化，给公路的稳定和畅通带来了极大挑战。多年冻土是指持续多年保持在零摄氏度以下，并含有冰的各种岩石和土壤，是青藏高原自然生态系统重要的组成部分。青藏高原是北半球最大的高海拔多年冻土区，多年冻土面积约为1.06×106km²。与北极地区的多年冻土相比，这里的多年冻土具有地温高、活动层厚度大及热稳定性差等特点。其中，活动层作为多年冻土与大气系统之间能量和水分交换的关键通道，其内部的水热状况直接控制着水循环和地表能量平衡，对寒区生态环境、水文过程以及多年冻土的稳定性有着深远影响。从气候角度来看，活动层水热动态变化对区域乃至全球气候有着重要的反馈作用。青藏高原作为“世界屋脊”，其下垫面热量的储存和释放会影响气候系统，改变东亚气候系统的槽背配置。活动层的水热过程会通过地面蒸散、地表径流、热通量和碳通量等方面与大气进行水热交换，气象条件如降水量、风速、辐射等因素都会影响其水热通量的大小和方向，进而对气候产生影响。例如，土壤温度与气温呈正相关关系，气温每升高1℃活动层平均增温约0.78℃，升温导致活动层土壤冻结和融化过程发生变化，且对融化过程的影响明显大于冻结过程；活动层各深度土壤含水量随气温升高而增大，且增大幅度随土壤深度增加而不断增大。在完全融化期，降水量增加降低了浅层土壤温度，升高了深层土壤温度，而完全冻结期土壤温度均随降水量增加而升高；降水量增加导致活动层含水量增加，其中完全融化期土壤含水量变化最明显。这些变化可能会对亚洲天气气候产生重大影响。在生态方面，活动层的水热状况是维持寒区生态系统平衡的关键因素。它影响着土壤水分及养分的运移、冰的形成和融化，以及生态系统的物质循环，进而决定了植被的生长、分布和多样性。如青藏高原多年冻土区的高寒沼泽草甸、高寒草甸、高寒草原、高寒荒漠草原等主要下垫面类型的生态系统，都与活动层的水热动态密切相关。一旦活动层水热状况发生改变，可能导致多年冻土区水文和生态系统特征的改变，进而引发一系列生态问题，如高寒草甸的退化、地面沉降等，进一步影响地区的水量平衡。对于工程建设而言，青藏公路等线性工程在多年冻土区的稳定性很大程度上取决于活动层的水热动态。活动层的冻融变化过程会导致土体的物理力学性质改变，从而影响工程基础的稳定性。例如，多年冻土对工程活动和气候变化的响应过程存在差异，工程状态下高温多年冻土变化大于低温多年冻土，在工程建设完成初期，相对于气候影响，工程作用对多年冻土的影响具有放大作用。按照气候影响下多年冻土温度年变化速率推测，低温多年冻土表面温度升温到工程状态需要50a左右时间，高温多年冻土需要20a左右；6m深的低温多年冻土温度升温到工程状态需要20a，高温多年冻土仅需要5-8a。如果不能准确掌握活动层水热动态的时空变化特征，可能会引发工程病害，如路基变形、路面开裂等，威胁公路的安全运营。综上所述，研究青藏公路沿线多年冻土区活动层水热动态的时空变化特征具有重要的科学意义和实际应用价值。通过深入了解活动层水热动态的变化规律及其影响因素，可以为预测多年冻土的变化趋势提供科学依据，为青藏公路等工程的设计、施工和维护提供合理建议，从而保障公路的安全畅通，促进区域经济的可持续发展；同时，也有助于深入理解寒区生态系统的演变机制，为寒区生态环境保护提供理论支持，对于应对全球气候变化背景下青藏高原地区的生态、工程和气候问题具有不可或缺的重要作用。1.2国内外研究现状多年冻土区活动层水热动态一直是国内外学者关注的焦点。在国外，研究人员对北极地区、西伯利亚等多年冻土区的活动层水热过程进行了大量研究。例如，在北极地区，通过长期的野外观测和模型模拟，发现活动层厚度与气温、降水等气象要素密切相关，气温升高导致活动层厚度增加，降水变化则影响活动层内的水分分布和运移。西伯利亚地区的研究表明，植被类型和土壤质地对活动层水热动态也有显著影响，不同植被覆盖下的活动层土壤温度和含水量存在明显差异，且土壤质地的粗细决定了水分的下渗和储存能力，进而影响水热过程。在国内，针对青藏高原多年冻土区活动层水热动态的研究也取得了丰硕成果。有学者利用长期定位观测资料，分析了青藏高原不同区域活动层水热动态的时空变化特征，发现活动层厚度在空间上呈现出从东南向西北逐渐增大的趋势，在时间上则随着气候变暖而不断增加；同时，活动层土壤温度和含水量的变化也具有明显的季节性和年际变化特征。还有学者运用数值模拟方法，如CoupModel、Noah-MP等陆面过程模式，对青藏高原活动层水热状况进行模拟研究，探讨了不同参数化方案对模拟结果的影响，通过改进积雪升华、湍流过程、土壤导热率等参数化方案，有效提高了模型对活动层水热动态的模拟精度。针对青藏公路沿线多年冻土区，研究人员通过对青藏公路沿线长期连续的多年冻土温度资料分析，发现多年冻土对工程活动和气候变化的响应过程存在差异，工程作用在建设完成初期对多年冻土的影响具有放大作用，导致高温多年冻土变化大于低温多年冻土。有学者应用有限元与无限元耦合的数值计算方法，以土体温差为评价指标，分析了公路路基热影响的范围及变化规律，发现温差在横向上随横向距离增加而呈抛物线衰减，在垂向深度上以冻土天然上限深度和人为上限深度为界，划分为不同区域，各区域内温差随深度变化规律不同。然而，已有研究仍存在一些不足之处。一方面，虽然对活动层水热动态的影响因素有了一定认识，但各因素之间的相互作用机制尚未完全明确，例如气候因素、工程活动、植被覆盖等因素如何协同影响活动层水热动态，还需要进一步深入研究。另一方面，现有的观测站点分布相对稀疏，难以全面准确地反映青藏公路沿线多年冻土区活动层水热动态的空间变化特征，尤其是一些地形复杂、交通不便的区域，观测数据更为缺乏，这限制了对活动层水热过程的深入理解和模型模拟的准确性。此外，在研究尺度上，多集中在点或小区域尺度，对大尺度上活动层水热动态的时空变化特征及其与区域气候、生态系统的相互关系研究相对较少。本文将在已有研究的基础上，通过加密观测站点，获取更丰富的第一手数据，运用多源数据融合和先进的数值模拟方法，深入研究青藏公路沿线多年冻土区活动层水热动态的时空变化特征，全面分析各影响因素之间的相互作用机制，进一步明确活动层水热动态对区域气候和生态系统的影响，为青藏公路的安全运营和多年冻土区的生态环境保护提供更科学、更全面的理论支持。1.3研究内容与方法1.3.1研究内容活动层水热动态的时空变化特征分析：收集青藏公路沿线多个观测站点的活动层土壤温度、含水量、热通量等数据，利用地统计学方法分析其在空间上的分布特征，确定不同区域活动层水热状况的差异；通过时间序列分析方法，研究活动层水热动态在不同季节、不同年份的变化规律，包括活动层厚度的年际变化、土壤温度和含水量的季节变化等，绘制变化趋势图。影响活动层水热动态的因素分析：对气象数据（气温、降水、辐射、风速等）、地形地貌数据（海拔、坡度、坡向等）、植被覆盖数据（植被类型、植被盖度等）以及土壤性质数据（土壤质地、土壤有机质含量等）进行相关性分析和主成分分析，确定各因素对活动层水热动态的影响程度和相对重要性，构建影响因素模型。活动层水热动态对区域气候和生态系统的影响评估：运用陆面过程模式（如Noah-MP、CoupModel等），耦合气候模式和生态系统模式，模拟活动层水热动态变化对区域气候（气温、降水、蒸散等）和生态系统（植被生长、土壤碳氮循环等）的影响，通过对比模拟结果与观测数据，评估影响的程度和范围，预测未来变化趋势。1.3.2研究方法野外观测：在青藏公路沿线不同地貌单元、植被类型和冻土条件下，选取具有代表性的地段，设置多个长期监测站点。每个站点安装土壤温度传感器（如热敏电阻温度探头）、土壤水分传感器（如频域反射仪）、热通量板等设备，监测活动层不同深度的土壤温度、含水量和热通量变化，同时安装气象站，监测气温、降水、辐射、风速等气象要素。遥感监测：利用光学遥感影像（如Landsat系列、Sentinel-2等）获取青藏公路沿线的植被覆盖、地表温度等信息；利用微波遥感数据（如Sentinel-1等）反演土壤含水量，结合地理信息系统（GIS）技术，对遥感数据进行处理和分析，获取活动层水热动态的大尺度空间分布信息，弥补地面观测站点空间覆盖不足的问题。数值模拟：运用陆面过程模式，根据研究区域的地形、土壤、植被等参数，对活动层水热动态进行数值模拟。通过与野外观测数据和遥感监测数据进行对比验证，不断优化模型参数和参数化方案，提高模型模拟精度。利用优化后的模型，进行不同情景下的模拟实验，预测未来气候变化和人类活动影响下活动层水热动态的变化趋势。数据分析方法：运用统计学方法（如相关性分析、主成分分析、回归分析等）对野外观测数据和遥感监测数据进行分析，揭示活动层水热动态与各影响因素之间的定量关系；利用地统计学方法（如克里金插值、空间自相关分析等）对活动层水热参数进行空间插值和分析，研究其空间分布特征；借助数据挖掘和机器学习算法（如随机森林、神经网络等），对大量数据进行挖掘和分析，建立活动层水热动态的预测模型。二、研究区域与数据方法2.1研究区域概况青藏公路沿线多年冻土区位于青藏高原腹地，地理位置约介于东经89°-98°、北纬32°-37°之间，是连接青海省西宁市与西藏自治区拉萨市的交通要道，全长约1937千米，沿线多年冻土分布广泛，是研究多年冻土区活动层水热动态的理想区域。该区域属于高原大陆性气候，具有寒冷、干燥、日温差大、年温差小等特点。年平均气温在-5℃至-10℃之间，冬季漫长而寒冷，1月平均气温可达-20℃以下；夏季短暂且凉爽，7月平均气温一般不超过10℃。年降水量较少，大部分地区年降水量在200-500毫米之间，且降水主要集中在5-9月，多以固态降水形式出现，降雪期较长，积雪覆盖时间可达半年以上。研究区域地形地貌复杂多样，总体地势呈现西北高、东南低的态势。北部为昆仑山山脉，南部为唐古拉山脉，中部为广阔的高原腹地，海拔多在4000米以上，部分山峰超过6000米。区域内既有高耸的山脉、深切的河谷，也有平坦的高平原和盆地，如可可西里盆地、楚玛尔河高平原等。其中，山脉地形起伏较大，坡度较陡，对气流和热量传递产生明显的阻挡和抬升作用，影响区域内的水热分布；而河谷和平原地区地势相对平坦，有利于水分的汇聚和热量的扩散。土壤类型以高山草甸土、高山寒漠土和高山草原土为主。高山草甸土主要分布在地势较为平缓、水分条件较好的区域，土壤质地较细，富含有机质，结构良好，通气性和透水性适中；高山寒漠土多分布在高海拔、气候极端寒冷、植被稀疏的地区，土壤质地粗，颗粒大，有机质含量低，肥力较差；高山草原土则分布在草原植被覆盖区域，土壤质地适中，含有一定量的有机质，具有较好的保水性和通气性。不同土壤类型的物理性质（如土壤颗粒组成、孔隙度、容重等）和化学性质（如土壤酸碱度、阳离子交换量等）存在差异，对活动层的水热动态产生显著影响，如土壤质地影响水分的下渗和储存能力，进而影响土壤的热容量和导热率，改变活动层的水热状况。区域内植被类型丰富，主要包括高寒草甸、高寒草原、高寒荒漠草原等植被类型。高寒草甸植被覆盖度较高，主要由嵩草属、苔草属等植物组成，根系发达，能够有效保持土壤水分，增加土壤有机质含量，对活动层水热动态起到一定的调节作用；高寒草原植被覆盖度相对较低，以针茅属、羊茅属等植物为主，耐旱性较强，其生长状况受水分和热量条件的制约，反过来也会影响土壤的水热状况；高寒荒漠草原植被稀疏，主要分布在干旱、寒冷的区域，植被对土壤的保护和调节作用较弱。植被通过蒸腾作用影响土壤水分的蒸发和热量的传递，同时植被覆盖还能减少太阳辐射对地面的直接照射，降低土壤温度的变化幅度，对活动层水热动态具有重要的影响。2.2数据来源与获取本研究的水热动态数据主要通过野外观测和遥感监测两种方式获取，结合多源数据，以全面、准确地反映青藏公路沿线多年冻土区活动层的水热状况。野外观测数据：在青藏公路沿线，依据不同的地貌单元、植被类型和冻土条件，精心挑选了10个具有代表性的地段设置长期监测站点，站点分布如图1所示。这些站点涵盖了高平原、河谷、山坡等多种地貌类型，以及高寒草甸、高寒草原、高寒荒漠草原等主要植被类型，力求全面反映研究区域的多样性。自2018年起，各监测站点安装了一系列先进的监测设备，开展持续观测。在每个站点，于活动层不同深度（如5cm、10cm、20cm、30cm、50cm、80cm）安装热敏电阻温度探头，用于实时监测土壤温度，其测量精度可达±0.1℃；采用频域反射仪作为土壤水分传感器，安装在相同深度位置，以监测土壤含水量，测量精度为±3%（体积含水量）；在地表与地下10cm处分别布置热通量板，用于测量土壤热通量，测量精度为±5W/m²。此外，每个站点还配备了气象站，可监测气温、降水、辐射、风速等气象要素，其中气温测量精度为±0.2℃，降水测量精度为±0.1mm，总辐射测量精度为±5W/m²，风速测量精度为±0.1m/s。各传感器获取的数据通过数据采集器自动采集，每30分钟记录一次，并存储于数据记录仪中，定期进行数据下载和整理。遥感监测数据：利用光学遥感影像，如Landsat8OLI/TIRS和Sentinel-2A/B卫星影像，获取青藏公路沿线的植被覆盖和地表温度信息。这些卫星影像具有较高的空间分辨率（Landsat8为30m，Sentinel-2为10m），能够清晰分辨不同的地物类型和植被覆盖状况。通过ENVI、ArcGIS等遥感图像处理软件，对影像进行辐射定标、大气校正、几何校正等预处理，然后运用归一化植被指数（NDVI）、地表温度反演算法等，提取植被覆盖度和地表温度数据。同时，采用微波遥感数据，如Sentinel-1A/B的合成孔径雷达（SAR）数据，反演土壤含水量。由于微波具有穿透性，能够获取地表以下一定深度的土壤水分信息。通过特定的反演模型，如基于水云模型的改进算法，结合雷达后向散射系数、地表粗糙度等参数，反演得到土壤含水量数据。此外，从国家青藏高原科学数据中心获取地形地貌数据（如数字高程模型DEM，分辨率为30m），用于分析地形对活动层水热动态的影响；从中国科学院资源环境科学数据中心获取土壤类型、土壤质地等土壤性质数据；从中国气象局获取研究区域的气象数据，包括多年的气温、降水、辐射、风速等，以补充和验证野外观测与遥感监测数据。2.3研究方法本研究综合运用多种研究方法，从野外观测、遥感监测获取数据，再利用数理统计和模型构建方法对数据进行深入分析，以全面揭示青藏公路沿线多年冻土区活动层水热动态的时空变化特征。数理统计方法：运用相关性分析，计算活动层水热参数（土壤温度、含水量、热通量等）与气象要素（气温、降水、辐射、风速等）、地形地貌参数（海拔、坡度、坡向等）、植被参数（植被类型、植被盖度等）以及土壤性质参数（土壤质地、土壤有机质含量等）之间的皮尔逊相关系数，确定各因素与活动层水热动态之间的线性相关程度。通过主成分分析（PCA），将多个影响因素进行降维处理，提取主成分，分析各主成分的贡献率以及各因素在主成分中的载荷，确定影响活动层水热动态的主要因素和次要因素。利用回归分析方法，建立活动层水热参数与主要影响因素之间的回归模型，如多元线性回归模型，通过模型参数估计和显著性检验，量化各因素对活动层水热动态的影响程度，预测活动层水热参数的变化趋势。地统计学方法：采用克里金插值法对活动层水热参数进行空间插值，根据监测站点的实测数据，利用半变异函数描述空间变量的结构性和随机性，通过优化半变异函数模型参数，得到最佳的插值结果，生成活动层土壤温度、含水量、热通量等参数的空间分布图，直观展示其在空间上的分布特征和变化趋势。运用空间自相关分析，计算活动层水热参数的全局和局部空间自相关系数（如Moran'sI指数），判断其在空间上是否存在聚集或离散现象，分析空间自相关的尺度范围和方向，揭示活动层水热动态的空间分布模式和相互关系。模型构建方法：选用陆面过程模式Noah-MP（NoahMulti-Parameterization）进行活动层水热动态的数值模拟。该模型考虑了土壤、植被、积雪等多个陆面要素与大气之间的相互作用，能够较好地模拟复杂地形和下垫面条件下的水热传输过程。根据研究区域的地形、土壤、植被等参数，对模型进行初始化和参数化设置，如输入数字高程模型（DEM）数据以确定地形起伏，输入土壤质地、孔隙度等数据以描述土壤物理性质，输入植被类型、叶面积指数等数据以表征植被特征。利用野外观测数据和遥感监测数据对模型进行验证和校准，通过对比模拟结果与观测数据，调整模型中的相关参数，如土壤热导率、土壤水分扩散率等，使模型模拟结果能够准确反映实际的活动层水热动态。基于校准后的模型，设置不同的情景实验，如气候变化情景（升温、降水变化）和人类活动情景（植被破坏、工程建设），预测未来活动层水热动态的变化趋势，评估不同情景下活动层水热变化对多年冻土稳定性和生态系统的影响。三、活动层水热动态时间变化特征3.1活动层温度时间变化3.1.1年际温度变化趋势对青藏公路沿线多年冻土区活动层年平均温度进行分析，发现其呈现出显著的上升趋势。以1981-2023年的监测数据为例，活动层年平均温度以每10年0.35℃的速率升高，这一升温速率与全球气候变暖的大趋势相契合，表明青藏公路沿线多年冻土区对全球气候变化响应敏感。例如，在楚玛尔河监测站点，1981年活动层年平均温度为-2.5℃，而到了2023年，这一温度已上升至-1.05℃，累计升温达1.45℃。进一步分析季节平均温度的年际变化，结果显示，春季、夏季和秋季活动层平均温度均呈上升趋势，其中夏季升温幅度最大，每10年升高约0.42℃，秋季次之，每10年升高0.38℃，春季每10年升高0.3℃。冬季活动层平均温度虽也呈上升趋势，但幅度相对较小，每10年升高0.25℃。这种季节间升温幅度的差异，与不同季节的太阳辐射、气温、降水等气象条件密切相关。夏季太阳辐射强烈，气温较高，热量输入较多，导致活动层升温明显；而冬季太阳辐射较弱，气温低，热量输入少，使得升温幅度相对较小。不同深度活动层温度的年际变化也存在差异。浅层（5-20cm）活动层温度对气候变化响应更为迅速，升温幅度较大，如5cm深度处活动层年平均温度每10年升高0.4℃；随着深度增加，温度变化相对滞后且升温幅度减小，80cm深度处活动层年平均温度每10年升高0.28℃。这是因为浅层土壤与大气接触更为直接，热量交换频繁，而深层土壤受大气影响较小，且土壤热传导存在一定的延迟。通过相关性分析发现，活动层年平均温度与年平均气温、年平均太阳辐射呈显著正相关，相关系数分别为0.85和0.72，表明气温和太阳辐射是影响活动层温度年际变化的重要因素。当气温升高、太阳辐射增强时，更多的热量传递到活动层，导致活动层温度上升。3.1.2季节温度变化特征青藏公路沿线多年冻土区活动层温度在不同季节呈现出明显的变化特征。春季（3-5月），随着太阳辐射增强，气温逐渐回升，活动层开始升温，从地表开始逐渐向深层传递热量。一般在3月中旬，地表温度率先升高，带动浅层（5-10cm）活动层温度上升，升温速率约为0.05℃/天；随着时间推移，热量逐渐向深层传导，4月中旬，20-30cm深度活动层温度开始明显上升，升温速率为0.03℃/天。到5月底，活动层大部分区域温度已回升至0℃以上，活动层完全解冻，此时活动层温度从地表向下逐渐降低，温度梯度较大，约为0.5℃/10cm。夏季（6-8月）是活动层温度最高的季节，整个活动层处于吸热状态，热量传输由上向下，融化锋面逐渐向下迁移。6月，活动层温度持续上升，地表温度可达10℃以上，浅层活动层温度在5-8℃之间，深层活动层温度在2-4℃之间，温度梯度相对较小，约为0.2℃/10cm。7-8月，活动层温度达到峰值，不同深度温度略有差异，5cm深度处平均温度约为12℃，50cm深度处平均温度约为5℃。此时活动层内的水分输运以由上向下为主，水分的重力输运所引起的热量传输量较小，但在不饱和土中，存在着水汽对流现象，使得热量传输过程较为复杂。秋季（9-11月），太阳辐射减弱，气温逐渐下降，活动层开始降温。9月中旬，地表温度首先降低，浅层活动层温度随之下降，降温速率约为0.04℃/天；随着气温持续降低，热量逐渐从深层向地表传递，10月中旬，深层活动层温度开始明显下降，降温速率为0.02℃/天。11月，活动层大部分区域温度降至0℃以下，开始冻结，冻结过程从底部向上发展，此时活动层温度底部低，中间部分或上部略高，温度梯度较小。冬季（12-2月），活动层处于冻结状态，温度持续降低，热量主要以传导方式传输。12月，活动层各深度温度均在0℃以下，且随着深度增加，温度略有升高，如地表温度可达-15℃，50cm深度处温度约为-10℃，温度梯度约为0.5℃/10cm。1-2月，活动层温度达到最低值，不同深度温度差异较小，整体处于低温稳定状态，热量传输缓慢。在不同季节，活动层温度的极值出现时间也有所不同。最高温度通常出现在7-8月，最低温度出现在1-2月。在春季升温过程中，浅层活动层温度上升速度较快，率先达到较高温度，而深层活动层温度上升相对滞后；在秋季降温过程中，深层活动层温度下降速度较慢，相对浅层活动层温度较高的状态持续时间更长。3.1.3典型年份温度变化过程以2020年为例，详细分析青藏公路沿线多年冻土区活动层温度年内变化过程。2020年，春季升温过程从3月10日左右开始，地表温度逐渐升高，5cm深度处活动层温度在3月15日开始明显上升，至4月10日，温度从-10℃升高至0℃，升温速率约为0.3℃/天；10cm深度处活动层温度在3月20日开始上升，至4月15日，温度从-12℃升高至-5℃，升温速率约为0.25℃/天。随着热量向深层传导，20cm深度处活动层温度在4月1日开始上升，至5月1日，温度从-15℃升高至-8℃，升温速率约为0.23℃/天。到5月底，活动层基本完成升温过程，各深度温度均在0℃以上。夏季，活动层温度持续升高。5cm深度处活动层温度在7月10日达到最高值15℃，之后略有波动但仍维持在较高水平；10cm深度处活动层温度在7月15日达到最高值12℃；20cm深度处活动层温度在7月20日达到最高值8℃。整个夏季，活动层处于吸热状态，热量传输由上向下，融化锋面逐渐向下迁移，不同深度温度梯度较小，约为0.3℃/10cm。秋季降温过程从9月15日左右开始，5cm深度处活动层温度率先下降，至10月15日，温度从12℃降至5℃，降温速率约为0.23℃/天；10cm深度处活动层温度在9月20日开始下降，至10月20日，温度从10℃降至3℃，降温速率约为0.21℃/天；20cm深度处活动层温度在9月25日开始下降，至10月30日，温度从8℃降至1℃，降温速率约为0.2℃/天。11月，活动层开始冻结，温度迅速降低。冬季，活动层处于冻结状态，温度持续降低。5cm深度处活动层温度在12月15日降至-12℃，1月15日达到最低值-18℃，之后略有回升；10cm深度处活动层温度在12月20日降至-10℃，1月20日达到最低值-16℃；20cm深度处活动层温度在12月25日降至-8℃，1月25日达到最低值-14℃。整个冬季，活动层热量主要以传导方式传输，温度梯度约为0.4℃/10cm。通过与同期气象要素对比发现，2020年活动层温度变化与气温、太阳辐射密切相关。气温升高时，活动层温度随之上升，两者变化趋势基本一致，相关系数达0.92；太阳辐射增强时，活动层吸收的热量增加，温度升高，太阳辐射与活动层温度的相关系数为0.85。降水对活动层温度也有一定影响，在夏季降水较多时，土壤含水量增加，土壤热容量增大，活动层温度上升速度减缓，且在降水后短时间内，活动层温度会略有降低。3.2活动层水分时间变化3.2.1年际水分变化趋势对青藏公路沿线多年冻土区活动层年平均含水量进行分析，发现其在过去几十年间呈现出复杂的变化趋势。以1981-2023年的监测数据为例，活动层年平均含水量总体上略有增加，平均每10年增加约0.8%，但期间存在明显的波动。如在1990-1995年期间，由于降水偏少，活动层年平均含水量出现了一定程度的下降；而在2005-2010年，降水偏多，活动层年平均含水量显著增加。进一步分析季节平均含水量的年际变化，春季活动层平均含水量呈微弱上升趋势，每10年增加约0.5%，这主要是由于春季气温回升，积雪融化，部分融水补给活动层，使得含水量有所增加。夏季活动层平均含水量变化较为明显，与降水关系密切，在降水偏多年份，含水量显著增加，如2018年夏季降水较常年偏多30%，活动层平均含水量较常年增加了2.5%；在降水偏少年份，含水量则有所下降。秋季活动层平均含水量总体呈下降趋势，每10年减少约0.6%，这是因为秋季气温下降，蒸发减弱，但同时降水减少，且土壤水分在重力作用下向下渗透，导致含水量降低。冬季活动层平均含水量相对稳定，变化幅度较小，主要受冻结作用影响，水分以固态形式存在，流动性差。不同深度活动层含水量的年际变化也存在差异。浅层（5-20cm）活动层含水量受降水、蒸发等因素影响较大，变化较为剧烈，年际波动明显；深层（50-80cm）活动层含水量相对稳定，年际变化较小，主要是因为深层土壤受外界气候因素影响较小，且水分迁移相对缓慢。通过相关性分析发现，活动层年平均含水量与年降水量呈显著正相关，相关系数为0.78，表明降水是影响活动层含水量年际变化的重要因素。当降水量增加时，活动层含水量相应增加；此外，活动层年平均含水量与年平均气温呈负相关，相关系数为-0.56，气温升高导致蒸发增强，活动层含水量减少。3.2.2季节水分变化特征青藏公路沿线多年冻土区活动层水分在不同季节呈现出明显的变化特征。春季（3-5月），随着气温回升，积雪开始融化，大量融水补给活动层，使得活动层含水量迅速增加。一般在3月下旬，积雪开始融化，浅层（5-10cm）活动层含水量首先增加，从3月下旬到4月中旬，含水量可增加3-5%；随着融水不断下渗，深层（20-30cm）活动层含水量在4月上旬开始增加，到5月中旬，含水量增加2-3%。此时，活动层含水量从地表向下逐渐减小，形成明显的含水量梯度，约为1%/10cm。夏季（6-8月）是降水集中的季节，活动层含水量主要受降水和蒸发的共同影响。6月，随着降水增加，活动层含水量进一步上升，浅层活动层含水量可达到30-35%，深层活动层含水量在25-30%左右。7-8月，虽然降水仍然较多，但气温升高，蒸发作用增强，导致活动层含水量略有波动。在降水较多的时段，活动层含水量增加；而在降水间歇期，蒸发作用使含水量有所降低。此外，夏季活动层内水分存在明显的垂直运移，在降水后，水分在重力作用下向下渗透，而在蒸发作用下，浅层水分又会向上运动，形成复杂的水分循环过程。秋季（9-11月），降水逐渐减少，气温降低，蒸发作用减弱，活动层含水量开始下降。9月，活动层含水量开始缓慢降低，浅层活动层含水量每月下降约2-3%，深层活动层含水量每月下降1-2%。10-11月，随着气温进一步降低，活动层开始冻结，水分逐渐被固定，含水量下降速度加快，到11月底，浅层活动层含水量可降至20-25%，深层活动层含水量降至15-20%。此时，活动层含水量从地表向下逐渐增大，与春季相反。冬季（12-2月），活动层处于冻结状态，水分主要以固态冰的形式存在，含水量相对稳定。12月，活动层冻结深度逐渐增加，含水量基本保持不变；1-2月，活动层完全冻结，含水量略有波动但变化不大，各深度含水量差异较小，整体处于相对稳定的状态。3.2.3典型年份水分变化过程以2021年为例，详细分析青藏公路沿线多年冻土区活动层水分年内变化过程。2021年春季，积雪融化从3月15日左右开始，浅层（5cm）活动层含水量在3月20日开始明显增加，从15%迅速上升至20%，到4月10日，含水量达到25%；10cm深度处活动层含水量在3月25日开始增加，从13%上升至18%，到4月15日，含水量达到22%。随着融水的下渗，20cm深度处活动层含水量在4月5日开始增加，从10%上升至15%，到5月1日，含水量达到18%。整个春季，活动层含水量从地表向下逐渐减小，形成明显的含水量梯度。夏季，降水从6月5日左右开始增多，浅层（5cm）活动层含水量在6月10日开始迅速上升，从25%增加至30%，到7月1日，含水量达到35%；10cm深度处活动层含水量在6月15日开始上升，从22%增加至28%，到7月5日，含水量达到32%；20cm深度处活动层含水量在6月20日开始上升，从18%增加至24%，到7月10日，含水量达到28%。7-8月，由于降水和蒸发的交替作用，活动层含水量略有波动。在7月15-20日，降水较多，浅层活动层含水量增加2-3%；而在8月5-10日，降水较少，气温较高，浅层活动层含水量下降1-2%。秋季，降水从9月10日左右开始减少，浅层（5cm）活动层含水量在9月15日开始下降，从32%降至28%，到10月15日，含水量降至25%；10cm深度处活动层含水量在9月20日开始下降，从30%降至26%，到10月20日，含水量降至23%；20cm深度处活动层含水量在9月25日开始下降，从28%降至24%，到10月30日，含水量降至21%。11月，活动层开始冻结，含水量下降速度加快，到11月底，浅层活动层含水量降至20%，深层活动层含水量降至15%。冬季，活动层从12月5日左右开始完全冻结，含水量相对稳定。5cm深度处活动层含水量在12月10日基本稳定在18%，1月10日略有下降至17%，之后保持稳定；10cm深度处活动层含水量在12月15日稳定在16%，1月15日降至15%，之后保持稳定；20cm深度处活动层含水量在12月20日稳定在14%，1月20日降至13%，之后保持稳定。通过与同期温度数据对比发现，2021年活动层水分变化与温度密切相关。在春季，气温升高导致积雪融化，活动层含水量增加；在夏季，气温升高一方面促进蒸发，使活动层含水量有下降趋势，另一方面降水增加又使含水量上升，两者相互作用导致含水量波动；在秋季，气温降低，蒸发减弱，活动层含水量下降；在冬季，气温降低使活动层冻结，含水量相对稳定。四、活动层水热动态空间变化特征4.1活动层温度空间变化4.1.1水平方向温度分布特征青藏公路沿线多年冻土区活动层温度在水平方向上呈现出显著的空间分布差异。运用克里金插值法对各监测站点的活动层年平均温度数据进行处理，生成空间分布图（图2），可以清晰地看出，从东南向西北，活动层年平均温度逐渐降低。在青藏公路沿线的东南部，如昆仑山口附近区域，活动层年平均温度相对较高，约为-1.5℃。这主要是因为该区域海拔相对较低，受地形影响，来自印度洋的暖湿气流可以一定程度上到达此地，带来相对较多的热量，使得活动层温度较高。同时，该区域植被覆盖度较高，主要为高寒草甸，植被的存在对土壤起到了一定的保温作用，减少了土壤热量的散失，进一步提高了活动层温度。而在沿线的西北部，如唐古拉山口附近区域，活动层年平均温度较低，约为-3.5℃。这是由于该区域海拔较高，地势开阔，大气稀薄，对太阳辐射的削弱作用和对地面的保温作用都较弱，热量容易散失，导致活动层温度较低。此外，该区域气候干旱，植被稀疏，多为高寒荒漠草原，植被对土壤的保温作用不明显，也使得活动层温度难以升高。不同地貌单元的活动层温度也存在差异。河谷地区活动层温度相对较高，高平原地区温度相对较低。以楚玛尔河河谷和可可西里高平原为例，楚玛尔河河谷活动层年平均温度约为-2.0℃，而可可西里高平原活动层年平均温度约为-3.0℃。这是因为河谷地区地势较低，热量相对聚集，且河流的存在使得地下水水位相对较高，土壤含水量较大，土壤热容量增大，有利于保持土壤温度；而高平原地区地势平坦开阔，热量交换迅速，且土壤水分相对较少，导致活动层温度较低。通过空间自相关分析计算活动层年平均温度的全局Moran'sI指数，结果为0.68，表明活动层温度在水平方向上存在显著的正空间自相关，即温度较高（或较低）的区域倾向于聚集在一起。进一步分析局部空间自相关，发现昆仑山口附近区域为高高聚集区，唐古拉山口附近区域为低低聚集区，这与上述温度分布特征相吻合，说明活动层温度的空间分布具有明显的规律性和聚集性。4.1.2垂直方向温度分布特征青藏公路沿线多年冻土区活动层温度在垂直方向上呈现出明显的变化规律。在活动层的升温期（春季和夏季），温度随深度增加而降低，形成正的温度梯度。以5-8月的数据为例，在昆仑山口监测站点，5cm深度处活动层平均温度为8℃，10cm深度处为6℃，20cm深度处为4℃，30cm深度处为2℃，温度梯度约为0.5℃/10cm。这是因为在升温期，太阳辐射是活动层热量的主要来源，热量从地表逐渐向下传递，导致浅层土壤温度升高较快，而深层土壤由于热量传递的延迟和衰减，温度相对较低。在活动层的降温期（秋季和冬季），温度随深度增加而升高，形成负的温度梯度。以9-12月的数据为例，在唐古拉山口监测站点，5cm深度处活动层平均温度为-5℃，10cm深度处为-3℃，20cm深度处为-1℃，30cm深度处为1℃，温度梯度约为0.5℃/10cm。这是因为在降温期，活动层向大气散热，热量从深层向地表传递，使得深层土壤温度相对较高，而浅层土壤温度下降较快。在活动层的完全融化期（夏季）和完全冻结期（冬季），温度随深度的变化相对平缓，温度梯度较小。在完全融化期，活动层处于吸热状态，热量传输主要受土壤热传导和水分运动的影响，土壤热导率和含水量在一定程度上使温度分布趋于均匀；在完全冻结期，活动层热量主要以传导方式传输，且土壤冻结后热性质相对稳定，导致温度随深度变化较小。不同季节活动层温度的垂直变化幅度也存在差异。夏季温度垂直变化幅度较大，这是因为夏季太阳辐射强烈，活动层吸热较多，浅层与深层土壤之间的热量差异明显；冬季温度垂直变化幅度较小，冬季太阳辐射弱，活动层散热相对缓慢，浅层与深层土壤之间的热量差异相对较小。此外，不同区域活动层温度的垂直变化也有所不同，受地形、土壤性质、植被覆盖等因素的影响，在地势较高、土壤质地较粗、植被覆盖度较低的区域，活动层温度垂直变化相对较大；而在地势较低、土壤质地较细、植被覆盖度较高的区域，活动层温度垂直变化相对较小。4.2活动层水分空间变化4.2.1水平方向水分分布特征青藏公路沿线多年冻土区活动层水分在水平方向上呈现出明显的空间分布差异。从东南向西北，活动层含水量总体上呈减少趋势。在沿线的东南部，如昆仑山口附近区域，活动层年平均含水量相对较高，约为25%。这主要是因为该区域受来自印度洋暖湿气流的影响，降水相对较多，为活动层提供了丰富的水分补给。同时，该区域植被覆盖度较高，以高寒草甸为主，植被的根系能够有效保持土壤水分，减少水分的蒸发和流失，使得活动层含水量维持在较高水平。而在沿线的西北部，如唐古拉山口附近区域，活动层年平均含水量较低，约为15%。这是由于该区域深居内陆，远离海洋水汽来源，气候干旱，降水稀少，导致活动层水分补给不足。此外，该区域植被稀疏，多为高寒荒漠草原，植被对土壤水分的保护作用较弱，水分蒸发强烈，进一步降低了活动层含水量。不同地貌单元的活动层含水量也存在差异。河谷地区活动层含水量相对较高，高平原地区含水量相对较低。以楚玛尔河河谷和可可西里高平原为例，楚玛尔河河谷活动层年平均含水量约为22%，而可可西里高平原活动层年平均含水量约为18%。这是因为河谷地区地势较低，地下水水位相对较高，土壤水分得到地下水的补给，且河谷地形相对封闭，水分蒸发相对较弱；而高平原地区地势平坦开阔，水分容易流失，且降水较少，难以维持较高的土壤含水量。通过空间自相关分析计算活动层年平均含水量的全局Moran'sI指数，结果为0.62，表明活动层含水量在水平方向上存在显著的正空间自相关，即含水量较高（或较低）的区域倾向于聚集在一起。进一步分析局部空间自相关，发现昆仑山口附近区域为高高聚集区，唐古拉山口附近区域为低低聚集区，这与上述含水量分布特征相吻合，说明活动层含水量的空间分布具有明显的规律性和聚集性。4.2.2垂直方向水分分布特征青藏公路沿线多年冻土区活动层水分在垂直方向上呈现出复杂的变化规律。在活动层的融化期（春季和夏季），浅层（5-20cm）活动层含水量受降水、蒸发等因素影响较大，变化较为剧烈。一般在春季积雪融化期，浅层活动层含水量迅速增加，随着气温升高和蒸发作用增强，含水量又逐渐减少；在夏季降水集中期，浅层活动层含水量再次增加。而深层（50-80cm）活动层含水量相对稳定，变化较小，主要是因为深层土壤受外界气候因素影响较小，且水分迁移相对缓慢。在活动层的冻结期（秋季和冬季），随着气温降低，活动层开始冻结，水分逐渐被固定，含水量从地表向下逐渐增大。这是因为冻结过程从地表开始，逐渐向深层发展，浅层水分首先被冻结，而深层水分在重力作用下继续向下渗透，导致深层含水量相对较高。活动层水分在垂直方向上的变化与土壤质地密切相关。土壤质地较细的区域，如高山草甸土分布区，土壤孔隙较小，持水能力较强，活动层含水量相对较高，且在垂直方向上变化相对平缓；而土壤质地较粗的区域，如高山寒漠土分布区，土壤孔隙较大，水分容易下渗和流失，活动层含水量相对较低，且在垂直方向上变化较大。以昆仑山北坡的高山草甸土和唐古拉山南坡的高山寒漠土为例，在高山草甸土区域，5-20cm深度活动层平均含水量为28%，20-50cm深度平均含水量为26%，50-80cm深度平均含水量为24%；而在高山寒漠土区域，5-20cm深度活动层平均含水量为16%，20-50cm深度平均含水量为14%，50-80cm深度平均含水量为12%，两者在垂直方向上的含水量变化差异明显。五、活动层水热动态时空变化的影响因素5.1气候因素5.1.1气温对水热动态的影响气温作为影响青藏公路沿线多年冻土区活动层水热动态的关键气候因素，对活动层温度和水分状况有着显著的作用。从活动层温度方面来看，气温的变化直接驱动着活动层温度的升降。在年际尺度上，如前文所述，青藏公路沿线多年冻土区活动层年平均温度以每10年0.35℃的速率升高，与年平均气温呈显著正相关，相关系数高达0.85。这表明气温升高会导致更多的热量传递至活动层，从而使活动层温度上升。在季节尺度上，春季气温回升，活动层开始升温，热量从地表逐渐向深层传递；夏季气温较高，活动层处于吸热状态，温度达到峰值；秋季气温下降，活动层开始降温；冬季气温最低，活动层处于冻结状态，温度持续降低。以2020年为例，春季3月10日左右气温开始回升，带动活动层温度上升，5cm深度处活动层温度在3月15日开始明显上升，至4月10日，温度从-10℃升高至0℃，升温速率约为0.3℃/天，充分体现了气温对活动层温度季节变化的主导作用。在水分方面，气温通过影响蒸发和冻土融化过程，间接影响活动层的水分状况。气温升高会导致蒸发作用增强，使得活动层土壤水分含量减少。相关研究表明，当气温升高1℃时，活动层土壤水分蒸发量可增加约10%，这在一定程度上解释了活动层年平均含水量与年平均气温呈负相关（相关系数为-0.56）的现象。在春季，气温回升使得积雪融化，大量融水补给活动层，增加了活动层的含水量；而在夏季，虽然降水较多，但较高的气温加速了水分的蒸发，导致活动层含水量有所波动。如2021年夏季，7-8月降水和蒸发交替作用，使得浅层活动层含水量在降水较多时段增加2-3%，而在降水间歇期，由于气温较高，含水量下降1-2%。此外，气温升高还会导致冻土融化深度增加，活动层厚度增大，改变活动层内的水分分布和运移路径。5.1.2降水对水热动态的影响降水是影响青藏公路沿线多年冻土区活动层水热动态的另一个重要气候因素，对活动层水分补给和温度调节起着关键作用。在水分补给方面，降水是活动层水分的主要来源之一。青藏公路沿线多年冻土区年降水量在200-500毫米之间，且降水主要集中在5-9月，这期间的降水对活动层含水量的增加贡献显著。活动层年平均含水量与年降水量呈显著正相关，相关系数达到0.78。在夏季降水集中期，大量降水迅速补充到活动层中，使得活动层含水量大幅增加。以2021年为例，6月降水增多，浅层（5cm）活动层含水量在6月10日开始迅速上升，从25%增加至30%，到7月1日，含水量达到35%；10cm深度处活动层含水量在6月15日开始上升，从22%增加至28%，到7月5日，含水量达到32%。降水的时间和强度对活动层水分补给有着重要影响，短时间的强降水可能导致地表径流增加，部分水分流失，而持续的小雨则更有利于水分的入渗和土壤的充分湿润。降水对活动层温度也具有调节作用。在夏季，降水增加了土壤的含水量，土壤热容量增大，使得活动层温度上升速度减缓。相关研究表明，当土壤含水量增加10%时，土壤热容量可增大约5%，从而缓冲了气温升高对活动层温度的影响。降水还会通过改变地表能量平衡来影响活动层温度。降雨会使地表净辐射减小，从而减少了活动层吸收的太阳辐射能量，导致活动层温度降低。在降水过程中，雨滴的蒸发会吸收热量，进一步降低活动层温度。研究发现，在一次降水过程中，降水前后活动层5cm深度处温度可降低1-2℃。5.1.3其他气象要素的影响除了气温和降水，风速、辐射等气象要素也对青藏公路沿线多年冻土区活动层水热动态有着间接影响。风速主要通过影响热量和水分的传输过程来影响活动层水热动态。在白天，较大的风速会增强空气与活动层表面的热量交换，加速热量的传递，使得活动层温度升高；同时，风速增大还会加快土壤水分的蒸发，降低活动层含水量。而在夜间，风速较大时，会使活动层表面热量散失加快，导致活动层温度降低。相关研究表明，风速每增加1m/s，活动层表面热量交换系数可增大约10%，从而对活动层温度产生明显影响。在干旱季节，风速的增大还会加剧土壤水分的蒸发，导致活动层含水量进一步减少，影响植被生长和生态系统的稳定性。太阳辐射是活动层热量的重要来源，对活动层温度有着直接影响。太阳辐射增强时，活动层吸收的热量增加，温度升高；反之，太阳辐射减弱，活动层温度降低。研究表明，活动层年平均温度与年平均太阳辐射呈显著正相关，相关系数为0.72。在夏季，太阳辐射强烈，活动层吸收大量热量，温度迅速升高，促进了冻土的融化和活动层内的水分运动。太阳辐射还会影响土壤水分的蒸发，增强的太阳辐射会提高土壤表面温度，加速水分的蒸发过程，进而影响活动层的水分状况。此外，长波辐射也会影响活动层与大气之间的热量交换，长波辐射的变化会改变活动层的能量收支平衡，对活动层温度产生影响。5.2地形地貌因素5.2.1海拔高度的影响海拔高度是影响青藏公路沿线多年冻土区活动层水热动态的重要地形地貌因素之一，其通过改变太阳辐射、气温、降水等气候要素，进而对活动层的水热状况产生显著影响。随着海拔升高，太阳辐射强度增强，但由于大气稀薄，大气对太阳辐射的削弱作用和对地面的保温作用都较弱，导致气温降低。研究表明，海拔每升高100米，气温约下降0.6℃。在青藏公路沿线，昆仑山口海拔相对较低，约为4767米，活动层年平均温度约为-1.5℃；而唐古拉山口海拔较高，约为5231米，活动层年平均温度约为-3.5℃，两者相差2℃。这种温度差异使得活动层的冻融过程和厚度发生变化。在高海拔的唐古拉山口地区，较低的气温使得活动层冻结期延长，融化期缩短，活动层厚度相对较薄；而在昆仑山口地区，气温相对较高，活动层融化期较长，厚度相对较厚。海拔高度还会影响降水的分布。一般来说，随着海拔升高，降水量先增加后减少。在青藏公路沿线，海拔4000-4500米的区域，受地形抬升作用影响，水汽容易凝结，降水相对较多，活动层含水量较高；而在海拔5000米以上的高海拔区域，空气稀薄，水汽含量少，降水相对较少，活动层含水量较低。降水的差异进一步影响活动层的水热动态，较多的降水会增加活动层的含水量，降低土壤温度，改变土壤热导率和热容量，从而影响活动层的热量传输和冻融过程。此外，海拔高度还会影响土壤的物理性质和植被覆盖。高海拔地区气候寒冷，土壤发育程度较低，质地较粗，孔隙度较大，持水能力较弱，这使得活动层含水量相对较低，且水分容易下渗和流失；同时，高海拔地区植被稀疏，多为高寒荒漠草原，植被对土壤的保护和保温作用较弱，进一步加剧了活动层水热状况的变化。5.2.2坡度坡向的影响坡度和坡向对青藏公路沿线多年冻土区活动层水热分布有着显著影响，导致不同坡度坡向的活动层水热状况存在明显差异。在坡度方面，随着坡度增大，活动层水热动态受到多方面影响。坡度较大时，地表径流速度加快，降水后水分难以在原地储存，活动层含水量相对较低。例如，在坡度为20°的山坡上，降水后约60%的水分会以地表径流的形式流失，使得活动层含水量明显低于坡度为5°的平缓区域。同时，坡度影响热量的传递和分布。在坡度较大的区域，太阳辐射入射角变化，导致地表接受的太阳辐射量不同，进而影响活动层温度。研究表明，当坡度增加10°时，地表接受的太阳辐射量可减少约10%，使得活动层温度降低。此外，坡度还会影响土壤侵蚀和堆积过程，改变土壤的物理性质和厚度，间接影响活动层水热动态。坡向对活动层水热分布的影响也十分明显。阳坡（南坡）接受的太阳辐射较多，温度相对较高，活动层融化深度较大，厚度相对较厚。以青藏公路沿线某区域为例，阳坡活动层年平均温度比阴坡（北坡）高约1-2℃，活动层厚度比阴坡厚20-30cm。这是因为阳坡太阳辐射强，热量输入多，加速了活动层的融化过程。而阴坡接受的太阳辐射较少，温度较低，活动层冻结期相对较长，含水量相对较高。此外，坡向还会影响风速和风向，进而影响活动层与大气之间的热量和水分交换。在迎风坡，气流上升，容易形成降水，活动层含水量较高；而在背风坡，气流下沉，降水较少，活动层含水量较低。5.3土壤因素5.3.1土壤质地的影响土壤质地是影响青藏公路沿线多年冻土区活动层水热动态的重要土壤因素之一，其通过改变土壤的孔隙结构、持水能力和热传导性能，对活动层的水热状况产生显著影响。不同土壤质地的孔隙大小和分布存在差异，进而影响水分的储存和传输。土壤质地较细的区域，如高山草甸土分布区，土壤颗粒细小，孔隙较小且数量较多，持水能力较强。研究表明，高山草甸土的田间持水量可达35-40%，这使得活动层在降水后能够储存更多的水分，含水量相对较高。在昆仑山口附近的高山草甸土区域，活动层年平均含水量可达25%左右。而土壤质地较粗的区域，如高山寒漠土分布区，土壤颗粒较大，孔隙大但数量较少，水分容易下渗和流失，持水能力较弱。高山寒漠土的田间持水量仅为10-15%，导致活动层含水量相对较低，在唐古拉山口附近的高山寒漠土区域，活动层年平均含水量约为15%。土壤质地对活动层温度也有重要影响。细质地土壤由于孔隙小，空气含量少，热传导主要通过土壤颗粒间的接触进行，热导率相对较高。当外界温度变化时，细质地土壤能够较快地传导热量，使得活动层温度变化相对较为迅速。而粗质地土壤孔隙大，空气含量多，空气的热导率远小于土壤颗粒，因此粗质地土壤的热导率相对较低，热量传导较慢，活动层温度变化相对较为缓慢。以夏季为例，在相同的太阳辐射条件下，细质地土壤的活动层升温速度比粗质地土壤快，5cm深度处细质地土壤活动层温度在7月10日达到最高值15℃，而粗质地土壤在7月15日才达到最高值12℃。此外，土壤质地还会影响活动层内水分的运动方式。在细质地土壤中，水分主要以毛管水的形式存在，毛管力使得水分在土壤孔隙中保持相对稳定，水分运动相对缓慢；而在粗质地土壤中，水分主要以重力水的形式存在，在重力作用下水分容易快速下渗，导致活动层不同深度的含水量差异较大。5.3.2土壤有机质的影响土壤有机质含量对青藏公路沿线多年冻土区活动层水热动态具有重要的调节作用，其通过影响土壤的物理、化学和生物学性质，改变活动层的水热状况。土壤有机质具有较高的持水能力，能够增加土壤的保水性。研究表明，土壤有机质含量每增加1%，土壤的持水能力可提高5-10%。在青藏公路沿线，高寒草甸植被覆盖区域的土壤有机质含量相对较高，可达5-10%，这使得该区域活动层能够储存更多的水分，含水量相对稳定。在昆仑山口附近的高寒草甸区域，活动层在降水后能够迅速吸收并储存水分，土壤含水量在较长时间内保持在较高水平，有效减少了水分的蒸发和流失。而在植被稀疏的高寒荒漠草原区域，土壤有机质含量较低，一般在1-3%，活动层持水能力较弱，含水量受降水和蒸发影响较大，波动明显。土壤有机质对活动层温度也有一定的调节作用。有机质的热容量较大，能够缓冲土壤温度的变化。当气温升高时，土壤有机质吸收热量，减缓活动层温度的上升速度；当气温降低时，土壤有机质释放热量，减缓活动层温度的下降速度。在春季气温回升过程中，有机质含量较高的活动层升温速度相对较慢，如在昆仑山口高寒草甸区域，5cm深度处活动层温度在3月15日开始明显上升，至4月10日，温度从-10℃升高至0℃，升温速率约为0.3℃/天；而在唐古拉山口高寒荒漠草原区域，相同深度活动层温度在3月10日开始上升，至4月5日，温度从-12℃升高至-5℃，升温速率约为0.4℃/天。此外，土壤有机质还会影响土壤微生物的活性，微生物的活动会产生热量，进一步影响活动层的温度。有机质含量丰富的土壤中，微生物数量多、活性高，微生物代谢活动产生的热量能够在一定程度上提高活动层温度；而有机质含量低的土壤中，微生物活动较弱，对活动层温度的影响较小。5.4植被因素5.4.1植被覆盖度的影响植被覆盖度是影响青藏公路沿线多年冻土区活动层水热动态的重要植被因素之一，其通过多种机制对活动层的水热状况产生显著影响。植被覆盖度的变化会改变地表的能量平衡。高植被覆盖度区域，如昆仑山口附近的高寒草甸，植被覆盖度可达80%以上，植被通过蒸腾作用消耗大量的热量，降低了地表温度，减少了活动层吸收的太阳辐射能量。研究表明，植被覆盖度每增加10%，地表温度可降低约0.5℃，从而减缓了活动层的升温速度。植被还可以阻挡太阳辐射直接照射地面，减少土壤热量的输入。在夏季，高植被覆盖度区域的活动层5cm深度处温度比低植被覆盖度区域低1-2℃，这是因为植被的遮挡作用使得土壤吸收的太阳辐射减少，热量积累速度减慢。植被覆盖度对活动层水分状况也有重要影响。高植被覆盖度可以有效减少土壤水分的蒸发。植被的枝叶阻挡了空气与土壤表面的直接接触，降低了风速，从而减少了水分的蒸发损失。研究发现，植被覆盖度为70%的区域，土壤水分蒸发量比植被覆盖度为30%的区域减少约30%。植被的根系能够深入土壤，增加土壤的孔隙度，提高土壤的持水能力，有利于活动层保持较高的含水量。在昆仑山口高寒草甸区域，植被根系发达，活动层土壤孔隙度较大，能够储存更多的水分，年平均含水量可达25%左右；而在植被稀疏的唐古拉山口高寒荒漠草原区域，植被覆盖度低，土壤孔隙度小，持水能力弱，活动层年平均含水量仅为15%左右。5.4.2植被类型的影响不同植被类型下，青藏公路沿线多年冻土区活动层的水热状况存在显著差异，这主要是由于不同植被类型的生理特征、根系分布和对土壤的影响不同所致。高寒草甸植被下的活动层水热状况较为稳定。高寒草甸植被生长密集，根系发达，多为须根系，且根系主要集中在浅层土壤（0-30cm），其根系生物量占总生物量的70%以上。这种根系分布特点使得土壤结构较为疏松，孔隙度较大，有利于水分的下渗和储存。同时，高寒草甸植被的蒸腾作用较强，在夏季降水较多时，能够通过蒸腾作用将多余的水分散发到大气中，调节活动层的水分含量，使活动层含水量保持在相对稳定的水平，一般在25-30%之间。在热量方面，高寒草甸植被的覆盖作用使得土壤温度变化较为平缓，夏季土壤温度相对较低，冬季土壤温度相对较高，有利于维持活动层的热稳定性。高寒草原植被下的活动层水热状况与高寒草甸有所不同。高寒草原植被覆盖度相对较低，一般在40-60%之间，根系相对较浅，主要分布在0-20cm的土层中。由于植被覆盖度较低，土壤水分蒸发相对较快，活动层含水量相对较低，一般在18-25%之间。在热量方面，高寒草原植被对土壤的保温作用较弱，夏季太阳辐射直接照射地面，土壤升温较快，活动层温度较高；冬季土壤热量散失较快，活动层温度较低，导致活动层温度的季节变化幅度较大。高寒荒漠草原植被下的活动层水热状况更为特殊。高寒荒漠草原植被稀疏，覆盖度通常在20%以下，植被矮小，根系不发达。这种植被类型下，土壤水分蒸发强烈，活动层含水量很低，一般在10-15%之间。由于植被对土壤的保护作用微弱，太阳辐射几乎无遮挡地照射地面，活动层温度受气温影响较大，变化剧烈，夏季温度高，冬季温度低，活动层的冻融循环较为频繁，对多年冻土的稳定性产生较大影响。六、活动层水热动态变化的生态与工程响应6.1生态系统响应6.1.1植被生长与分布变化青藏公路沿线多年冻土区活动层水热动态的变化对植被生长和分布产生了显著影响。随着活动层温度的升高和水分条件的改变，植被的生长状况和物种分布格局发生了明显变化。在植被生长方面，活动层温度升高使得植物的生长季延长。研究表明，过去几十年间，青藏公路沿线多年冻土区的植被生长季平均延长了约10-15天。这为植物提供了更多的生长时间，有利于植物进行光合作用和物质积累，从而促进植被生长。在昆仑山口附近区域，高寒草甸植被在生长季内的生物量明显增加，与活动层温度升高导致的生长季延长密切相关。水分条件的变化对植被生长也有着重要影响。活动层含水量的增加，为植被生长提供了更充足的水分，有利于植被的生长和发育。如在降水增多的年份，青藏公路沿线部分区域的高寒草原植被覆盖度和生物量显著提高。然而，当活动层含水量减少时，植被生长会受到抑制，甚至导致植被退化。在唐古拉山口附近的高寒荒漠草原区域，由于活动层含水量较低，植被生长受到限制，植被覆盖度和生物量相对较低。在植被分布方面，活动层水热动态变化导致植被群落结构发生改变，物种分布范围发生迁移。随着活动层温度升高，一些原本分布在较低海拔或较温暖地区的物种开始向高海拔或较寒冷地区迁移。例如，在青藏公路沿线，一些耐寒性较弱的草本植物逐渐向高海拔地区扩展，而一些原本分布在高海拔地区的耐寒物种则面临生存压力，分布范围可能缩小。水分条件的变化也会影响植被的分布，在活动层含水量增加的区域，一些喜湿植物的分布范围可能扩大；而在含水量减少的区域，耐旱植物的比例可能增加。6.1.2土壤微生物活动变化活动层水热变化对青藏公路沿线多年冻土区的土壤微生物群落结构和活性产生了重要影响。土壤微生物作为生态系统中的重要分解者，参与土壤有机质的分解、养分循环等过程，其群落结构和活性的改变会对生态系统的功能和稳定性产生深远影响。温度是影响土壤微生物活动的关键因素之一。随着活动层温度升高，土壤微生物的活性增强。研究表明，当活动层温度升高1℃时，土壤微生物的呼吸速率可增加约10-20%，这意味着土壤微生物对土壤有机质的分解作用增强，加速了土壤中碳、氮等养分的循环。在昆仑山口附近的高寒草甸区域，活动层温度升高使得土壤微生物的活性提高，土壤中有机碳的分解速率加快，导致土壤有机碳含量有所下降。水分条件对土壤微生物群落结构和活性也有着显著影响。适宜的土壤含水量有利于土壤微生物的生长和繁殖，当活动层含水量增加时，土壤微生物的数量和多样性通常会增加。在降水较多的年份，青藏公路沿线部分区域的土壤微生物数量明显增多，群落结构更加复杂。然而，当活动层含水量过高或过低时，都会对土壤微生物产生不利影响。含水量过高会导致土壤通气性变差，使一些好氧微生物的活动受到抑制；含水量过低则会使土壤微生物的生存环境恶化，影响其活性和数量。活动层水热变化还会导致土壤微生物群落结构的改变。不同种类的土壤微生物对水热条件的适应能力不同，在水热条件变化时，微生物群落中各物种的相对丰度会发生变化。研究发现，随着活动层温度升高和水分条件的改变，青藏公路沿线多年冻土区土壤微生物群落中，一些嗜热、耐旱的微生物种类相对丰度增加，而一些对水热条件要求较为严格的微生物种类相对丰度则可能减少。这种群落结构的改变会影响土壤生态系统的功能，如土壤的养分供应能力、抗干扰能力等。6.2工程设施响应6.2.1青藏公路路基稳定性青藏公路作为连接内地与西藏的重要交通要道，其路基稳定性直接关系到公路的安全运营和区域经济的发展。活动层水热动态变化对青藏公路路基稳定性产生了显著影响，引发了一系列路基病害。随着全球气候变暖，青藏公路沿线多年冻土区活动层温度升高，导致多年冻土上限下降，活动层厚度增加。活动层厚度的变化改变了路基下部土体的受力状态和热传导路径。当活动层厚度增大时，路基下部多年冻土承受的上部荷载增加，且由于活动层热阻的变化，热量传递至多年冻土的速率加快，使得多年冻土更容易发生融化，从而降低了路基的承载能力。研究表明，在青藏公路部分路段，活动层厚度每增加10cm，路基的沉降量可增加约5-8mm，严重威胁路基的稳定性。活动层含水量的变化对路基稳定性也有着重要影响。在夏季，降水增加和积雪融化导致活动层含水量升高，土体饱和度增大，抗剪强度降低。当含水量超过一定阈值时，土体容易发生软化和塑性变形，导致路基出现沉陷、坍塌等病害。在青藏公路的一些高填方路段，由于活动层含水量的季节性变化，每年夏季都会出现不同程度的路基下沉现象，影响路面平整度和行车安全。此外，冬季活动层冻结时，水分结冰膨胀，会产生冻胀力，对路基结构产生破坏作用，导致路面隆起、开裂等病害。活动层水热动态变化还会引发路基的融沉病害。融沉是指多年冻土融化后，土体中的冰变成水，体积减小，导致地面下沉的现象。在青藏公路沿线，由于活动层水热条件的改变，多年冻土融化加剧，融沉病害日益严重。融沉病害不仅会导致路基沉降，还会使路基边坡失稳，增加了公路养护和维修的难度和成本。据统计，青藏公路沿线约有30%的路段存在不同程度的融沉病害，其中严重融沉路段占比约为10%，给公路的正常运营带来了极大的安全隐患。6.2.2其他工程设施影响除了路基稳定性，青藏公路沿线的桥梁、涵洞等工程设施也受到活动层水热变化的显著影响。桥梁基础多位于多年冻土区，活动层水热动态变化会改变桥梁基础周围土体的物理力学性质，进而影响桥梁的稳定性。活动层温度升高导致多年冻土融化，使得桥梁基础周围土体的强度降低，承载能力下降。当桥梁基础周围土体发生融化时，基础可能会出现不均匀沉降，导致桥梁结构受力不均，产生裂缝、变形等病害。在青藏公路沿线的一些桥梁中，由于活动层水热变化，桥梁基础出现了不同程度的沉降，部分桥墩出现了倾斜现象，严重影响了桥梁的安全使用。涵洞作为公路排水系统的重要组成部分，其正常运行对于保护路基、防止水毁至关重要。活动层水热变化对涵洞的影响主要体现在两