Global spectra of energy and enstrophy and their fluxes during July 1979

Transient and stationary spectra of kinetic energy (KE), available potential energy (APE) and enstrophy (EN), and their spectral fluxes as a function of the two-dimensional wavenumbern were computed for July 1979. Triangular truncation at zonal wavenumber 42 was used for computation. The slopes of various spectra in the wavenumber range 14≤n≤25 were obtained by fitting a straight line in log-log scale by the least square method. The transientKE, APE andEN spectra in the lower (upper) troposphere had slopes −2·21 (−2·30), −2·65 (−2·64) and −0·36 (−0·46), respectively. The effect of stationary and divergent motion on the slope values was investigated. The possible correlation between the slope and percentage of transient component in the combined energy and enstrophy was examined to identify the transient motion of the atmosphere with the two-dimensional homogeneous isotropic turbulence. The vertically averaged slope of kinetic energy and enstrophy in the lower (upper) troposphere was close to the value at 700 (200) hPa level. The spectral fluxes of kinetic energy and enstrophy in the wavenumber range 14≤n≤25 satisfied, to a very rough approximation, the criteria of inertial subrange. The stationary fluxes were small. The estimated stationary-transient component of flux was larger, comparable and less than the corresponding transient flux of APE, KE and EN. Representative levels for computation of energy and enstrophy spectra and their fluxes in the lower and upper troposphere were identified.     

PRELIMINARY STUDY OF BAROTROPIC STOCHASTIC DIFFERENCE MODEL APPLIED TO WEATHER PREDICTION

In order to consider both the deterministic and the stochastic property of atmospheric motion simul-taneously,in this paper,the weather prediction is proposed as the problem of the evolution of meteorologicalfield.The historical viewpoint of atmospheric motion is emphasized here.Based on time series analysis te-chnique,a stochastic-dynamical model with multiple initial fields is derived.Thus,weather forecasting is sum-meal up as a problem of solving a set of stochastic difference equations.For the barotropic atmosphere,thenumerical solutions of the equations are obtained by using the method of empirical orthogonal functions(EOF),and examples of medium-range weather prediction are given here.Meanwhile,selecting the order oftime series,i.e.,determining the number of initial fields properly,is also discussed.     

北半球冬季遥相关型环流异常的若干性质

用１９５１－１９８５年间北半球冬季７００，５００，１００ｈＰａ月平均位势高度场，分析了Ｗａｌｌａｃｅ等提出的５个遥相关型所对应的环流异常和它们的性质，发现它们在对流层中下部是相当正压结构，在对流层上部减弱消失，某些遥相关型（如ＥＵ，ＰＮＡ）环流异常存在明显；主要遥相关型环流异常的持续时间约为一个月。     

大气中的包络Rossby孤立波——纬度效应的作用

本文进一步对旋转正压大气中的包络Rossby孤立波进行了研究，结果发现当Rossby波的波数m（m为纬向波数）满足1≤m≤2时，旋转正压大气中才存在包络Rossby孤立波，并且只有m＝2的包络Rossby孤立波才具有阻塞高压的结构。我们还对这种m＝２的包络Rossby孤立波随纬度变化的持续性进行了计算，得到了许多结果。     

COMPARISON OF ON-LINE COUPLED AND CONSTANT TRANSFER SIMULATION METHODS FOR DIRECT RADIATIVE FORCING OF ANTHROPOGENIC SULFATE

1 INTRODUCTION Of three main methods for studying the radiative forcing of anthropogenic sulfate and climatic response on the regional scale, the first is, with given rates for transforming SO2 to sulfate, converting actually released SO2 into sulfate and acquiring the distribution of sulfate by computing transfer equations in the climate model.     

Role of barotropic,baroclinic and combined barotropic-baroclinic instability for the growth of monsoon depressions and mid-tropospheric cyclones

A detailed barotropic, baroclinic and combined barotropic-baroclinic stability analysis has been carried out with mean monsoon zonal currents over western India, eastern India and S.E. Asia. The lower and middle tropospheric zonal wind profiles over western India are barotropically unstable. The structure and growth rate of these modes agree well with the observed features of the midtropospheric cyclones. Similar profiles over eastern India and S.E. Asia, however, are barotropically stable. This is attributed to weak horizontal shear, inherent to these profiles. The upper tropospheric profiles, on the other hand, are barotropically unstable throughout the whole region. The features of these unstable modes agree with those of observed easterly waves. The baroclinic and combined barotropic-baroclinic stability analyses show that the baroclinic effects are not important in tropics. Though the barotropic instability of the mean zonal current seems to be res ponsible for the initial growth of the mid-tropospheric cyclones, neither barotropic nor baroclinic instability of the mean zonal current seem to explain the observed features of the monsoon depressions.     

Two modes of the silk road pattern and their interannual variability simulated by LASG/IAP AGCM SAMIL2.0

In this study, two modes of the Silk Road pattern were investigated using NCEP2 reanalysis data and the simulation produced by Spectral Atmospheric Circulation Model of IAP LASG, Version 2 (SAMIL2.0) that was forced by SST observation data. The horizontal distribution of both modes were reasonably reproduced by the simulation, with a pattern correlation coefficient of 0.63 for the first mode and 0.62 for the second mode. The wave train was maintained by barotropic energy conversion (denoted as CK) and baroclinic energy conversion (denoted as CP) from the mean flow. The distribution of CK was dominated by its meridional component (CK _y ) in both modes. When integrated spatially, CK _y was more efficient than its zonal component (CK _x ) in the first mode but less in the second mode. The distribution and efficiency of CK were not captured well by SAMIL2.0. However, the model performed reasonably well at reproducing the distribution and efficiency of CP in both modes. Because CP is more efficient than CK, the spatial patterns of the Silk Road pattern were well reproduced. Interestingly, the temporal phase of the second mode was well captured by a single-member simulation. However, further analysis of other ensemble runs demonstrated that the successful reproduction of the temporal phase was a result of internal variability rather than a signal of SST forcing. The analysis shows that the observed temporal variations of both CP and CK were poorly reproduced, leading to the low accuracy of the temporal phase of the Silk Road pattern in the simulation.     

The Neutral, Barotropic Planetary Boundary Layer, Capped by a Low-Level Inversion

The presence of a low-level, capping inversion layer will affect the height and structure of the planetary boundary layer (PBL). Results from models of varying levels of sophistication, including analytical, turbulent kinetic energy (TKE), second-order closure (SOC), large-eddy simulation (LES) and direct numerical simulation (DNS) models, are used to investigate this influence for the neutral, barotropic PBL. Predicted and observed profiles of stress and geostrophic departure components, and integral measures, such as the parameters of Rossby-number similarity theory, are compared for the KONTUR, Marine Stratocumulus, JASIN, Leipzig, Pre-Wangara and Upavon field experiments.Analytical models of the equilibrium value of inversion height z_i, which depend on the surface friction velocity u_*, and both the Coriolis parameter f and the free-flow Brunt-Väisälä frequency N, are found to give reasonable estimates of the PBL height. They also indicate that only the KONTUR and Marine Stratocumulus experiments were strongly influenced by N. More quantitative comparisons would require larger, more comprehensive datasets. The effects of the presence of a capping inversion on the profile structure were found to be insignificant for h_* = |f|z_i/u_* > 0.15.The simple analytical model performed quite well over all values of h_*; it predicted the profiles of the longitudinal stress component (in the direction of the surface stress) better than the lateral component. The more advanced models performed well for small values of h_* (for flow over the sea), but systematically underestimated the cross-isobaric angle for flow over land. These models predicted the profiles of the lateral stress component better than the longitudinal component. The profiles of the analytical model agreed with those of the advanced models when the constant eddy viscosity of the outer layer was increased.Agreement with DNS was achieved by increasing the eddyviscosity of the analytical model by a factor of 5.Zilitinkevich and Esau(2002, Boundary-Layer Meteorology, 104, 371–379)suggest that the neutral, barotropic values of A and B of Rossby-numbersimilarity theory are not universal constants, but depend on the ratio N/|f|. The dependence for A and B is calculated using the analytical model and TKE models. Over the sea (h_* 0.1; N/|f| 100, where we have used the Zilitinkevich-Esau relation to convert between h_* and N/|f|) there is agreement between the model predictions and observations; however over land where the equilibrium boundary-layer height is greater (h_* 0.35; N/|f| 10) the inconsistency between the advanced model predictions (TKE, SOC, LES, and DNS) and observations, as noted previously by Hess and Garratt, still exists. We attribute this disagreement to violations of the strict assumptions of steady, horizontally homogeneous, neutral, barotropic conditions implicit in the observations. At small values of z_i and a strongly stable background stratification (h_* 0.04; N/|f| 1000) both the TKE and analytical models predict that A and B depend significantly on h_*, however observations are unavailable to confirm these predictions. Zilitinkevich and Esau call this case the `long-lived near-neutral PBL', and state that it is found in cold weather at high latitudes.     

近30a来北半球对流层大气月均环流的涡动减弱现象

利用1980-2009年NCEP／NCAR月平均再分析资料,研究了在全球变暖背景下北半球对流层大气涡动减弱现象。结果表明：北半球涡度拟能30a来整体呈减弱趋势,在北太平洋地区和极地减弱尤为显著,12．5～50°N为影响北半球大气涡动变化的关键区域。由于对流层200hPa以下大气的增暖,特别是中高纬地区显著增暖,减弱了经向温度梯度,使得副热带西风急流的强度亦呈减弱趋势。这与涡度拟能的变化有显著的正相关关系。分析了北半球正压Rossby波动诊断量E和热量经向涡动通量的变化,表明中纬度波能传播在不同地区有不同趋势,而热量的经向涡动输送与涡度拟能的变化也呈现显著的正相关关系,减弱明显。研究结果对深刻认识大气环流变化规律和理解全球变暖的可能影响具有重要意义。