Steady-state characteristics of inversions capping a well-mixed planetary boundary layer

The type of inversion discussed in this paper is essentially defined by subsidence, mixing due to thermally generated turbulent energy and a radiative flux difference at the inversion. A concept similar to that of Lilly (1968) is applied, assuming a well-mixed layer below the inversion and including advective and radiative processes.The characteristics of the inversion and of the whole PBL (e.g., height of inversion, height of cloud base, strength of inversion, flux-profiles) are investigated for their dependence on external parameters such as horizontal wind field divergence, advection, surface temperature excess, wind speed and surface temperature. This is done for steady-state conditions and gives considerable insight into the processes maintaining the type of inversion under consideration.A second goal is to present typical inversion structures, which can be found in certain climatic regions. The profiles of the state parameters and the energy-fluxes for the Trade-Wind region, the cold water area off the west-coast of California, the Norwegian Sea and the Arctic Ocean differ considerably.     

Development of a Single-Column Model in RegCM4 and Its Preliminary Application for Evaluating PBL Schemes in Simulating the Dry Convection Boundary Layer

A single-column model (SCM) is developed in the regional climate model RegCM4. The evolution of a dry convection boundary layer (DCBL) is used to evaluate this SCM. Moreover, four planetary boundary layer (PBL) schemes, namely the Holtslag-Boville scheme (HB), Yonsei University scheme (YSU), and two University of Washington schemes (UW01, Grenier-Bretherton-McCaa scheme and UW09, Bretherton-Park scheme), are compared by using the SCM approach. A large-eddy simulation (LES) of the DCBL is performed as a benchmark to examine how well a PBL parameterization scheme reproduces the LES results, and several diagnostic outputs are compared to evaluate the schemes. The results show that the SCM is proper constructed. In general, with the DCBL case, the YSU scheme performs best for reproducing the LES results, which include well-mixed features and vertical sensible heat fluxes; the simulated wind speed, turbulent kinetic energy, entrainment flux, and height of the entrainment zone are all underestimated in the UW09; the UW01 has all those biases of the UW09 but larger, and the simulated potential temperature is not well mixed; the HB is the least skillful scheme, by which the PBL height, entrainment flux, height of the entrainment zone, and the vertical gradients within the mixed layer are all overestimated, and a inversion layer near the top of the surface layer is wrongly simulated.Although more cases and further testing are required, these simulations show encouraging results towards the use of this SCM framework for evaluating the simulated physical processes by the RegCM4.     

A note on surface layer parameterizations in the weather research and forecast model

Two surface layer parameterization schemes along with five planetary boundary layer (PBL) schemes in the Weather Research and Forecasting model (WRF) are analyzed in order to evaluate the performance of the WRF model in simulating the surface variables and turbulent fluxes over an Indian sub-continent region. These surface layer schemes are based on the fifth-generation Pennsylvania State University—National Center for Atmospheric Research Mesoscale Model (MM5) parameterization; (a) Old MM5 scheme having Businger-Dyer similarity functions and (b) revised MM5 scheme utilizing the functions that are valid for full ranges of atmospheric stabilities. The study suggests that each PBL scheme can reproduce the diurnal variation of 2 m temperature, momentum flux and sensible heat flux irrespective of the surface layer scheme used for the simulations. However, a comparison of model-simulated values of surface variables and turbulent fluxes with observed values suggests that each PBL scheme is found to systematically over-estimate the nocturnal 2 m temperature and 10 m wind speed with both the revised and old schemes during stable conditions.     

The sensitivity of k-ɛ model computations of the neutral planetary boundary layer to baroclinicity

While the importance of baroclinicity in determining the structure of the planetary boundary layer (PBL) is well recognized, the actual effect of baroclinicity on the structure is not well understood. Results based on simulations obtained using the turbulent kinetic energy-dissipation rate of turbulent kinetic energy closure model of the turbulent flow in a neutral baroclinic PBL provide additional insight into the role of baroclinicity. The baroclinic PBL is characterized by significant shear production of turbulent kinetic energy throughout the complete boundary-layer depth. The turbulent mixing length is bounded by the presence of a stable temperature inversion layer indicating that the depth of the baroclinic PBL is determined by the inversion height. Significant turbulent shear stresses exist throughout the baroclinic PBL and the air is relatively well-mixed except in the surface layer.     

A Study On Atmospheric Boundary-Layer Characteristics At Anand, India Using Lsp Experimental Data Sets

An attempt is made to study the planetary boundary layer (PBL) characteristics during the winter period at Anand (22.4°N, 72.6°E), a semi-arid region, which is located in the western part of India. A one-dimensional turbulent kinetic energy (TKE) closure model is used for the study. The structure of the PBL,which consists of profiles of zonal and meridional components of wind,potential temperature and specific humidity, is simulated. A one-dimensional soil heat and moisture transport parameterization scheme is incorporated for the accurate representation of the energy exchange processes at the soil-atmosphere interface. The diurnal variation of fluxes of sensible heat, latent heat, shortwave radiation, net radiation and soil flux, soil temperature at different depths, Richardson number and TKE at the height of the constant flux layer is studied. The model predictions are compared with the available observations obtained from a special Land Surface Processes (LSP) experiment.     

Climatological aspects of the tropical convective boundary layer

The characteristics of a boundary layer depend both on conditions at the surface and in the interior of the medium. In the undisturbed tropics, the latter are largely determined by subsidence and by infrared radiational cooling. One-dimensional models are used to establish relationships between the inversion height, subsidence, upper-air humidity and sea-surface temperature. In particular, it is shown that a universally colder tropical ocean would probably be covered by more extensive clouds.Contribution No. 1700 Rosenstiel School of Marine and Atmospheric Science, University of Miami.     

A Case Study Of An On-Ice Air Flow Over The Arctic Marginal Sea-Ice Zone

A case study of warm air advection over the Arctic marginalsea-ice zone is presented, based on aircraft observations with direct flux measurements carriedout in early spring, 1998. A shallow atmospheric boundary layer (ABL) was observed, which wasgradually cooling with distance downwind of the ice edge. This process was mainly connected with astrong stable stratification and downward turbulent heat fluxes of about 10–20 W m^-2, but wasalso due to radiative cooling. Two mesoscale models, one hydrostatic and the other non-hydrostatic,having different turbulence closures, were applied. Despite these fundamental differences betweenthe models, the results of both agreed well with the observed data. Various closure assumptions had amore crucial influence on the results than the differences between the models.Such an assumption was, for example,the parameterization of the surface roughness for momentum (z₀) and heat (z_T). This stronglyaffected the wind and temperature fields not only close to the surface but also within and abovethe temperature inversion layer. The best results were achieved using a formulation for z₀ that took intoaccount the form drag effect of sea-ice ridges together withz_T = 0.1z₀. The stability within theelevated inversion strongly depended on the minimum eddy diffusivity K_min. A simple ad hocparameterization seems applicable, where K_min is calculated as 0.005 timesthe neutral eddy diffusivity. Although the longwave radiative cooling was largest within the ABL, theapplication of a radiation scheme was less important there than above the ABL. This was related to theinteraction of the turbulent and radiative fluxes. To reproduce the strong inversion, it wasnecessary to use vertical and horizontal resolutions higher than those applied in most regional andlarge-scale atmospheric models.     

北京地区夏季边界层结构日变化的高分辨模拟对比

使用WRF中尺度数值模式, 分别选用两种不同的边界层参数化方案 (MYJ, YSU) 和3种陆面参数化方案 (SLAB, Noah, RUC), 对2004年7月1日08:00—7月4日20:00 (北京时) 北京地区夏季边界层结构进行1 km的高分辨模拟。对比分析了近地面层风场、温度场以及边界层的日变化特征, 结果发现:WRF模式基本模拟出了北京夏季边界层的日变化特征; 在边界层方案中, MYJ方案描述的边界层结构较YSU方案合理; Noah陆面模式较好地反映了城市的热岛效应; 无降水时, 风速及边界层高度对于陆面过程不敏感, 而降水发生后, 陆面过程对于边界层结构的影响增大; 各方案模拟的城区风速明显偏大, 这是因为没有充分考虑城市建筑物的阻力作用。     

An efficient PBL model for global circulation models-design and validation

An efficient, pianetary boundary layer (PBL) model is developed and validated with empirical data for applications in general circulation models (GCMs). The purpose of this PBL model is to establish the turbulent surface fluxes as a function of the principal external PBL parameters in a numerically efficient way. It consists of a surface layer and a mixed layer matched together with the conditions of constant momentum and heat flux at the interface. An algebraic solution to the mean momentum equations describes the mixed-layer velocity profile and thus determines the surface wind vector. The velocity profile is globally valid by incorporating the effect of variable Coriolis force without becoming singular at the equator. Turbulent diffusion depends on atmospheric stability and is modeled in the surface layer by a drag law and with first-order closure in the mixed layer. Radiative cooling in the stably stratified PBL is considered in a simple manner. The coupled system is solved by an iterative method. In order to preserve the computational efficiency of the large-scale model, the PBL model is implemented into the GISS GCM by means of look-up tables with the bulk PBL Richardson number, PBL depth, neutral drag coefficient, and latitude as independent variables.A validation of the PBL model with observed data in the form of Rossby number similarity theory shows that the internal feedback mechanisms are represented correctly. The model, however, underpredicted the sensible heat-flux. A subsequent correction in the turbulence parameterization yields better agreement with the empirical data. The behavior of the principal internal PBL quantities is presented for a range of thermal stabilities and latitudes.     

Evidence of a Mixed-Layer Dynamics Contribution to the Entrainment Process

The internal thermal boundary layer developing over the Mediterranean during a cold-air outbreak associated with a Tramontane event has been studied by means of airborne lidar, in situ sensors, and a modelling approach that consisted of nesting the University of Washington (UW) planetary boundary-layer (PBL) model in an advective zero-order jump model. This approach bypasses some of the deficiencies associated with each model: the absence of the dynamics in the mixed layer for the zero-order jump model and the lack of an inversion at the PBL top for the UW PBL model. Particular attention is given to the parameterization of the entrainment flux at the PBL top. Values of the entrainment closure parameter derived with the model when matching PBL structure observations are much lower than those derived with standard zero-order jump models. They also are in good agreement with values measured in different meteorological situations by other studies. This improvement is a result of the introduction of turbulent kinetic energy production in the mixed layer.     

Surface energy balance,parameterizations of boundary-layer heights and the application of resistance laws near an Antarctic Ice Shelf front

A study of the surface energy balance with turbulent fluxes obtained by the Monin-Obukhov similarity theory and a comparison with results for resistance laws are presented for the strong baroclinic conditions in the vicinity of the Filchner/Ronne Ice Shelf front. The data are taken from a field experiment in the Antarctic summer season 1983/84. For the first time in the coastal Antarctic region, this data set comprises synchronous energy balance measurements over the polynya and the ice shelf together with soundings of the boundary layer, yielding vertical profiles of the wind velocity and temperature over the ice shelf, at the ice shelf front and over the polynya.Over the ice shelf, the radiation balance is the largest component of the energy fluxes and is mainly compensated by the subsurface energy flux and the turbulent heat flux in the daily mean. Over the polynya, turbulent fluxes of sensible and latent heat lead to large energy losses of the water surface in the night-time and in situations of very low air temperatures.Different parameterizations for boundary-layer height are compared using tethered sonde and energy balance measurements. With the height of the inversion base over the polynya and the height of the critical bulk Richardson number over the ice shelf, external parameters for the application of resistance laws were determined. The comparison of turbulent surface fluxes obtained by the energy balance measurements and by the resistance laws shows good agreement for the convective conditions over the polynya. For the stably stratified boundary layer over the ice shelf with small amounts of the turbulent heat flux, the deviation is large for the case of a cold air outflow with a superposed inertial oscillation.     

Helicopter-Borne Flux Measurements in the Nocturnal Boundary Layer Over Land – a Case Study

This case study introduces measurements of turbulent fluxes in a nocturnal boundary layer in North Germany with the new helicopter-borne turbulence measurement system HELIPOD, a detailed data analysis and examination in regard of systematic errors of the instrument, and some comparison with local similarity theory and experiments of the past, in order to confirm the occurrence of small vertical turbulent fluxes. The examined nocturnal boundary layer offered excellent conditions to analyse the quality of the measurement system. In this connection, a detailed look at a strong ground-based inversion disclosed small turbulent fluxes with a spectral maximum at ten metres wavelength or less, embedded in intermittent turbulence. For verification of these fluxes, the measurements were compared with well established results from past experiments. Local similarity theory was applied to calculate dimensionless variances of the turbulent quantities, which were found in good agreement with other observations. Since shear and stratification varied significantly on the horizontal flight legs due to global intermittency, a method was developed to determine vertical gradients on a horizontal flight pattern, by use of small fluctuations of the measurement height. With these locally determined gradients, gradient transport theory became applicable and the turbulent diffusivities for heat and momentum, the Richardson number, and the flux Richardson number were estimated within isolated strong turbulent outbursts. Within these outbursts the flux Richardson number was found between 0.1 and 0.2. The functional relationship between the gradient Richardson number and the turbulent Prandtl number agreed well with observations in past experiments and large eddy simulation. The impact of the stratification on the vertical turbulent exchange, as already described for the surface layer using Monin–Obukhov similarity, was analogously observed in the very stably stratified bulk flow when local scaling was applied.     

The Uncertainty of Tropical Cyclone Intensity and Structure Based on Different Parameterization Schemes of Planetary Boundary Layer

Based on different parameterization schemes of planetary boundary layer (PBL), the uncertainty of intensity and structure of the Super-strong Typhoon Rammasun (1409) is investigated using the WRF model (v3.4) with six PBL parameterization schemes. Results indicate that PBL uncertainty leads to the uncertainty in tropical cyclone (TC) prediction, which increases with forecast time. The uncertainty in TC prediction is mainly reflected in the uncertainty in TC intensity, with significant differences in the TC intensity forecasts using various PBL schemes. The uncertainty in TC prediction is also reflected in the uncertainty in TC structures. Greater intensity is accompanied by smaller vortex width, tighter vortex structure, stronger wind in the near-surface layer and middle and lower troposphere, stronger inflow (outflow) wind at the lower (upper) levels, stronger vertical upward wind, smaller thickness of the eye wall, smaller outward extension of the eye wall, and warmer warm core at the upper levels of eye. PBL height, surface upward heat flux and water vapor flux are important factors that cause the uncertainty in TC intensity and structure. The more surface upward heat flux and water vapor flux and the lower PBL height, the faster TC development and the stronger TC intensity.     

A bulk model for the atmospheric planetary boundary layer

The integrated momentum and thermodynamic equations through the planetary boundary layer (PBL) are solved numerically to predict the mean changes of wind and potential temperature from which surface fluxes are computed using bulk transfer coefficients of momentum and heat. The second part of the study involves a formulation and testing of a PBL height model based on the turbulent energy budget equation where turbulent fluxes of wind and heat are considered as the source of energy. The model exhibits capability of predicting the PBL height development for both stable and unstable regimes of observed conditions. Results of the model agree favourably with those of Deardorff's (1974a) and Tennekes' (1973) models in convective conditions.Contribution number 396.     

陆面面积平均通量的参数化问

目前,气候和大气边界层物理研究中一个十分重要的研究方向就是面积平均通量的参数化问题。本文对这一研究方向中存在的问题、可能的解决方法和目前的研究进展情况进行了阐述,简要介绍了非均匀陆面影响的高度及其尺度划分,阐述了混合高度、参考层高度（观测高度和模式第一层高度）、近地层高度、内边界层高度、平衡层高度、粗糙度副层和边界层高度等之间的关系及其在非均匀尺度划分中的作用,并且阐述了整体输送公式在不同尺度的非均匀陆面中存在的问题及相应的可能解决办法。同时还对中尺度的非均匀陆面驱动的一类非经典中尺度环流的参数化,即中尺度通量的参数化问题进行了评述。最后针对内蒙古草原实验和青藏高原实验等具体问题,提出了边界层观测和非均匀陆面参数化方法的几点问题。     

Impact of waves on air-sea exchange of sensible heat and momentum

The impact of sea waves on sensible heat and momentum fluxes is described. The approach is based on the conservation of heat and momentum in the marine atmospheric surface layer. The experimental fact that the drag coefficient above the sea increases considerably with increasing wind speed, while the exchange coefficient for sensible heat (Stanton number) remains virtually independent of wind speed, is explained by a different balance of the turbulent and the wave-induced parts in the total fluxes of momentum and sensible heat.Organised motions induced by waves support the wave-induced stress which dominates the surface momentum flux. These organised motions do not contribute to the vertical flux of heat. The heat flux above waves is determined, in part, by the influence of waves upon the turbulence diffusivity.The turbulence diffusivity is altered by waves in an indirect way. The wave-induced stress dominates the surface flux and decays rapidly with height. Therefore the turbulent stress above waves is no longer constant with height. That changes the balance of the turbulent kinetic energy and of the dissipation rate and, hence the diffusivity.The dependence of the exchange coefficient for heat on wind speed is usually parameterized in terms of a constant Stanton number. However, an increase of the exchange coefficient with wind speed is not ruled out by field measurements and could be parametrized in terms of a constant temperature roughness length. Because of the large scatter, field data do not allow us to establish the actual dependence. The exchange coefficient for sensible heat, calculated from the model, is virtually independent of wind speed in the range of 3–10 ms^-1. For wind speeds above 10 ms^-1 an increase of 10% is obtained, which is smaller than that following from the constant roughness length parameterization.The investigation was in part supported by the Netherlands Geosciences Foundation (GOA) with financial aid from the Netherlands Organization for Scientific Research (NWO).     

MM5模式中城市冠层参数化方案的设计及其数值试验

文中在综合国外一些较先进的中尺度模式城市作用参数化方案的基础上 ,从城市下垫面结构对城市边界层大气作用的物理机制及实际应用两方面出发 ,对城市下垫面结构和人为活动等因素对边界层结构的影响及中尺度模式中城市化作用的合理体现等问题进行了较全面的考虑 ,改进和设计出能够较全面、细致地描述城市结构对大气边界层动力、热力结构的影响 ,且适合中尺度模式结构特点的城市冠层参数化方案 (UCP) ,并实现了其与MM5模式的耦合。进行了耦合后的UCP方案及采用原城市作用方案的MM5模式对BECAPEX试验期间北京地区气象条件多重嵌套细尺度进行了模拟试验 ,并与观测结果对比 ,结果表明 :相比于MM 5模式中原有表示城市作用的参数化方案来讲 ,设计的UCP方案在很大程度上提高了MM 5模式对城市边界层热力和动力结构的模拟能力。     

Airborne Measurements in the Stable Boundary Layer over the Larsen Ice Shelf,Antarctica

We present aircraft measurements of boundary-layer structure and surface turbulent fluxes from a flight over the Larsen Ice Shelf, Antarctica. Warm advection, associated with föhn flow, led to the formation of a stable boundary layer over the ice shelf, with a well-defined low-level jet at the top of the surface inversion. The strong shear associated with the jet kept the gradient Richardson number small and maintained a turbulent boundary layer over a depth of at least 600 m. The net surface energy balance amounted to 52 Wm^?2, equivalent to a melt rate of 13 mm water per day, with net radiation (48 Wm^?2) making the largest contribution to melt. The contribution from the sensible heat flux (13 Wm^?2) was largely balanced by an upwards latent heat flux (?9 Wm^?2). These measurements provide insight into the processes that control surface melt rates in an area that has experienced recent rapid warming and deglaciation.     

Simulation of surface energy fluxes using high-resolution non-hydrostatic simulations and comparisons with measurements for the LITFASS-2003 experiment

Land-surface heterogeneity effects on the subgrid scale of regional climate and numerical weather prediction models are of vital interest for the energy and mass exchange between the surface and the atmospheric boundary layer. High-resolution numerical model simulations can be used to quantify these effects, and are a tool used to obtain area-averaged surface fluxes over heterogeneous land surfaces. We present high-resolution model simulations for the LITFASS area near Berlin during the LITFASS-2003 experiment, which were carried out using the non-hydrostatic model FOOT3DK of the University of Köln with horizontal resolutions of 1 km and 250 m. The LITFASS-2003 experimental dataset is used for comparison. The screen level quantities show good quality for the simulated pressure, temperature, humidity and wind speed and direction. Averaged over the four week experimental period, simulated surface energy fluxes at land stations show a small bias for the turbulent heat fluxes and an underestimation of the net radiation caused by excessive cloudiness in the simulations. For eight selected days with low cloud amounts, the net radiation bias is close to zero, but the sensible heat flux shows a strong positive bias. Large differences are found for latent heat fluxes over a lake, which are partly due to local effects on the measurements, but an additional problem seems to be the overestimation of the turbulent exchange under stable conditions in the daytime internal boundary layer over the lake. In the area average over the LITFASS area of 20 ×  20 km², again a strong positive bias of 70 W m^?2 for the sensible heat is present. For the low soil moisture conditions during June 2003, the simulation of the turbulent heat fluxes is sensitive to variations in the soil type and its hydrological properties. Under these conditions, the supply of ground water to the lowest soil layer should be accounted for. Different area-averaging methods are tested. The experimental set-up of the LITFASS-2003 experiment is found to be well suited for the computation of area-averaged turbulent heat fluxes.