冬季北半球风暴轴的协同变化及其与大气环流的相互作用

基于1971—2016年NCEP/NCAR的逐日、逐月再分析资料,研究冬季北半球西伯利亚风暴轴(Siberian Storm Track,SIST)、北太平洋风暴轴(Pacific Storm Track,PST)和北大西洋风暴轴(Atlantic Storm Track,AST)的协同变化特征及其与大气环流的关系。结果表明:(1)三大风暴轴不仅各自的位置与强度变化存在显著相关性,风暴轴之间也存在一定的协同变化且年代际尺度上比年际尺度上更紧密。年际尺度上,SIST与AST的经度变化呈显著负相关,而PST和AST的协同性较差;年代际尺度上,SIST与PST的经、纬度变化均呈弱的负相关,SIST与AST的经度和强度变化均呈显著正相关,PST与AST的经、纬度变化均呈显著负相关。(2)由联合EOF分析得到北半球风暴轴的协同变化时空特征:在年际尺度上,第一模态主要表现为SIST偏弱(强),PST主体偏弱(强)、东南偏强(弱),AST略偏北(南)偏强(弱)但不显著的协同变化。PC1为正位相时,对应的大气环流异常为:500 hPa高度场上为太平洋北美(Pacific North America,PNA)型和欧亚(Eurasian,EU)型的正位相,东亚急流偏强且偏南;第二模态主要表现为SIST偏强(弱)且偏东(西),PST中东部偏南(北)、西部强度偏强(弱),AST偏强(弱)的协同变化。PC2为正位相时,对应的大气环流异常为:500 hPa高度场上为PNA型和大西洋东部(East Atlantic,EA)型的正位相,北美急流减弱;在年代际尺度上,第一模态主要表现为SIST偏西(东)且偏弱(强),PST偏东(西)且偏弱(强),AST偏西(东)且偏弱(强)的协同变化。PC1为正位相时,对应的大气环流异常为:500 hPa高度场上为西大西洋(West Atlantic,WA)型和EU型的正位相。第二模态主要表现为SIST偏强(弱)且偏北(南),PST偏南(北)且偏弱(强),AST北抬(南压)的协同变化。PC2为正位相时,对应的大气环流异常为:500 hPa高度场上为EU型和WA型的正位相,东亚急流强度加强且偏南,北美急流强度减弱。     

Representation of Northern Hemisphere winter storm tracks in climate models

Northern Hemisphere winter storm tracks are a key element of the winter weather and climate at mid-latitudes. Before projections of climate change are made for these regions, it is necessary to be sure that climate models are able to reproduce the main features of observed storm tracks. The simulated storm tracks are assessed for a variety of Hadley Centre models and are shown to be well modelled on the whole. The atmosphere-only model with the semi-Lagrangian dynamical core produces generally more realistic storm tracks than the model with the Eulerian dynamical core, provided the horizontal resolution is high enough. The two models respond in different ways to changes in horizontal resolution: the model with the semi-Lagrangian dynamical core has much reduced frequency and strength of cyclonic features at lower resolution due to reduced transient eddy kinetic energy. The model with Eulerian dynamical core displays much smaller changes in frequency and strength of features with changes in horizontal resolution, but the location of the storm tracks as well as secondary development are sensitive to resolution. Coupling the atmosphere-only model (with semi-Lagrangian dynamical core) to an ocean model seems to affect the storm tracks largely via errors in the tropical representation. For instance a cold SST bias in the Pacific and a lack of ENSO variability lead to large changes in the Pacific storm track. Extratropical SST biases appear to have a more localised effect on the storm tracks.     

Quasi-stationarity of centennial Northern Hemisphere midlatitude winter storm tracks

The winter storm activity on the Northern Hemisphere during the last one thousand years in a global climate simulation was analyzed by determining all midlatitude storms and their tracks, then consecutively clustering them for hundred years’ segments. Storm track clusters with longest lifetime and largest deepening rates are found over the oceans. The numbers of extratropical winter storms exhibit notable yearly variability but hardly any variability on centennial time scales. The clusters of these storm tracks also show only small differences between the centuries. The numbers of members in neighboring oceanic clusters are negatively correlated. A linear relationship was found between the numbers of members per storm track clusters over the Pacific or Atlantic Ocean and seasonal mean atmospheric circulation patterns by a canonical correlation analysis.     

Ocean mixed layer processes in the Pacific Decadal Oscillation in coupled general circulation models

It is investigated how the Pacific Decadal Oscillation (PDO) is simulated differently among various coupled general circulation models (CGCMs), and how it is related to the heat budget of the simulated ocean mixed layer, which includes the surface heat flux and ocean heat transport. For this purpose the dataset of the climate of the 20th Century experiment (20C3M) from nine CGCMs reported to IPCC’s AR4 are used, while the MRI and MIROC models are examined in detail. Detailed analyses of these two CGCMs reveal that the PDO is mainly affected by ocean heat transport rather than surface heat flux, in particular in the MRI model which has a larger contribution of ocean heat transport to the heat budget. It is found that the ocean heat transport due to Ekman advection versus geostrophic advection contributes differently to the PDO in the western and central North Pacific. Specifically, the strength of PDO tends to be larger for CGCMs with a larger ocean heat transport in the region.     

北半球冬季风暴轴与ENSO循环的关系及其能量平衡特征

利用欧洲中期天气预报中心的再分析资料,分析了北半球冬季对流层上层300 hPa面上风暴轴与ENSO循环的关系以及风暴轴的能量平衡特征.研究表明,由于El Nino年北太平洋的哈德莱环流增强,导致北太平洋西风急流向赤道和向东伸展,进而引起了北太平洋风暴轴的增强并向赤道和向东伸展;而La Nina年事件期间情形正好相反.在北大西洋,El Nino年其西风急流中心最大值有所减小,但整个西风急流区域有所扩大并向西和向赤道伸展,相应北大西洋风暴轴强度在El Nino年也有所减弱并向西和向赤道方向伸展;La Nina事件的情形正好相反.能量分析表明,斜压转换的正值中心位于风暴轴及其上游区域.在El Nino年,北太平洋斜压转换的正值中心向赤道和向下游伸展,而北大西洋斜压转换的正值中心向赤道和向西伸展,这与北太平洋风暴轴和北大西洋风暴的变化趋势是一致的.行星尺度扰动对斜压转换的贡献比天气尺度扰动的贡献要小,而行星尺度-天气尺度扰动的相互作用项对斜压转换的贡献更小.研究还表明,能量正压转换的正中心化于风暴轴的上游,负中心位于风暴轴区域.相对于La Nina年,El Nino年北太平洋东部的正压转换负值中心偏北,北美的正值中心位置偏南,而北大西洋的负值中心强度减弱范围变小.这些变化与从La Nina年剑El Nino年北太平洋和北大西洋风暴轴的变化是一致的.在能量的正压转换过程中,天气尺度扰动的贡献要大于行星尺度扰动,而行星尺度-天气尺度扰动的相互作用项的贡献更小.     

欧亚土壤湿度异常对北半球大气环流的显著影响

用44 a ERA40再分析资料的土壤湿度和大气环流变量场, 研究持续性的欧亚大陆土壤湿度异常对后期北半球大尺度大气环流的反馈作用。首先,运用经验正交函数分解去除ENSO遥相关及趋势影响后,分析了欧亚大陆中高纬度土壤湿度变率主要模态的季节变化特征,及相对应主分量时间序列显示的土壤湿度异常的衰减时间,结果表明土壤湿度异常的主要模态在全年都表现出很好的连续性。其次,对不同季节的连续3个月的月平均土壤湿度和500 hPa高度场进行滞后最大协方差分析,研究欧亚地区中高纬度土壤湿度异常与北半球大气环流异常之间的线性耦合。第一最大协方差模态的结果表明:全年的主导信号是大气强迫土壤湿度的变化,但在冬季和夏季,大气中类似于负位相北极涛动的环流型与之前月份(最长达4个月)土壤湿度的持续变化显著相关。最后,基于土壤湿度变率中心的回归分析也证实了秋季和春季欧亚土壤湿度,特别是北非副热带,欧亚内陆和西伯利亚地区的土壤湿度异常,分别与其后的冬季和夏季的大气环流显著相关。欧亚大陆土壤湿度异常超前大气环流的信号,将有助于改善冬季和夏季北半球季节气候预报能力。     

冬季北半球风暴轴的多尺度能量变化特征及其可能机制

基于1971—2016年NCEP/NCAR(美国环境预报中心和国家大气研究中心)的逐日再分析资料及NCPC(美国国家海洋和大气管理局气候预报中心)的海温、大气环流及海洋指数等资料通过多尺度能量分析(MS-EVA)等方法,把冬季北半球风暴轴看做一整体,分析了风暴轴区域多尺度的能量变化特征及其可能机制。主要结论概括如下:(1)多年气候平均状态下,风暴轴的动能来源主要表现为在风暴轴中上游先由低频尺度向天气尺度输送有效位能,随后在风暴轴主体区再由天气尺度有效位能转换为天气尺度动能,其中风暴轴西端可直接由低频尺度向天气尺度输送动能。(2)北半球三大风暴轴联合EOF结果表明:第一模态下,主要体现了北西伯利亚风暴轴与北太平洋风暴轴强度的减弱(增强),同时伴随着北大西洋风暴轴位置北抬(南压);第二模态下,主要体现了北西伯利亚风暴轴强度减弱(增强),同时北太平洋风暴轴位置北抬(南压)中东部强度增强(减弱),而北大西洋风暴轴位置南压(北抬)。(3)回归分析表明:北半球风暴轴异常在不同模态下与低频尺度环流联系密切。低频尺度波动可通过海温及西风急流等异常变化先影响风暴轴区域多尺度间的能量转换,进而影响风暴轴整体的异常变化。     

Numerical simulation of Eurasian teleconnection pattern in atmospheric circulation during the Northern Hemisphere winter

In this paper, the anomaly of disturbance height field over Northern Hemisphere due to SST anomaly in the tropical Atlantic Ocean is simulated by using the general circulation model of IAP. A quasi-geostrophic, 34-level spherical coordinate model is also used to compute the distribution of atmospheric circulation anomaly when there is an anomaly of heat source over the tropical Atlantic. The computed results show that, owing to the heat source anomaly over the tropical Atlantic, the EU-pattern anomaly in the winter circulation may be caused. Namely, a ridge will be enhanced over the northwest Europe, a trough will be deepened over Siberia, but a positive anomaly of disturbance height field will be formed over the northeast China, Japan and other areas of East Asia. Moreover, the numerically simulated results show that the above-mentioned EU-pattern anomalies of circulation are due to the propagations of planetary wave train. About 15 days after an anomaly of the heat source over the tropical Atlantic is injected, the EU-pattern anomaly of atmospheric circulation is formed. This is in good agreement with the results analysed theoretically. On the leave from Geophysical Institute, Faculty of Science in Tokyo University, Japan.     

Global warming mode of atmospheric circulation

Variations of global zonal-mean atmospheric circulation are studied using the National Center of Environment Prediction/National Center of Atmosphere Research (NCEP/NCAR) reanalysis data set from January 1948 to February 2001 and CCM3.6 (Community Climate Model Version 3.6) simulations for the period 1900–1999. Empirical orthogonal function (EOF) analysis indicates that variations of zonal-mean geopotential height in the tropics are usually opposite to those over the subpolar zone in both hemispheres. The first mode of EOF analysis shows that height is higher than normal in the tropics when it is lower over the subpolar zone in both hemispheres with much stronger westerlies over mid latitudes in both hemispheres, and vice versa. This mode explains about 50% variance and is predominant in the whole troposphere. The time series of EOF1 has a sharp transition near about 1977 and the polarity changes from negative to positive. This mode is closely related to the variations of global mean surface air temperature. The detrended correlation coefficient between EOF1 time series and the surface air temperature is 0.74 in the boreal winter. Furthermore, the lowest correlation coefficient among the other three seasons, annual mean, seasonal mean and monthly mean is higher than 0.42 which indicates the fairly good relationship between this mode and the surface air temperature. This result has been verified using CCM3.6 simulations.     

Simple general atmospheric circulation and climate models with memory

This article examines some general atmospheric circulation and climate models in the context of the notion of “memory”. Two kinds of memories are defined： statistical memory and deterministic memory. The former is defined through the autocorrelation characteristic of the process if it is random （chaotic）, while for the latter, a special memory function is introduced. Three of the numerous existing models are selected as examples. For each of the models, asymptotic （at t →∞） expressions are derived. In this way, the transients are filtered out and that which remains concerns the final behaviour of the models.     

Modelled atmospheric response to changes in Northern Hemisphere snow cover

The surface boundary conditions are altered in a numerical simulation of January climate by prescribing (a) higher and (b) lower than average snow extent over Northern Hemisphere land masses. The anomalies in snow cover are shown to have quite a strong impact on the mean climatic state. Associated with an increase in the areal extent of the snow, there is a significant reduction in temperature throughout the lower troposphere. There are also large increases in sea-level pressure over most land areas. Significant responses in the mass field are also seen at 500 hPa where reductions in atmospheric thickness lead to significant negative anomalies in the height field. Responses are also seen non-locally, over both the North Pacific and North Atlantic basins. The impact of increased snow on cyclone tracks is also examined. A reduction in cyclones is noted over both continents and over the western sectors of both ocean basins. Over the North Atlantic basin this reduction extends across over Europe, significantly weakening the storm track. In the North Pacific, cyclone density is reduced in the west while in the east, there is actually a strengthening of the storm tracks. There are corresponding changes in the genesis of cyclones in both of these regions. The change in cyclogenesis, intensity and density is demonstrated to be associated with changes in baroclinicity between the two experiments. The anomalous snow boundary conditions lead to significant changes in the meridional temperature gradients over both ocean basins which impact on the baroclinic zones. Received: 5 January 1996 / Accepted: 4 May 1996     

冬季北太平洋风暴轴的年际型态及其与大气环流的关系

利用NCEP/NCAR再分析资料,运用31点带通数字滤波、线性相关和合成分析方法,研究了1961/1962—2010/2011年冬季北太平洋风暴轴西部、东部区域强度指数的年际演变特征,划分了风暴轴的典型型态,并进一步探讨了与同期北半球500 hPa位势高度场和SLP的关系。结果表明:风暴轴气候态的极大值区域位于中纬度北太平洋中西部,最大值点的频数集中区域和均方差分布的异常中心都有两个。风暴轴西部和东部区域强度指数(WI和EI)的年际演变具有独立性,典型型态可分为单、双中心型两类。WI(EI)指数与北半球500 hPa位势高度场的相关分布类似于WP(PNA)遥相关型;单中心型风暴轴偏强时,极涡南扩,平均槽加深;呈双中心型时,极涡明显偏西。WI(EI)指数与SLP的相关分布类似于NPO(NAO)遥相关型;单中心型风暴轴偏强(弱)时,SLP距平场呈AO遥相关型的正(负)异常位相。     

北方涛动同北半球温带大气环流的遥相关(一)——基本结构

本文利用1951—1980年逐季的平均值资料(共120个季)讨论了北方涛动和与其相联系的北太平洋海温与北半球海平面气压场、500hPa位势高度场遥相关的基本结构,并与南方涛动和赤道东太平洋海温的结果进行了对比分析.发现北太平洋Namias海区和加利福尼亚海流区海温的变化与北方涛动具有很密切的联系;北方涛动和这两个海区的海温同北半球中高纬度大气环流特别是PNA型和NAO型环流异常存在明显的遥相关关系;南方涛动和赤道太平洋海温同WP型或NPO型环流异常关系比较密切,而与PNA型和NAO型的关系不如北方涛动和Namias海区及加利福尼亚海流区海温的显著.     

北方涛动同北半球温带大气环流的遥相关(二)——季节变化

在以前工作的基础上,本文进一步分析了各个季节北方涛动和与其相联系的北太平洋海温同北半球温带大气环流遥相关的特征,并与南方涛动和赤道太平洋海温的结果进行了比较,发现它们同北半球温带大气环流的遥相关,不仅具有很强的地区性,而且具有明显的季节变化。 冬季,北方涛动和北太平洋海温与PNA型相关非常明显,且比南方涛动和赤道太平洋海温的更密切。除PNA型外,北方涛动还同500hPa高度的亚洲-北美(ANA)流型有联系。 夏季,PNA遥相关型的表现仍然非常清楚,但位置较冬季偏酉和偏北,并在北美西海岸具有特殊的分叉现象。     

Response of Southern Ocean circulation to global warming may enhance basal ice shelf melting around Antarctica

We investigate the large-scale oceanic features determining the future ice shelf–ocean interaction by analyzing global warming experiments in a coarse resolution climate model with a comprehensive ocean component. Heat and freshwater fluxes from basal ice shelf melting (ISM) are parameterized following Beckmann and Goosse [Ocean Model 5(2):157–170, 2003]. Melting sensitivities to the oceanic temperature outside of the ice shelf cavities are varied from linear to quadratic (Holland et al. in J Clim 21, 2008). In 1% per year CO₂-increase experiments the total freshwater flux from ISM triples to 0.09 Sv in the linear case and more than quadruples to 0.15 Sv in the quadratic case after 140 years at which 4 × 280 ppm = 1,120 ppm was reached. Due to the long response time of subsurface temperature anomalies, ISM thereafter increases drastically, if CO₂ concentrations are kept constant at 1,120 ppm. Varying strength of the Antarctic circumpolar current (ACC) is crucial for ISM increase, because southward advection of heat dominates the warming along the Antarctic coast. On centennial timescales the ACC accelerates due to deep ocean warming north of the current, caused by mixing of heat along isopycnals in the Southern Ocean (SO) outcropping regions. In contrast to previous studies we find an initial weakening of the ACC during the first 150 years of warming. This purely baroclinic effect is due to a freshening in the SO which is consistent with present observations. Comparison with simulations with diagnosed ISM but without its influence on the ocean circulation reveal a number of ISM-related feedbacks, of which a negative ISM-feedback, due to the ISM-related local oceanic cooling, is the dominant one.     

北半球大气环流能量循环的气候特征

利用1958—2011年NCEP/NCAR逐日再分析资料,根据Lorenz能量循环理论框架,分析了北半球大气能量循环的年变化特征,在此基础上给出了更具普适性的多年平均的大气能量循环框图.结果表明:北半球大气能量循环的年变化特征十分明显.大气能量及能量转换率均表现为冬季高、夏季低、春秋季过渡的演变特征;纬向平均有效位能、纬向平均动能和涡动动能中有少许能量在冬季时由南半球向北半球进行越赤道输送,夏季时则由北半球向南半球输送,而涡动有效位能的输送方向则与此相反;纬向平均有效位能的制造在秋季最大,涡动有效位能的制造在夏季最大;动能的耗散冬季最强,夏季最弱.就年平均而言,相较于能量转化过程,能量越赤道交换过程非常微弱.在经向上,纬向平均有效位能主要分布于高纬地区,纬向平均动能主要分布于中低纬地区,而涡动能量主要贮存在中纬和高纬地区;此外,能量转化过程一般在中纬度地区较活跃.     

7月不同形态南亚高压与北半球大气环流的关系

利用NCEP/NCAR月平均高度场和风场再分析资料,分析了1948—2013年7月南亚高压多中心特征及其与北半球大气环流的关系。结果表明:1)南亚高压存在1至5个中心不等,其中以双中心类和三中心类为主,占总样本数的82%,其次是单中心类,占总样本数的14%,四中心类和五中心类仅占总样本数的3%。2)根据高压中心个数、经向位置和环流特征,将不同类高压分成了不同型,其中单中心类分为Ⅰ1型和Ⅰ2型,分别占该类的44%和56%,双中心类分为Ⅱ1型、Ⅱ2型和Ⅱ3型,分别占该类的66.7%、18.5%和14.8%,三中心类仅考虑了Ⅲ1型,占该类的67%。3)Ⅰ1型高压中心在伊朗高原上空,Ⅰ2型高压中心在青藏高原上空,Ⅱ1型两高压中心分别在伊朗高原和青藏高原上空,Ⅱ2型两高压中心分别在伊朗高原和我国东部西太平洋上空,Ⅱ3型两高压中心分别在青藏高原和我国东部西太平洋上空,Ⅲ1型三个高压中心分别在伊朗高原、青藏高原和我国东部西太平洋上空。4)不同类型的高压中心所在地区高层位势高度场和对流层中上层温度场都表现为显著正异常,且不同区域温度场异常的维持机制不相同。     

Predictability of SST forced climate signals in two atmospheric general circulation models

 The predictability of atmospheric responses to global sea surface temperature (SST) anomalies is evaluated using ensemble simulations of two general circulation models (GCMs): the GENESIS version 1.5 (GEN) and the ECMWF cycle 36 (ECM). The integrations incorporate observed SST variations but start from different initial land and atmospheric states. Five GEN 1980–1992 and six ECM 1980–1988 realizations are compared with observations to distinguish predictable SST forced climate signals from internal variability. To facilitate the study, correlation analysis and significance evaluation techniques are developed on the basis of time series permutations. It is found that the annual mean global area with realistic signals is variable dependent and ranges from 3 to 20% in GEN and 6 to 28% in ECM. More than 95% of these signal areas occur between 35 °S–35 °N. Due to the existence of model biases, robust responses, which are independent of initial condition, are identified over broader areas. Both GCMs demonstrate that the sensitivity to initial conditions decreases and the predictability of SST forced responses increases, in order, from 850 hPa zonal wind, outgoing longwave radiation, 200 hPa zonal wind, sea-level pressure to 500 hPa height. The predictable signals are concentrated in the tropical and subtropical Pacific Ocean and are identified with typical El Ni?o/ Southern Oscillation phenomena that occur in response to SST and diabatic heating anomalies over the equatorial central Pacific. ECM is less sensitive to initial conditions and better predicts SST forced climate changes. This results from (1) a more realistic basic climatology, especially of the upper-level wind circulation, that produces more realistic interactions between the mean flow, stationary waves and tropical forcing; (2) a more vigorous hydrologic cycle that amplifies the tropical forcing signals, which can exceed internal variability and be more efficiently transported from the forcing region. Differences between the models and observations are identified. For GEN during El Ni?o, the convection does not carry energy to a sufficiently high altitude, while the spread of the tropospheric warming along the equator is slower and the anomaly magnitude smaller than observed. This impacts model ability to simulate realistic responses over Eurasia and the Indian Ocean. Similar biases exist in the ECM responses. In addition, the relationships between upper and lower tropospheric wind responses to SST forcing are not well reproduced by either model. The identification of these model biases leads to the conclusion that improvements in convective heat and momentum transport parametrizations and basic climate simulations could substantially increase predictive skill. Received: 25 April 1996 / Accepted: 9 December 1996