Spectral analysis method for a limited area using the eigenmodes of the spherical Laplacian operator

We propose a spectral analysis method using the eigenmodes of the spherical Laplacian operator on the limited area domain. Two numerical methods are considered for the horizontal discretization: One uses the half-ranged Fourier series for both longitudinal and latitudinal directions, and the other uses the Fourier finite-element method with piecewise linear basis functions for the latitudinal direction. The field variable for the two numerical algorithms is represented as linear combinations of the eigenvectors of the Laplacian operator on the limited area domain; we define the one-dimensional spectrum with the eigenvector coefficients as a function of the indices equivalent to the total wavenumbers of the Laplacian operator on the global domain. The spatial robustness of this method was verified through the self-consistency test comparing the spectra of isotropic Gaussian bells on the sphere. We used the method in the kinetic energy spectral analysis for a limited area with global atmospheric data, and compared the results for different seasons. The kinetic energy spectra represented the well-known characteristics with scale and different powers with season.     

ENERGY,ENERGY FLUX AND LAGRANGIAN OF ROSSBY WAVE

The energy flux derived from the barotropic vorticity equation differs from that obtained directly from the momentum equation.We re-study this problem raised in the early 1960s.The results show that if the momentum equation is rewritten in such a way that it contains the same conditions as that for the barotropic vorticity equation,then the same form of average energy flux can be obtained for the waves with constant amplitudes.With this new momentum equation,the potential energy of Rossby wave is derived and Lagrangian of nonlinear barotropic vorticity equation can be approximately found with this potential energy.This provides a physical basis for studying the dynamics of nonlinear Rossby wave with the approach of calculus of variation.     

DEVELOPMENT OF A TROPICAL CYCLONE IN BAROTROPIC ENVIRONMENTAL FLOWS

Development of a tropical cyclone in barotropic environmental flows is investigated with a shallow water model. It is found that the tropical cyclone develops only when it is embedded in an environmental now with a southward relative vorticity gradient. The related physical mechanism is explained by analyzing the kinetic energy conversion between the environmental flow and tropical cyclone circulation.     

Kinetic equations and stationary energy spectra of weakly nonlinear internal gravity waves

An ensemble of random-phase internal gravity waves is considered in the dynamical framework of the Euler–Boussinesq equations. For flows with zero mean potential vorticity, a kinetic equation for the mean spectral energy density of the waves is obtained under hypothesis of Gaussian statistics with zero correlation length. Stationary scaling solutions of this equation are found for almost vertically propagating waves. The resulting spectra are anisotropic in vertical and horizontal wave numbers. For flows with small but non-zero mean potential vorticity, under the same statistical hypothesis applied to the wave part of the flow, it is shown that the vortex part and the wave part decouple. The vortex part obeys a limiting slow dynamics equation exhibiting vertical collapse and layering which may contaminate the wave-part spectra. Relation of these results to the in situ atmospheric measurements and previous work on oceanic gravity waves is discussed.     

Mean dynamical characteristics of the Asian summer monsoon with a global analysis — Forecast system

Summary The dynamical characteristics of Asian summer monsoon are examined with a global spectral model. In addition to the seasonal circulation features, the large scale budgets of kinetic energy, vorticity and angular momentum are examined making use of mean analysis and forecast fields (upto day 5) for summer season comprising June, July and August, (JJA) 1994. Apart from elucidating the systematic errors over the monsoon region, the study expounded the influence of these errors on associated dynamics.The significant errors in the low level flow (850 hPa) evince (i) weakening of south easterly trades and cross equatorial flow into the Northern Hemisphere and (ii) weakening of westerly flow over the North Indian Ocean. Similarly, in the upper level flow (200 hPa) they indicate (i) weakening of Tibetan anti-cyclone and (ii) reduction of return flow into the Southern Hemisphere.The balance requirements in the mean analysis as well as forecast fields are fairly in agreement with the observed over the summer monsoon. The monsoon domain is discerned as the source region of kinetic energy and vorticity. Both are produced in the monsoon region and transported horizontally across. The model forecasts fail to retain the analyzed atmospheric variability in terms of mean as well as transient circulations which is revealed by kinetic energy budget. This is further corroborated by vorticity and angular momentum budgets. It is noticed that the forecasts in all ranges produce feeble monsoon circulation which weakens considerably with increase in the forecast period.With 16 Figures     

SYMMETRIC AND ASYMMETRIC MOTIONS IN THE BAROTROPIC FILTERED MODEL ATMOSPHERE

By expressing the stream function in terms of the symmetric and the antisymmetric components with respect to the equator,the barotropic vorticity equation can be separated into two equations,one describing the behaviours of the atmospheric symmetric motion,and the other those of the atmospheric antisymmetric motion.From them,for the entire global surface at the equivalent barotropic level,the conservations of several basic quantities,such as vorticity,angular momentum etc.in the symmetric and antisymmetric cases,have been discussed.In addition,the energy budget equations and the energy conservation expressions for the two kinds of motion are given.It can be seen from them that there are not only the conversion between the zonal mean kinetic energy and the disturbed kinetic energy in the atmosphere,but also the conversion between the kinetic energy generated by the symmetric and antisymmetric motions.In the case of including orography and horizontal diffusion into the vorticity equation,a mechanism of the generation of asymmetric behaviours of the atmosphere is proposed.The results show that the asymmetric distribution of orography and that of horizontal diffusion coefficients are likely the causes leading to the asymmetric motion.Finally,a theoretical comparison between the global model and the hemispheric one from the physical point of view is made.     

The effect of barotropic shear on baroclinic instability Part I: Normal mode problem

The effect of barotropic shear in the basic flow on baroclinic instability is investigated using a linear multilevel quasi-geostrophic β-plane channel model and a nonlinear spherical primitive equation model. Barotropic shear has a profound effect on baroclinic instability. It reduces the growth rates of normal modes by severely restricting their structure, confirming earlier results with a two-layer model. Dissipation, in the form of Ekman pumping and Newtonian cooling, does not change the main characteristics of the effect of the shear on normal mode instability.Barotropic shear in the basic state, characterized by large shear vorticity with small horizontal curvature, also effects the nonlinear development of baroclinic waves. The shear limits the energy conversion from the zonal available potential energy to eddy energy, reducing the maximum eddy kinetic energy level reached by baroclinic waves. Barotropic shear, which controls the level of eddy activity, is a major factor which should be considered when parameterizing the eddy temperature and momentum fluxes induced by baroclinic waves in a climate model.     

一次入海气旋快速发展的动力和热力学特征分析

利用扩展Zwack-Okossi方程(扩展Z-O方程)对2007年3月3-5日黄、渤海海上快速发展的温带气旋进行了诊断分析。结果表明,潜热释放和暖平流在气旋发展初期起主要作用,潜热释放是气旋接近爆发性发展的最主要强迫项,这与以往研究的主要受涡度平流强迫的中国沿海温带气旋不同,说明在中国近海快速发展的一部分温带气旋主要是由潜热释放项强迫的。绝对涡度平流在气旋发展后期才成为主要强迫项。进一步利用MM5模式模拟了潜热释放的作用以及其他一些影响气旋发展的因素,海表面热通量、水汽通量以及气旋路径上的地形对这次气旋的发展有积极贡献,但作用相对潜热释放较小。由于这种在中国近海接近爆发性发展的温带气旋容易造成风暴潮,因此利用动能方程对气旋发展中期的渤海大风进行了诊断分析。结果表明,动能水平平流是动能增加的主要作用项,但动能达到最大值后渤海海域的动能主要由有效位能转换和动量下传提供。     

Partitioning Eddy-Covariance Methane Fluxes from a Shallow Lake into Diffusive and Ebullitive Fluxes

Methane (\(\mathrm {CH}_{4}\)) is known to be emitted from lakes to the atmosphere via processes such as diffusion and ebullition (i.e., bubble emission). We developed a practical method for partitioning eddy-covariance \(\mathrm {CH}_{4}\) fluxes from a shallow lake into diffusive and ebullitive fluxes using a wavelet analysis based on local scalar similarity between the \(\mathrm {CH}_{4}\) concentration and other reference scalars, such as the air temperature or water vapour concentration, in the wavelet time-scale domain, with the hypothesis that similar and dissimilar fluctuation components are related to diffusive and ebullitive \(\mathrm {CH}_{4}\) fluxes, respectively. Our method is applied to approximately two weeks of data obtained at a shallow mid-latitude lake. The partitioned diffusive flux has a physically sound relationship with wind speed, supporting the validity of the method. The ratio of the diffusive flux to the total flux is typically 0.11 with flow from an area of steady bubble emission, but otherwise 0.36. Further validation is required using a larger dataset and data from other lakes. The proposed method can be easily applied to historical data because it requires only 10-Hz data of \(\mathrm {CH}_{4}\) concentration and other reference scalars, along with an empirical parameter.     

大气中动能的双向转化和单向转化

根据绝热无耗散大气总动能和总位涡拟能守恒的特性，本文讨论了二维平面上动能谱分布随时间变化的规律。结果表明，大气谱动能变化必然涉及至少三种不同的空间尺度；或遵从双向转化原则，或遵从单向转化原则。在讨论能量的双向转化和单向转化条件的基础上，本文还分别研究斜压大气和正压大气中能量转化的一般状况。     

ENERGY DISPERSION EFFECTS ON EXPLOSIVE DEVELOPMENT OF MARINE CYCLONES

A dominant role played by energy dispersion in the explosive development of extratropicalmarine cyclones over the Northwest Pacific has been revealed based on both the eddy energyequations and the energy flux vectors of nonlinear wave packet.At the initial and explosive time.the eddy energy from neighboring upstream systems beyond the radius of Rossby deformation isdispersed into the eddy energy center associated with the cyclone via ageostrophic geopotentialfluxes,and results in the rapid increase of eddy kinetic energy and the occurrence of explosivecyclogenesis.When the cyclone begins to decay,its corresponding eddy energy is exporteddownstream and hence triggers the growth of new perturbation downstream.Through generalizingthe energy flux vectors of quasigeostrophic wave packets to the nonlinear forms and making use ofthe relationship between the energy flux vectors and the total eddy energy,the approximationexpressions of the total group velocity and relative group velocity are derived,and then they areused to compute an explosive case.The normalized ageostrophic geopotential fluxes by dividing thevolume integral of ageostrophic geopotential fluxes by the integral of the total eddy energydetermine the relative group velocity at which the eddy energy is spreading out,and they can beused to evaluate the position of next new disturbance.The nonlinear advective fluxes influenceprimarily the phase speed and translation of the cyclones.The results in this paper facilitate toexpanding the mechanism research on explosive cyclones and have great significance for predictingthe explosive intensification and downstream disturbance growth.     

雷暴大风过程中对流层中低层动量通量和动能通量输送特征研究

2014年5月31日北京发生一次雷暴大风过程。以雷达资料同化结果为初始场,对此次过程进行高分辨率数值模拟。采用非静力平衡和非地转平衡的经向动量方程和质量权重动能方程,利用模拟资料,对雷暴大风过程中经向动量和质量权重动能进行收支分析,以此来研究雷暴过程中对流层中低层动量通量和动能通量输送特征,讨论地面大风的可能成因。分析结果表明,在对流层中低层,经向动量通量散度是影响经向动量局地变化的主要强迫项。雷暴系统后部的入流把中低层的经向动量倾斜向下输送,系统前部对流云区中低层的下沉气流也向下输送经向动量。这两支下传动量通量先后与近地面经向动量的水平通量汇合,向系统前沿输送经向动量。在北京西北部地形阻挡作用下,经向动量通量在系统前端近地面辐合,促进那里的经向动量局地增长,有利于增强那里的南风。质量动能收支的特征与经向动量收支类似,在近地面层质量动能的局地变化主要是由质量动能通量散度引起的。系统后部入流把中层质量动能向下传输到近地面层,然后与近地面质量动能的水平通量汇合,向系统前沿输送质量动能。相对来说,近地面层经向动量和质量动能的水平通量比下传通量更重要,这主要与低层较强的东南急流有关。     

Turbulent flux observations and modelling over a shallow lake and a wet grassland in the Nam Co basin,Tibetan Plateau

The Tibetan Plateau plays an important role in the global water cycle and is strongly influenced by climate change. While energy and matter fluxes have been more intensely studied over land surfaces, a large proportion of lakes have either been neglected or parameterised with simple bulk approaches. Therefore, turbulent fluxes were measured over wet grassland and a shallow lake with a single eddy-covariance complex at the shoreline in the Nam Co basin in summer 2009. Footprint analysis was used to split observations according to the underlying surface, and two sophisticated surface models were utilised to derive gap-free time series. Results were then compared with observations and simulations from a nearby eddy-covariance station over dry grassland, yielding pronounced differences. Observations and footprint integrated simulations compared well, even for situations with flux contributions including grassland and lake. The accessibility problem for EC measurements on lakes can be overcome by combining standard meteorological measurements at the shoreline with model simulations, only requiring representative estimates of lake surface temperature.     

Simulation of the Bare Soil Surface Energy Balance at the Tongyu Reference Site in Semiarid Area of North China

The performance of a 1-D soil model in a semiarid area of North China was investigated using observational data from a cropland station at the Tongyu reference site of the Coordinated Enhanced Observing Period (CEOP) during the non-growing period, when the ground surface was covered with bare soil. Comparisons between simulated and observed soil surface energy balance components as well as soil temperatures and water contents were conducted to validate the soil model. Results show that the soil model could produce good simulations of soil surface temperature, net radiation flux and sensible heat flux against observed values with the RMSE of 1.54oC, 7.71 W m-2 and 27.79 W m-2, respectively. The simulated volumetric soil water content is close to the observed values at various depths with the maximal difference between them being 0.03. Simulated latent heat and ground heat fluxes have relatively lower errors in relative to net radiation and sensible heat flux. In conclusion, the soil model has good capacity to simulate the bare soil surface energy balance at the Tongyu cropland station and needs to be further tested in longer period and at more sites in semiarid areas.     

边界层过程对暴雨影响的敏感性试验

用MM5v3模式对江淮暴雨进行了计入与不计入边界层和计入与不计入地面通量的敏感性试验,结果指出二者的影响是一致的,地面通量是边界层影响的主要部分。一般说,不计边界层影响,则湿位涡、水汽通量散度、涡度、低空急流等分布发生变化,其中心位置改变,强度也有改变(大部分减弱),造成总的降水减弱和降水位置的改变。边界层正是通过湿位涡、水汽通量散度、涡度、低空急流等因素的综合影响来影响暴雨,并影响暴雨的落区和强度。     

应用E-ε湍流动能闭合湖泊热力学过程模型对东太湖湖-气交换的模拟

采用考虑沉水植物影响的E-ε湍流动能闭合湖泊热力学过程模型,模拟2013年8月东太湖湖-气交换过程,并利用太湖的站点观测数据对模型进行了验证。太湖水温的模拟值与观测值吻合较好,模型计算的各层水温与观测值相比,均方根误差均未超过1℃。同时模型也较好地模拟出太湖表面感热通量和潜热通量,潜热通量的模拟值与观测值的标准差为54.7 W/m2。由于湖水较浅,太湖的水温层结会明显受到天气状况的影响。晴朗小风条件下的湖水呈现显著的热分层现象,当风速为0.8 m/s,高层和底层的温差达到7.9℃。大风天气条件驱动较强的水体湍流混合,水温的热分层消失,风速为12 m/s,湖泊上层与底层的水温差仅0.12℃。此外,模拟结果较好地呈现出了东太湖沉水植物的存在通过增大湖体消光系数,减小到达湖体内部的热量,并增加对湖水的阻力,影响湖体中湍流动能的分布,并进而影响湖水温度的分布。综上所述,该模型能够较好地模拟出浅水大湖湖-气交换的过程。      

An Updated Coupled Model for Land-Atmosphere Interaction. Part I： Simulations of Physical Processes

A new two-way land-atmosphere interaction model （R42_AVIM） is fulfilled by coupling the spectral atmospheric model （SAMIL_R42L9） developed at the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences （LASG/IAP/CAS） with the land surface model, Atmosphere-Vegetation-Interaction-Model （AVIM）. In this coupled model, physical and biological components of AVIM are both included. Climate base state and land surface physical fluxes simulated by R42_AVIM are analyzed and compared with the results of R42_SSIB [which is coupled by SAMIL_R42L9 and Simplified Simple Biosphere （SSIB） models]. The results show the performance of the new model is closer to the observations. It can basically guarantee that the land surface energy budget is balanced, and can simulate June-July-August （JJA） and December-January- February （DJF） land surface air temperature, sensible heat flux, latent heat flux, precipitation, sea level pressure and other variables reasonably well. Compared with R42_SSIB, there are obvious improvements in the JJA simulations of surface air temperature and surface fluxes. Thus, this land-atmosphere coupled model will offer a good experiment platform for land-atmosphere interaction research.     

一次梅雨暴雨的动能谱特征分析及收支诊断

利用WRF模式和NCEP再分析资料,对一次梅雨暴雨个例进行了数值模拟,基于高分辨率的模拟结果计算了水平动能谱、涡旋动能谱及辐散动能谱,诊断了动能收支谱方程。结果表明:在暴雨发展阶段,各个高度上都有中尺度动能增长,其显著增加始于中尺度低端,这导致了在中尺度波段出现谱转折特征,但在不同高度转折尺度不同;对流层高层涡旋动能大于辐散动能,平流层低层反之;不同降水阶段、不同高度和不同尺度上动能的来源不同,对流层高层,中尺度动能倾向由非线性项和气压项贡献;平流层低层,气压项的作用更为明显。     

北京北郊冬季大风过程湍流通量演变特征的分析研究

利用中国科学院大气物理研究所325 m气象观测塔1993年12月～1994年1月大气边界层实验资料, 计算分析了大风过境过程中47 m和120 m高度湍流通量演变特征及其影响因子, 以及与风速、 稳定度等参数的关系。结果表明: 大风过程对近地面层的物质能量输送有着重要影响, 大风之前出现短时间动量上传和热量下传; 大风过程中的湍流通量数值明显高于过境后, 水平方向湍流通量数值和能量增加幅度大于垂直方向; 当风速大于临界值5 m/s时, 湍流通量与风速、 湍流动能的相关迅速增大; 湍流谱特征表现为湍流能量的低频部分增加、 湍流谱曲线变宽; 大风能强烈影响近地面层的能量收支。     

Influence from Surrounding Land on the Turbulence Measurements Above a Lake

Turbulence measurements taken at a Swedish lake are analyzed. Although the measurements took place over a relatively large lake with several km of undisturbed fetch, the turbulence structure was found to be highly influenced by the surrounding land during daytime. Variance spectra of both horizontal velocity and scalars during both unstable and stable stratification displayed a low frequency peak. The energy at lower frequencies showed a daily variation, increasing in the morning and decreasing in the afternoon. This behaviour is explained by spectral lag, where the low frequency energy due to large eddies that originate from the convective boundary layer above the surrounding land. When the air is advected over the lake the small eddies rapidly equilibrate with the new surface forcing. However, the large eddies remain for an appreciable distance and influence the turbulence in the developing lake boundary layer. The variances of the horizontal velocity and scalars are increased by these large eddies, while the turbulent fluxes are mainly unaffected. The drag coefficient, Stanton number and Dalton number used to parametrize the momentum flux, heat flux and latent heat flux respectively all compare well with current parametrizations developed for open sea conditions. The diurnal cycle of the partial pressure of methane, $p\mathrm{CH}_{4}$ , observed at this site is closely related to the diurnal cycle of the lake-air methane flux. An idealized two-dimensional model simulation of the boundary layer at a lake site indicates that the strong response of $p\mathrm{CH}_{4}$ to the surface methane flux is due to the shallow internal boundary layer that develops above the lake, allowing methane to accumulate in a relatively small volume.