资料和方法-中国气象局国家气候中心任国玉课题组.doc

第二章 资料和方法 §2.1 资料 地面观测资料 来源于国家气象信息中心气候资料室，包括1961—2003年中国北方340个站的扬沙、沙尘暴、大风日数、降水量、温度、风速、湿度、水面蒸发量的逐年逐月资料；风速的逐日资料、记录有沙尘暴发生初始时刻的日值资料;1961—2001年逐月的积雪深度资料。其中，蒸发量是用小型蒸发皿观测的，可代表水面蒸发量或最大可能蒸发量，特别是在干燥地区与实际蒸发量存在较大的差异。考虑到蒸发资料1961—1963年中月缺测值比较多，2002年以后记录仪器型号的变动，选出了1964—2001年38年的蒸发量资料与沙尘日数资料进行相关分析。四季划分按气象季节划分方法，即前一年12月～当年2月为冬季，3月～5月为春季，6～8月为夏季，9～11月为秋季。对于序列中存在缺测的部分，根据需要进行了处理（用平均值取代）和剔除。 植被指数资料 来源于（）美国地球资源观测系统数据中心的探路者数据库(NOAA/NASA Pathfinder AVHRR Land Data Sets) 所提供的1982－2001年归一化植被指数（NDVI）格点数据集，这是迄今最长的卫星遥感植被观测数据，空间分辨率为1o×1o经纬度。文中将离台站最近的格点资料插值到台站作为该台站的植被指数数据资料。 北大西洋涛动（NAO）指数 NAO是北大西洋地区大气最显著的模态。很早人们就发现北大西洋地区的气压变化存在一种南北向跷跷 NAO的定量表征指标是NAO指数，来源于CRU（东Anglia大学气候研究机构,UK,P.D.Jones的逐年逐月的NAO指数值，时间为1961—2003年。 南方涛动指数（SOI） 沃克发现东南太平洋（塔希提岛以东）和印度洋（澳大利亚的达尔文站以西）气压之间有一种跷跷板式的关系，也就是说，其中的一个地区气压偏高，另一个地区的气压反过来就偏低。指塔希提与达尔文站海平面气压差。 与太平洋温度的周期变化类似，太平洋上空的气流也有近似的周期变化。“拉马德雷”是一种高空气压流，亦称为太平洋十年涛动 (PDO)。是存在于太平洋SST变率中的一种缓变空间状态，当北太平洋中部维持一个冷信号时，在北美西海岸则为暖信号。 来源于CPC (NOAA,气候预测中心)的逐年逐月的PDO指数标准化值，时间为1961—2003年。 北极涛动（AO） 指数 Thompson 和Wallace基于北半球大气环流的分析，提出了北极涛动(AO：Arctic Oscillation)的概念，Wallace最近指出，AO与NAO二者本质上是一致的，是同一事物在不同侧面的两种表现，实际上反映的都是中纬西风的强弱；只不过AO尺度更大，而NAO是其在北大西洋区域的表现。 来源于CPC (NOAA,气候预测中心)的逐年逐月的AO指数值，时间为1961—2003年. 太阳黑子数 用来表示太阳活动的强弱的指数，来源于比利时太阳影响资料分析中心（SIDC），选取1961—2003年内的逐年逐月数据。 蒙古国降水资料 来源于NOAA/NCEP地面降水量再分析格点场资料，空间分辨率为2.5o×2.5o经纬度. 以及国家气象中心气候资料室提供的蒙古国45个台站地面降水量的逐年逐月资料. 青藏高原及我国北方积雪资料 来源于国家气象中心气候资料室，包括青藏高原60个站点及我国北方340个站点的地面积雪日数的逐年逐月值。 §2.2 方法 文中主要采用了气象统计学[56]上的相关分析方法。 计算区域平均的方法 根据1996年Jones等[57]提供的用经纬度网格面积加权平均来计算某一气候要素区域平均值的方法，建立了各区多类多要素的平均时间序列，计算了多组相关性。这种方法的基本思路是：以逐年序列为例，首先将某区域按一定的经纬度划分为n个网格，再将每个格点内某年某气候要素值求算术平均，然后根据网格面积做加权平均，得到此年此气候要素的区域平均值，同样再求出其它年份的平均值，就可以建立起此气候要素区域平均值的时间序列。计算公式如下： 。 其中，即为要计算的某要素x的区域平均值， 为网格数，（i=1,2,…n）为每个网格的中心点纬度，（j=1,2,…）为区域中某一网格内共参与计算的个点的气候要素值。 文中使用了2.5°×2.5°经纬度的网格点。 （2）求线性趋势及检验趋势变化是否显著的方法 文中主要利用此方法计算了沙尘日数逐年序列的线性变化趋势。将某气候要素X的时间序列看成时间T的一元线性回归函数——X＝a + b T，用最小二乘法求得线性回归方程的斜率b： 其中， （i=1,2,…n）为某气候要素逐年序列；（i=1,2,…n）为对应的年份时间序列，、分别为两序列的平均值。b的大小和符号即反映了此气候要素随时间变化的快慢和总体趋势。当b0时，说明气候要素X随时间T的增加