Simulation of Atmospheric Boundary Layer Characteristics during Indian Summer Monsoon using Observations from Monsoon Trough Boundary Layer Experiment at Jodhpur, India

The diurnal structure of the boundary layer during Indian summer monsoon period is studied using a one-dimensional meteorological boundary layer model and the observations collected from the Monsoon Trough Boundary Layer Experiment conducted in 1990 at Jodhpur, India. The model was initialized with the observed temperature profiles at 0530 LST on 17 July, 1990 at Jodhpur and was run for 26 hours. The study is carried out with a geostrophic wind speed of 9.5 m s⁻¹ corresponding to the strong wind simulation. The mean thermodynamic and wind structure simulated by the model are in good agreement with those observed from 30 m tower. The computed surface layer characteristics such as the surface fluxes, TKE and standard deviations of velocity components are found to be reasonably in good agreement with those based on turbulence measurements. The shear and buoyancy budget computed from the model are also compared with the turbulence measurements. The integrated cooling budget in the nocturnal boundary layer is examined.     

Stable-boundary-layer regimes from the perspective of the low-level jet

This paper reviews results from two field studies of the nocturnal stable atmospheric boundary layer (SBL) over the Great Plains of the United States. Data from a scanning remote-sensing system, a High-Resolution Doppler Lidar (HRDL), provided measurements of mean and turbulent wind components at high spatial and temporal resolution through the lowest 500–1000 m of the atmosphere. This data set has allowed the characteristics of the low-level jet (LLJ) maximum (speed, height, direction) to be documented through entire nights. LLJs form after sunset and produce strong shear in the layer below the LLJ maximum or nose, which is a source of turbulence and mixing in the SBL. Simultaneous HRDL measurements of turbulence quantities related to turbulence kinetic energy (TKE) has allowed the turbulence in the subjet layer to be related to LLJ properties. Turbulence structure was found to be a function of the bulk stability of the subjet layer. For the strong-LLJ (> 15 m s⁻¹), weakly stable cases the strength of the turbulence is proportional to the strength of the LLJ. For these cases with nearly continuous turbulence in the subjet layer, low-level jet scaling, in which lengths are scaled by the LLJ height and velocity variables are scaled by the LLJ speed, was found to be appropriate. For the weak-wind (< 5 m s⁻¹ in the lowest 200 m), very stable boundary layer (vSBL), the boundary layer was found to be very shallow (sometimes < 10 m deep), and turbulent fluxes between the earth’s surface and the atmosphere were found to be essentially shut down. For more intermediate wind speeds and stabilities, the SBL shows varying degrees of intermittency due to various mechanisms, including shearinstability and other gravity waves, density currents, and other mesoscale disturbances.     

A Model Study of the Strong and Weak Wind,Stably Stratified Nocturnal Boundary Layer: Influence of Gentle Slopes

With the exception of intermittency and waves, a brief review of the observed and modeled mean structure of the nocturnal boundary layer (NBL) is presented. The effect of gentle slopes on strong and weak wind NBL was investigated here using a one-dimensional model, with a simple correction term to account for the slope effects, identical to the one used by Brost and Wyngaard (1978). The study indicates that the wind profiles, temperature profiles and surface layer turbulence characteristics are extremely sensitive to the imposed geostrophic wind when small slopes are present especially for light winds. This is due to the complex interaction between the buoyancy driven slope flow and the imposed geostrophic wind that in turn influence the shear generation of turbulence. Finally, the current issues in the modeling of weak wind boundary layer are discussed.     

斜压对流边界层中风廓线和湍流特征的研究

在美国Ｋａｎｓａｓ和Ｎｅｂｒａｓｋａ交界处于１９９２年进行的ＳＴＯＲＭ－ＦＥＳＴ（ＳＴｏｒｍ－ｓｃａｌｅ Ｏｐｅｒａｔｉｏｎａｌ　ａｎｄ　Ｒｅｓｅａｒｃｈ　Ｍｅｔｅｏｒｏｌｏｇｙ－Ｆｒｏｎｔｓ　Ｅｘｐｅｒｉｍｅｎｔ Ｓｙｓｔｅｍ　Ｔｅｓｔ）资料基础上,分析了 斜压对流边界层中温度、湿度和风速的垂直分布．分析表明,在对流边界层中温度和湿度常 显示很好的混合;但风速分量有时混合很好,有时则存在切变．计算了存在和不存在风切变 两种情况下的湍流动能收支．最后,讨论了对流边界层中风切变形成的可能原因．     

对流槽湍流涡旋结构特征的小波分析

为研究夹卷层的湍流特性,用对流槽模拟大气对流边界层并做湍流温度测量.能谱分析表明,混合层湍流谱幂律接近“-5/3”次律;而夹卷层湍流谱幂律有明显的分段特征.利用离散正交小波对其湍流脉动信号进行去噪和尺度分解,将原始信号分成均匀各向同性小尺度涡旋成分和大尺度涡旋成分,并对分解后的信号做能谱分析.结果表明,不论是在混合层还是在夹卷层,小涡都能很好地满足Kolmogorov的“-5/3”次律;而大尺度涡旋则在混合层和夹卷层中表现出不同的特征.     

Turbulence structure of the boundary layer below marine clouds in the SOFIA experiment

The SOFIA (Surface of the Ocean: Flux and Interaction with the Atmosphere) experiment, included in the ASTEX (Atlantic Stratocumulus Transition EXperiment) field program, was conducted in June 1992 in the Azores region in order to investigate air-sea exchanges, as well as the structure of the atmospheric boundary layer and its capping low-level cloud cover. We present an analysis of the vertical structure of the marine atmospheric boundary layer (MABL), and especially of its turbulence characteristics, deduced from the aircraft missions performed during SOFIA. The meteorological situations were characteristic of a temperate latitude under anticyclonic conditions, i.e., with weak to moderate winds, weak surface sensible heat flux, and broken capping low-altitude cloud cover topped by a strong trade inversion. We show that the mixed layer, driven by the surface fluxes, is decoupled from the above cloud layer. Although weak, the surface buoyancy flux, and the convective velocity scale deduced from it, are relevant for scaling the turbulence moments. The mixed layer then follows the behavior of a continental convective boundary layer, with the exception of the entrainment process, which is weak in the SOFIA data. These results are confirmed by conditional sampling analysis, which shows that the major turbulence source lies in the buoyant moist updrafts at the surface.     

西藏高原斜压对流边界层风、温、湿廓线特征

利用1998年夏季第2次青藏高原大气科学试验当雄观测站的边界层观测资料以及拉萨、改则和武汉等地探空资料,分析讨论西藏高原斜压对流边界层风、温、湿廓线的特征. 研究结果表明,高原地区白天对流边界层发展可高达2200m,显著超过中纬度平原地区和海面上对流边界层高度. 高原对流边界层中温度廓线具有较好的混合特征,湿度廓线有时在某一高度上出现湿度极大值. 高原对流边界层内热量和水汽收支分析表明,水平平流作用对边界层结构具有重要作用. 在对流边界层中平均风速垂直分布存在风切变现象. 水平温度梯度形成较强的斜压性是形成边界层风切变的主要原因.     

Turbulence in the presence of internal waves in the bottom boundary layer of the California inner shelf

Turbulence measurements were collected in the bottom boundary layer of the California inner shelf near Point Sal, CA, for 2 months during summer 2015. The water column at Point Sal is stratified by temperature, and internal bores propagate through the region regularly. We collected velocity, temperature, and turbulence data on the inner shelf at a 30-m deep site. We estimated the turbulent shear production (P), turbulent dissipation rate (ε), and vertical diffusive transport (T), to investigate the near-bed local turbulent kinetic energy (TKE) budget. We observed that the local TKE budget showed an approximate balance (P?≈?ε) during the observational period, and that buoyancy generally did not affect the TKE balance. On a finer resolution timescale, we explored the balance between dissipation and models for production and observed that internal waves did not affect the balance in TKE at this depth.     

MYJ和YSU方案对WRF边界层气象要素模拟的影响

研究新一代中尺度气象模式WRF中两种大气边界层方案(MYJ,YSU)对沈阳冬季大气边界层结构模拟的影响,重点分析温度层结、低层风场、边界层高度等对污染物扩散有重要影响的气象要素.和观测数据的比对表明WRF基本能够模拟出温度风速的日变化特征,但模拟风速偏大.YSU方案由于模拟的边界层顶卷挟和边界层内混合作用较强,夜间接地逆温强度低于MYJ方案,逆温维持时间比MYJ方案短4小时,同时模拟边界层高度也高于MYJ方案,有利于污染物垂直扩散.边界层高度的3种计算方法中,湍流动能方法计算的边界层高度最高,Richardson数方法次之,位温方法得到的高度最低.Richardson数方法对临界值的选取较敏感.     

边界层参数化方案对边界层热力和动力结构特征影响的比较

本文利用高分辨率中尺度WRF模式,通过改变边界层参数化方案进行多组试验,评估该模式对美国北部森林地区边界层结构的模拟能力,同时比较了五种不同边界层参数化方案模拟得出的边界层热力和动力结构.结果表明：除个别方案外,配合不同边界层方案的WRF模式都能成功模拟出白天对流边界层强湍流混合特征和夜间稳定边界层内强逆温、逆湿和低空急流等热力和动力结构.非局地YSU、ACM2方案在白天表现出强的湍流混合和卷夹,相比于局地MYJ、UW方案,模拟的对流边界层温度更高、湿度更低、混合层高度更高、感热通量更大,更接近实际观测,这表明在不稳定层结下考虑非局地大涡输送更为合理,但局地方案在风速和风向的预报上存在一定优势.TEMF方案得到的白天局地湍流混合强度为所有方案中最弱,混合层难以发展,无法体现对流边界层内气象要素垂直分布均匀的特点.对于夜间稳定边界层的模拟,不同参数化方案之间的差异较小,但是YSU方案在一定程度上高估了机械湍流,导致局地湍流混合偏强,从而影响了其对稳定边界层的模拟能力.     

Turbulent entrainment at an inversion

Three separate conditions controlling the structure of the inversion are discussed. The rate of change of boundary layer temperature and the vertical gradient of potential temperature in the overlying air set the maximum limit to possible rise rates. In transient conditions, the rise of the inversion layer is sometimes dominated by turbulent transport and erosion, however this cannot continue indefinitely. The horizontal advection above the inversion can also change the relative temperature increase at the inversion, and so may partly determine inversion history.In this paper the theory of the plume-connective boundary layer has been extended to allow the plumes to erode material from the overlying stable air, and so permit study of how the boundary layer could deepen by eroding upward into the capping inversion. Only clear-air processes without condensation are included in the model described here.When the mechanism of convection is examined, it is demonstrated that, because of energetic considerations, the rate of erosion, by buoyancy generated turbulence of the overlying air down across the inversion, will not usually determine the rate of rise of the inversion, except briefly. Other effects often dominate in determining the average erosion rate. Furthermore, because of the limitations set by the buoyancy of the entrained air, sharp inversions, exceeding about 1°C, cannot be formed without the action of latent heat from evaporating droplets. An intermediate buffer layer forms in the inversion zone whenever the erosion rate is high enough to produce a temperature increase above the boundary layer, sufficient to inhibit direct entrainment.     

A three‐dimensional numerical investigation of the idealized planetary boundary layer

Abstract

The nonlinear equations of motion are integrated numerically in time for a region of x‐y‐z space of volume 3h × h × h, where h turns out to be a height slightly above the level where the wind first attains the geostrophic flow direction. Only the ideal case is treated of a horizontal lower boundary, neutral stability, horizontal homogeneity of all dependent mean variables except the mean pressure, and statistically steady state. The resulting flow patterns are turbulent and the eddies transport required amounts of momentum vertically.

Topics which are investigated include the relative directions of stress, wind shear and wind; differences in Ekman wind spirals for the neutral numerical case and a stable atmospheric case; profiles of dimensionless turbulence statistics; effect of allowing the mean density to be either constant or to decrease with height; effect of the wind direction or latitude upon the turbulence intensities; and characteristic structure of the eddies in the planetary boundary layer.     

Determination of critical shear stress of non‐cohesive soils using submerged jet test and turbulent kinetic energy

Laboratory tests using Jet Erosion Testing (JET) apparatus, impinging normally on a horizontal boundary, were conducted to determine the critical shear stress (τ_c) of non‐cohesive soil samples. A three‐dimensional (3D) SonTek/YSI 16 MHz Micro‐Acoustic Doppler Velocimeter (MicroADV) was used to measure turbulent kinetic energy (TKE) at a radial limit of entrainment in the wall jet zone and the measurements were used to calculate τ_c of the samples. The results showed that TKE increases exponentially with increasing particle size. The τ_c from this study were comparable (R² = 0.8) to the theoretical τ_c from Shields diagram after bed roughness scale ratio (D/k_s), due to the non‐uniform bed conditions, was accounted for. This study demonstrated that JET and TKE can be used to determine τ_c of non‐cohesive soils. The use of JET and TKE was found to be faster and easier when compared to the conventional approach of using flumes. A relationship of TKE at the onset of incipient motion (TKE_c) and samples’ D₅₀ developed in this study can be used to predict τ_c of non‐cohesive soils under similar non‐uniform conditions. Copyright © 2015 John Wiley & Sons, Ltd.     

Near‐bed shear stress,turbulence production and dissipation in a shallow and narrow tidal channel

Near‐bed, highly resolved velocity profiles were measured in the lower 0.03 m of the water column using acoustic Doppler profiling velocimeters in narrow tidal channels in a salt marsh. The bed shear stress was estimated from the velocity profiles using three methods: the log‐law, Reynolds stress, and shear stress derived from the turbulent kinetic energy (TKE). Bed shear stresses were largest during ebbing tide, while near‐bed velocities were larger during flooding tide. The Reynolds stress and TKE method gave similar results, while the log‐law method resulted in smaller bed shear stress values during ebbing tide. Shear stresses and turbulent kinetic energy followed a similar trend with the largest peaks during ebbing tide. The maximum turbulent kinetic energy was on the order of 1 × 10^{? 2} m²/s². The fluid shear stress during flooding tide was approximately 30% of the fluid shear stress during ebbing tide. The maximum TKE‐derived shear stress was 0.7 N/m² and 2.7 N/m² during flooding and ebbing tide, respectively, and occurred around 0.02 m above the bed. Turbulence dissipation was estimated using the frequency spectrum and structure function methods. Turbulence dissipation estimates from both methods were maximum near the bed (~0.01 m). Both the structure function and the frequency spectrum methods resulted in maximum dissipation estimates on the order of 4 × 10^{? 3} m²/s³. Turbulence production exceeded turbulence dissipation at every phase of the tide, suggesting that advection and vertical diffusion are not negligible. However, turbulence production and dissipation were within a factor of 2 for 77% of the estimates. The turbulence production and dissipation decreased quickly away from the bed, suggesting that measurements higher in the water column cannot be translated directly to turbulence production and dissipation estimates near the bed. Copyright © 2015 John Wiley & Sons, Ltd.     

Simulation of airflow around rain gauges: Comparison of LES with RANS models

Wind is responsible for systematic errors that affect rain gauge measurements. The authors investigate the use of computational fluid dynamics (CFD) to calculate airflow around rain gauges by applying a high-resolution large eddy simulation (LES) model to determine the flow fields around a measuring system of two rain gauges. The simulated air flow field is characterized by the presence of massive separation which induces the formation and shedding of highly unsteady eddies in the detached shear layers and wakes. Parts of these detached structures occur over the orifice of the rain gauges and may substantially affect the dynamics of the raindrops in this critical region. Non-dissipative LES methods used with fine enough meshes can successfully predict these eddies and their associated fluctuations. The authors compare statistics from LES with steady-state Reynolds averaged Navier–Stokes (RANS) simulations using the k–ε and shear stress transport k–ω turbulence models. They find that both RANS and LES models predict similar mean velocity distributions around the rain gauges. However, they determine the distribution of the resolved turbulent kinetic energy (TKE) to be strongly dependent on the RANS model used. Neither RANS model predictions of TKE are close to those of LES. The authors conclude that the failure of RANS to predict TKE is an important limitation, as TKE is needed to scale the local velocity fluctuations in stochastic models used to calculate the motion of raindrops in the flow field.     

On total turbulent energy and the passive and active role of buoyancy in turbulent momentum and mass transfer

Measurements of turbulent fluctuations of horizontal and vertical components of velocity, salinity and suspended particulate matter are presented. Turbulent Prandtl numbers are found to increase with stratification and to become larger than 1. Consequently, the vertical turbulent mass transport is suppressed by buoyancy forces, before the turbulent kinetic energy (TKE) and vertical turbulent momentum exchange are inhibited. With increasing stratification, the buoyancy fluxes do not cease, instead they become countergradient. We find that buoyantly driven motions play an active role in the transfer of mass. This is in agreement with trends derived from Monin–Obukhov scaling. For positive Richardson flux numbers (Ri _f ), the log velocity profile in the near-bed layer requires correction with a drag reduction. For negative Ri _f , the log velocity profile should be corrected with a drag increase, with increasing |Ri _f |. This highlights the active role played by buoyancy in momentum transfer and the production of TKE. However, the data do not appear to entirely follow Monin–Obukhov scaling. This is consistent with the notion that the turbulence field is not in equilibrium. The large stratification results in the decay of turbulence and countergradient buoyancy fluxes act to restore equilibrium in the energy budget. This implies that there is a finite adjustment timescale of the turbulence field to changes in velocity shear and density stratification. The energy transfers associated with the source and sink function of the buoyancy flux can be modeled with the concept of total turbulent energy.     

The influence of riparian vegetation on near‐bank turbulence: a flume experiment

Measurements from a fixed‐bed, Froude‐scaled hydraulic model of a stream in northeastern Vermont demonstrate the importance of forested riparian vegetation effects on near‐bank turbulence during overbank flows. Sections of the prototype stream, a tributary to Sleepers River, have increased in channel width within the last 40 years in response to passive reforestation of its riparian zone. Previous research found that reaches of small streams with forested riparian zones are commonly wider than adjacent reaches with non‐forested, or grassy, vegetation; however, driving mechanisms for this morphologic difference are not fully explained. Flume experiments were performed with a 1:5 scale, simplified model of half a channel and its floodplain, mimicking the typical non‐forested channel size. Two types of riparian vegetation were placed on the constructed floodplain: non‐forested, with synthetic grass carpeting; and forested, where rigid, randomly distributed, wooden dowels were added. Three‐dimensional velocities were measured with an acoustic Doppler velocimeter at 41 locations within the channel and floodplain at near‐bed and 0·6‐depth elevations. Observations of velocity components and calculations of turbulent kinetic energy (TKE), Reynolds shear stress and boundary shear stress showed significant differences between forested and non‐forested runs. Generally, forested runs exhibited a narrow band of high turbulence between the floodplain and main channel, where TKE was roughly two times greater than TKE in non‐forested runs. Compared to non‐forested runs, the hydraulic characteristics of forested runs appear to create an environment with higher erosion potential. Given that sediment entrainment and transport can be amplified in flows with high turbulence intensity and given that mature forested stream reaches are wider than comparable non‐forested reaches, our results demonstrated a possible driving mechanism for channel widening during overbank flow events in stream reaches with recently reforested riparian zones. Copyright © 2007 John Wiley & Sons, Ltd.     

白天混合层顶部夹卷层厚度的特征研究

本文首先对Deardorff的一阶模型给予解释，在此基础上分析对流边界层湍流动能方程，分析机械湍流和对流湍流对边界层发展的贡献，提出一个新的速度尺度，混合层顶速度尺度，定义了全理查森数，给出夹卷层厚度的参数化方案，并用Boers和Elotanta的雷达观测数据进行验证。参数化方案与实验数据符合得很好。当夹卷层厚度表示为夹卷速度或夹卷理查森数的函数时，该函数曲线随边界层发展通常表现为磁滞回线形状现象，利用本文的理论进行了解释。     

森林冠层和森林边界层大涡模拟

在采用各向异性湍流动能闭合方案和3阶Runge Kutta时间积分方案的大涡模式中,引入由森林冠层粗糙元造成的动量拖曳项、热量输入项和TKE耗散项,以模拟森林冠层和森林边界层的气象场. 通过中性和不稳定层结条件下不同叶面积指数算例的模拟分析及与已有观测结果的比较表明,本文所建大涡模式对森林冠层和森林边界层有较好的模拟效果. 进一步研究表明：不稳定层结条件下较稠密的森林冠层中特有的Kinking & Pairing湍涡结构与森林边界层中湍流的大涡运动相互作用,形成了森林冠层附近的温度斜坡型结构.