2. 国家气候中心,中国气象局气候研究开放实验室,北京 100081;
3. 南京信息工程大学气象灾害预报预警与评估协同创新中心,南京 210044
2. Laboratory for Climate Studies, National Climate Centre, CMA, Beijing 100081;
3. Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters,Nanjing University of Information Science and Technology, Nanjing 210044
对东亚及中国夏季气候的年际、年代际变化已有大量研究,多是针对降水。对于气温,关注的季节侧重于冬季(孙林海和赵振国,2004;李维京等, 2013;王东阡等,2013;2015;伍红雨等,2014;司东等,2014; 2016;江琪等,2016),其中李维京等(2013)全面研究了我国冬季气温与北极海冰、冬季风、西伯利亚高压等影响因子的关系,指出其相关关系发生了显著的年代际变化,要在冷、暖期分时段建立冬季气温预测模型。而对于夏季气温,关注重点在于变化趋势、区域气候特征及其影响系统。在全球变暖背景下,中国气候变化存在着明显的季节性和区域性特征,如“冬暖夏凉”和“北暖南凉”等(丁裕国和江志红,1999;秦爱民等,2005;唐国利和丁一汇,2006)。区域气温的重点是东北低温,已有研究探讨了积雪、马斯克林高压等下垫面外强迫信号及大气环流异常对东北夏季气温的影响(孙建奇和王会军,2006;张茜和李栋梁,2012;曾宇星等,2014;韩荣青等,2014;Chen et al,2016)。还有研究对南方地区夏季低温和阶段性高温特征及其成因做出分析(蔡佳熙和管兆勇,2007;唐恬等,2014;崔童等,2015),尤其是夏季变冷趋势最明显的长江流域,从年际、年代际变化角度,解析了长江流域夏季气温变化型,重点研究其与西太平洋副热带高压等环流异常的重要联系(管兆勇等,2010;蔡佳熙等, 2011;蔡佳熙和管兆勇,2011)。此外,黄嘉佑(1991)还发现夏季气温与降水场在长江中下游、华南及华北具有较好的反相关,还有研究分别从北半球雪盖(王光宇和曾群柱,1994)、印度洋及中东太平洋海温(杨霞等,2007)、澳大利亚高压(张蓬勃等,2010)等海洋及环流因子出发,探讨与气温、降水的关系及可能联系。实际上,夏季气温和降水关系较复杂,很多年份并没有表现为反相关。例如1998年夏季,我国大部分地区降水异常偏多,松嫩流域和长江流域发生流域性洪涝,但全国大部气温仍然明显偏高。说明我国夏季气温的变化有其特殊性,值得作为独立的变量深入分析。此外,过去的研究侧重于整个夏季平均或者某月气温的变化特征或者某个区域的变化特征,没有全面考虑我国夏季各月气温距平的持续性空间模态演变特征及其对应的环流型。因此本文试图从我国夏季季内气温变化的主要模态及造成各模态空间分布的环流特征角度,探讨夏季气温型演变和发展的规律,为进一步诊断和预测全国夏季气温异常分布型提供参考。
1 资料与方法本文所用资料有:中国气象局气象信息中心整编的全国284站逐月气温资料;NCEP/NCAR再分析资料(Kalnay et al, 1996),包括月平均位势高度场、U/V风场,水平分辨率为2.5°×2.5°,场资料的垂直方向分为17层。所用资料时段均为1961年1月至2014年12月共54年;文中气候均值采用1981—2010年共30年平均。
主要采用经验正交函数分解(empirical orthogo-nal function,EOF)(Ghil et al,2002;魏凤英,2007)、合成分析、相关分析、t检验(魏凤英,2007)等方法。其中t检验是构造检验两个总体的均值有无显著性差异的统计量,其中x和y代表样本均值,s1和s2代表样本标准差,n1和n2代表样本量,在确定性水平α之后,遵从自由度t=n1+n2-2的分布。若t的绝对值不小于tα,则拒绝原假设,认为在α这个显著性水平上,两组样本值有显著性差异。
2 夏季气温变化的分型及环流特征选取全国284站作为研究站点,为认识6、7、8月这3个月气温在时间上持续的空间分布特征和时间演变规律,以全国站点同年6、7、8月气温距平场联合扩展成为该年夏季的气温距平场,将1961—2014年历年构建的夏季逐月气温距平场进行EOF分解,所得到的空间模态即包括6、7、8月这3个月的空间分布结构。经North检验(North et al,1982),前3个特征值可分离,前3个特征模态的方差贡献率分别为22.0%、9.3%和7.0%,同时对前3个模态的时间系数进行标准化处理。为客观地判断与各空间模态相似的年份,因此将历年夏季6、7、8月气温联合构造的实况距平分布场与EOF前3个模态分别进行空间相似分析,以相似系数是否超过+/-0.4作为是否达到显著性水平检验的标准,同时结合标准化后的时间系数PC是否通过+/-1个标准差来作为挑选异常年份的标准。
2.1 第一模态(UM型)及其环流特征第一模态逐月的空间分布(图 1)表现出全国同位相变化特征,简称为UM型(uniform mode)。从空间型来看,异常大值主要位于北方,包括西北大部、内蒙古大部、华北及东北大部;从时间来看,全国一致的异常自6月持续至8月。UM模态对应的时间系数为正(负)异常时,我国夏季逐月气温自北向南一致偏高(低),其中北方偏高(低)的变率更为明显。UM型对应的标准化时间系数序列(PC1)看出,UM型存在明显的线性趋势,同时具有明显的年代际变化特征(图 2a)。从20世纪60年代到21世纪,PC1由负转正,其中在90年代中期以前PC1多为负值,意味着全国为一致偏冷型;而90年代中期以后多为正值,为全国一致持续偏高型。该演变与全球变暖的大背景以及80年代后期气温突变的特征吻合(向辽元和陈星,2006)。尤其是西北地区东部、内蒙古、华北、东北西南部等地增暖的趋势明显,长江中下游流域在6、7月增暖明显,8月增暖相对较弱。由历年夏季各月气温场与第一模态空间分布场的相关系数(图 2b)可知,近10年一致增暖型的特征比较典型。
依据异常年的挑选标准,得到正位相年有2001、2002、2005、2011年,负位相年有1969、1970、1976及1979年,可见这两类异常年处于完全不同的年代际阶段。
夏季逐月气温EOF分解的不同模态对应着不同的环流形势。图 3给出了UM型对应的时间系数PC1回归的6、7、8月500 hPa高度场。自6—8月东亚大部多为正值,其中6月东亚通过显著性水平检验的区域主要有贝加尔湖以东至日本群岛、以南至西太平洋的广大区域(图 3a)。7月显著的区域为贝加尔湖以南至华南的广大地区(图 3b),8月回归分布与7月类似,只是显著区域南扩至南亚和东南亚地区(图 3c)。
以6、7、8月都通过显著性水平检验的东亚地区作为重点分析对象,图 4a为PC1回归的(20°~55°N、90°~120°E)范围内高度距平的气压-时间剖面图。可以看出,UM型为正位相时,从对流层低层到平流层低层均为正高度距平,并且自2月持续到10月,在夏季持续的特征尤其显著,表明控制中国大范围的高度场异常为深厚正压系统,大陆高压较为明显,对应的我国大部气温偏高。PC1回归的500 hPa上90°~120°E平均高度距平的纬度-时间剖面图(图 4b)及30°~60°N平均高度距平的时间-经度剖面图(图 4c)显示,从春季到夏季,东亚区域的距平异常存在三个特征,一是由高纬向中纬传播,二是夏季时段内从西向东传播,三是在夏季时段内基本维持正(负)异常不变。
图 5是第二模态的空间分布,从时间上来看,6—8月都呈现出中部与南、北反相的三极子型分布特征,简称TM型(tripole mode)。TM处于正位相时,对应时间系数为正,长江流域不仅是该型的负异常中心,也是该型异常值最大的区域。8月异常值最大,且中间负异常范围最广,气温偏低;对应着东北、内蒙古、西北大部以及江南、华南气温为正异常,即气温偏高。TM型为负位相时则相反。TM型的时间系数PC2变化(图 6a)显示,TM型同样具有明显的年代际变化特征,PC2在20世纪80年代以前多为负值,以后则多为正值。PC2的5年滑动平均表现出这种年代际转折的特征,尤其在近几年振幅加大,表明TM型年际变化特征在增强。相比图 2b中UM型明显的年代际变化特征,TM型的相关系数与历年夏季各月气温场的相关表现出很强的年际变化特征(图 6b)。即80年代末期以后,虽然全球存在着明显的增暖,但是对中国夏季气温而言,这种增暖更突出的是体现在中国黄河以北地区和华南大部,而以江淮为中心的中部地区自6—8月存在持续的变冷现象。
依据异常年挑选标准,得到TM正位相年有1980、1989、1993、1999、2005、2014年,TM负位相年有1964、1966、1967、1973、1978、2013年。2013年夏季,长江流域出现持续高温天气;而2014年夏季长江流域为凉夏,其中6月气温偏低的区域位于长江上游东部、中游及长江下游北部,7月位于长江中下游,而8月低温范围相对7月低温区有所扩展,且长江中下游气温偏低程度加强,在流域范围内110°E以东的平均气温居于历史倒数第二位。这也说明了在气候变暖背景下,气温的年际变化差异仍然较大,尤其是TM型在近些年出现频繁,年际变化明显(图 6b)。
图 7为TM型正/负位相年的6、7、8月气温差值t检验。6月显著的正差异区域包括新疆西北部、内蒙古大部、西南南部及华南大部,显著的负差异区域为江汉地区;7—8月与6月相似,但显著的负差异中心随时间的推移范围明显扩大,7月主要包括长江中下游地区,8月扩至长江上游中东部及黄淮地区。
第二模态TM型对应的时间系数PC2回归的500 hPa高度距平场(图 8)及TM正(负)位相年850 hPa矢量风距平差值t检验场(图 9)反映出:TM正位相的6月,副热带地区的高度场整体偏高,在贝加尔湖以南为显著的异常正值中心(图 8a),对应850 hPa为反气旋式距平环流(图 9a),长江上中游位于其南侧,受反气旋东侧南下的偏北气流影响,气温易随之下降,内蒙古等北方地区受反气旋环流控制,气温偏高,华南沿海一带位于副热带高压的北侧,受副高南侧偏暖气流影响,气温也容易偏高。
至7月,以内蒙古为中心的正异常增强并东扩至鄂霍次克海,北扩至贝加尔湖,中心较6月略向东移动(图 8b),低层该反气旋环流控制范围相对6月明显扩大、东移,反气旋最东侧延伸至鄂霍次克海,偏东北气流经日本群岛、渤海至黄淮及长江流域,偏北气流控制的位置相对东移,主要影响长江中下游,此外我国东北地区的东南部也受偏北气流影响,气温下降。长江中下游至日本群岛呈现出明显的负高度异常,低层为气旋式异常环流。副热带地区维持正异常,低层在菲律宾以东为反气旋式异常环流。总体来看,7月自西太平洋、东亚沿岸、蒙古地区存在“A(反气旋)”—“C(气旋)”—“A”的异常环流(图 9b)波列,这种分布减弱了热带季风槽,导致西太平洋副热带高压偏南,有利于冷空气沿东路南下影响长江流域,造成该地区低温的出现。
至8月,中高纬的正异常范围略北缩,长江中下游至日本群岛的负异常增强且范围向南、向北扩展,受此影响,西太平洋副热带高压略有南退(图 8c)。对应的低层气旋式距平环流加强,偏北气流增强,夏季风明显偏弱,使得低温范围相对7月扩大,低温程度加深(图 9c)。总体说来,8月东亚地区的“A”—“C”—“A”异常环流波列仍然存在,且经向度加大。
PC2回归的200 hPa高度场同图 8中的分布较为一致,但经日本群岛至江淮地区显著的负异常中心范围更广,西扩至新疆和高原北部。鉴于东亚区域的这一负异常(对应着低层的气旋性环流)在东亚地区南北向“A”-“C”-“A”异常环流型中起着重要作用,下面将重点分析该区域的时空演变特征。
从(30°~40°N、110°~140°E)的时间-高度剖面图及在110°~140°E上的高度-纬度剖面图(图 10)可以看出,自3月开始在平流层中上层已经有高度异常出现,并向对流层下传,至夏季,虽然6月的强度变化不显著,但7—8月的深厚特征非常清晰。6月高度差值的减弱是因为负差值区在6月相对偏北,至7月负差值区在热带地区的平流层及30°~60°N平流层底层至对流层,通过显著性水平检验的区域主要是位于80°~90°N的高纬地区对流层的正异常、30°~60°N平流层底层至对流层的负异常及10°~20°N的对流层中下层的正异常(图 10b),这种自高纬到低纬的“+-+”差值分布与图 8、图 9的环流型结论较为一致。此外,通过这种高度-时间差值剖面分布图,可知前期(30°~40°N、110°~140°E)区域春季平流层至对流层高层高度场的高低异常可作为后期夏季各月持续南北高(低)、中间低(高)环流模态的前兆因子。
通过以上分析说明中国夏季6、7、8月气温的TM型对应500 hPa高度环流场从南向北也为三极子分布型,即气温TM正(负)位相时,东亚上空环流自上至下为深厚的相当正压结构,配合蒙古地区的阻塞形势、西太平洋的副热带高压强弱一起构成“+-+”(“-+-”)的东亚-太平洋遥相关型波列分布。
2.3 第三模态(DM型)及其环流特征图 11是第三模态的空间分布, 从6—8月较为一致呈现出我国南北地区反位相的偶极型分布,简称DM型(tripole mode)。从空间分布来看,DM为正位相时,在6月仅东北和内蒙古、华北北部气温为负异常,其他大部为正异常;至7和8月,负异常区南扩至长江,以长江为界,南正北负,尤其是在7月华南为正异常极大值中心。DM为负异常时,特征相反。DM型的时间系数PC3变化及5年滑动平均(图 12)显示,在20世纪70年代末以前时间系数为负,80—90年代转为正值,90年代至2000年转为负值,而在近些年又处于正值,经历了四个阶段的变化(图 12a)。时间系数的变率显示在70年代以前年际变化小,而在90年代至2010年期间振荡明显,年际变化大。气温实况场与第三模态空间分布场的场相关系数图(图 12b)显示,通过显著性水平检验的年份多在20世纪80年代以后。依据异常年挑选标准,得到DM正位相年有1979、1993、2003、2009年,DM负位相年有1974、1994、1997、1999、2000、2001年。下面分析其对应的环流特征。
从第三模态DM型对应的时间系数PC3回归的500 hPa高度距平场可看出(图 13),6月,我国东北地区为明显的负高度距平控制,我国除东北以外的其他地区高度场虽不显著但为正值,这种分布下东北冷涡活跃,东北、内蒙古东部容易形成低温而我国其他大部气温偏高;至7月,东北地区的负高度异常控制范围略有南移,我国江南至南海的正高度异常开始显著;8月与7月的情况较为类似,但东北负值的显著性程度减弱,中心值略有北移。DM型正位相对应的东亚地区为北低南高的环流特征,主要的影响系统是东北冷涡的位置和强度以及西北太平洋高度场(副热带高压)的强度。
根据逐月气温资料、NCEP/NCAR再分析资料,利用EOF分解方法、场相关分析、合成分析及回归分析等统计分析方法,将中国夏季气温可分为三个主要的模态:
第一模态为全国一致型(UM),解释方差为约为22%,UM正位相,对应时间系数为正,6—8月全国气温一致偏高,具有明显的线性趋势,与全球变化背景一致,其中以西北地区东部、内蒙古、华北、东北西南部变暖趋势明显,同时还具有显著的年代际变化特征,在20世纪90年代中期以前位一致偏冷型,90年代中期以后位一致偏暖型。UM正位相对应同期环流为大陆高压较强,在东亚大范围地区从对流层低层到平流层低层的高度场均为正异常,呈深厚正压系统,该异常从春季到夏季时段由高纬向中纬传播,并在夏季时段内从西向东传播。该特征可作为一致型模态的前期信号。
第二模态为全国南、北与中间反位相的三极子型(TM),解释方差约为9%,TM正位相,对应时间系数为正,6—8月东北、内蒙古、西北大部以及江南、华南气温偏高,而长江流域、黄淮等地气温偏低,其中以8月低温范围最大。TM型具有明显的年代际和年际变化特征,20世纪80年代以后以正位相特征为主,但近些年际信号较强且振荡明显,即全球变暖背景下,中国黄河以北地区和华南大部是增暖的,但以江淮为中心的中部地区自6—8月存在持续的变冷现象。TM正位相时同期环流特征为西太平洋、东亚沿岸、蒙古地区的“A”—“C”—“A”异常环流波列,并且(30°~40°N、110°~140°E)区域的平流层至对流层高层高度场的异常在初春的平流层中上层已经出现,由春到夏向对流层下传,可作为夏季各月持续南北高(低)、中间低(高)环流模态的前兆因子。
第三模态为我国南北反位相的偶极型(DM),解释方差约为7%,DM正位相时,对应时间系数为正,6月东北地区、内蒙古及华北北部气温偏低,其他大部偏高,至7、8月低温范围明显南扩,但以长江为界北低南高。同期环流特征,6月我国东北大部为明显负高度异常,随时间发展至7、8月负异常略有扩展和南移,而江南、华南地区为正高度异常。这种北低南高的环流分布显示东北冷涡活跃、西太平洋副热带高压偏强,东北、内蒙古东部容易形成低温而我国其他大部气温偏高。
本文对夏季季内月气温前3个模态及对应的环流特征、前兆信号进行了初步分析,未来还需继续深入研究前兆信号的传播和影响机理及预测方法。
蔡佳熙, 管兆勇, 2007. 我国南方地区夏季低温变化特征及其成因[J]. 南京气象学院学报, 30(6): 799-806. |
蔡佳熙, 管兆勇, 2011. 长江流域夏季气温变化型及其成因Ⅰ:年际变化与遥相关[J]. 气象学报, 69(1): 99-111. DOI:10.11676/qxxb2011.009 |
蔡佳熙, 管兆勇, 于田田, 等, 2011. 长江流域夏季气温变化型及其成因Ⅱ:年代际变化[J]. 气象学报, 69(1): 112-124. DOI:10.11676/qxxb2011.010 |
崔童, 王东阡, 李多, 等, 2015. 2014年夏季我国气候异常及成因简析[J]. 气象, 41(1): 121-125. DOI:10.7519/j.issn.1000-0526.2015.01.015 |
丁裕国, 江志红, 1999. 中国近50年严冬和冷夏演变趋势与区划[J]. 应用气象学报, 10(S1): 88-96. |
管兆勇, 蔡佳熙, 唐卫亚, 等, 2010. 长江中下游夏季气温变化型与西太平洋副高活动异常的联系[J]. 气象科学, 30(5): 666-675. |
韩荣青, 高辉, 李维京, 2014. 旋转经验正交函数分解回归方法在东北夏季气温季节预测和成因诊断中的应用[J]. 气象学报, 72(2): 291-305. DOI:10.11676/qxxb2014.015 |
黄嘉佑, 1991. 我国夏季气温、降水场的时空特征分析[J]. 大气科学, 15(3): 124-132. |
江琪, 马学款, 王飞, 2016. 2016年1月大气环流和天气分析[J]. 气象, 42(4): 514-520. DOI:10.7519/j.issn.1000-0526.2016.04.016 |
李维京, 李怡, 陈丽娟, 等, 2013. 我国冬季气温与影响因子关系的年代际变化[J]. 应用气象学报, 24(4): 385-396. DOI:10.11898/1001-7313.20130401 |
秦爱民, 钱维宏, 蔡亲波, 2005. 1960—2000年中国不同季节的气温分区及趋势[J]. 气象科学, 25(4): 338-345. |
司东, 李清泉, 柳艳菊, 等, 2014. 2013/2014年东亚冬季风异常偏弱的可能成因[J]. 气象, 40(7): 891-897. DOI:10.7519/j.issn.1000-0526.2014.07.014 |
司东, 马丽娟, 王朋岭, 等, 2016. 2015/2016年冬季北极涛动异常活动及其对我国气温的影响[J]. 气象, 42(7): 892-897. DOI:10.7519/j.issn.1000-0526.2016.07.013 |
孙建奇, 王会军, 2006. 东北夏季气温变异的区域差异及其与大气环流和海表温度的关系[J]. 地球物理学报, 49(3): 622-671. |
孙林海, 赵振国, 2004. 我国暖冬气候及其成因分析[J]. 气象, 30(12): 57-60. |
唐国利, 丁一汇, 2006. 近44年南京温度变化的特征及其可能原因的分析[J]. 大气科学, 30(1): 56-68. |
唐恬, 金荣花, 彭相瑜, 等, 2014. 2013年夏季我国南方区域性高温天气的极端性分析[J]. 气象, 40(10): 1207-1215. DOI:10.7519/j.issn.1000-0526.2014.10.005 |
王东阡, 崔童, 司东, 等, 2015. 2014/2015年东亚冬季风活动特征及其可能成因分析[J]. 气象, 41(7): 907-914. DOI:10.7519/j.issn.1000-0526.2015.07.013 |
王东阡, 周兵, 孙丞虎, 等, 2013. 2012/2013年东亚冬季风活动特征及其可能成因分析[J]. 气象, 39(7): 930-937. DOI:10.7519/j.issn.1000-0526.2013.07.014 |
王光宇, 曾群柱, 1994. 北半球雪盖与我国夏季降水的相关分析[J]. 冰川冻土, 16(2): 181-184. |
魏凤英, 2007. 现代气候统计诊断与预测技术: 2版[M]. 北京: 气象出版社.
|
伍红雨, 潘蔚娟, 王婷, 2014. 华南冬季气温异常与ENSO的关系[J]. 气象, 40(10): 1230-1239. DOI:10.7519/j.issn.1000-0526.2014.10.007 |
向辽元, 陈星, 2006. 近55年中国大陆气温突变的区域特征和季节特征[J]. 气象, 32(6): 44-47. DOI:10.7519/j.issn.1000-0526.2006.06.007 |
杨霞, 管兆勇, 朱保林, 2007. IOD对ENSO影响中国夏季降水和气温的干扰作用[J]. 南京气象学院学报, 30(2): 170-177. |
曾宇星, 封国林, 赵俊虎, 等, 2014. 春季马斯克林高压对中国东北地区夏季气温的影响[J]. 高原气象, 33(5): 1374-1382. DOI:10.7522/j.issn.1000-0534.2013.00078 |
张蓬勃, 管兆勇, 蔡佳熙, 2010. SVD分析揭示的澳大利亚高压年际变化对中国夏季气温的可能影响[J]. 大气科学学报, 33(1): 58-66. |
张茜, 李栋梁, 2012. 欧亚和我国东北冬春季积雪对东北夏季气温的影响[J]. 冰川冻土, 34(2): 284-295. |
Chen Wei, Hong Xiaowei, Lu Riyu, et al, 2016. Variation in summer surface air temperature over Northeast Asia and its associated circulation anomalies[J]. Adv Atmos Sci, 33(1): 1-9. DOI:10.1007/s00376-015-5056-0 |
Ghil M, Allen M R, Dettinger M D, et al, 2002. Advanced spectral methods for climatic time series[J]. Rev Geophys, 40(1): 1-41. |
Kalnay E, Kanamitsu M, Kistler R, et al, 1996. The NCEP/NCAR 40-year reanalysis project[J]. Bull Am Meteor Soc, 77(3): 437-471. DOI:10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2 |
North G R, Bell T L, Cahalan R F, et al, 1982. Sampling errors in the estimation of empirical orthogonal functions[J]. Mon Wea Rev, 110(7): 699-706. DOI:10.1175/1520-0493(1982)110<0699:SEITEO>2.0.CO;2 |