On drag coefficient parameterization with post processed direct fluxes measurements over the ocean

This study presents an evaluation of the atmospheric factors influencing the post-processing for fast-response data of horizontal momentum, vertical wind component, temperature, and water vapor to measure turbulent fluxes. They are observed at the Ieodo ocean research station over the Yellow Sea during the period of October 2004 to February 2008. The post process methods employed here are composed of quality control and tilt correction for turbulent flux measurement. The present result of quality control on the fast-response data shows that total removal ratio of the data generally depends on the factors such as a wind speed, relative humidity, significant wave height, visibility, and stability parameter (z/L). Especially, the removal ratio of water vapor data is significantly increased on light wind and strong stability conditions. The results show that the total removal ratio of water vapor data increases when wind speed is less than 3 m s^?1 and wave height is less than 1 m. The total removal ratio of water vapor data also increases with the value of the stability parameter. Three different algorithms of tilt correction methods (double rotation, triple rotation, and planar fit) are applied to correct the tilt of the sonic anemometer used in the observation. Friction velocities in near neutral state are greater than friction velocity in other states. Drag coefficients are categorized in terms of stabilities and seasons.     

The geostrophic drag coefficient

A power-law formulation for the geostrophic drag coefficient in a neutral atmosphere is presented.Deceased, 27 December, 1973.     

On the drag coefficient over shallow water

It is shown that if the wind-wave spectrum in shallow water is approximately independent of wind speed due to the combined effects of white-capping and bottom friction, then the wave-induced drag coefficient has a maximum value when the wind speed is twice the maximum wave-speed; as the wind speed increases further, the drag coefficient slowly decreases. This result is consistent with the observations of Hicks et al. (1974).     

On the parameterization of drag over small-scale topography in neutrally-stratified boundary-layer flow

Predictions of the surface drag in turbulent boundary-layer flow over two-dimensional sinusoidal topography from various numerical models are compared. For simple 2D terrain, the model results show that the drag increases associated with topography are essentially proportional to (slope)² up to the steepness at which the flow separates. For the purposes of boundary-layer parameterisation within larger-scale models, we propose a representation of the effects of simple 2D topography via an effective roughness length, z ₀ ^eff. The form of the varation of z ₀ ^eff with terrain slope and topographic wavelength is established for small slopes from the model results and a semi-empirical formula is proposed.     

Momentum budget diagnosis and the parameterization of subgrid-scale orographic drag in global GRAPES

The initial tendency approach is used to diagnose systematic errors in global GRAPES (Global/Regional Assimilation Prediction System), including overly strong westerlies in the northern midlatitudes, cold/warm bias dipoles in the vicinity of the tropopause, and excessively strong southerlies in downstream regions of the Tibetan Plateau. This approach, involving the use of the assimilation system, focuses on the first few time steps of numerical weather forecasts to identify the deficiencies in diabatic forcing. The results show that there is insufficient diabatic dissipation in the upper troposphere and lower stratosphere of the northern midlatitudes and the lower troposphere of most latitudes, which results from the absence of a parameterization of subgrid orographic drag in global GRAPES. A scheme to parameterize the effects of these drags is therefore tested and the experiments indicate that the newly introduced scheme reduces zonal momentum budget residuals, weakens the northern midlatitude westerlies and southerlies in the downstream regions of the Tibetan Plateau, decreases the cold/warm bias dipoles, and leads to improved objective verification scores.     

GRAPES中地形重力波拖曳物理过程的引进和应用试验

在中国新一代全球中期/区域中尺度同化与预报系统(GRAPES)模式中引进了ECMWF地形重力波拖曳物理过程,填补了GRAPES全球中期数值预报系统中物理过程的空白。重新计算了地形重力波过程需要的地形静态资料数据,并与原ECMWF模式的地形静态参数进行了对比分析,验证了模式地形参数的正确性。利用GRAPES模式,进行了地形重力波拖曳物理过程影响的敏感性数值试验;结果表明:引进地形重力波拖曳过程以后,在存在大地形的区域,风场会发生变化,当纬向风遇到青藏高原时,一部分气流会产生爬坡效应而越过高原,使高原上空的西风气流减弱;另一部分气流会绕过高原,在高原的南侧产生绕流;随着模式积分时间的延长,风场变化会越来越明显,地形越复杂,风场的变化也越复杂;连续的模式积分试验结果显示,引进地形重力波过程,可以延长GRAPES模式的可用预报时效,提高了全球形势预报的准确率。通过对一次降水过程的模拟,对地形重力波过程影响降水预报的原因进行了简单分析。结果显示:引进地形重力波拖曳过程后,改变了大气流场的分布,使预报的流场更接近于大气真实状态,从而提高了降水预报的准确率。     

Test of a convectively forced gravity wave drag parameterization in a general circulation model

The influence of gravity wave drag induced by cumulus convection (GWDC) on a simulated boreal summer climate was evaluated in a general circulation model. For this, the GWDC scheme developed by Chun and Baik was implemented into a version of the National Centers for Environmental Prediction (NCEP) global spectral model (GSM). Ensemble simulations with the two different convection schemes, the simplified Arakawa-Schubert (SAS) scheme and Community Climate Model (CCM) convection scheme, were conducted for the boreal summer of 1996. A cloud factor to modulate the stress intensity with respect to the cloud type was introduced in this study, in order to prevent unrealistic behaviors of the GWDC scheme in GSM. The effect of gravity wave drag on the zonal mean of wind and temperature fields was focused. On the whole, the effect of GWDC in this study is positive on the simulated seasonal climate. It is evident that biases in temperature in the polar region as well as in the zonal and meridional winds in the upper atmosphere are reduced. The percentage of reduction of the bias in zonal winds is about 10–20%. Such a response of the GWDC forcing widely appears not only in tropical regions but also in mid-latitude regions. These characteristics are prominent in the case of the SAS scheme, which is due to the various convective cloud types. The magnitude of GWDC forcing is generally small, but still positive, in the case of the CCM scheme, which is due to rather homogeneous cloud types. It is also found that the role of a particular GWDC forcing depends upon the inherent systematic biases of a particular model. It is concluded that incorporation of the GWDC parameterization in GCMs should be taken into account to improve the seasonal prediction.     

Aircraft estimates of the geostrophic drag coefficient and the rossby similarity functions A and B over the sea

Estimates of the geostrophic drag coefficient and the Rossby similarity functions, A and B obtained from data collected by an instrumented aircraft over the sea are presented. The average value of the geostrophic drag coefficient is 0.027 and is independent of the geostrophic windspeed. The dependence of the similarity functions A and B on boundary-layer parameters is investigated. The function A is found to depend on baroclinicity parameters, while B depends on the parameter u _*/fh (where u _* is the surface friction velocity, f is the Coriolis parameter, and h is the boundary-layer depth). Using the geostrophic drag coefficient found here and the results of surface drag coefficient studies, a relationship between geostrophic windspeed and surface windspeed is obtained which shows good agreement with empirical data.     

Impact of terrain complexity on the turbulence drag coefficient: A case study from the Indian Himalayan region

The aerodynamic drag coefficient (C_D) is conjectured to change (or remains almost uniform) with the horizontal wind speed (U) over a flexible (or fixed) surface element, which is represented with the surface roughness (z₀). This conjecture is tested for the near neutral atmospheric turbulence (i.e. when surface stability z/L is almost equal to 0, where z is the measurement height and L is Obukhov length) of monsoon and winter season at an on-slope and a ridge-top site in the Indian Himalaya, wherein the ridge-top site is associated with a higher degree of sensitivity to the roughness element and terrain attributes. This hypothesis is successfully verified for two conditions, (i) the monsoon period observations of ridge-top site are found to have higher z₀ due to vegetative growth than the winter period for flows having similar terrain signature, and (ii) the monsoon and winter period observations of on-slope site are noted to have similar z₀ for flows having signature of steep terrain. Subsequently, constants (i.e. a and b) of the power-law relationships between C_D and U (i.e. C_D = aU^b), as a function of z₀, are optimized. It is noted that the relationship between C_D and U has higher sensitivity towards the terrain slope than the vegetative growth.     

Matching classical boundary-layer solutions toward a geostrophic drag coefficient relation

The relationship between geostrophic, Ekman-layer, and surface-layer flow is explored. Formalized limit solutions for each layer are developed; together they form a continuous solution for the semi-infinite flow over a surface. These classical solutions are re-derived and related, to prepare for a two-layer composite boundarylayer solution with patching criteria. The patching process yields a resistance law which indicates a correspondence between the stress and the geostrophic flow deviation involving a similarity relation with two arbitrary parameters.     

Study on diurnal variation of bulk drag coefficient over different landsurfaces

Summary The magnitude and diurnal change of turbulent bulk drag coefficients over land have been analysed using mean velocity and temperature gradient data of the planetary boundary layer. The turbulent drag coefficients can be about an order of magnitude larger over the rough land surface than over the sea surface. We computed these coefficients by the same method for three typical underlying surfaces represented by urban, grassland and Gobi desert. The results show that there are significant differences in the turbulent transfer among the three typical underlying surfaces.With 6 Figures     

A parameterization of the stable atmospheric boundary layer

Two formulations of the stable atmospheric boundary layer are proposed for use in weather forecasting or climate models. They feature the log-linear profile near the surface, but are free from the associated critical Richardson number. The diffusion coefficients in the Ekman layer are a natural extension of the surface layer. They are locally determined using wind shear in one case and turbulent kinetic energy in the other. The parameterizations are tested in a one-dimensional model simulating the evolution of the nocturnal boundary layer with and without radiative cooling. Both formulations give very similar results, except near the top of the boundary layer where the transition to the free atmosphere is smoother with the wind shear formulation. A distinctive feature of these schemes is that they retain their simulating skill when resolution is reduced. This is verified for a wide range of situations. In practice, this means that there is no need for a large-scale model to have a level below 50 m or so.     

A subgrid parameterization of orographic precipitation

Summary Estimates of the impact of global climate change on land surface hydrology require climate information on spatial scales far smaller than those explicitly resolved by global climate models of today and the foreseeable future. To bridge the gap between what is required and what is resolved, we propose a subgrid-scale parameterization of the influence of topography on clouds, precipitation, and land surface hydrology. The parameterization represents subgrid variations in surface elevation in terms of probability distributions of discrete elevation classes. Separate cloud, radiative, and surface processes are calculated for each elevation class. Rainshadow effects are not treated by the parameterization; they have to be explicitly resolved by the host model. The simulated surface temperature, precipitation, and snow cover for each elevation class are distributed to different geographical locations according to the spatial distribution of surface elevation within each grid cell.The subgrid parameterization has been implemented in the Pacific Northwest Laboratory's climate version of the Penn State/NCAR Mesoscale Model. The scheme is evaluated by driving the regional climate model with observed lateral boundary conditions for the Pacific Northwest and comparing simulated fields with surface observations. The method yields more realistic spatial distributions of precipitation and snow cover in mountainous areas and is considerably more computationally efficient than achieving high resolution by the use of nesting in the regional climate model.With 17 Figures     

The variation of the GPS dropwindsonde drag coefficient and its impact on the wind retrieval from dropwindsonde measurements

The wind finding equations currently used to retrieve horizontal winds from Global Positioning System dropwindsonde wind measurements are derived based on a point object model in which the drag coefficient of the dropwindsonde is assumed to be a constant. The wind tunnel tests performed as part of this study showed, however, that the dropwindsonde aerodynamic coefficients vary appreciably with angles of attack. To investigate the impact of this finding, the dropwindsonde motion in a pseudo-stochastic wind field has been simulated using a motion model more sophisticated than the point object model. The results showed that, although the constant drag coefficient assumption is not supported by the wind tunnel test results, the wind finding equations still correctly calculate both the mean and the turbulence intensity profiles. In addition, a revised method to calculate the vertical wind was proposed based on the derivation of the improved motion model, which enhanced the accuracy of vertical wind estimates by including the real-time dropwindsonde drag coefficient and the dropwindsonde vertical acceleration into calculation.     

A parameterization of the boundary layer of a tropical cyclone

In an earlier paper by one of the authors (Smith, 1968), a momentum integral method was developed to parameterize the gross constraint imposed by the surface boundary layer of a steady, axisymmetric, tropical cyclone on the meridional circulation within the vortex itself. Specifically, the method provides an effective means of estimating the radial variation of mean upflow/downflow induced by the boundary layer, compatible with a prescribed radial variation of azimuthal velocity just above the boundary layer,V _gr. However, it relies on a judicious choice of vertical profiles of radial and azimuthal velocity components within the boundary layer. An especially suitable set of profiles is discussed herein; these are Ekman-like profiles in which turbulent mixing is characterized by a vertically constant eddy diffusivityK _M , matched to a constant stress sublayer just above the sea surface. An attractive feature of the formulation is that a suitable value forK _M as a function of radius, which is extremely difficult to extract from observational data, can be calculated when the state of the sea surface, described by a roughness lengthZ ₀, is prescribed. Although observations ofZ ₀ at high wind speeds are not yet available, the effect of radial variations in sea surface roughness can be assessed and it is shown that these affect the upflow to a significant degree.     

Assessment of surface drag coefficient parametrizations based on observations and simulations using the Weather Research and Forecasting model

海气界面动量交换对台风发展起着重要作用,其估算公式依赖于拖曳系数。本文利用观测和台风模拟试验对八种拖曳系数参数化方案进行系统的评估。结果表明,相对于早期的参数化方案,近些年提出的参数化方案对拖曳系数的描述有了显著的改进,主要体现在高风速条件下。不同的拖曳系数参数化方案对超强台风"海燕"的模拟路径几乎没有影响,但是对其强度、尺寸和结构的模拟结果有显著影响。     

Aerodynamic drag coefficient and roughness length for three seasons over a tropical western Indian station

Land Surface Processes Experiment (LASPEX) was conducted over semi-arid region of western India in 1997. As a part of this program, wind and temperature observations were taken using slow as well as fast response sensors over a semi-arid station Anand (22°35′N, 72°55′E) situated in Gujarat state of India. Turbulent parameters such as drag coefficient and sensible heat flux were estimated using eddy correlation method and aerodynamic roughness length was estimated using wind profiles. The analysis has been carried out for the data representing summer, monsoon and winter seasons. It was found that the wind speed does not exceed 5 ms^− 1 during the observational period considered in this study. Relationship of aerodynamic drag coefficient and roughness length with wind speed and stability has been investigated. Aerodynamic roughness length was greater in the stable conditions when the wind speed was low and it reduced drastically during convective conditions. The resulting values of aerodynamic roughness length and drag coefficient for the monsoon period agree well with values reported in literature over Indian subcontinent for homogeneous grass covered surfaces.     

The drag coefficient as determined by the dissipation method and its relation to intermittent convection in the surface layer

The influence of intermittent convection on surface-layer stress estimates during the GARP Atlantic Tropical Experiment (GATE) is described. A negative correlation between the drag coefficient (C _D) and the wind speed (% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabmyvayaara% aaaa!36DE!\[\bar U\]) is found when short averaging periods are used. Well-defined, discrete events produce this negative correlation, and these events are shown to correspond to the passage of convective plumes. Constraints on averaging times necessary to obtain reasonable stress estimates using the bulk method are discussed.Conditional sampling is used to produce average values of dissipation (% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafqyTduMbae% baaaa!37AB!\[\bar \varepsilon \]), wind speed (% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabmyvayaara% aaaa!36DE!\[\bar U\]), and virtual temperature (% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabmivayaara% WaaSbaaSqaaiaaiw8aaeqaaaaa!385B!\[\bar T_\upsilon \]) for each high turbulent intensity event, and for the quiescent periods in between. Such statistics indicate that the highly turbulent states coincide with the presence of plumes and account for the negative correlation between C _D and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabmyvayaara% aaaa!36DE!\[\bar U\]. Some of these statistics are also stability dependent.The probability distributions of the dissipation rate are bimodally log-normal which suggests that turbulence generated at two different heights is being sampled. This, along with other results of this paper, support a picture of a boundary layer which is dominated by vertical exchange.Contribution Number 409, Department of Atmospheric Sciences, University of Washington.     

一次局地强雹暴天气过程分析

通过对2004-06-06延安市强雹暴天气过程分析,得出该过程前期大气层结较为稳定,高空无垂直风切变,在天气形势向不稳定发展的趋势下可形成弱的垂直风切变;雷暴单体是在接近对流层中层中小尺度槽、切变线附近发展旺盛的,当地面形成中小尺度辐合场时有利于低层暖湿空气抬升,是雷暴单体产生降雹时段。当低层有强的干冷空气入侵,0oC层高度与常年同期比较明显偏低时,造成雷暴单体反射率因子0 dB z高度偏低。     

一次罕见的冬季冰雹天气过程分析

对河池市2002-12-18冰雹天气的天气系统和物理量场进行分析,结果表明,此次冰雹天气过程是高空槽和冷空气过境引发位势不稳定大气中的能量释放的结果,卫星云图在这次强对流天气中也起到很好的预报作用。