一维大气边界层二阶闭合的有限元数值模式 I. 对流边界层模拟

利用有限元法用于二阶湍流闭合大气对流边界层模式,并模拟了Wangara试验33天的对流边界层的发展过程,计算结果表明,该模式能真实的模拟出对流边界层的发展过程,并且准确地反映对流边界层中湍流扩散和输送过程。     

Laboratory Simulation Of Vertical Plume Dispersion Within A Convective Boundary Layer

Mesoscale measurements of the vertical dispersion coefficient ₂ by using a composite turbulence water tank were validated through a comparison with CONDORS (Convective Diffusion Observed with Remote Sensors) field data, and were analysed with respect to the intensity of the thermal flux, mechanical turbulence, and plume release height.It seems possible to correct the plume _z values for different release heights below 0.5z_i (z_i is the mixing height) by applying an equation expressing the height dependency of turbulence intensity. The downwind distance where the plume's mass centre height approaches its final level was also analysed with respect to the above three parameters, and an empirical equation to estimate the downwinddistance derived.     

大气对流边界层发展的模拟研究

室内水槽模拟是大气边界层研究的一种重要手段。利用室内模拟水槽对大气边界层的发展进行了模拟,通过处理平均温度廓线和光斑图像得到了对流边界层顶部位置h2和边界层高度zi。结果表明,不同测量方法得到的结果一致性很好,与实际大气的边界层发展情况也较为接近。同时,根据试验情况确定初始条件和边界条件,使用边界层参数化模型进行了数值模拟,其结果与室内模拟的结果也较吻合。     

Fluid Modelling Of Atmospheric Dispersion In The Convective Boundary Layer

A laboratory convection tank has been established following thepioneering work of Willis and Deardorff, but with many improvements and enhancements thattake advantage of modern technology. The main emphasis in the current design was toprovide the ability to conduct a virtually unlimited number of realizations under essentiallyidentical conditions in order to obtain reliable statistics on the dispersion of plumes and puffsreleased within the simulated atmospheric convective boundary layer. Described herein is the tankitself and its auxiliary systems, including a laser-induced-fluorescence and video-imaging system for makingnon-intrusive, full-field measurements of concentrations, and the interfacing of varioussubsystems with a master controller that automates essentially all operation and measurement functions.The current system provides unprecedented resolution, control, and data volumes. Exampleresults are presented from two types of releases: continuous plumes and instantaneous puffs.These data sets clearly show penetration of the highly buoyant plumes and puffs into theinversion above the convective boundary layer, gravity spreading within the inversion, andrapid diffusion within the mixed layer. They also show extreme `spottiness' in the instantaneousconcentration cross-sections.     

Dissimilarity of Scalar Transport in the Convective Boundary Layer in Inhomogeneous Landscapes

A land-surface model (LSM) is coupled with a large-eddy simulation (LES) model to investigate the vegetation-atmosphere exchange of heat, water vapour, and carbon dioxide (CO₂) in heterogeneous landscapes. The dissimilarity of scalar transport in the lower convective boundary layer is quantified in several ways: eddy diffusivity, spatial structure of the scalar fields, and spatial and temporal variations in the surface fluxes of these scalars. The results show that eddy diffusivities differ among the three scalars, by up to 10–12%, in the surface layer; the difference is partly attributed to the influence of top-down diffusion. The turbulence-organized structures of CO₂ bear more resemblance to those of water vapour than those of the potential temperature. The surface fluxes when coupled with the flow aloft show large spatial variations even with perfectly homogeneous surface conditions and constant solar radiation forcing across the horizontal simulation domain. In general, the surface sensible heat flux shows the greatest spatial and temporal variations, and the CO₂ flux the least. Furthermore, our results show that the one-dimensional land-surface model scheme underestimates the surface heat flux by 3–8% and overestimates the water vapour and CO₂ fluxes by 2–8% and 1–9%, respectively, as compared to the flux simulated with the coupled LES-LSM.     

Analytical Model of Roll Vortices in the Convective Boundary Layer

Analytical solutions of convective waves in the convective boundary layer (CBL) were obtained with two-layer linearized atmospheric equations including Rayleigh friction, which represents the turbulent viscosity in the lower CBL. The analytical model shows that the interaction between the convection in the lower layer and gravity waves in the upper layer is one of the causes for the formation of convective bands. The flow and temperature fields obtained by the analytical model present the main characteristics of convective bands found in field observations. We have also investigated the influences of atmospheric conditions on the characteristics of the bands. Results accord with previous knowledge about these phenomena.     

Doppler Lidar Measurements of Vertical Velocity Spectra in the Convective Planetary Boundary Layer

We utilized a Doppler lidar to measure spectra of vertical velocity w from 390m above the surface to the top of the daytime convective boundary layer (CBL). The high resolution 2μm wavelength Doppler lidar developed by the NOAA Environmental Technology Laboratory was used to detect the mean radial velocity of aerosol particles. It operated continuously during the daytime in the zenith-pointing mode for several days in summer 1996 during the Lidars-in-Flat-Terrain experiment over level farmland in central Illinois, U.S.A. The temporal resolution of the lidar was about 1 s, and the range-gate resolution was about 30m. The vertical cross-sections were used to calculate spectra as a function of height with unprecedented vertical resolution throughout much of the CBL, and, in general, we find continuity of the spectral peaks throughout the depth of the CBL. We compare the observed spectra with previous formulations based on both measurements and numerical simulations, and discuss the considerable differences, both on an averaged and a case-by-case basis. We fit the observed spectra to a model that takes into account the wavelength of the spectral peak and the curvature of the spectra across the transition from low wavenumbers to the inertial subrange. The curvature generally is as large or larger than the von Kármán spectra. There is large case-to-case variability, some of which can be linked to the mean structure of the CBL, especially the mean wind and the convective instability. We also find a large case-to-case variability in our estimates of normalized turbulent kinetic energy dissipation deduced from the spectra, likely due for the most part to a varying ratio of entrainment flux to surface flux. Finally, we find a relatively larger contribution to the low wavenumber region of the spectra in cases with smaller shear across the capping inversion, and suggest that this may be due partly to gravity waves in the inversion and overlying free atmosphere.     

An Alternative Eddy-Viscosity Model for the Horizontally Uniform Atmospheric Boundary Layer

This paper explores the utility of specifying the eddy viscosity for the horizontally uniform boundary layer as the product of the variance of vertical velocity and an empirical time scale τ _w , as opposed to the more usual formulation where k is the turbulent kinetic energy (TKE), λ _k is a length scale and α is a dimensionless coefficient. Simulations were compared with the observations on Day 33 of the Wangara experiment, and with a plausible specification of τ _w (or λ _k ) each model simulated convective boundary-layer development reasonably well, although the closure produced a more realistic width for the entrainment layer. Under the light winds of Day 33, and with the onset of evening cooling, an excessively shallow and strongly-stratified nocturnal inversion developed, and limited its own further deepening. Boundary-layer models that neglect radiative heat transport and parametrize convective transport by eddy viscosity closure are prone to this runaway (unstable) feedback when forced by a negative (i.e. downward) surface flux of sensible heat.     

A Similarity Model of Subfilter-Scale Energy for Large-Eddy Simulations of the Atmospheric Boundary Layer

A scale-similarity model to estimate the subfilter-scale energy using the trace of the Leonard stress tensor is proposed and evaluated for large-eddy simulations of the atmospheric boundary layer (ABL). The model is derived from a stability-dependent model of the energy spectrum in the ABL, which accounts for the effects of buoyancy and mean shear as a function of z/L, the Monin–Obukhov stability variable. An a priori test using ABL turbulence data demonstrates that the model has accurate performance for dimensionless filter widths of Δ/z = 2, 1, and 0.5 for stabilities of −1 ≤ z/L ≤ 0.5, and improves considerably upon a similar model that is derived using an infinite κ ^−5/3 spectrum. This improvement is especially significant in the first several grid points near the surface in large-eddy simulations of the ABL, where Δ/z is necessarily large. The modelling procedure is then extended to develop a similarity model for the subfilter-scale scalar variance; it is shown to have robust performance for temperature.     

Single-Column Model Intercomparison for a Stably Stratified Atmospheric Boundary Layer

The parameterization of the stably stratified atmospheric boundary layer is a difficult issue, having a significant impact on medium-range weather forecasts and climate integrations. To pursue this further, a moderately stratified Arctic case is simulated by nineteen single-column turbulence schemes. Statistics from a large-eddy simulation intercomparison made for the same case by eleven different models are used as a guiding reference. The single-column parameterizations include research and operational schemes from major forecast and climate research centres. Results from first-order schemes, a large number of turbulence kinetic energy closures, and other models were used. There is a large spread in the results; in general, the operational schemes mix over a deeper layer than the research schemes, and the turbulence kinetic energy and other higher-order closures give results closer to the statistics obtained from the large-eddy simulations. The sensitivities of the schemes to the parameters of their turbulence closures are partially explored.     

A Parameterization of Third-order Moments for the Dry Convective Boundary Layer

We describe one-dimensional (1D) simulations of the countergradient zone of mean potential temperature observed in the convective boundary layer (CBL). The method takes into account the third-order moments (TOMs) in a turbulent scheme of relatively low order, using the turbulent kinetic energy equation but without prognostic equations for other second-order moments. The countergradient term is formally linked to the third-order moments and , and a simple parameterization of these TOMs is proposed. It is validated for several cases of a dry CBL, using large-eddy simulations that have been realized from the MESO-NH model. The analysis of the simulations shows that TOMs are responsible for the inversion of the sign of in the higher part of the CBL, and budget analysis shows that the main terms responsible for turbulent fluxes and variances are now well reproduced.     

Coherence and Scale of Vertical Velocity in the Convective Boundary Layer from a Doppler Lidar

We utilized a Doppler lidar to measure integral scale and coherence of vertical velocity w in the daytime convective boundary layer (CBL). The high resolution 2 μm wavelength Doppler lidar developed by the NOAA Environmental Technology Laboratory was used to detect the mean radial velocity of aerosol particles. It operated continuously in the zenith-pointing mode for several days in the summer 1996 during the “Lidars in Flat Terrain” experiment over level farmland in central Illinois. We calculated profiles of w integral scales in both the alongwind and vertical directions from about 390 m height to the CBL top. In the middle of the mixed layer we found, from the ratio of the w integral scale in the vertical to that in the horizontal direction, that the w eddies are squashed by a factor of about 0.65 as compared to what would be the case for isotropic turbulence. Furthermore, there is a significant decrease of the vertical integral scale with height. The integral scale profiles and vertical coherence show that vertical velocity fluctuations in the CBL have a predictable anisotropic structure. We found no significant tilt of the thermal structures with height in the middle part of the CBL.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

A Simple Model Of The Convective Internal Boundary Layer And Its Application To Surface Heat Flux Estimates Within Polynyas

A simple model of the convective (thermal) internalboundary layer has been developed for climatologicalstudies of air-sea-ice interaction, where in situobservations are scarce and first-order estimates ofsurface heat fluxes are required. It is amixed-layer slab model, based on a steady-statesolution of the conservation of potentialtemperature equation, assuming a balance betweenadvection and turbulent heat-flux convergence. Boththe potential temperature and the surface heat fluxare allowed to vary with fetch, so the subsequentboundary-layer modification alters the fluxconvergence and thus the boundary-layer growth rate.For simplicity, microphysical and radiativeprocesses are neglected.The model is validated using several case studies.For a clear-sky cold-air outbreak over a coastalpolynya the observed boundary-layer heights,mixed-layer potential temperatures and surface heatfluxes are all well reproduced. In other cases,where clouds are present, the model still capturesmost of the observed boundary-layer modification,although there are increasing discrepancies withfetch, due to the neglected microphysical andradiative processes. The application of the model toclimatological studies of air-sea interaction withincoastal polynyas is discussed.     

Vertical Structure in the Marine Atmospheric Boundary Layer and its Implication for the Inertial Dissipation Method

The structure of the marine atmospheric boundarylayer and the validity ofMonin–Obukhov similarity theory over the seahave been investigated using longterm measurements. Three levels of turbulencemeasurements (at 10 m, 18 mand 26 m) at Östergarnsholm in themiddle of the Baltic Sea have beenanalysed. The results show that turbulentparameters have a strong dependenceon the actual height due to wave influence.The wind profile and thus thenormalised wind gradient are very sensitiveto wave state. The lower part of theboundary layer can be divided into three heightlayers, a wave influenced layerclose to the surface, a transition layer andan undisturbed ordinary surfacelayer; the depth of the layers is determinedby the wave state. This heightstructure can, however, not be found for thenormalised dissipation, which is onlya function of the stability, except duringpronounced swell where the actualheight also has to be accounted for. Theresults have implications for the heightvariation of the turbulent kinetic energy(TKE) budget. Thus, the imbalancebetween production and dissipation willalso vary with height according to thevariation of wave state. This, in turn,will of course have strong implicationsfor the inertial dissipation method, inwhich a parameterisation of the TKEbudget is used.     

Equilibrium Evaporation and the Convective Boundary Layer

A theory is developed for surface energy exchanges in well-mixed, partlyopen systems, embracing fully open and fully closed systems as limits.Conservation equations for entropy and water vapour are converted intoan exact rate equation for the potential saturation deficit D in a well-mixed, partly open region. The main contributions to changes in D arise from (1) the flux of D at the surface, dependent on a conductance g_q that is a weighted sum of the bulk aerodynamic and surface conductances; and (2) the exchange flux of D with the external environment by entrainment or advection, dependent on a conductance g_e that is identifiable with the entrainment velocity when the partly open region is a growing convective boundary layer (CBL). The system is fully open when g_e/g_q , and fully closed when g_e/g_q 0. The equations determine the steady state surface energy balance (SEB) in a partly open system, the associated steady-state deficit, and the settling time scale needed to reach the steady state. The general result for the steady-state SEB corresponds to the equations of conventional combination theory for the SEB of a vegetated surface, with the surface-layer deficit replaced by the external deficit and with g_q replaced by the series sum (g_q ^-1 + g_e ^-1)^-1. In the fully open limit D is entirely externally prescribed, while in the fully closed limit, D is internally determined and the SEB approaches thermodynamic equilibrium energy partition. In the case of the CBL, the conductances g_q and g_e are themselves functions of D through short-term feedbacks, induced by entrainment in the case of g_e and by both physiological and aerodynamic (thermal stability) processes in the case of g_q. The effects of these feedbacks are evaluated. It is found that a steady-state CBL is physically achievable only over surfaces with at least moderate moisture availability; that entrainment has a significant accelerating effect on equilibration; that the settling time scale is well approximated by h/(g_q + g_e), where h is the CBL depth; and that this scale is short enough to allow a steady state to evolve within a semi-diurnal time scale only when h is around 500 m or less.     

A Variable Mesh Spacing for Large-Eddy Simulation Models in the Convective Boundary Layer

A variable vertical mesh spacing for large-eddy simulation (LES) models in a convective boundary layer (CBL) is proposed. The argument is based on the fact that in the vertical direction the turbulence near the surface in a CBL is inhomogeneous and therefore the subfilter-scale effects depend on the relative location between the spectral peak of the vertical velocity and the filter cut-off wavelength. From the physical point of view, this lack of homogeneity makes the vertical mesh spacing the principal length scale and, as a consequence, the LES filter cut-off wavenumber is expressed in terms of this characteristic length scale. Assuming that the inertial subrange initial frequency is equal to the LES filter cut-off frequency and employing fitting expressions that describe the observed convective turbulent energy one-dimensional spectra, it is feasible to derive a relation to calculate the variable vertical mesh spacing. The incorporation of this variable vertical grid within a LES model shows that both the mean quantities (and their gradients) and the turbulent statistics quantities are well described near to the ground level, where the LES predictions are known to be a challenging task.     

Impacts of Aerosol Shortwave Radiation Absorption on the Dynamics of an Idealized Convective Atmospheric Boundary Layer

We investigated the impact of aerosol heat absorption on convective atmospheric boundary-layer (CBL) dynamics. Numerical experiments using a large-eddy simulation model enabled us to study the changes in the structure of a dry and shearless CBL in depth-equilibrium for different vertical profiles of aerosol heating rates. Our results indicated that aerosol heat absorption decreased the depth of the CBL due to a combination of factors: (i) surface shadowing, reducing the sensible heat flux at the surface and, (ii) the development of a deeper inversion layer, stabilizing the upper CBL depending on the vertical aerosol distribution. Steady-state analytical solutions for CBL depth and potential temperature jump, derived using zero-order mixed-layer theory, agreed well with the large-eddy simulations. An analysis of the entrainment zone heat budget showed that, although the entrainment flux was controlled by the reduction in surface flux, the entrainment zone became deeper and less stably stratified. Therefore, the vertical profile of the aerosol heating rate promoted changes in both the structure and evolution of the CBL. More specifically, when absorbing aerosols were present only at the top of the CBL, we found that stratification at lower levels was the mechanism responsible for a reduction in the vertical velocity and a steeper decay of the turbulent kinetic energy throughout the CBL. The increase in the depth of the inversion layer also modified the potential temperature variance. When aerosols were present we observed that the potential temperature variance became significant already around $0.7z_i$ (where $z_i$ is the CBL height) but less intense at the entrainment zone due to the smoother potential temperature vertical gradient.     

Improving Non-local Parameterization of the Convective Boundary Layer

We propose improvements in the “non-local” parameterization scheme of the convective boundary layer. The countergradient terms for components of the momentum fluxes are introduced in a form analogous to those for other scalars. The scheme also includes explicit expressions for entrainment fluxes of momentum, temperature, and humidity. A simplified procedure for calculating the boundary-layer height is proposed, consisting of two steps: the evaluation of the convection level, followed by the assessment of the depth of the interfacial layer.     

Structure Inclination Angles in the Convective Atmospheric Surface Layer

Two-point correlations of the fluctuating streamwise velocity are examined in the atmospheric surface layer over the salt flats of Utah’s western desert, and corresponding structure inclination angles are obtained for neutral, stable and unstable conditions. The neutral surface-layer results supplement evidence for the invariance of the inclination angle given in Marusic and Heuer (Phys Rev Lett 99:114504, 2007). In an extension of those results it is found that the inclination angle changes drastically under different stability conditions in the surface layer, varying systematically with the Monin–Obukhov stability parameter in the unstable regime. The variation is parametrized and subsequently can be used to improve existing near-wall models in the large-eddy simulation of the atmospheric surface layer.