从湍流经典理论到大气湍流非平衡态热力学理论

湍流是日常生活中一种普遍的自然现象,也是经典物理学仍未完全解决的难题。湍流更是大气运动的最基本特征。本文系统地回顾了大气湍流经典理论发展简史,进一步详细介绍了大气湍流非平衡态热力学理论。大气湍流非平衡态热力学理论在熵平衡方程中引入动力过程,进而统一推导出大气湍流输送的Fourier定律、Flick定律和Newton定律,证明了Dufour效应、Soret效应、可逆动力过程与热力不可逆湍流输送过程之间的交叉耦合效应,以及湍流强度定理。这些定律和定理中得到了观测的事实验证,同时它们的唯象系数也由观测资料所确定。湍流强度定理揭示,湍流发展的宏观原因是速度和温度的剪切效应,Reynolds湍流和Rayleigh-Bénard湍流共存于大气湍流中。热力过程和动力过程间耦合效应现象的发现突破了传统湍流输送理论,即Fourier定律、Flick定律和Newton定律的观点——一个宏观量的输送通量等价于这个宏观量的梯度湍流输送通量。热力和动力过程间的耦合原理认为,一个宏观量的输送通量包括这个量的梯度湍流输送通量和速度耦合输送通量两部分。因此,能量和物质的垂直输送通量除了相应物理量梯度造成的湍流输送外,还应包括垂直速度耦合效应,即辐散或辐合运动造成的耦合效应。在一个很宽的尺度范围内,地表面的空间特征是非均匀的。下垫面非均匀性造成的对流运动将引起大气的辐散或辐合运动。这可能是导致地表能量收支不平衡的重要原因之一。垂直速度对垂直湍流输送的交叉耦合效应为非均匀下垫面大气边界层理论的发展,并为克服地表能量收支不平衡问题及非均匀下垫面大气边界层参数化遇到的困难提供了可能的线索。     

大气边界层湍流能量交换特征研究进展

从一维湍流能量平衡方程出发,回顾了近几十年湍流能量平衡方程中的各项以及Karman常数k的研究成果,总结了大气边界层湍流能量交换特征的研究概况和热点问题,并对今后发展趋势做了展望。实验研究表明,湍流能量平衡方程中的各项在不同条件下有不同的形式;传统的能量产生和耗散的局地平衡假设存在不足,特别是在不稳定条件下,垂直的湍流输运和压力脉动对湍流能量收支起了非常重要的作用。Karman常数与Rossby数和Reynolds数无关,在比较光滑的下垫面上k近似等于040±001。     

大气边界层气象学研究综述

文中回顾了大气边界层气象学的发展历史，总结了目前大气边界层气象学的主要进展，并指出国内外在未来大气边界层气象学研究方面面临的一些主要科学问题，以及对未来大气边界层气象学的发展方向提出若干建议，同时还指出了大气边界层气象学在思想上和方法上应该注意的一些相关问题。     

试探大气边界层湍流运动的若干问题

对湍流运动的几个基本理论问题作了一定程度的探讨，并根据湍流动能方程分析了大气边界层湍流产生和发展的能量机制，以此为判据推导了滴运动的凡主范围和物理含意较为清楚的混合长表达式。     

局地环流与大气边界层湍流相互作用的数值研究

本文建立了一个二维定常大气中尺度数值模式，并用该模式讨论了由下垫面粗糙分布非均匀（局地动力强迫）和温度分布非均匀（局地热力强迫）而产生的局地大气环流与大气边界层湍流的相互作用。结果指出：下垫面存在温度分布非均匀时运动方程中的湍流交换项与水平气压梯项一样可促使局地热力环境的形成，由粗糙度分布非均匀强迫产生的局地环流与由下垫面温度分布非均匀强迫产和的局地热力环流间的非线笥相互作用是通过湍流交换实现的。     

线性热力学对大气系统的应用(I)大气系统的线性唯象关系和热力学性质(英文)

从一般的热力学原理或其它自然原理对唯象关系所强加的限制,能够演绎出大气系统的一系列热力学性质。利用非平衡态线性热力学导出了湍流K闭合理论中湍流交换系数同唯象系数的关系,从理论上证明大气系统热量湍流输送同水汽之间存在交叉耦合,还导出了湍流强度同速度和位温梯度的关系,从而证明速度和位温空间分布的非均匀性是湍流之源。并证明湍流强度定理,不可压缩气体和各向同性湍流大气中,湍流强度正比于速度与位温梯度的标积。进而证明大气涡旋定理,位温的切变将导致涡旋运动或各种环流运动,速度涡度等于速度同位温相对梯度的矢积。展现了线性热力学在大气系统的应用前景。     

非均匀下垫面大气边界层研究进展

总结了水平非均匀下垫面的大气边界层研究,指出干旱区的绿洲具有明显的“冷湿岛效应”,并可形成边缘逆湿和区域环流,而植被和水域的分布则显著改变了大气边界层结构．阐述了复杂地形下垫面对局地风场和小尺度环流产生的影响、城市热岛效应的产生及特点和城市中绿地和水体的微气候效应．比较分析了不同陆面参数化方案以及边界层气候模式,提出未来在热力与动力非均匀的结合、湿地下垫面的数值模拟、城市复杂下垫面参数化以及气候模式的耦合等方面尚需进行深入的研究．     

非均匀网格的过渡湍流理论及其大气边界层数值模拟中的应用

本文将Ｓｕｌｌ提出的均匀网格下的过渡湍流理论推广到非均匀网格情形，推广的非均匀网格的过渡湍流理论满足Ｓｔｕｌｌ提出的对过渡矩阵系数的要求并具有清晰的物理意义。然后，将非均匀网络的过滤湍流理论应用于一维大气边界层数值模式中，对Ｗａｎｇａｒａ资料进行了模拟，并与均匀网格情形进行了对比。     

城市边界层湍流和下垫面空气动力学参数观测研究

城市大气边界层中的湍流交换和下垫面空气动力学参数的观测研究是城市环境气象研究的重要内容。文章概述了城市大气边界层中湍流的观测方法和计算城市下垫面空气动力学参数的可行性方案。对城市湍流的基本特性和城市边界层的垂直结构、城市边界层观测的历史、城市下垫面空气动力学参数的估算以及存在的问题进行了讨论。     

北京冬季大风过程大气边界层特征分析

利用北京中国科学院大气物理研究所325 m气象观测塔的气象梯度资料和湍流资料,分析了2014年11月29日至12月5日北京两次大风过程中气象要素和湍流输送特征的变化。第一次大风过程的强度和持续时间均高于第二次大风过程。强烈的风速垂直切变主要集中在距地面100 m高度范围内,最强风速垂直切变达到0.31 s~(-1)。大风过程中,阵风系数呈现随高度减小的趋势,越接近地面,阵风系数愈大。阵风强度的变化与阵风系数相似,100 m以下高度时,阵风强度随高度增大而减小。大风过程自上而下改变边界层结构,平均动能、湍流动能和摩擦速度最先从上层(280 m)发生变化且迅速增加。近地层由于风速垂直梯度的显著差异,近地层垂直方向的湍流强度最大。大风时各功率谱在低频区(0.01 s~(-1))达到峰值,大风过后各高度的能量都有所下降。     

天津城市边界层湍流统计特征

利用2008年2月1～28日在天津气象塔上观测的超声风、温资料和常规风、温及湿梯度资料,计算了天津城市边界层无量纲湍流速度方差、湍流温度方差、感热通量、动量通量和湍流动能。结果表明,在不稳定层结条件下,40m高度上无量纲湍流速度方差和湍流温度方差遵循莫宁-奥布克霍夫相似理论,在220m高度上只有垂直方向上的湍流速度方差...     

局地相似性理论在林地及林木湿地下垫面塔层中的应用研究

对局地相似性理论在“粗糙副层”范围的林地及在塔层范围的林木湿地下垫面大气边界层的适用性进行了研究。结果表明，对于二阶或更高阶的统计参数，局地相似性理论在林地及林木湿地是有效的。无量纲速度、垂直速度三阶矩、温度、湿度仅是无量纲长度ζ=z-d/L1的函数。其中，在对σT/T*l的分析中发现，σT/T*l按感热通量分层明显，感热通量大的数据点能更好地符合相似性关系。类似地，潜热通量大的数据点σT/T*l相似性关系也很好。     

北京城市化发展对大气边界层特性的影响

利用中国科学院大气物理研究所大气边界层物理与大气化学国家重点实验室的北京325 m气象塔1993年～2003年夏季 (7月～9月) 的观测资料, 统计分析了各年的风速与温度廓线分布特征。统计分析结果表明, 随着城市化的发展, 相对风速有逐年减小的趋势, 并且越靠近地面, 相对风速的减小越明显, 这反映了城市建筑对近地面层空气流动的摩擦作用。对风速廓线进行线性拟合得到风速随高度的垂直递增率, 发现无论是100 m以下的近地面层还是较高层, 风速的垂直递增率都随城市化发展存在逐年增大的趋势, 表明粗糙下垫面的影响已经向高层扩展。根据温度廓线计算了各年的温度垂直递减率, 发现其有增大的趋势, 这表明城市化发展对边界层热力结构同样有显著影响。本文还依据统计整理得到的近中性层结下的风速廓线资料, 利用莫宁-奥布霍夫相似理论计算了下垫面的空气动力学参数, 结果表明, 地表粗糙度、 零平面位移随着城市化发展皆有明显增加的趋势。同时, 分析了各空气动力学参数与平均风速及无量纲风速的关系。其中, 摩擦速度和平均风速二者基本成正相关, 且摩擦速度随平均风速的增大而增大的趋势越发明显。本文研究结果对研究城市化发展对区域大气边界层结构、 气候和环境影响有参考意义, 可为城市大气边界层模式和区域气候模式提供参数化依据。     

Wavelet Analysis of Intermittent Turbulent Transport in the Atmospheric Surface Layer over a Monsoon Trough Region

The structure of the turbulence in the atmospheric surface layer over a monsoon trough region has been studied using structural analysis based on wavelet transform. The observational site is located at the eastern (wet) end of the monsoon trough region, characterized by high moisture in the atmospheric surface layer. On the average relative humidity varied from 70% to 100% during the experiment. The wind and temperature data, collected at Kharagpur (22°25' N, 87°18' E) at six observational hours of a day in June 1990 during the Monsoon Trough Boundary Layer Experiment (MONTBLEX), have been utilized in the study. The wind and instantaneous momentum flux time series were decomposed into 12 scales using the Haar wavelet transform. The eddies exhibited a large temporal variability generating intermittency in the energy and flux distributions. A criterion based on the isotropy has been suggested for separating the large eddies from the small eddies. At the separation scale the isotropy coefficient drops sharply. It is shown that the intermittency in the small eddies resulted from the spatial variation of energy, and deviation of velocity statistics from the Gaussian distribution known as flatness. The deviation from the -5/3 power law has been attributed to the increased mean values of, (i) the coefficient of variation of energy, and (ii) the flatness factor, in the inertial subrange. The decomposition of the instantaneous momentum flux time series reveals that the major contribution to the total flux arises from the large eddies. The quadrant analysis of the momentum flux shows that ejections and sweeps account for a substantial part of the total flux, and quantifies the relative importance of the various spatial scales that contribute to the transport of momentum.     

On the Possible Impact of a Following-Swell on the Atmospheric Boundary Layer

A simple model of the atmospheric boundary layer over the ocean where the swell impact on the atmosphere is explicitly accounted for is suggested. The model is based on Ekman’s equations, where the stress in the wave boundary layer is split into two parts: the turbulent and wave-induced stress. The turbulent stress is parameterized traditionally via the eddy viscosity proportional to the generalized mixing length. The wave-induced stress directed upward (from swell to the atmosphere) is parameterized using the formalism of the wind-over-waves coupling theory. The model can be seen as an extension of the model by Kudryavtsev and Makin (J Phys Oceanogr 34:934–949, 2004) to the scale of the entire atmospheric boundary layer by including the Coriolis force into the momentum conservation equation and generalizing the definition of the mixing length. The regime of low winds for swell propagating along the wind direction is studied. It is shown that the impact of swell on the atmosphere is governed mainly by the swell parameter—the coupling parameter that is the product of the swell steepness and the growth rate coefficient. When the coupling parameter drops below − 1 the impact of swell becomes significant and affects the entire atmospheric boundary layer. The turbulent stress is enhanced near the surface as compared to the no-swell case, and becomes negative above the height of the inner region. The wind profile is characterized by a positive gradient near the surface and a negative gradient above the height of the inner region forming a characteristic bump at the height of the inner region. Results of the model agree at least qualitatively with observations performed in the atmosphere in presence of swell.     

Determination of the Atmospheric Boundary Layer Height from Radiosonde and Lidar Backscatter

The height of the atmospheric boundary layer is derived with the help of two different measuring systems and methods. From radiosoundings the boundary layer height is determined by the parcel method and by temperature and humidity gradients. From lidar backscatter measurements a combination of the averaging variance method and the high-resolution gradient method is used to determine boundary layer heights. In this paper lidar-derived boundary layer heights on a 10 min basis are presented. Datasets from four experiments – two over land and two over the sea – are used to compare boundary layer heights from both methods. Only the daytime boundary layer is investigated because the height of the nighttime stable boundary layer is below the range of the lidar. In many situations the boundary layer heights from both systems coincide within ±200 m. This corresponds to the standard deviation of lidar-derived 10-min values within a 1-h interval and is due to the time and space variability of the boundary layer height. Deviations appear for certain situations and depend on which radiosonde method is applied. The parcel method fails over land surfaces in the afternoon when the boundary layer stabilizes and over the ocean when the boundary layer is slightly stable. An automatic radiosonde gradient method sometimes fails when multiple layers are present, e.g. a residual layer above the growing convective boundary layer. The lidar method has the advantage of continuous tracing and thus avoids confusion with elevated layers. On the other hand, it mostly fails in situations with boundary layer clouds     

Large-Eddy Simulation of the Stable Atmospheric Boundary Layer using Dynamic Models with Different Averaging Schemes

Large-eddy simulation (LES) of a stable atmospheric boundary layer is performed using recently developed dynamic subgrid-scale (SGS) models. These models not only calculate the Smagorinsky coefficient and SGS Prandtl number dynamically based on the smallest resolved motions in the flow, they also allow for scale dependence of those coefficients. This dynamic calculation requires statistical averaging for numerical stability. Here, we evaluate three commonly used averaging schemes in stable atmospheric boundary-layer simulations: averaging over horizontal planes, over adjacent grid points, and following fluid particle trajectories. Particular attention is focused on assessing the effect of the different averaging methods on resolved flow statistics and SGS model coefficients. Our results indicate that averaging schemes that allow the coefficients to fluctuate locally give results that are in better agreement with boundary-layer similarity theory and previous LES studies. Even among models that are local, the averaging method is found to affect model coefficient probability density function distributions and turbulent spectra of the resolved velocity and temperature fields. Overall, averaging along fluid pathlines is found to produce the best combination of self consistent model coefficients, first- and second-order flow statistics and insensitivity to grid resolution.     

The Adaptation of the Atmospheric Boundary Layer to a Change in Surface Roughness

The adaptation of the atmospheric boundary layer to a change in the underlying surface roughness is an interesting problem and hence much research, theoretical, experimental, and numerical, has been undertaken. Within the atmospheric boundary layer an accurate numerical model for the turbulent properties of the atmospheric boundary layer needs to be implemented if physically realistic results are to be obtained. Here, the adaptation of the atmospheric boundary layer to a change in surface roughness is investigated using a first-order turbulence closure model, a one-and-a-half-order turbulence closure model and a second-order turbulence closure model. Perturbations to the geostrophic wind and the pressure gradients are included and it is shown that the second-order turbulence closure model, namely the standard k - model, is inferior to a lower-order closure model if a modification to limit the turbulent eddy size within the atmospheric boundary layer is not included within the model.