青藏高原的热量变化特征

滨江学院学年论文题目：青藏高原热量的变化特征学生姓名刘凡学号20112316011院系大气与遥感专业大气科学指导教师张春莹二Ｏ一四年十二月二十九日论文目录摘要......................................................11.引言....................................................12.资料选取和方法介绍......................................22.1资料选取.............................................22.2方法介绍.............................................33.高原感热潜热总辐射的时空分布............................43.1感热的时空分布特征...................................43.2潜热的时空分布特征...................................53.3总辐射的时空分布特征.................................64.结论....................................................74.1感热.................................................84.2潜热.................................................84.3总辐射...............................................85.参考文献................................................9青藏高原热量的变化特征刘凡南京信息工程大学滨江学院大气科学专业，南京210044摘要：利用NCEP/NCAR地面感热资料再分析格点资料，分析了1951-2010年青藏高原（下称高原）地区地面感热通量的基本气候特征,年纪与年代变化及其空间分布，采用滑动t检验和小波分析研究了高原年平均感热通量的变化的突变特征。结果表明，就全年平均而言，高原感热通量大部分地区为正值，说明高原为热源；冬季是全年感热通量最小的季节，为负值；其余季节感热均为正值。高原地表热源表现出明显的非均匀性，高原东西部热源呈反向变化。近60年高原的感热通量出现了不同程度的减少。利用热带测雨卫星（tropicalrainfallmeasuringmission.TRMM）上微波成像仪（TRMMmicrowaveimage,TM1）的观测资料，利用卫星遥感对夏季高原地区的凝结潜热(下称潜热)水平分布形式，潜热垂直结构及其变化特征进行了分析，潜热水平分布形式与基于NCEP/NCAR降水资料的基本结果一致，表明TM1的资料在高原地区的可用性。结果显示了夏季高原上三个比较稳定的潜热加热中心，显示了显著地高空热源。利用高原及其附近的22个日射站将近40a的总辐射及百分率资料确定了Angstrom-Prescott（APM）系数，在结合高原及其邻近地区的资料估算了将近40a的总辐射。结果表明：高原主体光照充沛，空间分布上，高原西部为高值区：年代变化上，高原总体减小，但波动较小，较稳定。关键词：青藏高原感热通量凝结潜热总辐射时空分布1引言青藏高原（下称高原）作为世界上最高和最大的高原，具有独特的地理面貌特征，是全球气候变化研究的重要区域，也是“全球变化的驱动机与放大器”【1】。作为世界的“第三极”，高原已经成为气候变化研究领域继南极与北极之后的另一热点。高原大气对大气环流，降水以及海气相互作用过程有着极为重要的热力作用，吴国雄等【2】指出，高原的热力作用就像一个巨大的“感热气泵”调节着周边乃至全球的大气运动。大气有3种热量来源，即地面感热，大气净辐射以及凝结潜热【3】。高原对大气的热力影响在干季以感热为主，在雨季则由于降水量的增加，凝结潜热可以与感热达到同量级，成为雨季影响大气环流的重要能量之一。感热方面：20世纪后半叶高原地区气候发生了显著的变暖【4】。高原本身温度的升高伴随着其热力作用变化。而热力作用变化的具体体现是看其冷、热源作用过程是否变化?如感热通量，那么影响它变化的因素是什么及其变化的天气气候意义如何?数值模拟试验表明，高原气温升高引起大气加热能够加强东亚亚热带锋而雨的爆发【5】。李栋梁等【6】利用高原气象观测站，通过经验方法计算了高原地而感热通量并进行了异常诊断分析。Duan等【7】也用传统的经验方法估算了高原的感热通量，并指出感热通量在减弱，尤其在春季。Yang等【8】利用气象观测资料，通过对比分析两类计算方法(物理方法和经验方法)获得高原地表感热通量的变化趋势，结果表明，由物理方法估计的年均感热通量以每10年2%的速率减弱，经验方法不能很好地反映感热的变化趋势，并指出在气候变化研究中如果采用经验的计算方案(即使应用很广的方案)，结果的偏差和确定性可能非常大，甚至可能与实际完全相反。既然常规气象观测中是没有地表通量的直接观测，更不可能有区域尺度范围的格点值。而且中国气象台站在高原的分布相当不均匀，中东部较密集，西部很少，低海拔地区台站较多，高海拔地区台站很少，没有海拔超过4800m的台站【9】。因此，高原地表感热估算一直存在争议。NCEP/NCAR提供了连续性好、时间序列长的全球地面通量等加热场的再分析格点资料【10】。宋敏红等【11】研究表明高原地区的NCEP/NCAR地面和大气热通量再分析资料能较好地反映该区热源强度的年及年际变化特征。不同学者，【12】使用NCEP/NCAR通量资料分别研究了高原视热源的年际变化、感热通量、潜热通量等的气候状况及其与初夏东亚大气环流之间的关系。潜热方面：长期以来，有许多关于凝结潜热方面的研究，有基于地面观测站资料【13】和FGGE资料【14】，其分析结果与前人的基本一致。徐国昌等【15】和李栋梁等【6】等应用OLR资料对凝结潜热进行估算，认为高原上OLR与降水有很好的相关性，都可以反映出凝结潜热的分布特征。但是，迄今为止，因全球气候变暖，高原大气潜热状况在时间和空间上如何变化，整个凝结潜热总量如何变化等问题还没有系统的结论。总辐射方面：目前，地面接收到的短波辐射能的计算大体有三类方案:其一是利用辐射传输模式计算辐射通量【16】，此类方案由于计算模式相对复杂，模式需要输入的实际大气参数难以获得以及需要较多的时机，因而该类方案在实际应用中受到了一定的限制；其二是遥感反演【17】，此类方案在区域及全球尺度太阳辐射研究有较好的应用前景，但受相应卫星过境时段的影响，选用不同卫星的资料，其反演结果缺乏对比性，高原地区地面辐射观测站点本来就很少；其三是辐射计算的经验模型，此类方案利用日照时数、气温、降水、相对湿度以及总云量等较容易获得的地面气象数据与太阳辐射实测值之间建立回归关系式来计算辐射值。此类方案由于简单、所用的资料容易获取，因而在国内外辐射气候的研究中得到了较广泛的推广应用。总辐射的气候计算方法研究是其中比较早也比较重要的，从20世纪20年代开始，国外先后有多位学者给出了一些总辐射的计算式，并不断演变、发展完善【18】。相对于国外研究，中国的总辐射气候计算方法研究起步较晚，主要是从1957年开始建立全国日射站网后，才真正发展起来的.尹宏【19】，萧文俊【20】分别利用当时有限的实测资料，验证Angstrom-Prescott模型(APM)在我国的适用性.1963年左大康等【21】最先系统地研究了我国总辐射的气候计算问题，提出了以Angstrom式为基础的全国统一的总辐射气候计算式，并首次绘制出了具有气候学意义的中国年、月总辐射分布图，讨论了其分布规律.这一工作对推动我国整个辐射气候研究起到了重大作用.之后许多的研究结果表【22】，APM模型可用来估算中国地区总辐射APM模型由于误差较小，模型输入参数较少，且相应资料易于获取等优点，近年来己在相关研究中得到了较为广泛的应用【23】。本文在感热上利用NCEP/NCAR再分析资料中的地表感热通量分析高原地区的感热通量季节变化特征，研究高原感热通量的年际、年代际变化及其空间分布，并用突变检验和小波分析研究了年平均感热通量变化的突变特征，从而对高原地而感热通量异常变化的空间结构和时间演变趋势作了较系统的诊断研究，利用不同资料分析了影响感热通量变化的因素，最后探讨了高原感热通量的变化对东亚、南亚夏季风的影响。在潜热上本文针对地形比较复杂的青藏高原东部。采用最新获取的BEM资料，基于GIS软件技术对高原东部46a(1961-2006)以来的降水和凝结潜热进行估算，并给出不同分区不同时间的凝结潜热总量，力求为高原大气热量源汇的研究探索出一种新的方法。在总辐射上本文利用高原及其周边22个日射观测站的实测资料确定APM模型参数，以此模型估算了高原及周边116个气象参站1961-2000时段40a的总辐射，研究了近40a青藏高原地区总辐射的分布变化特征。２资料选取和方法介绍2.1资料选取感热的资料：所用的气温和降水资料来自英国EastAnglia大学的1951年1月一2009年12月的高分辨率全球逐月格点数据集(CRUTS3.1,)，水平网格距为0.50×0.50。所用的10m处月平均风速和地温以及感热通量资料均来自NCEP/NCAR再分析资料，水平分辨率为高斯网格(192×94个格点)，研究时段为1951一2010年。潜热的资料：本文研究区域包括青海省、西藏大部、四川大部以及甘肃大部分地区。本文研究所需的DEM(数字高程模型)是美国NASA(NationalAeronauticsandSpaceAdminis-tration)和NGA(NationalGeospatial-IntelligenceAgency)于2000年2月历经11d所获得的高分辨率高程数据，使用的仪器是TheShuttleRadarTo-pographyMissionCSRTM)探测器，DEM数据分辨率为3，数据来源于马里兰大学网站。总辐射的资料：高原及周边22个辐射站的总辐射资料源于国家气象中心出版的中国地而辐射资料年册。本文选取APM模型估算总辐射，该模型输入参数少，使用简单方便，是目前较为广泛采用总辐射估算模型。2.2方法介绍1、感热通量NCEP/NCAR资料中感热通量的计算公式采用整体输送法，即其中:H为感热通量；为大气密度；为定压比热；V为10m风速；为地表温度；为2m气温；为整体交换系数，是一个经验值。若交换系数己知，则感热通量可以直接由式（1）计算。对高原年平均感热通量用滑动，检验、小波分析【24】进行突变检验、趋势变化和周期分析。t检验统计量:其中：；n1，n2分别为两段时间的年数；和分别为两段子系列的平均值；和分别为方差。一维信号函数f(t）的小波变换为其中：a为频率参数；b为时间参数，表示波动在时平移；R为实数域。小波函数有多种函数形式，本文选用高斯型小波，也叫墨西哥帽小波，即:2、潜热通量高原雨季降水具有明显的垂直分布特征，根据降水量随海拔高度的上升成三阶多项式的变化特点，将研究区域分为三类，即h1400m,1400mh3600m,h3600m，估算方法见文献【25】。凝结潜热指大气中水汽凝结所释放的热量，降水凝结潜热QW的计算公式为其中:QW必为凝结潜热，量纲为焦耳（J）P为网格区平均降水量(cm)；S为小区面积(cm2)；P为水的密度(1g·Cm-3)；Lv为凝结潜热系数，即Lv=2497J·g-1。3、总辐射青藏高原地区APM可写成如下形式：（4）模型估算效果由平均偏移误差(MBE)、平均绝对误差(MAE)、均方根误差(RMSE)及平均相对误差(MRE)等量进行评估，相关细节见文献【26】。