BOUNDARY LAYER GROWTH AND LAPSE RATE CHANGES DETERMINED BY LIDAR AND SURFACE HEAT FLUX IN SOFIA

In this study the results from a boundary layer experiment,conducted in autumn 1991 over a flat,build-up urbanarea in Southeast Sofia,together with some models for mixed layer growth rates are used to investigate the layered struc-ture of the vertical atmospheric stability distribution in the Sofia Valley.Lidar measurements of aerosol layer heightsand morning boundary layer development are combined with surface eddy correlation measurements of kinematic heatand moisture fluxes,profiles of temperature and humidity,wind speed and wind direction.A diagnostic method is pres-ented for determining vertical lapse rates using surface meteorological measurements and lidar returns observed duringthe transition from nighttime stable stratification to daytime convective boundary layer after the sunrise.     

Observing and Simulating Wind-Turbine Wakes During the Evening Transition

Wind-turbine-wake evolution during the evening transition introduces variability to wind-farm power production at a time of day typically characterized by high electricity demand. During the evening transition, the atmosphere evolves from an unstable to a stable regime, and vertical stratification of the wind profile develops as the residual planetary boundary layer decouples from the surface layer. The evolution of wind-turbine wakes during the evening transition is examined from two perspectives: wake observations from single turbines, and simulations of multiple turbine wakes using the mesoscale Weather Research and Forecasting (WRF) model. Throughout the evening transition, the wake’s wind-speed deficit and turbulence enhancement are confined within the rotor layer when the atmospheric stability changes from unstable to stable. The height variations of maximum upwind-downwind differences of wind speed and turbulence intensity gradually decrease during the evening transition. After verifying the WRF-model-simulated upwind wind speed, wind direction and turbulent kinetic energy profiles with observations, the wind-farm-scale wake evolution during the evening transition is investigated using the WRF-model wind-farm parametrization scheme. As the evening progresses, due to the presence of the wind farm, the modelled hub-height wind-speed deficit monotonically increases, the relative turbulence enhancement at hub height grows by 50%, and the downwind surface sensible heat flux increases, reducing surface cooling. Overall, the intensifying wakes from upwind turbines respond to the evolving atmospheric boundary layer during the evening transition, and undermine the power production of downwind turbines in the evening.     

BOUNDARY LAYER GROWTH AND LAPSE RATE CHANGES DETERMINED BY LIDAR AND SURFACE HEAT FLUX IN SOFIA*

In this study the results from a boundary layer experiment,conducted in autumn 1991 over a flat,build-up urban area in Southeast Sofia,together with some models for mixed layer growth rates are used to investigate the layered structure of the vertical atmospheric stability distribution in the Sofia Valley.Lidar measurements of aerosol layer heights and morning boundary layer development are combined with surface eddy correlation measurements of kinematic heat and moisture fluxes,profiles of temperature and humidity,wind speed and wind direction.A diagnostic method is presented for determining vertical lapse rates using surface meteorological measurements and lidar returns observed during the transition from nighttime stable stratification to daytime convective boundary layer after the sunrise.     

Water vapour flux profiles in the convective boundary layer

Summary Water vapour flux profiles in the atmospheric boundary layer have been derived from measurements of water vapour density fluctuations by a ground-based Differential Absorption Lidar (DIAL) and of vertical wind fluctuations by a ground-based Doppler lidar. The data were collected during the field experiment LITFASS-2003 in May/June 2003 in the area of Lindenberg, Germany. The eddy-correlation method was applied, and error estimates of ±50 W/m² for latent heat flux were found. Since the sampling error dominates the overall measurement accuracy, time intervals between 60 and 120 min were required for a reliable flux calculation, depending on wind speed. Rather large errors may occur with low wind speed because the diurnal cycle restricts the useful interval length. In the lower height range, these measurements are compared with DIAL/radar-RASS fluxes. The agreement is good when comparing covariance and error values. The lidar flux profiles are well complemented by tower measurements at 50 and 90 m above ground and by area-averaged near surface fluxes from a network of micrometeorological stations. Water vapour flux profiles in the convective boundary layer exhibit different structures mainly depending on the magnitude of the entrainment flux. In situations with dry air above the boundary layer a positive entrainment flux is observed which can even exceed the surface flux. Flux profiles which linearly increase from the surface to the top of the boundary layer are observed as well as profiles which decrease in the lower part and increase in the upper part of the boundary layer. In situations with humid air above the boundary layer the entrainment flux is about zero in the upper part of the boundary layer and the profiles in most cases show a linear decrease.     

Near-Surface Vertical Flux Divergence in the Stable Boundary Layer

Flow in the stable boundary layer is examined at four contrasting sites with greater upwind surface roughness. The surface heterogeneity is disorganized and in some cases weak as commonly occurs. With low wind speeds, the vertical divergence (or convergence) of the momentum and heat fluxes can be large near the surface in what is normally assumed to be the surface layer where such divergence is neglected. For the two most heterogeneous sites, a shallow “new” boundary layer is captured by the tower observations, analogous to an internal boundary layer but more complex. Above the new boundary layer, the magnitudes of the downward fluxes of heat and momentum increase with height in a transition layer, reach a maximum, and then decrease with height in an overlying regional boundary layer. Similar structure is observed at the site with rolling terrain where the shallow new boundary layer at the surface is identified as cold-air drainage generated by the local slope above which the flow undergoes transition to an overlying regional flow. Significant flux divergence near the surface is generated even over an ice floe for low wind speeds and in a shallow Ekman layer that forms during the polar night. For higher wind speeds, the magnitude of the downward fluxes decreases gradually with height at all levels as in a traditional boundary layer.     

Sensitivity of simulated wintertime Arctic atmosphere to vertical resolution in the ARPEGE/IFS model

The current state-of-the-art general circulation models, including several of those used by the IPCC, show considerable biases in the simulated present day high-latitude climate compared to observations and reanalysis data. These biases are most pronounced during the winter season. We here employ ideal vertical profiles of temperature and wind from turbulence-resolving simulations to perform a priori studies of the first-order eddy-viscosity closure scheme employed in the ARPEGE/IFS model. This reveals that the coarse vertical resolution (31 layers) of the model cannot be expected to realistically resolve the Arctic stable boundary layer. The curvature of the Arctic inversion and thus also the vertical turbulent-exchange processes cannot be reproduced by the coarse vertical mesh employed. To investigate how turbulent vertical exchange processes in the Arctic boundary layer are represented by the model parameterization, a simulation with high vertical resolution (90 layers in total) in the lower troposphere is performed. Results from the model simulations are validated against data from the ERA-40 reanalysis. The dependence of the surface air temperature on surface winds, surface energy fluxes, free atmosphere stability and boundary layer height is investigated. The coarse-resolution run reveals considerable biases in these parameters, and in their physical relations to surface air temperature. In the simulation with fine vertical resolution, these biases are clearly reduced. The physical relation between governing parameters for the vertical turbulent-exchange processes improves in comparison with ERA-40 data.     

The Morning Transition of the Nocturnal Boundary Layer

Six years of observations from a surface instrument site have been analysed to determine timings and factors influencing developmental changes in the near-surface wind and turbulent heat fluxes during the morning heating of the atmospheric boundary layer. A simple relationship has been found between near-surface wind speed and screen temperature, together with a predictive equation for the morning transition air temperature. Profile measurements from a probe mounted on a tethered balloon have beenused to supplement the surface data and study the processes underlying these surface relationships. The results have confirmed earlier work and have shown that both before and immediately after morning transition, almost all heating in the surface layer is due to turbulent diffusion from above. In order to explain the mechanisms involved in the relationships, a simple finite difference model has been run and validated against the profile data. The model predictions are compared with observations during both the morning and evening and the differences related to the different temperature profiles. Numerical forecasting rules for the surface wind speed and transition temperature are derived from the results.     

Mean Profiles of Moisture Fluxes in Snow-Filled Boundary Layers

Profiles of moisture fluxes have been examined for convective boundary layers containing clouds and snow, using data derived from aircraft measurements taken on four dates during the 1983/1984 University of Chicago lake-effect snow project. Flux profiles were derived from vertical stacks of aircraft cross-wind flight legs taken at various heights over Lake Michigan near the downwind shore. It was found that, if ice processes are taken into account, profiles of potential temperature and water content were very similar to those presented in past studies of convective boundary layers strongly heated from below. Profiles of total water content and equivalent potential temperature adjusted for ice were nearly invariant with height, except very near the top of the boundary layer, suggesting that internal boundary-layer mixing processes were rapid relative to the rates at which heat and vapour were transported into the boundary layer through entrainment and surface fluxes. Ice was found to play a significant, measurable role in boundary-layer moisture fluxes. It was estimated that 40 to 57% of the upward vapour flux was returned to the surface in the form of snow, converting about 45 to 64% of the surface latent heat flux into sensible heat in the snow-producing process. Assuming advective fluxes are relatively small (thought to be appropriate after the first few tens of km over the lake as suggested by past studies), the boundary layer was found to warm at a rate faster than could be explained by surface heat fluxes and latent heat releases alone, the remainder of the heating presumably coming from radiational processes and entrainment. Discussions of moisture phase change processes throughout the boundary layer and estimates of errors of these flux measurements are presented.     

Radiative Processes in the Stable Boundary Layer: Part II. The Development of the Nocturnal Boundary Layer

The interaction between radiation and turbulence in the stable boundary layer over land is explored using an idealized model, with a focus on the surface layer after the evening transition. It is shown that finer vertical resolution is required in transitional boundary layers than in developed ones. In very light winds radiative cooling determines the temperature profile, even if similarity functions without a critical Richardson number are used; standard surface similarity theory applied over thick layers then yields poor forecasts of near-surface air temperatures. These points are illustrated with field data. Simulations of the developing nocturnal boundary layer are used to explore the wider role of radiation. Comparatively, radiation is less significant within the developed stable boundary layer than during the transition; although, as previous studies have found, it remains important towards the top of the stable layer and in the residual layer. Near the ground, reducing the surface emissivity below one is found to yield modest relative radiative warming rather than intense cooling, which reduces the potential importance of radiation in the developed surface layer. The profile of the radiative heating rate may be strongly dependent on other processes, leading to quite varied behaviour.     

Elevated stratified layers observed with sodar during VTMX

Summary A combination of low frequency sodar, radar wind profiler and in-situ balloon-borne measurements of temperature and water vapor have been used to investigate the structure of elevated stratified layers within the transition layer above the nocturnal boundary layer during the Vertical Transport and Mixing Field Campaign in Salt Lake City Utah, during October, 2000. Elevated layers determined from sodar and radar vertical time sections were penetrated with a balloon-born instrument package to determine the fine scale temperature and moisture structure of the layers. As expected a potential temperature increase was found in the upper half of the layers; however the magnitude was considerably smaller than found above the daytime well-mixed layer and the vertical distance of the increase was quite variable. Mixing ratio, in the mean was found to have a relative maximum in the lower portion of the layers. It was found that the potential temperature within the layers decreased with time relative to background values, regardless of whether the layer descended or ascended.     

Representing the atmospheric boundary layer in climate models of intermediate complexity

In this study the role of atmospheric boundary layer schemes in climate models is investigated. Including a boundary layer scheme in an Earth system model of intermediate complexity (EMIC) produces only minor differences in the estimated global distribution of sensible and latent heat fluxes over land (upto about 15% of the net radiation at the surface). However, neglecting of boundary layer processes, such as the development of a well-mixed layer over land or the impact of stability on the exchange coefficient in the surface layer, leads to erroneous surface temperatures, especially in convective conditions with low wind speeds. As these conditions occur frequently, introducing a boundary layer scheme in an EMIC gives reductions in June-July-August averaged surface temperature of 1–2 °C in wet areas, to 5–7 °C in desert areas. Even a relatively simple boundary layer scheme provides reasonable estimates of the surface fluxes and surface temperatures. Detailed schemes that solve explicitly the turbulent fluxes within the boundary layer are only required when vertical profiles of potential temperature are needed.     

Resolved Versus Parametrized Boundary-Layer Plumes. Part I: A Parametrization-Oriented Conditional Sampling in Large-Eddy Simulations

A conditional sampling based on the combination of a passive tracer emitted at the surface and thermodynamic variables is proposed to characterise organized structures in large-eddy simulations of cloud-free and cloudy boundary layers. The sampling is evaluated against more traditional sampling of dry thermals or clouds. It enables the characterization of convective updrafts from the surface to the top of the boundary layer (or the top of cumulus clouds), describing in particular the transition from the sub-cloud to the cloud layer, and retrieves plume characteristics, entrainment and detrainment rates, variances and fluxes. This sampling is used to analyze the contribution of boundary-layer thermals to vertical fluxes and variances.     

地表热力非均匀性对近地层相似理论适用性影响的大涡模拟

利用非均匀地表加热的大涡模拟试验,研究了不稳定条件下地表热力非均匀性对近地层相似理论适用性的影响。结果发现,边界层的平均廓线基本不受地表热力非均匀性的影响。进一步分析发现,较大尺度的地表非均匀加热可以激发出有组织的大尺度次级环流,冷暖斑块的通量直到边界层上部才混合均匀;而当地表非均匀尺度较小时,次级环流难以形成有组织的结构,冷暖斑块的通量很快就可以混合均匀。然而,不管是哪种尺度的非均匀地表,非均匀斑块间的平流都对各斑块近地层结构产生重要影响,进而斑块近地层通量—梯度关系与相似理论产生偏差,其中风速梯度关系的偏差更为明显。最后,对目前大气模式中常用的基于相似理论的次网格非均匀地表通量参数化方法——Mosaic方法提出了改进思路。     

A vertical profile of Ozone concentration in the atmosphericboundary layer over central Taiwan

Summary  In the central region of Taiwan, ozone episodes occur most often during autumn. Two field experiments were conducted during the autumns of 1998 and 1999 to analyze the vertical profile of the boundary layer and determine its effects on ozone concentration over the region. The vertical virtual potential temperature and wind profiles were derived from tethersonde data. The NO_x, NMHC and O₃ concentration vertical profiles were monitored up to a height of 500 meters using black-covered Teflon tedler sampling bags. During the experimental periods, nighttime terrestrial long wave radiation could cause the inversion height to reach 500 meters by the following morning. It was shown that these types of synoptic structures suppress the vertical diffusion of NO_x, NMHC and O₃. During the daytime, measurements indicate that pollutants were well mixed in the upper portion of the mixing layer. At night, the ground level ozone concentration was on the decrease but increased with altitude to a height of 500 m. The NO_x decreased with altitude whereas the NMHC showed no significant variations. Received April 13, 2000 Revised July 24, 2000     

Verification of RASS-measured temperature profiles using a radio-controlled model glider

Potential sonic temperature profiles from a continuous electromagnetic, pulsed acoustic, radio acoustic sounding system (RASS) were compared with profiles recorded by a commercially available temperature, relative humidity and pressure recorder mounted in a radio-controlled model motor-glider. The glider profiles covered the period of the morning transition in the lowest 200 m of an initially stably stratified boundary layer. Comparison of the profiles shows that the sonic temperature can be calculated based on the average speed of sound in the boundary layer; this removes the need to correct for vertical velocity in each temperature profile, thus avoiding the possibility of contaminating the temperatures with measurement noise from the vertical velocity profiles. Further, the low-level cold bias that occurs with the spatially separated transmit and receive antennas of a bistatic RASS system was not significant for the present measurements as the separation between the antennas was minimised. The comparison of RASS and glider temperatures gives confidence in the use of RASS-derived temperatures for investigating the performance of boundary layer models.     

Wind profiles,momentum fluxes and roughness lengths at Cabauw revisited

We describe the results of an experiment focusing on wind speed and momentum fluxes in the atmospheric boundary layer up to 200 m. The measurements were conducted in 1996 at the Cabauw site in the Netherlands. Momentum fluxes are measured using the K-Gill Propeller Vane. Estimates of the roughness length are derived using various techniques from the wind speed and flux measurements, and the observed differences are explained by considering the source area of the meteorological parameters. A clear rough-to-smooth transition is found in the wind speed profiles at Cabauw. The internal boundary layer reaches the lowest k-vane (20 m) only in the south-west direction where the obstacle-free fetch is about 2 km. The internal boundary layer is also reflected in the roughness lengths derived from the wind speed profiles. The lower part of the profile (< 40 m) is not in equilibrium and no reliable roughness analysis can be given. The upper part of the profile can be linked to a large-scale roughness length. Roughness lengths derived from the horizontal wind speed variance and gustiness have large footprints and therefore represent a large-scale average roughness. The drag coefficient is more locally determined but still represents a large-scale roughness length when it is measured above the local internal boundary layer. The roughness length at inhomogeneous sites can therefore be determined best from drag coefficient measurements just above the local internal boundary layers directly, or indirectly from horizontal wind speed variance or gustiness. In addition, the momentum and heat fluxes along the tower are analysed and these show significant variation with height related to stability and possibly surface heterogeneity. It appears that the dimensionless wind speed gradients scale well with local fluxes for the variety of conditions considered, including the unstable cases.     

Three-dimensional numerical study of the height and mean structure of a heated planetary boundary layer

The heated boundary layer for DAY 33 of the Wangara data of southeast Australia (Clarke et al., 1971) is studied numerically with a three-dimensional model using 64000 grid points within a volume 5 km on a side and 2 km deep. Subgrid-scale transport equations were utilized in place of eddy-coefficient formulations. The rate of growth of the mixed layer is examined and parameterized, and the vertical profiles of heat flux, moisture flux and momentum fluxes are examined. The momentum boundary layer is found to coincide essentially with the mixed layer, and to grow with the latter during the hours of solar heating of the surface.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Boundedness of Turbulent Temperature Probability Distributions,and their Relation to the Vertical Profile in the Convective Boundary Layer

Higher-order moments, minima and maxima of turbulent temperature and water vapour mixing ratio probability density functions measured with an eddy-covariance system near the ground were related to each other and to vertical boundary-layer profiles of the same scalars obtained through airborne soundings. The dependence of kurtosis on squared skewness showed a kurtosis intercept below the Gaussian expectation, suggesting a compression of the probability density function by the presence of natural boundaries. This hypothesis was corroborated by comparing actual minima and maxima of turbulent fluctuations to estimates obtained from the first four sample moments by fitting a four-parameter beta distribution. The most sharply defined boundaries were found for the minima of temperature datasets during the day, indicating that negative temperature fluctuations at the sensor are limited by the availability of lower temperatures in the boundary layer. By comparison to vertical profiles, it could be verified that the turbulent minimum of temperature near the ground is close to the minimum of potential temperature in the boundary layer. The turbulent minimum of water vapour mixing ratio was found to be equal to the mixing ratio at a height above the minimum of the temperature profile. This height roughly agrees with the top of the non-local unstable domain according to bulk Richardson number profiles. We conclude that turbulence statistics measured near the surface cannot be solely explained by local effects, but contain information about the whole boundary layer including the entrainment zone.     

乌鲁木齐大气边界层结构的四季特征初探

利用乌鲁木齐市2011～2012年08时、20时L波段(1型)雷达探测的高空资料建立了乌鲁木齐大气边界层气象要素数据库,分析了乌鲁木齐边界层内气温、风向、风速和相对湿度的垂直分布及其时间变化特征。结果表明：边界层内温度廓线的日变化和季节变化比较显著,各月均有逆温出现,且08时较20时更易出现逆温,冬季08时逆温层厚度较厚且强度最大。边界层内夏、冬两季风速随高度变化波动较大,春、秋两季变化较小。近地层春、夏、秋三季08时盛行西南偏南风,冬季盛行偏东风和西南风;20时春季盛行东北风,夏秋盛行偏北风和西北风,冬季则盛行东风和东北偏东风。08时、20时风向均随高度的增加呈明显的向右偏转趋势,且日风向的变化具有明显的“山谷风”特点。08、20时的相对湿度冬季最大,夏季最小,且随高度增加,春、夏两季08、20时相对湿度的变化较大。     

Numerical simulations of the tropical air-sea planetary boundary layer

A one-dimensional numerical model of the planetary boundary layer was used to investigate thermal and kinetic energy budgets. The simulation experiments were based on two sets of data. The first set was based on a ‘typical’ June with climatological data extracted for the oceanic region slightly northeast of Barbados. The second set used data from the third phase of project BOMEX, for approximately the same area and time of year as the first set. Comparison with observations of three simulated elements (viz., sea surface temperature and wind and humidity at 6 m) which are important in determining the near-interface energy transports shows that:

1. the model is capable of realistic simulations of both ‘typical’ conditions, and conditions for a specific four-day period;
2. the model is capable of realistically simulating the differences between prevailing values of these parameters in the two cases (‘typical’ and specific four-day period).

The simulated interface fluxes are those of incoming and outgoing short- and long-wave radiation; transmitted radiation at -0.5 m in the ocean, sensible heat transfer into the ocean and air, and latent heat flux of evaporation. Comparison with observational analyses shows that the diurnal variations in net radiation and heat storage in the mixed layer are realistically simulated. The simulated values of evaporation are consistent with other estimates for both ‘typical’ conditions and specific conditions during this four-day period. The rate of heat storage varies between +51 and -37 percent of the diurnal maximum incoming radiation, and the evaporation varies between +16% and -13% of this term. The non-dimensional transfer coefficients (C _D, C_T, C_q) computed from the model show general agreement with the coefficients calculated from observations in the simulated region (Pondet al., 1971). The simulated vertical profiles of temperature are in general agreement with observed profiles, except in the uppermost portions of the atmospheric boundary layer where deviations of approximately 1.5C occur. Simulated vertical profiles of wind speed are generally consistent with observed profiles, with the largest deviations appearing to be of the order of 0.5 m s^-1. Simulated vertical profiles of the eddy fluxes of sensible heat, water vapor, and momentum are generally consistent with Bunker's (1970) aircraft-based measurements of these quantities. The time averages of these simulated profiles show regular decreases with height, while simulated profiles for specific hours of the day show intermediate maxima and minima, which are also seen in the measured profiles. The vertically integrated kinetic energy budgets of the modelled atmospheric layer are presented through the four terms of the kinetic energy budget, viz., the upper and the lower boundary drags, dissipation, and potential-to-kinetic conversion. The dominant terms in the atmospheric energy budgets are the production and dissipation terms, with kinetic energy being exported both to the overlying atmospheric layer and to the underlying oceanic layer at rates of about 2 to 6% of the production, respectively. Comparisons between the climatological and BOMEX simulations are presented. The vertically integrated humidity budgets are presented for the two simulation experiments. Under ‘typical’ conditions, the humidity budget reveals an upper boundary flux of about +29% of the lower boundary flux with the vertically integrated advective flux being -59% of the lower flux. For the specific four-day simulation, the upper boundary flux and advection are about +28 and -70%, respectively, of the lower boundary flux.