对流边界层顶部特性的对流槽实验模拟研究

利用对流槽研究对流边界层的顶部特性.实验结果表明,产生于混合层的上冲热泡可在对流边界层顶上的覆盖逆温层中激发出重力波;夹卷层的湍流结构表现出各向异性,水平尺度大于垂直尺度,与混合层中的湍流结构明显不同;利用可视图像直观显示了夹卷过程,夹卷层的温度谱表现出独特的结构特征,湍流能谱有明显的分区现象,谱幂率与现有理论分析结果有较大偏离.在实验结果的基础上,提出了关于夹卷速率的新参数化方案,夹卷速率可由地面热通量、混合层高度和覆盖逆温强度确定,方案中的系数C由实验测量数据拟合得出:C=0.95,与Deardroff等的对流槽实验数据的拟合结果(C=1.11)非常接近.     

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

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

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

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

云覆盖对流边界层顶部湍流结构参数的研究

应用飞机探测资料分析研究云覆盖对流边界层顶部温度和湿度湍流结构,在考虑对流边界层顶部夹卷过程的基础上得到计算温度和湿度结构参数的公式。应用实际观测资料分析了云覆盖对流边界层顶部的湍流特征．资料分析表明,云外晴空区温度和湿度结构函数值明显高于云内的值．云顶边界清晰,通过界面温度和湿度具有明显的跃变特征．应用观测资料检验了温度和湿度结构参数计算公式,计算结果与观测结果符合较好．     

夜间稳定边界层中小尺度地形激发的形式阻力和波动阻力

通过求解含有摩擦耗散的线性化大气动力学方程组,得到了在夜间稳定大气边界层中小尺度地形产生的波动阻力和形式阻力的解析解.结果表明边界层中的稳定度、风速和湍流状态、边界层厚度、上部残余层中的稳定度和风速以及地形高度和坡度,都会影响波动阻力和形式阻力的大小,应在数值模式的参数化方案中给予考虑.分析还表明,当地形坡度减到一定程度时,形式阻力可以忽略不计.     

大气边界层研究进展

大气边界层对云和对流的发展、演变有重要作用.本文回顾了在大气边界层高度计算方法,边界层的时空分布特征、结构和发展机理,以及边界层参数化方案等方面的主要研究进展.大气边界层高度计算方法主要分为基于大气廓线观测数据计算和基于模式参数化方案计算两大类;大气边界层高度频率分布形态具有明显的日变化特征,并且稳定、中性和对流边界层...     

边界层参数化对海南岛山地环流结构和湍流特征模拟的影响

本文基于WRF-ARW(V4.0)中尺度数值模式,选用耦合同一近地层参数化方案(Eta)的五种边界层参数化方案(MYJ、MYNN2、MYNN3、BouLac、UW),对2020年5月1—2日海南岛一次典型山地环流个例进行模拟,对比分析了这五种方案所模拟的山地环流结构和湍流特征的差异,旨在为模式应用于山地环流研究和模式改进提供一定的科学依据.研究结果表明,这五种湍流动能边界层参数化方案均能模拟出山谷风特征,对环流结构和湍流特征的影响表现为谷风时段大于山风时段.对于山地环流水平结构的模拟,因平原风推进距离的不同,五种方案模拟的近地面风速差异可达4 m·s^-1以上,其中MYJ方案模拟的谷风最强盛,而MYNN3方案最弱,山区多为静风或小风.对于山地环流垂直结构的模拟,MYNN2、UW方案模拟的谷风环流较强,表现为谷风厚度较厚,推进距离较远,同时由于模拟的谷风环流可越过山顶,从而模拟的高海拔地区上升区的覆盖范围和强度均较大;MYJ、BouLac方案模拟的谷风环流均未能越过山顶,且其中BouLac方案的平原风环流未能与上坡风环流耦合;而MYNN3方案模拟的环流结构最不明显.湍...     

北京城市大气边界层低层垂直动力和热力特征及其与污染物浓度关系的研究

利用2001年3月19~29日和2003年8月11~25日中国科学院大气物理研究所325m大气边界层观测塔资料,分析研讨了北京城市大气边界层低层的垂直动力结构特征及其与污染物浓度分布的关系.对比分析了北京城市大气边界层低层不同高度的风、温度和湿度梯度资料、大气湍流和大气化学观测系统资料,综合分析获取了无因次速度、温度湍流方差和湍流通量、湍能分布特征及其与污染物浓度空间分布的关系,同时分析了北京地区沙尘天气过程中城市边界层低层垂直结构特征及其污染物浓度的分布与变化特征.分析结果表明,在不稳定层结条件下,47和120m高度上无因次速度湍流方差和温度湍流方差遵循莫宁-奥布克霍夫相似规律,并给出相应的拟合公式.稳定大气边界层可按层结参数z′/L分成二分区,z′/L<0.1为弱稳定区,此时相似规律可适用,z′/L>0.1为强稳定区,在此区内无因次速度方差随稳定度增大有增大的趋势,而无因次温度方差则保持不变.白天近地层包含了47和120m,而280m已在近地层之上.对2001年3月北京地区一次沙尘天气过程的城市边界层资料分析发现,320m高度上总悬浮颗粒物浓度最高达到913.3μg/m3,在边界层内大气颗粒物从上层向低层输送,这与锋面过境时低空急流从上层向下发展过程并伴随的强下沉运动有关.     

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

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

黄土高原复杂地形上边界层低空急流对近地层湍流的影响

利用中尺度气象数值模式（Weather Research and Forecasting Model,WRF）模拟风场,结合兰州大学半干旱气候与环境观测站（Semi-Arid Climate and Environment Observatory of Lanzhou University,SACOL）湍流观测资料,分析了黄土高原复杂地形上稳定边界层低空急流对近地层湍流活动的影响.黄土高原复杂地形上稳定边界层低空急流的形成与地形作用引发的局地环流有关.低空急流对近地层湍流活动有强烈影响,剪切作用使小尺度湍涡活动加剧,湍动能增大,同时非平稳运动被压制.低空急流发生时,观测数据有87.3%是弱稳定情形（梯度理查森数小于0.25）;而无低空急流时,对应时段的观测表明65.4%属于强稳定层结（梯度理查森数大于0.3）,非平稳运动造成湍流功率谱在低频端迅速增大.与无低空急流和弱低空急流情形相比,强低空急流发生时,近地层湍动能增大1倍,湍动能在垂直方向上的传递增大1个量级,且方向向下,约为-3 × 10^-3 m³·s^-3,湍流在上层产生并向下传递.     

Effects of shear in the convective boundary layer: analysis of the turbulent kinetic energy budget

Effects of convective and mechanical turbulence at the entrainment zone are studied through the use of systematic Large-Eddy Simulation (LES) experiments. Five LES experiments with different shear characteristics in the quasi-steady barotropic boundary layer were conducted by increasing the value of the constant geostrophic wind by 5 m s^-1 until the geostrophic wind was equal to 20 m s^-1. The main result of this sensitivity analysis is that the convective boundary layer deepens with increasing wind speed due to the enhancement of the entrainment heat flux by the presence of shear. Regarding the evolution of the turbulence kinetic energy (TKE) budget for the studied cases, the following conclusions are drawn: (i) dissipation increases with shear, (ii) the transport and pressure terms decrease with increasing shear and can become a destruction term at the entrainment zone, and (iii) the time tendency of TKE remains small in all analyzed cases. Convective and local scaling arguments are applied to parameterize the TKE budget terms. Depending on the physical properties of each TKE budget contribution, two types of scaling parameters have been identified. For the processes influenced by mixed-layer properties, boundary layer depth and convective velocity have been used as scaling variables. On the contrary, if the physical processes are restricted to the entrainment zone, the inversion layer depth, the modulus of the horizontal velocity jump and the momentum fluxes at the inversion appear to be the natural choices for scaling these processes. A good fit of the TKE budget terms is obtained with the scaling, especially for shear contribution.     

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.     

Equivalent time-average and effective medium for periodic layers

The propagation of acoustic waves through a periodic layered medium is analyzed by an eigenvalue decomposition of the propagator matrix. This reveals how the velocity and attenuation of the layered medium vary as function of the periodic structure, material parameters and frequency. There are two important parameters which control the wave propagation in the periodic medium: the reflection coefficient and the ratio between one‐way traveltimes of the two parts of the cyclic layered medium. For low frequencies (large values of wavelength to layer thickness), the layered structure behaves as an effective medium, then there is a transition zone, and for higher frequencies (small values of wavelength to layer thickness) the medium is described by the time‐average velocity. In this paper we mostly concentrate on the transition zone between an effective medium and time‐average medium regimes. The width of the transition zone increases with larger values of the reflection coefficient. The transition zone corresponds to a blocking regime for which the transmission response of the layered structure is close to zero. For even higher frequencies, the time‐average medium is replaced by a new transition zone, and then again a time‐average medium. This pattern is periodically repeated with higher frequencies. For small values of the reflection coefficient, the transition between effective medium and time‐average medium occurs around a value of wavelength to layer thickness equal to 4.     

Bed-load transport rate based on the entrainment probabilities of sediment grains by rolling and lifting

A function for the bed-load sediment transport rate is derived. This function is obtained by using the entrainment probabilities of the rolling and lifted sediment grains, and by introducing two travel lengths, respectively. The predictions from the new bed-load function agree well with experimental results over the entire experimental range and show significant improvement over the commonly used formula for the bed-load transport rate. The new function shows that, in terms of contributing to the bed-load transport rate, the total entrainment probability of the sediment grains is a weighted summation of those for the lifted and rolling grains, rather than a simple addition of the two. The function is also used to predict the total entrainment probability, saltation length, and the bed layer thickness at a high bed-load transport rate. These predictions all agree well with the experimental results. It is found that, on average, the travel length for the rolling sand grains is about an order of magnitude less than that of the lifted grains.     

Influence of permeability on hydrothermal circulation in the sediment-buried oceanic crust

1 Introduction Hydrothermal circulation is the key process of hydrothermal activity. Modern seafloor hydrothermal circulation can be divided into three parts: convective cells in the oceanic curst, interface between seafloor and ocean and hydrothermal plume. Hydrothermal convection in the crust is the dominant part of the whole seafloor hydrothermal circulation. The distribu-tion and nature of hydrothermal system in the oceanic crust are controlled by crust thermal structures and permeability …     

OBSERVATIONS ON ESTUARINE FLUID MUD ENTRAINMENT

I. INTRODUCTIONPrediction of mud bed erosion by forcing due to tidal currents usually requires a numerical solution of the advection--dispersion equation for sediment mass transport. Key role is of course played inthis by the bottom boundary conditions defining erosion and deposition fluxes. The issue of erosion isbriefly considered here. noting that it is customary to calculate the rate of erosion as a function of thebed shear stress in excess of the erosion shear strength of the bed (Me…     

The structure of thermal convection and geotherms beneath a continent and ocean

The formation of the thermal cross section of the lithosphere and mantle upon the interaction between the mantle convection and the immobile continent surrounded by the oceanic lithosphere is studied by numerical modeling. The convective temperature and velocity fields and then the averaged geotherms for subcontinental and suboceanic regions up to the boundary with the core are calculated from the solution of convection equations with a jump in viscosity in the continental zone. Using the experimental data on the solidus temperature in the rocks of the upper mantle, the average thickness of the continental and oceanic lithosphere is estimated at 190 and 30 km, respectively. The effect of a hot spot formed in the subcontinental upper mantle at a depth of 250–500 km, which has not been previously noted, is revealed. Although the temperature in this zone is typically assumed to be close to adiabatic, the calculations show that it is actually higher than adiabatic by up to 200°C. The physical mechanism responsible for this effect is associated with the accumulation of convective heat beneath the thermally insulating layer of the continental lithosphere. The revealed anomalies can be important in studying the phase and mineral transformations at the base of the lithosphere and in the regional geodynamical reconstructions.     

Impact of soil moisture heterogeneity length scale and gradients on daytime coupled land‐cloudy boundary layer interactions

Using a coupled large‐eddy simulation–land surface model framework, the impact of two‐dimensional soil moisture heterogeneity on the cloudy boundary layer under varied free‐atmosphere stabilities is investigated. Specifically, the impacts of soil moisture heterogeneity length scale and heterogeneity in terms of soil moisture gradients on micrometeorological states, surface fluxes, boundary layer characteristics, and cloud development are examined. The results show that mesoscale circulations due to surface heterogeneity in soil moisture play an important role in transferring water vapour within the boundary layer and in regulating cloud distribution at the entrainment zone, which, in turn, provides feedbacks on boundary layer/surface energy budgets. The initial domain‐averaged soil moisture is identical for all homogenous and heterogeneous cases; however, the soil moisture heterogeneity in gradient and length scale between dry and wet regions has a significant impact on the estimates of near‐surface micrometeorological properties and surface fluxes, which further affect the boundary layer states and characteristics. Both liquid water potential temperature and liquid water mixing ratio increase with an increasing soil moisture gradient, whereas the amount of specific humidity decreases. Heterogeneity length scale and free atmosphere stability also amplify these impacts on the boundary layer structure and cloud formation. In a low atmospheric stability condition that potentially allows for a deeper boundary layer and a higher entrainment rate, cloud base height and cloud thickness significantly increase as the soil moisture gradient and length scale increase. Analysis to differentiate the influences of surface heterogeneity type (i.e. length scale vs gradient) shows that in general soil moisture gradient provides a larger impact than heterogeneity length scale, although the heterogeneity length scale is large enough to initiate circulation features responsible for differences in the coupled system between homogeneous and heterogeneous soil moisture cases. Copyright © 2012 John Wiley & Sons, Ltd.