THE EFFECTS OF TOPOGRAPHY ON THE SUMMER ATMOS-PHERIC ENERGETICS OF THE NORTHERN HEMISPHERE IN A LOW-RESOLUTION GLOBAL SPECTRAL MODEL

An analysis is made of the effects of topography on the summer atmospheric energetics of the Northern Hemisphere in a low-resolution global spectral model. The numerical model is a global, spectral, primitive equation model with five equally spaced sigma levels in the vertical and triangular truncation at wavenumber 10 in the horizontal. The model includes comparatively full physical processes.Each term of the energy budget equations is calculated in four specific latitudinal belts (81.11°S-11.53°S; 11.53°S-11.53°N; 11.53°N-46.24°N; 46.24°N-81.11°N) from a five-year simulation with mountains and a one-year simulation without mountains, respectively. Differences between them are compared and statistically tested. The results show that synoptical scale waves transport available potential energy and kinetic energy to long waves and increase conversion from available potential energy of the zonal flow to eddy’s and from the eddy kinetic energy to the zonal kinetic energy in region 3 (11.53°N-46.24°N) due to mountains; topography intensifies the atmospheric baroclinity in region 3, consequently the baroclinic conversion of atmosphere energy is increased. The seasonal characteristics associated with the summer atmospheric energy source in region 3 are caused by seasonal variation of the solar radiation and the land-ocean contrasts and independent of topographic effects. The mechanism of topographic effects on the increase of long wave kinetic energy is also discussed.     

The effects of topography on the summer atmospheric energetics of the Northern Hemisphere in a low-resolution global spectral model

An analysis is made of the effects of topography on the summer atmospheric energetics of the Northern Hemisphere in a low-resolution global spectral model. The numerical mode! is a global, spectral, primitive equation model with five equally spaced sigma levels in the vertical and triangular truncation at wavenumber 10 in the horizontal. The model includes comparatively full physical processes. Each term of the energy budget equations is calculated in four specific latitudinal belts (81.11°S–11.53°S; 11.53°S–11.53°N; 11.53°N–46.24°N; 46.24°N–81.11°N) from a five-year simulation with mountains and a one-year simulation without mountains, respectively. Differences between them are compared and statistically tested. The results show that synoptical scale waves transport available potential energy and kinetic energy to long waves and increase conversion from available potential energy of the zonal flow to eddy's and from the eddy kinetic energy to the zonal kinetic energy in region 3 (11.53°N-46.24°N) due to mountains; topography intensifies the atmospheric baroclinity in region 3, consequently the baroclinic conversion of atmosphere energy is increased. The seasonal characteristics associated with the summer atmospheric energy source in region 3 are caused by seasonal variation of the solar radiation and the land-ocean contrasts and independent of topographic effects. The mechanism of topographic effects on the increase of long wave kinetic energy is also discussed.     

Changes in the normal mode energetics of the general atmospheric circulation in a warmer climate

Changes in the normal mode energetics of the general atmospheric circulation are assessed for the northern winter season (DJF) in a warmer climate, using the outputs of four climate models from the Coupled Model Intercomparison Project, Phase 3. The energetics changes are characterized by significant increases in both the zonal mean and eddy components for the barotropic and the deeper baroclinic modes, whereas for the shallower baroclinic modes both the zonal mean and eddy components decrease. Significant increases are predominant in the large-scale eddies, both barotropic and baroclinic, while the opposite is found in eddies of smaller scales. While the generation rate of zonal mean available potential energy has globally increased in the barotropic component, leading to an overall strengthening in the barotropic energetics terms, it has decreased in the baroclinic component, leading to a general weakening in the baroclinic energetics counterpart. These global changes, which indicate a strengthening of the energetics in the upper troposphere and lower stratosphere (UTLS), sustained by enhanced baroclinic eddies of large horizontal scales, and a weakening below, mostly driven by weaker baroclinic eddies of intermediate to small scales, appear together with an increased transfer rate of kinetic energy from the eddies to the zonal mean flow and a significant increase in the barotropic zonal mean kinetic energy. The conversion rates between available potential energy and kinetic energy, C, were further decomposed into the contributions by the rotational (Rossby) and divergent (gravity) components of the circulation field. The eddy component of C is due to the conversion of potential energy of the rotational adjusted mass field into kinetic energy by the work realized in the eddy divergent motion. The zonal mean component of C is accomplished by two terms which nearly cancel each other out. One is related to the Hadley cell and involves the divergent component of both wind and geopotential, while the other is associated to the Ferrel cell and incorporates the divergent wind with the rotationally adjusted mass field. Global magnitude increases were found in the zonal mean components of these two terms for the warmer climate, which could be the result of a strengthening and/or widening of both meridional cells. On the other hand, the results suggest a strengthening of these conversion rates in the UTLS and a weakening below, that is consistent with the rising of the tropopause in response to global warming.     

Optimisation of simplified GCMs using circulation indices and maximum entropy production

Two kinds of objective functions for parameter optimisation in simplified general circulation models (SGCMs) are introduced and tested with an SGCM employing linear parameterisations for diabatic heating, surface friction and horizontal diffusion. (a) A set of circulation indices is introduced to characterise the zonal mean primary and secondary circulation and the global energetics. The objective function is then given by the distance between the modelled and a reference (e.g. observed) circulation in a state space spanned by these indices. (b) The global and time mean entropy production and kinetic energy dissipation are introduced as additional objective functions, following the maximum entropy production principle. It is found that both methods lead to optimal parameter values close to the standard configuration of the model, though the method of the second kind is restricted to those model parameters associated with internal processes such as heat and momentum fluxes.     

An approximate energy cycle for inter-member variability in ensemble simulations of a regional climate model

The presence of internal variability (IV) in ensembles of nested regional climate model (RCM) simulations is now widely acknowledged in the community working on dynamical downscaling. IV is defined as the inter-member spread between members in an ensemble of simulations performed by a given RCM driven by identical lateral boundary conditions (LBC), where different members are being initialised at different times. The physical mechanisms responsible for the time variations and structure of such IV have only recently begun to receive attention. Recent studies have shown empirical evidence of a close parallel between the energy conversions associated with the time fluctuations of IV in ensemble simulations of RCM and the energy conversions taking place in weather systems. Inspired by the classical work on global energetics of weather systems, we sought a formulation of an energy cycle for IV that would be applicable for limited-area domain. We develop here a novel formalism based on local energetics that can be applied to further our understanding IV. Prognostic equations for ensemble-mean kinetic energy and available enthalpy are decomposed into contributions due to ensemble-mean variables (EM) and those due to deviations from the ensemble mean (IV). Together these equations constitute an energy cycle for IV in ensemble simulations of RCM. Although the energy cycle for IV was developed in a context entirely different from that of energetics of weather systems, the exchange terms between the various reservoirs have a rather similar mathematical form, which facilitates some interpretations of their physical meaning.     

An application of the E-ε turbulence model for studying coupled air-sea boundary-layer structure

An E- turbulence model is used to study air-sea interaction characteristics and turbulence structure using a coupled model for air-sea boundary layers. The E- turbulence model consists of equations for the turbulent kinetic energy, the energy-dissipation, and for the turbulent exchange coefficient expressed in terms of turbulent kinetic energy and energy-dissipation. The energy-dissipation equations for the air-sea interface are solved analytically to obtain boundary conditions for energy-dissipation at the interface. The air-sea interaction and turbulence characteristics of the E- model are compared with those of the mixing-length model and with available observations.The simulations demonstrate that the air-sea interaction parameters obtained by the E- model agree well with observations. The numerical studies also show that the E- turbulence model with appropriate constants can give good results in modeling coupled air-sea boundary-layer flows.     

On the vorticity and energy budgets of the cold vortex in Northeast China: a case study

This study examines the vorticity budgets, turbulent extended exergy and kinetic energy evolution equations to investigate the major dynamical and energy conversion processes contributing to the initiation and intensification of the cold vortex over Northeast China that occurred during June 19–22, 2009. The results show that the cyclonic vorticity was initiated in the lower troposphere due to the intense convergence of horizontal winds. The growth of cyclonic vorticity in the middle troposphere is mainly due to the vertical transportation of the vorticity, yet the increase of cyclonic vorticity in the upper troposphere primarily results from the horizontal advection of vorticity. Of special interest in this study is the evaluation of the role of thermal advections in the baroclinic development of the cold vortex. The results indicate that the rising of the air over relatively warm areas and the sinking of the air in relatively cold regions are favorable for releasing turbulent extended exergy $ \left( {e_{\text{t}} } \right) $ , which is later converted to turbulent kinetic energy $ \left( {k_{\text{t}} } \right) $ , and this process occurs during the initiation and intensification of the cold vortex. In addition, barotropic energy conversion is another important process that contributes to the growth of k _t, and it strengthens gradually after the initiation of the cold vortex. Other than frictional consumption, the flux of k _t in the vertical direction also depletes some of k _t. The fluxes of e _t, baroclinic energy conversions and diabatic generations are favorable factors for the growth of e _t, whereas it decreases with time as a result of a large amount of e _t that is released. Most of the energy conversion processes, including the baroclinic and the barotropic energy transformations and the energy conversions from e _t to k _t, as well as the fluxes of e _t, are stronger in the lower troposphere than the other areas during the formation of the cold vortex. This accounts for the initiation of the cyclonic vorticity in the lower troposphere. Finally, the fact that the turbulent extended exergy releases primarily in the middle troposphere through the vertical thermal circulation is consistent with our understanding based on the vorticity budget analyses.     

在大气运动中考虑重力加速度空间变化的问题

给出了在旋转坐标系中考虑重力加速度g的空间变化的大气控制方程组,并证明了大气总质量和总能量的守恒性,以及和取常值g时得到的结果相一致的动、位能和动、内能之间的转换关系.还讨论了在球坐标系中应用方程组时可能出现的困难,给出了在高度近似下在该系中考虑g的空间变化的方案,它可以用来建立完全弹性的非静力模式.     

Localized multiscale energy and vorticity analysis: I. Fundamentals

A new methodology, multiscale energy and vorticity analysis (MS-EVA), is developed to investigate the inference of fundamental processes from oceanic or atmospheric data for complex dynamics which are nonlinear, time and space intermittent, and involve multiscale interactions. Based on a localized orthogonal complementary subspace decomposition through the multiscale window transform (MWT), MS-EVA is real problem-oriented and objective in nature. The development begins with an introduction of the concepts of scale and scale window and the decomposition of variables on scale windows. We then derive the evolution equations for multiscale kinetic and available potential energies and enstrophy. The phase oscillation reflected on the horizontal maps from Galilean transformation is removed with a 2D large-scale window synthesis. The resulting energetic terms are analyzed and interpreted. These terms, after being carefully classified, provide four types of processes: transport, transfer, conversion, and dissipation/diffusion. The key to this classification is the transfer–transport separation, which is made possible by looking for a special type of transfer, the so-called perfect transfer. The intricate energy source information involved in perfect transfers is differentiated through an interaction analysis.The transfer, transport, and conversion processes form the basis of dynamical interpretation for GFD problems. They redistribute energy in the phase space, physical space, and space of energy types. These processes are all referred to in a context local in space and time, and therefore can be easily applied to real ocean problems. When the dynamics of interest is on a global or duration scale, MS-EVA is reduced to a classical Reynolds-type energetics formalism.     

On the effect of thermal forcing on the global atmospheric angular momentum and the general circulation

The effect of latitudinal differential heating on the atmospheric general circulation is studied using a simple general circulation model driven by different heating rates. It is found that an increase in differential heating leads to a strengthening of the general circulation as characterized by an increase in global available potential energy, kinetic energy and atmospheric angular momentum. The strengthening of the circulation results in three circulation regimes characterized by different eddy activities. One is a zonally symmetric Hadley regime with no eddy activity while the other two are Rossby regimes dominated by eddies with intermediate and low wave numbers, respectively. Relative to other global indices, the global relative atmospheric angular momentum is superior in detecting the transition of circulation regimes. The regime changes in mid- and high-latitudes resemble the response of the atmosphere to large changes in Earths rotation rate.     

Theoretical studies of the parameterization of the non-neutral surface boundary layer

The parameterization of the non-neutral atmospheric surface layer has been reexamined using the basic principles of small-scale energetics and thermodynamics. On the basis of this more complete treatment, theK-parameterization has been reformulated. It is found that the linear regression laws between fluxes and driving gradient forces of the turbulent heat and humidity exchanges in the surface layer can be derived in a much more comprehensive manner than by using the commonly used K-theory. With respect to stationary conditions and in the context of similarity concepts, a system of algebraic equations has been formulated which provides reasonable estimates of the distributions of the dimension less rates of viscous energy dissipation as well as turbulent kinetic and thermal-diffusive energy transport as functions of the variablez/L. Quantitative calculations have been performed using the scaling height formulations of Takeuchi and Yokoyama, Prandtl, and von Karman as closure conditions of the equations.     

Turbulence Reynolds number and the turbulent kinetic energy balance in the atmospheric surface layer

The relation between the turbulence Reynolds numberR and a Reynolds numberz* based on the friction velocity and height from the ground is established using direct measurements of the r.m.s. longitudinal velocity and turbulent energy dissipation in the atmospheric surface layer. Measurements of the relative magnitude of components of the turbulent kinetic energy budget in the stability range 0 >z/L 0.4 indicate that local balance between production and dissipation is maintained. Approximate expressions, in terms of readily measured micrometeorological quantities, are proposed for the Taylor microscale and the Kolmogorov length scale .     

Comparative energetics analysis of CCM2 with different horizontal resolutions

 Comprehensive global energetics analysis is carried out for the NCAR CCM2 with different horizontal resolutions of R15, T42, T63, and T106 to assess the effect of various model truncations on the global energetics characteristics in climate models. Both the energy levels and energy transformations are examined over the zonal wave number domain during a northern winter and summer. In addition to the simulated atmosphere, the ECMWF global analysis during 1986 to 1990 is analyzed for comparison using the same diagnostic scheme. Previous studies have revealed that zonal kinetic energy is supplied by synoptic disturbances in terms of the zonal-wave interactions of kinetic energy. According to our result, however, such an energy flow from eddies to zonal motions is valid only for zonal wave numbers up to about 30. We find that the zonal-wave interactions of kinetic energy change sign beyond wave number 30 where the energy is transformed from zonal to eddies for both the ECMWF and CCM2-T106. The large-scale zonal motions are diffusive against the short waves beyond wave number 30, which may well be parameterized by various forms of the diffusion schemes. We suggest from this result that the atmospheric disturbances with wave numbers lower than 30 are necessary to represent accurately the two-way interactions between zonal and eddy motions, because these waves can actively influence the behavior of the zonal motions. Based on this finding, we suggest that the model resolution of R15 is inadequate for climate studies from the energetics point of view, and that resolution of T42 is the minimum requirement to represent the general circulation adequately. Some other discrepancies are discussed in detail for the coarse resolution climate models. Received: 15 July 1996/Accepted: 3 January 1997     

A global available potential energy‐kinetic energy budget in terms of the two‐dimensional wavenumber for the FGGE year

Abstract

A global vertically integrated available potential energy‐kinetic energy budget in terms of the two‐dimensional wavenumber is formulated using spherical harmonics. Results of the budget equations applied to the four mid‐season months of the FGGE year are given.     

Global diagnostic energetics of five state-of-the-art climate models

In this study, the global energy cycle of five state-of-the-art climate models is evaluated in the wave number domain for all seasons. The energy cycle estimates are based on 30?years of 6-hourly data obtained at pressure levels of all models. The models energetics are compared to those obtained from three reanalysis datasets (ERA-40, JRA-25 and NCEP-R2). The results show that the distributions of the energetics integrands and the shape of the various wave number spectra are reasonably well simulated. Many important features can be found in most models, namely both the upscale and downscale energy cascade for the wave?Cwave interactions of kinetic energy, the downscale energy cascade for the wave?Cwave interactions of available potential energy and the downscale energy transfer for the zonal?Cwave interactions of kinetic energy. However, the magnitude in the integrands distributions is generally excessive, yielding too much energy and an overactive energy cycle in the models. Accordingly, this energy excess is also reflected in the various spectra, specially but not exclusively, at the synoptic scale wave numbers for the energy conversion/transfer rates. The well known cold pole bias and the too strong tropospheric jets, along with their dislocation in some cases, still persist in the climate models. These are some of the deficiencies in the models directly implicated in the energy cycle. Apparently, simply increasing the horizontal and vertical resolutions is not enough to eliminate these deficiencies due to somewhat opposite effects achieved by refining both spatial resolutions. Therefore, more accurate physics parameterisations as well as improved numerical schemes and resolution dependence of parameterisations seem to be essential for a significant improvement in the models energetics. Moreover, efforts should be made to improve the physical processes controlling the generation of zonal available potential energy and dissipation of eddy kinetic energy, in which the synoptic scale should be fundamental, as inferred from the excessive energy conversion/transfer rates in the models spectra.     

Sensitivity analysis of the equations for a convective mixed layer

Jump or slab models are frequently used to calculate the depth of the convectively mixed layer and its potential temperature during the course of a clear day. Much attention has been paid theoretically to the parameterization of the budget for turbulent kinetic energy that is required in these models. However, for practical applications the sensitivity of the solutions of the model equations to variations in the entrainment formulation and in the initial and boundary conditions is also very important. We analyzed this sensitivity on the basis of an analytical solution for the model which uses the well-known constant heat flux ratio. The initial conditions for the mixed-layer height (h) and potential temperature ( _m ) quickly lose their influence. Only the initial temperature deficit is important. The mixed-layer temperature at noon on convective days is insensitive to the entrainment coefficient c. It is governed by the integral of the heat input and by the stable lapse rate. A change in c from 0.2 to 0.5 leads to a variation of 20% in h. This is not very much considering the accuracy in the determination of h from actual observations.     

Eta model forecasts of tropical cyclones from Australian Monsoon Experiment: Dynamical adjustment of initial conditions

Summary The University of Belgrade/National Meteorological Centre, Washington (UB/NMC) limited area Eta Model predicted the development, structure, associated precipitation and tracks of the Australian Monsoon Experiment (AMEX) (10 January through 15 February 1987) tropical cyclonesConnie, Irma, Damien andJason. The initial positions and intensities of the tropical cyclone vortices from the global European Centre for Medium-Range Weather Forecasts (ECMWF) analyses, which are used as initial data, do not quite agree with the observations. These disagreements produce additional erros in predicting the tropical cyclone tracks.To improve the initial position of the vortex, the flow is split into the small and the large scale motions, and during the first two hours of the integration, the small-scale part is forced in small steps towards its observed position. The adjustment is performaed with the reduced model dynamics (adjustment processes only) and no physics.With the adjustment during the first two hours of the integration, the model successfully adjusts to the new position of the initial vortex. After the completion of the adjustment stage, the model runs normally, i.e., without any modification. The tracks of the 48-h forecasts with the adjusted initial vortices are parallel to the tracks obtained in the control forecasts without the adjustment. However, e.g., the mean absolute error of the positions during 48-h forecast of the tropical cycloneConnie was reduced from 174 km in the control case to 129 km in the case with adjustment of the position.The latent heat, the thermal energy, the kinetic energy and the total energy of the extracted small scale vortices are calculated every three hours of the integration time. These small-scale energies obtained in the 48-h control forecast are compared to those of the rund with the initial vortex adjustment to monitor the spin up of the model.With 7 Figures     

全球大气能量的时空特征及变化趋势分析

大气能量学是大气科学重要的组成部分，了解大气能量的时空分布和变化特征，能够为大气科学研究，尤其是气候变化研究提供新的思路和手段。本文基于1948～2016年NCEP逐月再分析资料，从大气的总能量及其内能、位能、潜热和动能的分布、变化趋势和主模态变化等方面阐释了全球大气能量变化的整体特征。主要结论如下:（1）除高海拔地区外，总能量呈现从赤道向两极逐渐递减的分布，且全球大部分地区呈增加趋势，内能和位能的分布和变化与总能量较为接近；潜热能的极大值区和显著变化区均位于赤道及低纬地区；动能的极大值区分布在中纬度长波槽和西风急流出口区，其在南半球双西风急流区的变化最为显著。（2）总能量呈现出不连续的阶段性跳跃式增长特征；北半球的总能量多于南半球，而增速却慢于南半球，即两半球间的能量呈趋同趋势；海洋上空的总能量多于陆地，且海陆间差额有增大趋势；火山爆发事件可能对大气能量在年际尺度上的减少有重要作用。（3）大气各能量第一模态的空间特征与其各自变化趋势分布非常相似，并先后在1975年左右发生了年代际突变。就第二模态而言，大气的总能量、内能和位能从整体上反映出南北极与其它地区呈反向变化的特征；部分低纬度地区的潜热能与其它地区呈反向变化；动能主要呈现从热带太平洋向南北两极的经向波列分布；它们的时间系数均有一定的多年代际变化特征，可能与气候系统的内部变率有关。     

中纬度对流层上部的能量变化

能量在大气中是一个主要的因素,大气中一切的现象之所以能发生,都是因为大气有适当的能量的供给和变化。大气中能量的式样很多,不同式样的能力有不同的重要性,研究各种式样的能力的多寡和它的变化是很有意义的,因为大气环流