地形影响的飞机颠簸及其数值仿真实验

首先利用中尺度数值模式ARPS模拟气流过山生成飞行数值仿真所需要的风场，然后利用飞机载荷因数变量方程进行飞机飞行的数值仿真试验。研究结果表明，气流过山产生的山脉重力波由于风切变临界层破碎，一方面能在对流层产生较强的湍流引起晴空飞机颠簸，另一方面也能在山脉背风面产生强烈下坡风，背风转子环流及低空湍流，影响飞机的起飞和着陆。揭示了飞机在过山脉地形背风面所产生的大气湍流中飞行时引起飞机颠簸的物理机制，有助于增强飞机颠簸的预测能力和飞行气象保障能力。     

中南地区一次飞机颠簸的物理量场分析

利用ERA-interim再分析资料和中南地区航空器空中报告,分析2017年2月4日20:00前后的颠簸报告位置与物理量场分布特征,结果表明:与颠簸位置对应最好的是散度场的强梯度区和温度平流的等零度线两侧,靠近涡度大值中心的强梯度区和垂直速度的大值区或大梯度区附近也容易发生颠簸,若满足以上条件的物理量场达3个或以上的区...     

Boundary-layer evolution and nocturnal inversion dispersal—Part II

Observations from an instrumented aircraft are used to study the small-scale structure of turbulence and convection in well-mixed boundary layers and the erosion processes in the nocturnally-formed inversions above them. The ways in which turbulence statistics for temperature, humidity and vertical velocity scale with height in the mixed layer are compared with the results of a three-dimensional model by Deardorff (1974a, b), and agreement is found in many aspects. Conditional sampling enables the statistics of thermals and their environment to be considered separately and, in particular, shows that the mode of the vertical velocity in thermals markedly decreases with height in the upper half of the mixed layer. Thermals may be recognized equally readily by either their excess of temperature or humidity. Transfers of heat and moisture through the nocturnal inversions influence the structure of the upper region of the mixed layer and there is strong evidence that these transfer processes are turbulent and not organized on scales similar to convective thermals.     

Analysis of aircraft measurements of momentum flux in the subcloud layer over the tropical atlantic ocean during gate

The momentum flux data obtained by the gust probe aboard the NOAA DC-6 aircraft during GATE are analyzed. Vertical profiles are obtained for Phases I and III and correlated with vertical wind velocity profiles using the geostrophic departure method. Reasonable agreement is obtained using the horizontal equations of motion with negligible advective acceleration. The vertical profiles of momentum flux and wind speed variance compare well with the numerical model results of Deardorff (1972) and Wyngaard et al. (1974). Vertical distributions of power spectra for vertical eddy motion and cospectra corresponding to the momentum flux components are obtained along with the height variation of the dominant length scales of vertical eddy motion and the dissipation rate of turbulence kinetic energy. When normalized by mixed-layer similarity, these results agree well with previous determinations in the boundary layer over tropical oceans and over land.     

Multi-beam techniques for deriving wind fields from airborne doppler radars

Summary Two techniques for deriving horizontal and vertical air motions using vertically scanning airborne Doppler radar data are presented and discussed. These techniques make use of the scanning ability of the NOAA P-3 tail-mounted radar antenna to view a region of space from at least two vantage points during a straight-line flight track. The scanning methodology is termed the Fore/Aft Scanning Technique or FAST because the antenna is alternately scanning forward and then aft of the flight track. The major advantages of FAST over flying two quasi-orthogonal flight tracks with the antenna scanning normal to the flight track are that the data are collected in roughly half the time and the aircraft does not have to execute a right-angle turn. However, accuracy of the resulting wind field is compromised slightly because the beam intersection angle is reduced from 90° to about 50°. The reduction of area covered because of large drift angles is also discussed.A three-dimensional wind field can be constructed using the dual-Doppler equations from FAST data using the two radial velocity estimates and vertical integration of the continuity equation with a boundary condition of no vertical motion at cloud top and the Earth's surface. To keep errors in the calculated winds acceptably small, the elevation angles are typically restricted to ±45° from the horizontal to minimize contamination of the horizontal wind by terminal fallspeeds.A different, and perhaps more believable vertical velocity, can be derived using a second technique that utilizes two (or more) airborne Doppler radar equipped aircraft each using FAST to observe the echo-top vertical velocity at common point (e.g., two aircraft flying parallel flight, paths, or by using an L-shaped flight track with a single aircraft). This technique results in 4 (or more) radial velocity estimates at each point (hence is called the quad-Doppler technique). Horizontal winds can be derived using either an overdetermined three-equation solution or an overdetermined dual-Doppler solution, whichever is more accurate. For the calculation of vertical velocity a new approach is proposed that utilizes the overdetermined triple-Doppler solution for vertical particle motion near cloud top, minus an estimate of terminal fallspeeds, as a top boundary condition for the downward vertical divergence integration to derive vertical air velocity elsewhere in the domain. In addition, this approach allows measurements at steep elevation angles allowing for more depth of coverage for a given range.To show the utility of the method, analyses of data collected using FAST are compared to conventional dual-Doppler-derived wind fields constructed from data collected simultaneously by S-band ground-based Doppler radars. An example of the quad-Doppler technique is also presented from the recently completed Tropical Oceans/Global Atmospheres Coupled Ocean/Atmosphere Response Experiment (TOGA/COARE). Comparisons of quad-Doppler vertical velocity are made with in-situ derived vertical air motions collected by the NASA DC-8 to judge the quality of the approach.With 9 Figures     

Cross‐shore variations of near‐surface wind velocity and atmospheric turbulence at the land‐sea boundary during CASP

Abstract

Airborne measurements of mean wind velocity and turbulence in the atmospheric boundary layer under wintertime conditions of cold offshore advection suggest that at a height of 50 m the mean wind speed increases with offshore distance by roughly 20% over a horizontal scale of order 10 km. Similarly, the vertical gust velocity and turbulent kinetic energy decay on scales of order 3.5 km by factors of 1.5 and 3.2, respectively. The scale of cross‐shore variations in the vertical fluxes of heat and downwind momentum is also 10 km, and the momentum flux is found to be roughly constant to 300 m, whereas the heat flux decreases with height. The stability parameter, z/L (where z = 50 m and L is the local Monin‐Obukhov length), is generally small over land but may reach order one over the warm ocean. The magnitude and horizontal length scales associated with the offshore variations in wind speed and turbulence are reasonably consistent with model results for a simple roughness change, but a more sophisticated model is required to interpret the combined effects of surface roughness and heat flux contrasts between land and sea.

Comparisons between aircraft and profile‐adjusted surface measurements of wind speed indicate that Doppler biases of 1–2 m s^?1 in the aircraft data caused by surface motions must be accounted for. In addition, the wind direction measurements of the Minimet anemometer buoy deployed in CASP are found to be in error by 25 ± 5°, possibly due to a misalignment of the anemometer vane. The vertical fluxes of heat and momentum show reasonably good agreement with surface estimates based on the Minimet data.     

一种飞机空中积冰潜势诊断指数与检验评估

飞机积冰是影响飞行安全的重要危险天气,准确的积冰诊断对民航运行与安全具有重大意义,特别是对国产大飞机适航取证具有重要作用。本文参考改进的CIP(Current Icing Potential)指数,使用垂直速度和云中液态水含量对积冰潜势诊断指数进行优化得到SCIP(Simplified Current Icing Potential)指数,该指数直接基于大气层结资料给出较为准确地飞机积冰潜势,具有命中率高、实用性强的特点。用2021年2月28日至12月31日民航航空器报告对原CIP指数和SCIP指数进行对比验证评估,发现相比于原CIP指数,SCIP指数具有更优异的ROC曲线(Receiver Operating Characteristic Curve)和TSS评分(True Skill Score)。但是对于不同程度的积冰,原CIP指数和SCIP指数均表现出较一致诊断效果,表明SCIP指数对诊断飞机积冰潜势具有重要作用,但无法判断积冰强度。根据2020年3月16日陕西区域飞机结冰个例验证表明增加垂直速度和云中液态水含量的影响对飞机积冰的描述更加细致,对提升积冰潜势诊断效果具有重要意义...     

Mechanisms Controlling Turbulence Development Across A Forest Edge

In this paper we discuss the development of turbulence back from the transition fromopen moorland to a forest. Data from a field study and a wind-tunnel experiment arepresented. These show that the variance in the streamwise velocity begins to adjust tothe new surface between 2 to 4 tree heights downwind of the transition. This is soonerthan either the vertical velocity variance or the shear stress, both of which begin to adjust in a zone 3 to 5 tree heights downwind of the edge. Key terms in the prognostic equations for streamwise and vertical velocity variance are evaluated in order to explain these differences. The flow distortion caused by the forest edge, which extends to 4 tree heights downwind of the forest edge, is shown to be crucial in the delayed turbulence development. Initially the shear production term, which is the dominant source for the streamwise velocity variance, is counteracted by a sink in the vertical advection term. After the flow levels out the pressure redistribution (return-to-isotropy) term becomes the main sink of streamwisevelocity variance and feeds energy into the vertical velocity component. Therefore, thedevelopment of the vertical velocity variance and shear stress cannot begin until afterdevelopment of an increase in the streamwise velocity variance. Results are comparedwith other experiments, including the flow across shelterbelts, and large-eddy simulations of forest flow.     

Observations of planview flux patterns within convective structures of the marine atmospheric surface layer

Air/sea flux variability on horizontal scales from 50 m to several km results, in part, from the presence of coherent convective structures within the atmospheric boundary layer. The horizontal distribution of fluxes within these convective updrafts and downdrafts is, therefore, central to studies of air/sea interaction and remote sensing of sea surface wind and wave fields. This study derives these flux patterns from observations of the Marine Atmospheric Surface Layer (MASL).Research aircraft flights through the MASL provide an optimal means for sampling large numbers of the above-mentioned coherent structures. The NCAR Electra flew numerous legs through the MASL at a height of 50 m during the 1987 stratocumulus phase of Project FIRE (First ISSCP (International Satellite Cloud Climatology Program) Regional Experiment).In situ measurements from these legs serve as the dataset for this paper. The data are processed in such a way as to retain only the turbulence fluctuations. Conditional sampling, based on the vertical velocity field, results in the isolation of convective updrafts and downdrafts. Compositing of the data for these two classes of convective drafts results in horizontal planviews of the vertical fluxes of buoyancy, absolute humidity, along-meanwind component of momentum, and vertical velocity. To ensure dynamical similarity, these horizontal planviews are oriented in a coordinate system aligned with the mean wind.     

Etude de la couche limite convective sahelienne en presence de brumes seches (Expérience ECLATS)

The ECLATS experiment was conducted in order to investigate the influence of radiative processes on the dynamics of the atmospheric boundary layer during its diurnal evolution. This experiment was carried out over Niger, near Niamey, by measuring continuously the energy balance at ground level and by using an instrumented aircraft for turbulence, radiative fluxes and aerosol measurements in the boundary layer during dusty conditions (brumes sèches). This paper is restricted to an analysis of the turbulent structure in the homogeneous and stationary convective boundary layer. The turbulence moments for kinetic energy and the spectral characteristics of the vertical velocity are discussed. These results are compared with a set of data obtained for clear convective boundary layers. The differences observed are quite important and seem, at least in part, due to radiative processes (infrared radiative divergence in the surface layer and absorption of solar radiation in the boundary layer).     

风廓线雷达对塔克拉玛干沙漠晴天边界层的探测分析

利用塔克拉玛干沙漠大气环境观测试验站2010年8和10月边界层风廓线雷达资料,分析了沙漠晴天边界层湍流、大气温度、水平风速风向以及垂直速度的发展演变特征和日变化规律。研究表明：（1）大气折射率结构常数（C_n~2）能较好地反映晴空湍流对电磁波的后向散射能力,可以详细刻画湍流发展旺盛区域的高度、强度及其演变特征;沙漠夏季白天湍流发展剧烈,旺盛区域顶部可达4000 m高度左右。（2）RASS系统对沙漠边界层大气温度的探测具有较好的可信度,其近地边界层温度符合一般的日变化规律,昼夜温差显著,白天高温维持时间长,升温过程相对滞后于近地表气温。（3）风廓线雷达对大气风场的探测结果与地面风速风向一致,沙漠晴天主要受东风和东北风控制,风速较小,平均在2.0～6.0 m·s~（-1）范围变化。（4）沙漠腹地大气垂直速度变化符合静力平衡理论,铅直方向运动很弱,一般在-1.0～1.0 m·s~（-1）范围波动。     

Wave-Modified Flux and Plume Dispersion in the Stable Boundary Layer

The effects of a pressure jump and a following internal gravity wave on turbulence and plume diffusion in the stable planetary boundary layer are examined. The pressure jump was accompanied by a sudden increase in turbulence and plume dispersion. The effects of wave perturbations on turbulence statistics are analysed by calculating fluxes and variances with and without the wave signal for averaging times ranging from 1 to 30 min. The wave signals are obtained using a band-pass filter. It is shown that second-order turbulence quantities calculated without first subtracting the wave perturbations from the time are greater than those calculated when the wave signal is separated from the turbulence. Estimates of the vertical dispersion of an elevated tracer plume in the stable boundary layer are made using an elastic backscatter lidar. Plume dispersion observed 25 m downwind of the source increases rapidly with the arrival of the flow disturbances. Measured plume dispersion and plume centreline height correlate with the standard deviation of the vertical velocity but not with the wave signal.     

A Homogeneous Langevin Equation Model, Part i: Simulation of Particle Trajectories in Turbulence with a Skewed Velocity Distribution

A Lagrangian stochastic model for the time evolution of the velocity of a fluid particle is presented. This model is based on a one-dimensional generalized Langevin equation, and assumes the velocity probability distribution of the turbulent fluid is skewed and spatially homogeneous. This has been shown to be an effective approach to simulating vertical dispersion in the convective boundary layer. We use a form of the Langevin equation that has a linear (in velocity) deterministic acceleration and a random acceleration that is a non-Gaussian, skewed process. For the case of homogeneous fluid velocity statistics, this 'linear-skewed' Langevin equation can be integrated explicitly, resulting in an efficient numerical simulation method. Model simulations were tested using cases for which exact, analytic statistical properties of particle velocity are known. Results of these tests show that, for homogeneous turbulence, a linear-skewed Langevin equation model can overcome the difficulties encountered in applying a Langevin equation with a skewed random acceleration. The linear-skewed Langevin equation model results are compared to results of a 'nonlinear-Gaussian' Langevin equation model, and show that the linear-skewed model is significantly more efficient.     

Measurements and Langevin simulations of mean tracer concentration fields downwind from a circular line source inside an alfalfa canopy

Two Langevin simulations of trajectories of marked fluid elements in inhomogenous turbulence, where the Lagrangian length and vertical velocity scales are height dependent, were compared with field data. A CO₂ tracer was released from a circular line source and the concentration profiles were measured for diffusion distances of 50 and 100 cm inside and above an alfalfa canopy.One of the simulations, suggested by Wilson et al. (1983), biases the vertical velocities by adding a mean upward drift. The second simulation proposed here by-passes this difficulty by reflecting marked particles according to a probability calculated from the gradient in vertical velocity variance between the beginning and the end of each step. This simulation also makes use of a constant time-scale within the canopy, following preliminary results from a turbulence experiment within a forest (Leclerc, 1987).Comparing the results of these simulations with the field data shows that the simulation proposed by Wilson et al. (1983) does not correctly reproduce the difusion for the larger fetch in systems exhibiting strong gradients in vertical velocity variance. Instead, the modelled plumes exhibit a bulge at the source height whereas the field data show smooth profiles. In addition, the modelled plumes overestimate the vertical spread of the plumes, which is possibly due to the inadequacy of the approach in severely inhomogeneous systems. In contrast, the results from the tracer experiments indicate that the diffusion can be better reproduced with the use of a reflection probability calculated at each step. The discrepancies between the experimental results and the simulation using a reflection probability are attributed to stability effects.     

A Homogeneous Langevin Equation Model, Part ii: Simulation of Dispersion in the Convective Boundary Layer

We present a Lagrangian stochastic model of vertical dispersion in the convective boundary layer (CBL). This model is based on a generalized Langevin equation that uses the simplifying assumption that the skewed vertical velocity probability distribution is spatially homogeneous. This approach has been shown to account for two key properties of CBL turbulence associated with large-scale coherent turbulent structures: skewed vertical velocity distributions and long velocity correlation time. A 'linear-skewed' form of the generalized Langevin equation is used, which has a linear (in velocity) deterministic acceleration and a skewed random acceleration. 'Reflection' boundary conditions for selecting a new velocity for a particle that encounters a boundary were investigated, including alternatives to the standard assumption that the magnitudes of the particle incident and reflected velocities are positively correlated. Model simulations were tested using cases for which exact, analytic statistical properties of particle velocity and position are known, i.e., well-mixed spatial and velocity distributions. Simulations of laboratory experiments of CBL dispersion show that (1) the homogeneous linear-skewed Langevin equation model (as well as an alternative 'nonlinear-Gaussian' Langevin equation model) can simulate the important aspects of dispersion in the CBL, and (2) a negatively-correlated-speed reflection boundary condition simulates the observed dispersion of material near the surface in the CBL significantly better than alternative reflection boundary conditions. The homogeneous linear-skewed Langevin equation model has the advantage that it is computationally more efficient than the homogeneous nonlinear-Gaussian Langevin equation model, and considerably more efficient than inhomogeneous Langevin equation models.     

黄土高原复杂地形受中尺度运动影响的稳定边界层湍流特征

利用兰州大学半干旱气候与环境观测站(Semi-Arid Climate and Environment Observatory of Lanzhou University,简称SACOL)2008年12月观测资料,研究了稳定边界层湍流特征.使用涡动相关资料研究湍流通量时,定义湍流的平均时间τ内的中尺度运动是造成湍流统计量变化范围大的主要原因,稳定情形? τ取几十秒至几分钟.对梯度理查森数大于0.3的强稳定情形的湍流尺度分解(MRD)谱分析表明,感热通量在112.4～449.9 s存在谱隙,尺度大于谱隙的中尺度运动造成了通量观测资料离散性大,甚至有支配性影响.动量通量的谱隙在112.4～224.9 s之间.弱风时,中尺度运动的影响更大,垂直风速标准差以0.1的比率随中尺度风速变化;垂直风速标准差同广义风速表现出很好的相关性,并随着广义风速消失而消失.三维风速标准差与摩擦速度呈很好的线性关系,垂直、水平、横风风速的无量纲标准差分别为1.35、2.54、2.21.对湍流动能的研究发现,在梯度理查森数大于0.3的条件下,仍然存在连续的湍流.以湍动能为依据,分析了湍流的平稳时间长度,其长度随稳定度变化而变化,2008年12月7～11日从133.5 s变化到856.2 s,湍流平稳时间长度反映了中尺度运动的发生频率.     

Increased Light,Moderate,and Severe Clear-Air Turbulence in Response to Climate Change

Anthropogenic climate change is expected to strengthen the vertical wind shears at aircraft cruising altitudes within the atmospheric jet streams. Such a strengthening would increase the prevalence of the shear instabilities that generate clear-air turbulence. Climate modelling studies have indicated that the amount of moderate-or-greater clear-air turbulence on transatlantic flight routes in winter will increase significantly in future as the climate changes. However,the individual responses of light,moderate,and severe clear-air turbulence have not previously been studied,despite their importance for aircraft operations. Here,we use climate model simulations to analyse the transatlantic wintertime clear-air turbulence response to climate change in five aviation-relevant turbulence strength categories. We find that the probability distributions for an ensemble of 21 clear-air turbulence diagnostics generally gain probability in their right-hand tails when the atmospheric carbon dioxide concentration is doubled. By converting the diagnostics into eddy dissipation rates,we find that the ensembleaverage airspace volume containing light clear-air turbulence increases by 59%(with an intra-ensemble range of 43%–68%),light-to-moderate by 75%(39%–96%),moderate by 94%(37%–118%),moderate-to-severe by 127%(30%–170%),and severe by 149%(36%–188%). These results suggest that the prevalence of transatlantic wintertime clear-air turbulence will increase significantly in all aviation-relevant strength categories as the climate changes.     

Large-eddy simulation of turbulent flow above and within a forest

A large-eddy simulation has been performed of an atmospheric surface layer in which the lower third of the domain is occupied by a drag layer and heat sources to represent a forest. Subgridscale processes are treated using second-order closure techniques. Lateral boundaries are periodic, while the upper boundary is a frictionless fixed lid. Mean vertical profiles of wind velocity derived from the output are realistic in their shape and response to forest density. Similarly, vertical profiles of Reynolds stress, turbulent kinetic energy and velocity skewness match observations, at least in a qualitative sense. The limited vertical extent of the domain and the artificial upper boundary, however, cause some departures from measured turbulence profiles in real forests. Instantaneous turbulent velocity and scalar fields are presented which show some of the features obtained by tower instrumentation in the field and in wind tunnels, such as the vertical coherence of vertical velocity and the slope of structures revealed by temperature patterns.     

高风速相干结构对通量输送影响的实验研究

切变湍流的相干结构是湍流研究中的重大发现,它表明湍流运动并非完全随机,其中具有可检测的有序结构.本文通过处理南京浦口地区大气边界层观测数据,来分析不稳定层结中高风速相干结构特征.本次观测项目包括对场地中央的气象铁塔上2 m和40 m高度上超声风速仪的脉动速度、温度测量以及风廓线雷达对边界层风速廓线的测量.对超声水平风速时间序列数据进行小波变换 (时间尺度400 s),通过阈值来识别这种高风速相干结构.与多普勒风廓线雷达测量结果对比后发现,这种方法确定的相干结构符合常规的认识,具有较长的时间尺度和较大的垂直尺度 (接近边界层厚度).分析三天相干结构特性得到无量纲空间间隔约为6,即每隔6个边界层厚度的水平位置出现一个高速相干结构.通过与垂直风速小波系数的比较,发现高风速相干结构与向下垂直风速之间有较好相关,这与湍流中 “阵风” 现象的研究结论相似.使用四象限分析方法分类得到两种动量通量输送为负的运动:较小水平风速的上扬 (ejection) 运动 (简称为上扬运动) 和较大水平风速的下扫 (sweep) 运动 (简称为下扫运动),这两种运动在整个湍流活动中处于主导地位.高风速相干结构通过促进下扫运动和抑制上扬运动来影响动量通量的输送.