摘要:青藏高原是地气相互作用相对活跃的地区，深入了解青藏高原土壤热扩散率的变化，才能正确计算地表能量平衡进而准确认识青藏高原对全球和区域气候变化的影响。根据青藏高原1980—2001年39个观测站点实测的0.8m和3.2m土壤温度资料，利用热传导对流法结合最小二乘法拟合求得各站点的土壤热扩散率，并分析了土壤热扩散率的时空变化规律。结果表明，1980—2001年期间青藏高原土壤热扩散率在20世纪90年代以前波动较大，20世纪90年代以后波动较小。青藏高原东部地区的深层土壤热扩散率从春季至夏季增大，夏季至秋季减小，秋季至冬季减小；夏季最大值出现在青、川、甘三省的交界处，土壤热扩散率的值为8×10-6m2s-1，冬季最大值为5.1×10-7m2s-1；而除东部以外的青藏高原其他地区的土壤热扩散率，春季至夏季减小，夏季至秋季增加，秋季至冬季减小，该区域土壤热扩散率的变化范围为1.2×10-7m2s-1～9.2×10-7m2s-1。土壤热扩散率的多年平均最大值出现在青海省和甘肃省西南部以及四川西部的青藏高原东部地区，土壤热扩散率的极值为6.4×10-6m2s-1。最小值出现在祁连山地区，土壤热扩散率为1.2×10-7m2s-1，中部地区为相对高值区。
关键词:青藏高原,土壤热扩散率,空间分布,季节差异,土壤温度,热传导对流法
青藏高原在亚洲季风中扮演了重要角色，青藏高原的气候变化不仅直接驱动中国东部和西南部气候的变化，而且对整个亚洲和北半球也有巨大影响，甚至对全球气候变化也具有明显的敏感性、超前性和调节性。此外，青藏高原的地气物理过程对全球气候异常、东亚大气环流及中国灾害性天气的发生发展有重大影响。青藏高原的热力作用是由地气相互作用决定的。地气之间的能量和物质交换是决定气候状态和大气环流特征的最主要物理过程之一，而定量描述地面的热力过程、土壤热传导和土壤热输送是数值预报和气候模式物理中陆面过程参数化的任务之一。土壤的热力学性质影响着地表和大气之间的物质和能量交换，对地球大气边界层物理过程、大气环流和区域气候产生重要影响，正确地模拟土壤温度和热通量是陆面物理过程模式、大气环流模式、区域气候模式数值模拟研究的重要环节之一。因此，深入研究土壤的有关物理参数（土壤热扩散率、土壤热传导率、土壤体积热容等）有利于更好地分析能量分配，从而为利用一维、中尺度和大尺度模式模拟地气相互作用和天气过程提供参数化方案。土壤热扩散率是计算土壤热通量的输入项，土壤中热量的传输是陆面过程的一部分，也是影响地表能量平衡的重要因素之一，土壤中热量的吸收和损失均与气候变化息息相关。因此，对土壤热扩散率时空分布的研究，有利于深刻认识土壤热通量的物理过程，为各种天气、气候及生态模型的建立和检验提供重要参数和理论支持。
由于受观测技术和观测资料的限制，国内外对土壤热扩散率的观测研究并不多。因为土壤的热扩散率取决于诸多因素，如土壤质地、矿物组成、含盐量和土壤含水量等，使得热参数的估计非常困难。许多研究人员试图通过在实验室条件下使用土壤含水量、密度和形状因子来参数化某些土壤的热扩散率，然而，由于尺度效应和土壤形状、颗粒大小以及包装等的复杂影响，当实验室结果外推至田间实验时，通常会产生误差。为了避免这些误差，开发了使用土壤温度观测数据估计土壤热扩散率的方法。董文杰和汤懋苍用相位差法计算了多年平均土壤热扩散率，分析了我国土壤热流场与深层大地热流场的关系。1993年王宝灵和汤懋苍用谐波法计算了土壤的导温率，并分析了我国土壤热流的空间分布。范新岗和汤懋苍利用同时考虑了土壤热传输和土壤水分热对流的热传导对流法计算了土壤的热扩散率，得出了传导热通量与对流热通量具有相同的量级，且分析了热通量与地震的关系。高志球等用不同方法求出了土壤热传导方程的解析解，并对那曲地区土壤热扩散率的变化进行了分析。缪育聪等用谐波法、振幅法、相位法和热传导对流法计算了5～20cm的土壤热扩散率，发现热传导对流法的拟合效果好于振幅法和相位法，在理论上能够更好地反映土壤中的热传输过程。
以上学者的研究基本均是通过理论分析，或者对浅层土壤或者局限于某一区域的土壤热扩散率的分析研究，对整个青藏高原的研究并不多。在诸多模式研究中也将不同土壤类型的土壤热扩散率取为定值，但实际上相同土壤类型的土壤热扩散率也有其年际变化和空间变化。目前尚无对青藏高原土壤热扩散率空间分布的研究。本研究利用实测土壤温度资料和热传导对流法分析了青藏高原深层土壤热扩散率的时空分布特征。
1材料与方法
1.1 土壤温度资料来源
1.2 本研究利用青藏高原39个气象站（如图1）实测的0.8m和3.2m深的月平均土壤温度，资料年限为1980—2001年。本资料来自中国气象局，资料选择的依据是：（1）此39个站点土壤温度资料相对齐全、无缺测；（2）80cm以上各层地温因受大气的强烈影响，难以说明深层变化，而且0.8m和3.2m两层深度是中国气象局地面观测的规范深度，包括了0.8m～1.6m和1.6m～3.2m两个深度。通过计算，得到土壤热扩散率的月和季平均值，对青藏高原土壤热扩散率的时空分布特征进行分析。
1.2土壤热扩散率的计算——热传导对流法
计算土壤热扩散率的方法有很多，如振幅法、相位法、谐波法和热传导对流法等。本文采用热传导对流法，因为土壤的热传输过程包含了热传导和多孔介质热对流两个过程。Gao等指出土壤温度的日变化不仅受土壤热传导作用的影响，还与土壤中液态水热对流作用有关，给出了耦合热传导与热对流的一维热传导对流方程，并利用一阶谐波模型推导出了热传导对流方程的解析解和土壤热扩散率的计算公式。Gao等利用青藏高原和黄土高原等地的土壤温度资料，验证了热传导对流法的可靠性，认为热传导对流法能更好地描述土壤中的热传输物理过程。在以上理论的基础上发展了一套包括热传导和热对流的计算土壤热通量方法，该方法可以利用气象站地温资料反演获得地热流场信息，从热传导方程和热对流方程出发推导出了包括传导和对流的线性传热方程：
2结果与讨论
利用热传导对流法计算了青藏高原39个气象站各站点1980—2001年的多年月平均土壤热扩散率和多年季节平均土壤热扩散率。在季节划分上，本文仍以3月—5月为春季，6月—8月为夏季，9月—11月为秋季，12月—次年2月为冬季。
由于土壤特性和气候特征的差异，青藏高原地区土壤热扩散率在时空上表现为非均匀性。因此，本文根据各站土壤热扩散率的多年月平均值和多年季平均值分析了青藏高原土壤热扩散率的时空分布特征。
2.1青藏高原土壤热扩散率年际变化特征
用热传导对流法计算的青藏高原土壤热扩散率K的年际变化如图2所示，可以看出，土壤热扩散率在20世纪90年代以前波动较大，在1991—2001年有小幅度波动。K的年际变化范围为3.4×10-7m2s-1～8.9×10-7m2s-1。根据杨玮等的研究，1980—2001年青藏高原的长期降水变化也呈现出1990年之前波动较大，90年代之后波动较小的趋势，说明土壤热扩散率的年际变化与降水的年际变化较为相似。青藏高原土壤热通量的这种变化趋势还与张雯等研究的青藏高原土壤湿度的年际变化特征一致。根据邵明安等的研究，土壤热扩散率与土壤含水量有关，开始会随着土壤含水量的增加而增加，达到一定程度之后反而随着土壤含水量的增加而减小。此外，图2给出了土壤热扩散率标准差的变化，标准差反映的是离散程度大小。可以看出，每年的标准差各不相同，说明土壤热扩散率的变化不同。1993年标准差较小，说明1993年土壤热扩散率变化小，即空间差异小。而1985年土壤热扩散率的标准差相对较大，说明土壤热扩散率变化大，即空间差异大，出现这种现象的原因可能与该年青藏高原各地区降水量的不同有关，这种关系值得进一步研究。
2.2青藏高原不同季节土壤热扩散率空间分布特征
图3为青藏高原不同季节多年平均土壤热扩散率的空间分布。从图中可以看出，春季（图3(a)）西南部、中部和东部小部分地区的土壤热扩散率为相对高值区，最大值出现在藏北高原那曲东部、怒江上游的南姜塘大湖盆区，达6×10-6m2s-1。其余部分为相对低值区，土壤热扩散率的范围为1.3×10-7m2s-1～9.2×10-7m2s-1。夏季（图3(b)）青藏高原东部小部分仍为相对高值区，最大值出现在青、川、甘三省的交界处，为8×10-6m2s-1，土壤热扩散率相对春季有所增大；其余大部分地区的土壤热扩散率减小，范围为1.21×10-7m2s-1～7.4×10-7m2s-1，特别是青藏高原北部地区土壤热扩散率减小明显。该地区位于长江河源高平原连续多年冻土区。根据土壤的不同冻融过程，在春季处于升温期、夏季处于融化期，冰的导热率几乎是水导热率的4倍，因而冻土的导热率大于融土的导热率，而土壤体积热容的取值与季节无关，是由土壤类型决定的。青藏高原北部地区主要为粉砂土，土壤的体积热容相对固定，因此会出现夏季导温率较冬季导温率小。位于青海省和甘肃省西南部、四川西部的青藏高原东部地区为高原亚寒带湿润气候区，海拔在3000m以上，地势复杂，受高原季风影响年降水量400～780mm，并主要集中在夏季。而土壤热扩散率的大小主要受质地和含水量的影响，含水量在夏季增大，因此，青藏高原东部土壤热扩散率的最大值出现在夏季。
秋季（图3(c)）东部小部分地区的土壤热扩散率开始减小，但仍为相对高值区，最大值为6.2×10-6m2s-1；中部地区也为相对高值区，其他地区的土壤热扩散率相对夏季有所增大，土壤热扩散率的变化范围为1.3×10-7m2s-1～8.9×10-7m2s-1。冬季（图3(d)），青藏高原大部分地区的土壤热扩散率均减小，西南部、中部和东部小部分地区为相对高值区，最大值为5.1×10-6m2s-1。从秋季至冬季寒冷干燥、降水量减小，因而东部小部分地区土壤热扩散率减小；而其余大部分地区的土壤热扩散率减小是因为从秋季至冬季土壤温度减小至零度以下。根据Inaba的实验，在负温度范围内，土壤有效导热率的测量值强烈依赖于土壤温度、含水量及粒度，对于细粒土壤而言，导热率随负温的降低而减小。并且，在冬季，未冻结的水含量较小，所以土壤热扩散率减小。
此外，从冬季和夏季分布图可以看出，冬季高原中北部地区的土壤热扩散率较夏季大，因为此地区位于巴颜喀拉山脉南麓，该地区积雪覆盖时间高达40d左右，积雪的出现会阻挡热量的交换，使深层温度下降小，而在此期间土壤热通量有小幅度增大，导致土壤热扩散率增大。
2.3青藏高原多年平均土壤热扩散率空间分布特征
图4给出1980—2001年青藏高原土壤热扩散率的多年平均值空间分布，可以看出，青藏高原大部分地区土壤热扩散率的量级为10-7m2s-1，中部和东部小部分地区为相对高值区，量级为10-6m2s-1。
受地理位置、土壤特性、含水量和积雪日数等因素的影响，在西藏那曲东北部的青藏高原中部地区和东部地区为高原亚寒带湿润气候区，由于降水较多，土壤含水量较大，为土壤热扩散率的高值区。祁连山地区为低值区，该地区为高原高寒带半干旱气候，为多年冻土区，含水量少，土壤热扩散率范围在1×10-7～4.8×10-7m2s-1。
3结论
利用实测土壤温度资料分析了青藏高原土壤热扩散率在时空上的分布规律，对土壤热扩散率变化特征进行了研究，得到如下结论：土壤热扩散率在20世纪90年代前期波动较大，在1991—2001年有小幅度波动；受气候要素、土壤特性和含水量的影响，土壤热扩散率在空间上的分布是不同的。青藏高原东部地区，从春季至夏季土壤热扩散率增大，夏季至冬季土壤热扩散率逐渐减小；除青藏高原东部以外的其他地区，从春季至夏季土壤热扩散率减小，夏季至秋季土壤热扩散率增大，秋季至冬季土壤热扩散率减小。青藏高原土壤热扩散率月平均最大值出现在青海省和甘肃省西南部以及四川西部的青藏高原东部地区。青藏高原中部的那曲东北部为相对高值区，其他区域为相对低值区，土壤热扩散率的最低值分布在祁连山脉地区。青藏高原大部分地区土壤热扩散率的量级为10-7m2s-1，中部和东部小部分地区为相对高值区，量级为10-6m2s-1。
立即询价
