CISK理论中的若干问题

在本文中,对Charney与Eliassen等人所提出的CISK理论所包含的低层风速不随高度变化的假定作了分析与改进,并讨论了它对结果的影响。     

Large-Scale Motions in the Marine Atmospheric Surface Layer

Multi-level turbulent wind data from the Risø Air-Sea Experiments (RASEX) were used to examine the structure of large-scale motions in the marine atmospheric surface layer. The quadrant technique was used to identify flux events (ejections/sweeps). Ejections, which appear to occur in groups, are seen to occur first at the upper level, moving successively to lower levels with small time delays. A strong correlation between events at different heights suggests that they may all be part of a single large structure. Cross-correlation between velocity signals was used to estimate orientation of the structure using Taylor's hypothesis. The inclination of this structure is shallow ( 15°) near the surface and increases with height. Spatial representations of the fluctuating wind vectors show a structure that is strikingly similar to conceptual models of transverse vortices and shear layers seen in laboratory flows and direct numerical simulation (DNS) of low Reynolds number flows. Spatial visualization of velocity fluctuations during other time periods and conditions clearly shows the existence of shear layers, transverse vortices, plumes, and downdrafts of various sizes and strengths. A quantitative analysis shows an increase in the frequency of shear related events with increasing wind speed.     

湖南省油菜春季涝渍过程灾变判别指标

以湖南省油菜为研究对象,综合洪涝与连阴雨的致灾因子,利用历史灾情数据反演涝渍过程样本,构建临灾、受灾状态样本,采用均值t检验方法,定量分析不同降水累积衰减系数w情况下,降水因子差异的显著性,构建有效累积降水量P_E。基于Fisher判别准则,计算涝渍灾害临灾与受灾的临界线,构建逐日动态的油菜春季涝渍过程灾变判别指标y,并进行涝渍过程样本反演及独立样本验证。分析春季涝渍受灾频率与相对气象产量的关系,构建灾害影响评价模型。结果表明:当w=0.899时,受灾与临灾样本P_E的差异最显著。灾变判别指标y可逐日动态监测涝渍灾变过程及灾害后续影响,为动态监测区域油菜涝渍过程提供了数据支撑;指标指示的灾变范围与实际受灾情况基本一致,为区域防灾减灾提供了科学参考。开花期及结荚期湖南油菜涝渍受灾频率较高,成熟期较低。湖南省油菜春季涝渍灾害影响指数呈东南高西北低,长沙东部、株洲中部及北部、湘潭、永州及郴州中部油菜产量受涝渍灾害影响最重。     

A Wind-Tunnel Study of Airflow in Waving Wheat: An EOF Analysis of the Structure of the Large-Eddy Motion

We have conducted an Empirical Orthogonal Functionanalysis (EOF) of a three dimensional, 2-pointvelocity covariance field, measured in a wind tunnel. The rate of convergence of the EOF sequence was usedas an objective test for the presence of distinctlarge turbulent structures. We found that in theroughness sublayer (2h > z > 0) the sequence convergedmuch more rapidly than in the lower surface layer(6h > z > 0), 75% of the total velocity variancebeing captured by the first three of 42 eigenmodes;h is the canopy height. The analysis was extended to three dimensions, whereover 50% of the variance and most of the spatialstructure of the covariance fields were captured by aneven smaller fraction of the total number ofeigenmodes. With some relatively weak additionalassumptions we were able to construct the velocityfield of a characteristic eddy or large coherentstructure. This consisted of a pair ofcounter-rotating streamwise vortices centred above thecanopy. The sense of rotation of the vortex pair wasopposite to that found in the wall region of boundarylayers but matched that found in plane mixing layers. A strong gust or sweep motion generated between thevortices was responsible for most of the shear stresscarried by the large eddies. The region of significanttransport of streamwise momentum by the characteristiceddy is much smaller than the region of coherence of theeddy's velocity field.     

Turbulent Intensities and Velocity Spectra for Bare and Forested Gentle Hills: Flume Experiments

To investigate how velocity variances and spectra are modified by the simultaneous action of topography and canopy, two flume experiments were carried out on a train of gentle cosine hills differing in surface cover. The first experiment was conducted above a bare surface while the second experiment was conducted within and above a densely arrayed rod canopy. The velocity variances and spectra from these two experiments were compared in the middle, inner, and near-surface layers. In the middle layer, and for the canopy surface, longitudinal and vertical velocity variances () were in phase with the hill-induced spatial mean velocity perturbation (Δu) around the so-called background state (taken here as the longitudinal mean at a given height) as predicted by rapid distortion theory (RDT). However, for the bare surface case, and remained out of phase with Δu by about L/2, where L is the hill half-length. In the canopy layer, wake production was a significant source of turbulent energy for , and its action was to re-align velocity variances with Δu in those layers, a mechanism completely absent for the bare surface case. Such a lower ‘boundary condition’ resulted in longitudinal variations of to be nearly in phase with Δu above the canopy surface. In the inner and middle layers, the spectral distortions by the hill remained significant for the background state of the bare surface case but not for the canopy surface case. In particular, in the inner and middle layers of the bare surface case, the effective exponents derived from the locally measured power spectra diverged from their expected  − 5/3 value for inertial subrange scales. These departures spatially correlated with the hill surface. However, for the canopy surface case, the spectral exponents were near  − 5/3 above the canopy though the minor differences from  − 5/3 were also correlated with the hill surface. Inside the canopy, wake production and energy short-circuiting resulted in significant departures from  − 5/3. These departures from  − 5/3 also appeared correlated with the hill surface through the wake production contribution and its alignment with Δu. Moreover, scales commensurate with Von Karman street vorticies well described wake production scales inside the canopy, confirming the important role of the mean flow in producing wakes. The spectra inside the canopy on the lee side of the hill, where a negative mean flow delineated a recirculation zone, suggested that the wake production scales there were ‘broader’ when compared to their counterpart outside the recirculation zone. Inside the recirculation zone, there was significantly more energy at higher frequencies when compared to regions outside the recirculation zone.     

Stratocumulus-capped mixed layers derived from a three-dimensional model

Results of a three-dimensional numerical model are analysed in a study of turbulence and entrainment within mixed layers containing stratocumulus with or without parameterized cloud-top radiative cooling. The model eliminates most of the assumptions invoked in theories of cloud-capped mixed layers, but suffers disadvantages which include poor resolution and large truncation errors in and above the capping inversion.For relatively thick mixed layers with relatively thick capping inversions, the cloud-top radiative cooling is found to be lodged mostly within the capping inversion when the cooling is confined locally to the upper 50 m or less of the cloud. It does not then contribute substantially towards increased buoyancy flux and turbulence within the well mixed layer just below.The optimal means of correlating the entrainment rate, or mixed-layer growth rate, for mixed layers of variable amounts of stratocumulus is found to be through functional dependence upon an overall jump Richardson number, utilizing as scaling velocity the standard deviation of vertical velocity existing at the top of the mixed layer (near the center of the capping inversion). This velocity is found to be a fraction of the generalized convective velocity for the mixed layer as a whole which is greater for cloud-capped mixed layers than for clear mixed layers.     

Radar and sodar probing of waves and turbulence in statically stable clear-air layers

This paper reviews the remote sensing of waves and turbulence in statically stable atmospheric layers, utilizing sodar and microwave radar echoes from the small-scale inhomogeneities in gaseous refractive index caused by localized fluctuations in temperature, humidity, and velocity. Scattering theory and sounding methodology are reviewed briefly, and the relative performance of typical radar and sodar systems compared. The main section of the paper takes the form of a summary and discussion of experimental progress since 1969, showing how the echo patterns obtained may be applied to the interpretation of multiple layering, gravity waves, internal fronts and the details of dynamic instability and the genesis of turbulence in stably stratified shear layers. In addition, methods for the measurement of the intensity of the small-scale ( /2) variability of wind, temperature and water vapor from the observed radar or sodar echo intensities, and the use of Doppler techniques for the measurement of mean velocity and turbulence are discussed.SODAR from SOund Detection And Ranging in analogy to RAdio Detection And Ranging.     

2004年7月两次强对流天气过程的对比分析

对两次强对流天气过程的影响系统、物理量场和雷达回波特征等进行了对比分析.分析表明：两次强对流过程有许多相似之处,高、低空西南急流适宜配置而引起的高低层动力场的耦合是强对流发生的动力条件,强对流发生在低层露点锋区及低层水汽通量和散度辐合、高层则相应辐散的大值区域.高低层相对散度正值越大,其所产生的雨量也就越大.多普勒雷达反射率因子图上弓状回波的后侧出现下沉入流急流,相应速度图上的速度值越大,回波移动速度越快,预示着所造成的地面风力也就越大.     

天气尺度系统稳定性结构、相路特征的分岔及热源时间振荡影响问题的探讨

本文研究了有关天气系统突变的非线性因子的影响效应,指出环境大尺度基本流场的平流输送及不同尺度系统的能量转换可能对天气尺度系统异变起了重要作用;由大气扰动偏离与趋于平衡态的过程,看出扰动发展与衰减的速率与平流输送、不同尺度能量转换的非线性作用有关,还研究了扰动稳定性结构与相路特征的分岔问题和热源对扰动解的影响效应,揭示出热源时间振荡与大气扰动存在着非线性的相关,可能导致扰动的突变。     

2011年冬季贵州低温雨雪冰冻天气的成因分析

利用常规观测资料和NCEP 1°×1°及2.5°×2.5°格点再分析资料,对贵州2011年初的低温雨雪冰冻过程的天气成因进行了初步分析。结果表明：东亚地区500 hPa高度距平场＂北高南低＂分布以及中高纬阻塞高压的稳定维持,有利于引导冷空气频繁南下影响贵州,是持续低温的主要原因。滇黔准静止锋也是重要的影响系统,它的长期存在有利于阴雨天气维持;而大范围的冻雨（雪）天气则是伴随着东移的南支槽以及向上伸展的水汽辐合。由于西太平洋副高较弱和中亚地区低值系统不活跃,没有稳定持续的水汽向贵州输送,所以过程表现出短暂的间歇性,尤其是距离滇黔准静止锋区相对较远的贵州东部地区,受灾程度远不如2008年严重。因为有偏北路径的冷空气入侵,贵州中西部地区冰冻强度超过2008年;垂直结构的差异是影响范围最大强度最强的两次过程（第1和第4次）降水相态不同的重要原因;冰冻的形成与增长是多种气象因子综合影响的结果,在一定的降水条件下,日平均气温是影响冰冻强度的重要因子。     

Theory And Measurements For Turbulence Spectra And Variances In The Atmospheric Neutral Surface Layer

Predictions from a new theory for high Reynolds number turbulent boundary layers during near-neutral conditions are shown to agree well with measurements of atmospheric surface-layer variances and spectra. The theory suggests surface-layer turbulence is determined by detached eddies that largely originate in the shearing motion immediately above the surface layer; as they descend into this layer, they are strongly distorted by the local shear and impinge onto the surface. Because the origin of these eddies is non-local, they are similar to those described in previous studies as `inactive' turbulence. However, they are, in fact, dynamically highly active, supplying the major mechanism for the momentum transport, including upward bursting on the time scale of the larger eddies. The vertical velocity results show that the variance and the low frequency parts of spectra increase with height in the surface layer, while in the self similar (k₁ ^-1) range the streamwise low frequency components are approximately constant with height. These large-scale longitudinal eddies extend to a length _s, which is equal to the boundary-layer height near the surface andincreases linearly to a maximum of about three times the boundary-layer height at roughly 15 m and decreases in the upper parts of the surface layer. This lower part of the surface layer, the eddy surface layer, is the region in which the eddies impinging from layers above are strongly distorted. This new result for the atmospheric boundary layer has practical application for calculating fluctuating wind loads on structures and lateral dispersion of pollution from local sources.     

Influence of foliar density and thermal stability on profiles of Reynolds stress and turbulence intensity in a deciduous forest

Observations were made of turbulence in an extensive deciduous forest on level terrain using a vertical array of seven three-dimensional sonic anemometer/thermometers within and above the canopy. Data were collected through the period of leaf fall and over a range of thermal stabilities. A bulk canopy drag coefficient was nearly independent of the density of the forest but decreased greatly with the onset of nocturnal stability. The depth of penetration of momentum into the forest increased with leaf fall but, again, was greatly curtailed by stable conditions. Turbulent velocities decreased with increasing depth in the forest but relative turbulence intensities increased to mid-canopy levels. Leaf density influenced turbulence levels but not as strongly as did thermal stability. Thermal effects were adequately described by the single parameter h/L, where h is the canopy height and L is the Monin-Obukhov length. The longitudinal and vertical velocity correlation coefficient was larger in magnitude than expected in the upper layers of the forest but decreased to a small value in the lowest layers where the Reynolds stress was small. The ratio _w /u ^*, where u ^* is the local friction velocity, reflected changes in the uw correlation, becoming smaller than usual in the upper canopy layers. It is believed that these effects result from the intermittent, spatially coherent structures that are responsible for a large fraction of the momentum flux to the forest.     

A new parameterization of eddy diffusivities for nocturnal boundary-layer modeling

Exchange coefficients and mixing lengths under stable stratification have been studied through measurements of mean wind velocity and temperature in the nocturnal boundary layer. For values of the gradient Richardson number lower than 0.15, our measurements fit well the relation of Delage (1974). Beyond Ri = 0.15, the decrease of mixing length is much slower. So a new parameterization of turbulent exchanges is suggested. When introduced in a model of the nocturnal boundary layer, it results in a thickening of the turbulent and inversion layers.     

A numerical study of the influence of an air temperature-inversion layer and a seawater density-jump layer on the structure of interacting boundary layers

The effects of an air-temperature inversion in the atmosphere and a seawater density jump in the ocean on the structure of the atmospheric and oceanic boundary layers are studied by use of a coupled model. The numerical model consists of a closed system of equations for velocities, turbulent kinetic energy, turbulent exchange coefficient, local turbulent length scale, and stratification expressions for both air and sea boundary layers. The effects of the temperature inversion and the density jump are incorporated into the equations of turbulent kinetic energy of the atmosphere and ocean by a parameterization. A series of numerical experiments was conducted to determine the effects of various strengths of the inversion layer and surface heat fluxes in the atmosphere and of the density-jump layer in the ocean on the structure of the interacting boundary layers.The numerical results show that the temperature inversion in the atmosphere and density jump in the ocean have strong influences on turbulent structure [especially on the turbulent exchange coefficient (TEC) and turbulent kinetic energy (TKE)] and on air-sea interaction characteristics. Maxima of TKE and TEC strongly decrease with increasing strength of the inversion layer, and they disappear for strong inversions in the atmosphere. Certain strengths (density differences between the upper and the lower layers) of the density-jump layer in the ocean (₂ 0.1 g/cm³) produce double maxima in TEC-profiles and TKE-profiles in the ocean. The magnitudes of air-sea interaction characteristics such as geostrophic drag coefficient, and surface drift current increase with increasing strength of the density-jump layer in the ocean. The density-jump layer plays the role of a barrier that limits vertical mixing in the ocean. The numerical results agree well with available observed data and accepted quantitive understanding of the influences of a temperature inversion layer and a density-jump layer on the interacting atmospheric and oceanic boundary layers.     

The wind in the baroclinic boundary layer with three sublayers incorporating the weak non-linear effect

In considering the weak non-linear effect, and using the small parameter expansion method, the analyt-ical expressions of the wind distribution within PBL (planetary boundary layer) and the vertical velocity at the top of the PBL are obtained when the PBL is divided into three layers and different eddy transfer coefficients K are adopted for the three layers. The conditions of barotropy and neutrality for the PBL are extended to that of baroclinity and non-neutral stratification. An example of a steady circular vortex is used to display the characteristics of the horizontal wind within the PBL and the vertical velocity at the top of the PBL. Some new results have been obtained, indicating that the magnitude of the speed in the lower height calculated by the present model is larger than that by the model in which k is a constant within the whole boundary layer, for example, in the classical Ekman boundary layer model and the model by Wu (1984). The angle between the wind at the top of the PBL and the wind near the surface calculated by the present model is less than that calculated by the single K model. These results are in agreement with the observations.     

A Note On The Contribution Of Dispersive Fluxes To Momentum Transfer Within Canopies

Dispersive flux terms are formed when the time-averaged meanmomentum equation is spatially averaged within the canopy volume.These fluxes represent a contribution to momentum transfer arisingfrom spatial correlations of the time-averaged velocity componentswithin a horizontal plane embedded in the canopy sublayer (CSL).Their relative importance to CSL momentum transfer is commonlyneglected in model calculations and in nearly all fieldmeasurement interpretations. Recent wind-tunnel studies suggestthat these fluxes may be important in the lower layers of thecanopy; however, no one study considered their importance acrossall regions of the canopy and for a wide range of canopy roughnessdensities. Using detailed laser Doppler anemometry measurementsconducted in a model canopy composed of cylinders within a largeflume, we demonstrate that the dispersive fluxes are onlysignificant (i.e., >10%) for sparse canopies. These fluxes arein the same direction as the turbulent flux in the lower layers ofthe canopy but in the opposite direction near the canopy top. Fordense canopies, we show that the dispersive fluxes are <5% atall heights. These results appear to be insensitive to theReynolds number (at high Reynolds numbers).     

春季一次暴雪过程的多普勒雷达动力学诊断

利用太原C波段多普勒天气雷达基数据资料和自动雨量站资料,对2006年4月11日发生在山西省的一次区域性暴雪过程进行了分析。应用改进的EVAD技术,定量计算垂直高度层的平均散度和平均垂直速度,并分析平均散度和平均垂直速度随时间和高度的变化以及与降雪的对应关系。结果表明：2．5km以下始终维持一个较强的上升气流,是强降雪维持的基本动力条件;整层出现辐合上升运动,且强上升速度中心的高度随时间的演变逐渐下降,同时中高层出现强辐散下沉气流与之配合时,未来2小时将出现强降雪;降雪的强度与雷达探测范围内各高度层的辐舍、辐散有着密切的关系,且辐合、辐散的增强与减弱时间早于降雪强度的增大和减小时间,这对预报降雪的生消、雪强的增大与减小提供了一定的理论依据。     

Assessing mitigation-adaptation scenarios for reducing catastrophic climate risk

Countries can use both mitigation and adaptation strategies to protect their citizens from catastrophic risk posed by climate change (e.g., shift in the jet stream). A nation can mitigate by reducing CO₂ emissions, which reduces the probability of a catastrophic event; it can adapt by altering the infrastructure so that damages can be reduced in the event a catastrophe is realized. Herein we add to the current literature by extending the endogenous risk framework into a dynamic framework permitting analysis of both mitigation and adaptation while allowing for the dynamic process of global climate change. Our results suggest adaptation to catastrophe is a small fraction of the national climate protection budget relative to mitigation when nations cooperate fully, when damages are both continuous and catastrophic, and when nations have a short planning horizon. Adaptation becomes more important relative to mitigation when nations are unlikely to cooperate, when damages are mainly catastrophic, or when the nation’s planning horizon increases.     

二维多路多普勒雷达风场自动退模糊算法应用研究

本文选取Zhang Jian的二维多路自动退模糊算法做雷达径向速度退模糊应用研究,并做适当改进.它首先通过在一个仰角扫描数据中搜寻速度线附近的弱风区来确定初始参考点集,并对弱风区内三条邻近的径线进行连续性检测和模糊订正,然后从顺时针和逆时针两个方向进行两轮二维退模糊处理.改进后,在第一轮和第二轮二维退模糊过程中,先不在切向上对可疑模糊点进行订正,而是根据空间连续性原则,在切向上确定更多的非模糊点作为参考速度,然后在从径向上进行退模糊处理,这样可以有效避免了在真实的切变(辐合)区造成不恰当的模糊订正.对台风、龙卷、飑线等系统中的速度模糊订正试验结果表明,该算法退模糊的准确率高,在雷达观测资料无效点较多的情况下,该技术退模糊效果比WSR-98D更优,更加有利于实际业务运行和科研分析.     

Temperature distribution and current system in a convectively mixed lake

During spring and autumn, many lakes in temperate latitudes experience intensive convective mixing in the vertical, which leads to almost isothermal conditions with depth. Thus the regime of turbulence appears to be similar with that characteristic of convective boundary layers in the atmosphere. In the present paper a simple analytical approach, based on boundary-layer theory, is applied to convective conditions in lakes. The aims of the paper are firstly to analyze in detail the temperature distribution during these periods, and secondly to investigate the current system, created by the horizontal temperature gradient and wind action. For these purposes, simple analytical solutions for the current velocities are derived under the assumption of depth-constant temperatures. The density-induced current velocities are shown to be small, in the order of a few mm/sec. The analytical model of wind-driven currents is compared with field data. The solution is in good qualitative agreement with observed current velocities under the condition that the wind field is steady for a relatively long time and that residual effects from former wind events are negligible.The effect of the current system on an approximately depth-constant temperature distribution is then checked by using the obtained current velocity fields in the heat transfer equation and deriving an analytical solution for the corrected temperature field. These temperature corrections are shown to be small, which indicates that it is reasonable to describe the temperature distribution with vertical isotherms.Notation T temperature - t time - x, y, z cartesian coordinates - molecular viscosity - _h , _v horizontal and vertical turbulent viscosity - K _h ,K _v horizontal and vertical turbulent conductivity - Q heat flux through the water surface - D depth - u, v, w average current velocity components inx, y andz directions - f Coriolis parameter - p pressure - density - g gravity acceleration - a constant in the freshwater state equation - h _s deviation from the average water surface elevation - L _*,H _* length and depth scale - U _*,W _* horizontal and vertical velocity scale - T temperature difference scale - bottom slope - u _* friction velocity at the water surface - von Karman constant - L Monin-Obukhov length scale - buoyancy parameter - l turbulence length scale - C ₁,C ₂,C ₃ dimensionless constants in the expressions for the vertical turbulent viscosity - , dimensionless vertical coordinate and dimensionless local depth - angle between surface stress direction andx-axis - T _bx ,T _by bottom stress components - c bottom drag coefficient