华北西部地区月气候变化与水资源关系的插值分析.docxVIP

华北西部地区月气候变化与水资源关系的插值分析 1 气候变化与高效降水 随着世界经济的进程，中国北方平原的干旱问题变得越来越明显。如何解决这一地区的缺水问题是政府面前的一项重要任务。要解决这一问题，必须了解中国北方水资源未来变化的趋势，水资源变化与气候变化密切相关。以往的研究表明,80年代华北地区平均气温比50年代升高0.54 ℃,降水是近百年最少的10年,这种暖与干的搭配使得河川径流量、湖泊蓄水量、土壤含水量、地下水量都明显减少,引发了华北地区水资源的危机。对华北地区未来气候变化的预测成了国家是否建设南水北调工程的主要因素,80年代末至90年代初,许多科学家对华北地区气候变化做了研究,基于大气中温室气体含量逐年上升的现实,对气温变化的认识比较统一,认为华北地区继续升温的可能性较大,至2030年,华北地区温度比80年代可能增加2.0±1.0 ℃。华北降水的多少取决于季风的强弱,季风增强将使华北进入多雨期。近年来,在全球进一步变暖的形势下,华北地区的降水变化如何呢?本文将利用气温、降水和河川流量资料进行分析研究。 2 河流月流量测定 选用1980～1996年9个水文站点(响水堡、石匣里、张家坟、下会、滦县、小觉、北河店、观台和楚望)的月流量资料,同时选取该区域42个气象站的月降水和气温资料来进行分析研究,站点的分布如图1所示。观台和楚望两站缺1982年1～12月的流量值,观台还缺1992年12月、1993年5月的流量值,张家坟缺1983年10月、滦县缺1983年4月的流量值。由于获得的流量资料个别年份和个别月份缺测,这给我们的分析研究带来困难,为进一步工作的需要,必须首先对缺测资料进行插补。 以上分析表明,河流月流量不仅与同期降水、气温,而且与前期的降水、气温有一定的相关关系。利用相关分析插补或延长水文资料,相关系数的绝对值不能小于0.8。表1中的相关系数都达不到相关分析插值的要求,不宜用相关分析法对月流量缺测值进行插补。 3 缺测值的计算 EOF作为一种研究手段,综合了时空变化的信息,已得到十分广泛地应用,张邦林等用实际资料对不同季节和不同年份分别作EOF展开,证实了生成的经验正交函数对大气要素的展开精度是稳定的。董文杰将影响降水的各种因素列于同一EOF分解中,获得了有意义的成果。根据EOF的特性,我们将月降水、气温和月流量置于同一EOF中,既考虑它们的时空演变,又考虑到相互间的影响,来获取缺测值,具体方案如下: (1) 对1980～1996年各月降水、气温和流量分别求其年际标准化距平,缺测站月的年际标准化距平暂取为0。小觉站没有缺测值,假定小觉站1982年1～12月的值也为缺测值,相应地取年际标准化距平为0。然后将标准化距平值按1980年1月～1996年12月共204个月(i=1,204)的次序排列,作为时间维,将42个站的降水、气温、9个站的流量共93个站(j=1,93)作为空间维,进行EOF分解。 (2) 根据EOF的收敛性,选择累积方差占总方差一定比例的前n个特征向量场和时间系数,按下式对缺测变量重新组合。 Xij=∑k=1nak1Ekj.(1)Xij=∑k=1na1kEjk.(1) (3) 将缺测站月的年际标准化距平用(1)式计算结果代替,再进行EOF分解。 (4) 重复步骤(2)、(3), 将由(1)式计算得到的小觉站1982年1～12月的年际标准化距平还原为月流量值,与实测值相减,再除以实测值,得相对误差,求各月相对误差绝对值的平均值(称为平均相对误差)。当再进行一次EOF分解与组合所得的月流量值与实测值比较,平均相对误差比前一步所得到的要大时,停止步骤(3),将所得缺测站月的年际标准化距平还原为月流量值,这些月流量值就是我们得到的插值。 4 月流量的插值 EOF分解将最主要的成份排列在前,而将小扰动排列在后,我们选取了不同的截断数n及其前n个主成份占总方差的比值B,按照上述插补方案进行试验。结果表明,65%截断(截断数为6)且循环4次的平均相对误差最小,也就是说,前6个特征向量的重新组合达到最佳效果。选取65%截断且循环4次的前6个特征向量重新组合所得的标准化距平,将其还原为月流量值,作为缺测站月流量的插值结果。 图3给出了小觉站1982年1～12月流量的EOF插值(实线)和实测值(虚线),总的来说,二者已相当接近。由于受诸多气候因素的影响,月流量的年际变幅比较大,不宜用简单的线性插值对缺测值进行插补。为了比较,这里我们给出了小觉站1982年1～12月流量的线性插值(图3中点划线)的结果,即缺测月前一年和后一年该月流量的平均值。很明显,线性插值的效果较差。 图4为观台(实线)和楚望(虚线)1982年1～12月的流量EOF插补值,其年变特征是很明显的。 5 日变化时的特征向量场 影响水资源的气候因素主要是气温、降