西双版纳地区雾的数值模拟研究

建立了一个适用于复杂地形上的三维非定常雾模式,用于研究西双版纳地区雾的生消规律。模式详细考虑了湍流、长短波辐射、凝结、蒸发、重力沉降等影响因子,特别注意了植被和气溶胶辐射效应对雾的生消的影响。利用本模式对山谷风及雾的生消过程进行了实况模拟,结果表明,模式在一定程度上反映了实际情况。     

Do stratocumulus clouds detrain? FIRE I data revisited

Analyses of aircraft observations of the stratocumulus-topped boundary layer during the First ISCCP (International Satellite Cloud Climatology Project) Regional Experiment (FIRE I) show the frequent presence of clear, but relatively moist, air patches near the stratocumulus cloud-top interface. A conditional sampling of measurements in these clear air patches shows that their thermodynamic properties do more resemble boundary-layer air characteristics than those of free troposphere air. From an aircraft leg through cloud tops it is demonstrated that turbulent mixing across the cloud-top interface can lead to the local dissipation of the cloud top. Analogous to the terminology used for shallow cumulus parameterizations this process can be considered as detrainment, with which we mean that after a mixing event across the cloud-top boundaries, mixed unsaturated parcels become part of the clear environment of the cloud.     

Turbulence structure in temperature inversions and in convection fields as observed by doppler sodar

The structure of turbulence in an inversion layer and in an homogeneous convective field of the planetary boundary layer is described. In the first part of the paper, we validate the sodar estimates of turbulent dissipation, by using measurements with an hot-wire anemometric system in situ. Limitations of an ε measurement technique using structure function calculations are given, taking account of atmospheric properties and acoustic Doppler instrumental effects. By comparison between isopleths of backscattering intensity and of turbulent dissipation rates, we observe that in the early morning, turbulence is advected by mechanical turbulence generated by wind shear. The same mechanism seems to be operating in the case of an inversion layer capping thermal instability, when the convective activity is not too greatly developed. A turbulent kinetic energy budget is examined using aircraft, sodar, and tower measurements. This indicates a constant turbulent dissipation profile through a deep convective layer.     

Local Imbalance of Turbulent Kinetic Energy in the Surface Layer

We utilize experimental data collected in 2002 over an open field in Hanford, Washington, USA, to investigate the turbulent kinetic energy (TKE) budget in the atmospheric surface layer. The von Kármán constant was determined from the near-neutral wind profiles to be 0.36 ± 0.02 rather than the classical value of 0.4. The TKE budget was normalized and all terms were parameterized as functions of a stability parameter z/L, where z is the distance from the ground and L is the Obukhov length. The shear production followed the Businger–Dyer relation for −2 < z/L < 1. Contrary to the traditional Monin–Obukhov similarity theory (MOST), the shear, buoyancy and dissipation terms were found to be imbalanced due to a non-zero vertical transport over all stabilities. Motivated by this local imbalance, modified parameterizations of the dissipation and the turbulent transport were attempted and generated good agreement with the experimental data. Assuming stationarity and horizontal homogeneity, the pressure transport was estimated from the residual of the TKE budget.     

A parameterization of the stable atmospheric boundary layer

Two formulations of the stable atmospheric boundary layer are proposed for use in weather forecasting or climate models. They feature the log-linear profile near the surface, but are free from the associated critical Richardson number. The diffusion coefficients in the Ekman layer are a natural extension of the surface layer. They are locally determined using wind shear in one case and turbulent kinetic energy in the other. The parameterizations are tested in a one-dimensional model simulating the evolution of the nocturnal boundary layer with and without radiative cooling. Both formulations give very similar results, except near the top of the boundary layer where the transition to the free atmosphere is smoother with the wind shear formulation. A distinctive feature of these schemes is that they retain their simulating skill when resolution is reduced. This is verified for a wide range of situations. In practice, this means that there is no need for a large-scale model to have a level below 50 m or so.     

Implementation of a Stable PBL Turbulence Parameterization for the Mesoscale Model MM5: Nocturnal Flow in Complex Terrain

The difficulties associated with the parameterization of turbulence in the stable nocturnal planetary boundary layer (PBL) have been a great challenge for the nighttime predictions from mesoscale meteorological models such as MM5. As such, there is a general consensus on the need for better stable boundary-layer parameterizations. To this end, two new turbulence parameterizations based on the measurements of the Vertical Transport and Mixing (VTMX) field campaign were implemented and evaluated in MM5. A unique aspect of this parameterization is the use of a stability-dependent turbulent Prandtl number that allows momentum to be transported by the internal waves, while heat diffusion is impeded by the stratification. This improvement alleviates the problem of over-prediction of heat diffusion under stable conditions, which is a characteristic of conventional atmospheric boundary-layer schemes, such as the Medium Range Forecast (MRF) and Blackadar schemes employed in MM5. The predictions made with the new PBL scheme for the complex terrain airshed of Salt Lake City were compared with those made with a default scheme of MM5, and with observations made during the VTMX campaign. The new schemes showed an improvement in predictions, particularly for the nocturnal near-surface temperature. Surface wind predictions also improved slightly, but not to the extent of temperature predictions. The default MRF scheme showed a significantly higher surface temperature than observed, which could be attributed to the enhanced vertical heat exchange brought about by its turbulence parameterization. The modified parameterizations reduced the surface sensible heat flux, thus enhancing the strength of the near-surface inversion and lowering the temperature towards the observed values.     

Turbulent dissipation of cold air lake in a basin

Summary Problems of turbulent dissipation of a cold air lake (CAL) and the inversion layer bordering CAL on the upper boundary are presented and studied with the compound model. In wintertime such cold air lakes can persist for days even if rather strong winds are blowing above them. The required conditions for CAL dissipation are removed processes of its formation or maintenance, as well as a sufficiently strong invasion of turbulence in the inversion layer from above down-wards. By this, the inversion layer at first becomes stronger and dissipation is stopped, until the increase of turbulent kinetic energy of the upper flow enables further dissipation. Such turbulent dissipation process is shown by the model for typical conditions and for different initial values of the relevent variables.With 6 Figures     

白天谷地逆温层生消过程的数值模拟

本文用二维数值模式模拟了谷地白天边界层的发展过程.通过求解大气动力.热力学方程组.模拟出了谷地上空的温度场、流场以及在两种不同层结条件下的湍流动能q~2的分布特征.从而对白天谷地上稳定层结(或逆温层)的生消机制进行了初步探讨.模拟出的温度分布廓线等与有关实测资料基本一致.     

Planetary boundary-layer turbulence studies from acoustic echo sounder and in-situ measurements

Clear-air plume and wave-like structures are revealed in the atmospheric boundary layer by combined acoustic remote sensing and meteorological tower measurements. The magnitude of turbulent production and dissipation plus properties of velocity and temperature spectra determined from the tower measurements are well correlated with phenomena indicated by the acoustic sounder. Interpretation of either set of records is greatly enhanced by the other. For example, the onset of a sudden burst of turbulent production from the tower measurements may correspond to plume passage or breaking of stable waves recorded by the acoustic echo sounder.     

大气运动的动能在湍流边界层中的耗散

在大气能量学中,大气平均运动的动能通过湍流摩擦力转化为脉动运动的能量并最终转化为内能的“变性能量”是大气动能的“汇”,正确估计其大小对研究动能的收支,以至天气系统的生衰,发展均有一定意义。此项由于与湍流通量有关,因而难于从大尺度气象参数计算。以往的研究中,或是将动能平衡方程中其余各项算出,而将此项作为余项估计出,或用地转风由粗浅的边界层模式算出。     

Molecular dissipation of turbulent fluctuations in the convective mixed layer part II: Height variations of characteristic time scales and experimental test of molecular dissipation models

Extensive turbulence measurements from the Limagne and Beauce experiments were used to compute a characteristic time scale of the turbulence field (Τ = second moment/dissipation rate) for turbulent kinetic energy, temperature and humidity variances, and temperature-humidity covariance. The height variations of these time scales were analysed. The characteristic half-time scale Τ/2 of the turbulent velocity field was found, as expected, to be of the same order of magnitude as the large-eddy time scale Τ _L = Z_i/w*, showing that the turbulence structure is controlled by large eddies in the bulk of the mixed layer. The increase of Τ/2 above z/Z _i ～- 0.7 implies, however, that this time scale is no longer relevant to destruction of turbulent kinetic energy in the statically stable region with negative heat fluxes. An effective time scale Τ_eff, introduced by Zeman (1975), has been computed and its behaviour discussed. The scales for θ′ ², q′², and θ′q′ were found to be much shorter than Τ. Furthermore, a significant difference in behaviour was also revealed between the characteristic time scales of temperature and humidity fields in the stable layer. By using these experimental estimates, we tested some of the models for molecular dissipations, which are currently in use in higher order closure atmospheric boundary-layer models. The parameterized dissipation rates for θ′ ², and q′ ² agree well qualitatively with experimental estimates in the bulk of the mixed layer. In the stable layer, however, the parameterized dissipation rate ε _θ tends to become larger than the experimental ones although the parameterized dissipation rate ε _q still agrees with the experimental ones. For the molecular dissipation of θ′q′, this current model becomes physically inconsistent in the middle part of the mixed layer, because this term may become a production term for temperature-humidity covariance.     

On Monin–Obukhov Similarity In The Stable Atmospheric Boundary Layer

Atmospheric measurements from several field experiments have been combined to develop a better understanding of the turbulence structure of the stable atmospheric boundary layer. Fast response wind velocity and temperature data have been recorded using 3-dimensional sonic anemometers, placed at severalheights (1 m to 4.3 m) above the ground. The measurements wereused to calculate the standard deviations of the three components of the windvelocity, temperature, turbulent kinetic energy (TKE) dissipation andtemperature variance dissipation. These data were normalized and plottedaccording to Monin–Obukhov similarity theory. The non-dimensional turbulencestatistics have been computed, in part, to investigate the generalapplicability of the concept of z-less stratification for stable conditions. From the analysis of a data set covering almost five orders ofmagnitude in the stability parameter = z/L (from near-neutral tovery stable atmospheric stability), it was found that this concept does nothold in general. It was only for the non-dimensional standard deviation oftemperature and the average dissipation rate of turbulent kinetic energythat z-less behaviour has been found. The other variables studied here(non-dimensional standard deviations of u, v, and w velocity components and dissipation of temperature variance) did not follow the concept of z-less stratification for the very stable atmospheric boundary layer. An imbalance between production and dissipation of TKE was found for the near-neutral limit approached from the stable regime, which matches with previous results for near-neutral stability approached from the unstable regime.     

Modelling the arctic convective boundary-layer with different turbulence parameterizations

Different parameterizations of subgrid-scale fluxes are utilized in a nonhydrostatic and anelastic mesoscale model to study their influence on simulated Arctic cold air outbreaks. A local closure, a profile closure and two nonlocal closure schemes are applied, including an improved scheme, which is based on other nonlocal closures. It accounts for continuous subgrid-scale fluxes at the top of the surface layer and a continuous Prandtl number with respect to stratification. In the limit of neutral stratification the improved scheme gives eddy diffusivities similar to other parameterizations, whereas for strong unstable stratifications they become much larger and thus turbulent transports are more efficient. It is shown by comparison of model results with observations that the application of simple nonlocal closure schemes results in a more realistic simulation of a convective boundary layer than that of a local or a profile closure scheme. Improvements are due to the nonlocal formulation of the eddy diffusivities and to the inclusion of heat transport, which is independent of local gradients (countergradient transport).     

Molecular dissipation of turbulent fluctuations in the convective mixed layer part I: Height variations of dissipation rates

During the Limagne and Beauce experiments, the INAG-IGN Aerocommander FL 280 aircraft made extensive ‘in situ’ measurements of turbulent fluctuations in diurnally evolving convective boundary layers. In this paper, these measurements were used to investigate characteristics of the molecular dissipation of turbulent fluctuations through the mixed layer and well into the overlying stable layer. The dimensionless dissipation rates of turbulent kinetic energy, temperature and humidity variances, and temperature-humidity covariance (ψ, ψ_θ, ψ _qand ψ _θq) were computed and their height variations analysed. The behaviour of the dissipation rate ψ was found to differ significantly from those observed for the other rates. In the lowest region of the mixed layer, ψ does not obey the local free convection prediction. Instead, it follows practically a relationship similar to the one established in the surface layer by Wyngaard et al. (1971). The dissipation rate ψ remains fairly constant in the bulk of the mixed layer (0.3 ≤ z/Z _i≤ 0.8) and shows a very rapid decrease above the inversion. These results confirm those reported previously from the Minnesota and Ashchurch data by Kaimal et al. (1976), Caughey and Palmer (1979), etc. The height variations for the other dissipation rates were found to obey, as expected, the (z/Z _i)^-4/3 decrease predicted under the local free convection similarity hypothesis in the lowest region of the mixed layer. This region extends to the height z/Z _i- 0.4, 0.1, and 0.3, respectively, for ψ_θ, ψ_q, and ψ_θq. Above these levels, the dissipation rates ψ_θ and ψ_q show, on average, a slight increase to reach peak-values near the mixed-layer top, while the ‘dissipation’ rate ψ _θqchanges sign from positive to negative around the height z/Z _i, - 0.7. These characteristics confirm the fact that the structures of temperature and humidity fluctuations are considerably affected by their entrainment-induced fluctuations. Therefore, an attempt has been made to non-dimensionalize the dissipation rates near the mixed-layer top with the interfacial scaling factors.     

Note on the viscous heating term in the temperature variance equation

Experimental results in the atmospheric surface layer indicate that the viscous heating term in the equation for the mean-square value of the turbulent temperature fluctuation is small, but the correlation between temperature and turbulent energy dissipation fluctuations is significant.     

Parameterization and Explicit Modeling of Cloud Microphysics: Approaches,Challenges, and Future Directions

Cloud microphysical processes occur at the smallest end of scales among cloud-related processes and thus must be parameterized not only in large-scale global circulation models(GCMs) but also in various higher-resolution limited-area models such as cloud-resolving models(CRMs) and large-eddy simulation(LES) models. Instead of giving a comprehensive review of existing microphysical parameterizations that have been developed over the years, this study concentrates purposely on several topics that ...     

Direct near-surface measurements of sensible heat fluxes in the Arctic tundra applying eddy covariance and laser scintillometry��the Arctic Turbulence Experiment 2006 on Svalbard (ARCTEX-2006)

Recent climate warming in the Arctic requires improvements in permafrost and carbon cycle monitoring, accomplished here by setting up long-term observation sites with high-quality in situ measurements of turbulent atmospheric energy fluxes applying the eddy covariance method and/or laser scintillometry in Arctic landscapes. Accurate quantification and well-adapted parameterizations of turbulent energy fluxes, e.g., during neutral to stable stratified conditions, are a fundamental problem in soil?Csnow?Cice?Cvegetation?Catmosphere interaction studies. We present results from the Arctic Turbulence Experiment (ARCTEX-2006) performed on the island of Svalbard, Norway, during the winter/spring transition 2006 that focus on data correction and quality assessment, on synoptic weather conditions, as well as site-specific micrometeorological features. A quality assessment and data correction adapted to the environmental conditions of polar regions demonstrates that specific measurement errors common at a high Arctic landscape could be minimized. We discuss the role of the intermittency of the turbulent atmospheric fluctuation of momentum and scalars, the existence of a disturbed vertical temperature profile (sharp inversion layer) close to the surface, and the relevance of possible free convection events for the snow or ice melt in the Arctic spring at Svalbard. Recommendations and improvements regarding the interpretation of eddy flux data as well as the arrangement of the instrumentation under polar distinct exchange conditions and (extreme) weather situations are presented.     

E − ε − 〈θ 2〉 turbulence closure model for an atmospheric boundary layer including the urban canopy

A modified three-parameter model of turbulence for a thermally stratified atmospheric boundary layer (ABL) is presented. The model is based on tensor-invariant parameterizations for the pressure–strain and pressure–temperature correlations that are more complete than the parameterizations used in the Mellor–Yamada model of level 3.0. The turbulent momentum and heat fluxes are calculated with explicit algebraic models obtained with the aid of symbol algebra from the transport equations for momentum and heat fluxes in the approximation of weakly equilibrium turbulence. The turbulent transport of heat and momentum fluxes is assumed to be negligibly small in this approximation. The three-parameter $E - \varepsilon - {\left\langle {\theta ^{2} } \right\rangle }$ model of thermally stratified turbulence is employed to obtain closed-form algebraic expressions for the fluxes. A computational test of a 24-h ABL evolution is implemented for an idealized two-dimensional region. Comparison of the computed results with the available observational data and other numerical models shows that the proposed model is able to reproduce both the most important structural features of the turbulence in an urban canopy layer near the urbanized ABL surface and the effect of urban roughness on a global structure of the fields of wind and temperature over a city. The results of the computational test for the new model indicate that the motion of air in the urban canopy layer is strongly influenced by mechanical factors (buildings) and thermal stratification.