A simple model of the atmospheric boundary layer; sensitivity to surface evaporation

A simple formulation of the boundary layer is developed for use in large-scale models and other situations where simplicity is required. The formulation is suited for use in models where some resolution is possible within the boundary layer, but where the resolution is insufficient for resolving the detailed boundary-layer structure and overlying capping inversion. Surface fluxes are represented in terms of similarity theory while turbulent diffusivities above the surface layer are formulated in terms of bulk similarity considerations and matching conditions at the top of the surface layer. The boundary-layer depth is expressed in terms of a bulk Richardson number which is modified to include the influence of thermals. Attention is devoted to the interrelationship between predicted boundary-layer growth, the turbulent diffusivity profile, countergradient heat flux and truncation errors.The model predicts growth of the convectively mixed layer reasonably well and is well-behaved in cases of weak surface heat flux and transitions between stable and unstable cases. The evolution of the modelled boundary layer is studied for different ratios of surface evaporation to potential evaporation. Typical variations of surface evaporation result in a much greater variation in boundary-layer depth than that caused by the choice of the boundary-layer depth formulation.     

A model of the planetary boundary layer over a snow surface

A model of the planetary boundary layer over a snow surface has been developed. It contains the vertical heat exchange processes due to radiation, conduction, and atmospheric turbulence. Parametrization of the boundary layer is based on similarity functions developed by Hoffert and Sud (1976), which involve a dimensionless variable, ζ, dependent on boundary-layer height and a localized Monin-Obukhov length. The model also contains the atmospheric surface layer and the snowpack itself, where snowmelt and snow evaporation are calculated. The results indicate a strong dependence of surface temperatures, especially at night, on the bursts of turbulence which result from the frictional damping of surface-layer winds during periods of high stability, as described by Businger (1973). The model also shows the cooling and drying effect of the snow on the atmosphere, which may be the mechanism for air mass transformation in sub-Arctic regions.     

A Model of the Planetary Boundary Layer Over a Snow Surface

A model of the planetary boundary layer over a snow surface has been developed. It contains the vertical heat exchange processes due to radiation, conduction, and atmospheric turbulence. Parametrization of the boundary layer is based on similarity functions developed by Hoffert and Sud (1976), which involve a dimensionless variable, ζ, dependent on boundary-layer height and a localized Monin-Obukhov length. The model also contains the atmospheric surface layer and the snowpack itself, where snowmelt and snow evaporation are calculated. The results indicate a strong dependence of surface temperatures, especially at night, on the bursts of turbulence which result from the frictional damping of surface-layer winds during periods of high stability, as described by Businger (1973). The model also shows the cooling and drying effect of the snow on the atmosphere, which may be the mechanism for air mass transformation in sub-Arctic regions.     

Turbulent transfer at the ground: On verification of a simple predictive model

Observations on heat transfer from ground-based plates and evaporation from free water surfaces in the laboratory and in the field are compared with predictions from a simple model. The model relates the convective transfer coefficient (or boundary-layer resistance) at any point on a surface to the momentum transfer (friction velocity) in the boundary layer immediately above it and should be applicable to practically any soil surface, open or vegetated.Heat-transfer data showed a standard deviation of 25%; between predictions and observations. Evaporation data provided only order-of-magnitude confirmation of the model because of uncertainty in effective water vapor density above small free-water surfaces.     

Characteristics and Numerical Simulations of Extremely Large Atmospheric Boundary-layer Heights over an Arid Region in North-west China

Over arid regions in north-west China, the atmospheric boundary layer can be extremely high during daytime in late spring and summer. For instance, the depth of the observed convective boundary layer can exceed 3,000 m or even be up to 4,000 m at some stations. In order to characterize the atmospheric boundary-layer (ABL) conditions and to understand the mechanisms that produce such an extreme boundary-layer height, an advanced research version of the community weather research and forecasting numerical model (WRF) is employed to simulate observed extreme boundary-layer heights in May 2000. The ability of the WRF model in simulating the atmospheric boundary layer over arid areas is evaluated. Several key parameters that contribute to the extremely deep boundary layer are identified through sensitivity experiments, and it is found that the WRF model is able to capture characteristics of the observed deep atmospheric boundary layer. Results demonstrate the influence of soil moisture and surface albedo on the simulation of the extremely deep boundary layer. In addition, the choice of land-surface model and forecast lead times also plays a role in the accurate numerical simulation of the ABL height.     

A bulk model of the well-mixed boundary layer

A bulk boundary-layer model is developed to predict surface fluxes and conditions in the well-mixed layer between the surface and the lower troposphere. The model includes the effects of all the dominant processes, including advection, in a dry boundary layer. The numerical model is compared with theoretical predictions for the growth of an internal boundary layer, and it is used to simulate the generation of a sea breeze by the diurnal cycle of radiative heating.     

SENSITIVITY TESTS OF INTERACTION BETWEEN LAND SURFACE PHYSICAL PROCESS AND ATMOSPHERIC BOUNDARY LAYER

In this paper,an interactive model between land surface physical process and atmosphereboundary layer is established,and is used to simulate the features of soil environmental physics,surface heat fluxes,evaporation from soil and evapotranspiration from vegetation and structures ofatmosphere boundary layer over grassland underlying.The sensitivity experiments are engaged inprimary physics parameters.The results show that this model can obtain reasonable simulation fordiurnal variations of heat balance,soil volumetric water content,resistance of vegetationevaporation,flux of surface moisture,and profiles of turbulent exchange coefficient,turbulentmomentum,potential temperature,and specific humidity.The model developed can be used tostudy the interaction between land surface processes and atmospheric boundary layer in cityregions,and can also be used in the simulation of regional climate incorporating a mesoscalemodel.     

SENSITIVITY TESTS OF INTERACTION BETWEEN LAND SURFACE PHYSICAL PROCESS AND ATMOSPHERIC BOUNDARY LAYER*

In this paper,an interactive model between land surface physical process and atmosphere boundary layer is established,and is used to simulate the features of soil environmental physics,surface heat fluxes,evaporation from soil and evapotranspiration from vegetation and structures of atmosphere boundary layer over grassland underlying.The sensitivity experiments are engaged in primary physics parameters.The results show that this model can obtain reasonable simulation for diurnal variations of heat balance,soil volumetric water content,resistance of vegetation evaporation,flux of surface moisture,and profiles of turbulent exchange coefficient,turbulent momentum,potential temperature,and specific humidity.The model developed can be used to study the interaction between land surface processes and atmospheric boundary layer in city regions,and can also be used in the simulation of regional climate incorporating a mesoscale model.     

陆面过程和大气边界层相互作用敏感性实验

文中建立了一个研究陆面物理过程与大气边界层相互作用的模式。模拟了草原下垫面的土壤 环境物理、地面热量通量、蒸发、蒸散及大气边界层结构特征。并对主要的环境物理参数进 行了敏感性实验。结果表明，本模式能合理地模拟地表热量平衡、土壤体积含水量、植被蒸 发阻抗、地表水汽通量日变化和湍流交换系数、湍流动能、位温和比湿廓线等。该模式还可 进一步应用于研究城市陆面物理过程与大气边界层相互作用机制，及与中尺度大气模式耦合用于区域气候的研究。     

Interpretation of the positive low-cloud feedback predicted by a climate model under global warming

The response of low-level clouds to climate change has been identified as a major contributor to the uncertainty in climate sensitivity estimates among climate models. By analyzing the behaviour of low-level clouds in a hierarchy of models (coupled ocean-atmosphere model, atmospheric general circulation model, aqua-planet model, single-column model) using the same physical parameterizations, this study proposes an interpretation of the strong positive low-cloud feedback predicted by the IPSL-CM5A climate model under climate change. In a warmer climate, the model predicts an enhanced clear-sky radiative cooling, stronger surface turbulent fluxes, a deepening and a drying of the planetary boundary layer, and a decrease of tropical low-clouds in regimes of weak subsidence. We show that the decrease of low-level clouds critically depends on the change in the vertical advection of moist static energy from the free troposphere to the boundary-layer. This change is dominated by variations in the vertical gradient of moist static energy between the surface and the free troposphere just above the boundary-layer. In a warmer climate, the thermodynamical relationship of Clausius-Clapeyron increases this vertical gradient, and then the import by large-scale subsidence of low moist static energy and dry air into the boundary layer. This results in a decrease of the low-level cloudiness and in a weakening of the radiative cooling of the boundary layer by low-level clouds. The energetic framework proposed in this study might help to interpret inter-model differences in low-cloud feedbacks under climate change.     

Boundary-Layer Entrainment Estimation Through Assimilation of Radiosonde and Micrometeorological Data into a Mixed-Layer Model

In this study we estimate boundary-layer growth and entrainment bycombining radiosonde and micrometeorological observations with asimple coupled boundary-layer and land surface model. A variational(smoothing) approach is used to find the optimal estimate ofentrainment over the daytime window. This method is appealingbecause it accounts for the uncertainty in the various data streams,while enforcing the dynamics of the model (i.e., water and energybudgets in the boundary layer, mixed-layer growth, etc.). Thetraditional variational framework was modified in this study toinclude an ensemble approach, which not only yields a (mean) estimateof entrainment, but a measure of its uncertainty as well. Themethodology is applied to a field experiment site in Kansas. Resultsfrom this study indicate a much larger ratio of entrainment to surfacefluxes compared to early literature values from other sites. However,our results are consistent with recent estimates at the site usingindependent estimation methods. In tests where radiosonde data werewithheld, reasonable skill in entrainment estimation was still shown,suggesting the potential for more widespread applications where onlymicrometeorological data are available. Finally, the data assimilationframework presented here has not traditionally been used inatmospheric boundary-layer studies, and may provide a useful approachfor studying other aspects of the boundary layer in the future.     

Effect of surface layer thickness on simultaneous transport of heat and water in a bare soil and its implications for land surface schemes

Abstract

A physically‐based multi‐layer numerical model is developed to determine the coupled transport of heat and water in the soil and in the soil‐atmosphere boundary layer. Using inputs of standard weather data and initial soil conditions the model is capable of predicting the surface energy balance components as well as water content and temperature profiles in the soil. It is used to predict these variables for a bare silt loam soil under two tillage treatments, viz. culti‐packed and left loose after disc‐harrowing, and the predicted results are compared with measurements. Very good agreement between the model predictions and measured evaporation and heat fluxes and soil water and temperatures for a ten‐day period shows that the model is capable of simulating the coupled transport of soil heat and soil water and their transfer across the soil surface‐atmosphere interface adequately.

Model predictions were compared with those of CLASS (Canadian Land Surface Scheme). It is shown that CLASS, version 2.6, provides good estimates of evaporation and hence the latent heat flux density, Q_E, under wetter soil conditions, but overestimates Q_E at moderately wet soil conditions and underestimates it under dry soil conditions. Under dry to moderately wet soil conditions the calculation of evaporation from bare soil is very sensitive to the thickness of the top layer particularly as the thickness approaches 10 cm.     

UNCERTAINTY IN THE SPECIFICATION OF SURFACE CHARACTERISTICS: A STUDY OF PREDICTION ERRORS IN THE BOUNDARY LAYER

The effects of uncertainty in the specification of surface characteristics on simulated atmospheric boundary layer (ABL) processes and structure were investigated using a one-dimensional soil-vegetation-boundary layer model. Observational data from the First International Satellite Land Surface Climatology Project Field Experiment were selected to quantify prediction errors in simulated boundary-layer parameters. Several numerical 12-hour simulations were performed to simulate the convective boundary-layer structure, starting at 0700 LT 6 June 1987.In the control simulation, measured surface parameters and atmospheric data were used to simulate observed boundary-layer processes. In the remaining simulations, five surface parameters – soil texture, initial soil moisture, minimum stomatal resistance, leaf area index, and vegetation cover – were varied systematically to study how uncertainty in the specification of these surface parameters affects simulated boundary-layer processes.The simulated uncertainty in the specification of these five surface parameters resulted in a wide range of errors in the prediction of turbulent fluxes, mean thermodynamic structure, and the depth of the ABL. Under certain conditions uncertainty in the specifications of soil texture and minimum stomatal resistance had the greatest influence on the boundary-layer structure. A lesser but still moderately strong effect on the simulated ABL resulted from (1) a small decrease (4%) in the observed initial soil moisture (although a large increase [40%] had only a marginal effect), and (2) a large reduction (66%) in the observed vegetation cover. High uncertainty in the specification of leaf area index had only a marginal impact on the simulated ABL. It was also found that the variations in these five surface parameters had a negligible effect on the simulated horizontal wind fields. On the other hand, these variations had a significant effect on the vertical distribution of turbulent heat fluxes, and on the predicted maximum boundary-layer depth, which varied from about 1400–2300 m across the 11 simulations. Thus, uncertainties in the specification of surface parameters can significantly affect the simulated boundary-layer structure in terms of meteorological and air quality model predictions.     

A coupled surface boundary-layer-quasigeostrophic model

A three-dimensional model of the mesoscale surface boundary layer of the open ocean is developed through scale analysis of the primitive equations with mixing included. A set of surface boundary-layer equations appropriate for a broad range of oceanic and atmospheric scales is thereby derived. The essential basis of the model is a coupling between quasigeostrophic dynamics away from the boundary layer and arbitrary mixing models within the mixed layer. The coupling consists of advection of the boundary layer by the horizontal and vertical components of the interior quasigeostrophic flow and forcing of the interior by the boundary layer in the form of divergence within the boundary layer which leads to vortex stretching/compression in the interior. The divergence is generalized for mesoscale wind-driven flows and includes nonlinear interaction between the directly wind-driven boundary-layer flow and the interior flow in the form of interior relative vorticity advection by the wind-driven flow. The nature of the equations leads us to apply a numerical algorithm to their solution. This algorithm is calibrated through application to idealized problems to determine the temporal and spatial grid requirements. The model is initialized with a realistic ocean flow having the properties of the Gulf Stream.     

Sea-Breeze Modification Of The Growth Of A Marine Internal Boundary Layer

A numerical mesoscale model is used to make a high-resolutionsimulation of the marine boundary layer in the Persian Gulf, duringconditions of offshore flow from Saudi Arabia. A marine internal boundary layer(MIBL) and a sea-breeze circulation (SBC) are found to co-exist. The sea breeze develops in the mid-afternoon, at which time its frontis displaced several tens of kilometres offshore. Between the coastand the sea-breeze system, the MIBL that occurs is consistent with a picture described in the existing literature. However, the MIBL isperturbed by the SBC, the boundary layer deepening significantly seaward of the sea-breeze front. Our analysis suggests that thisstrong, localized deepening is not a direct consequence offrontal uplift, but rather that the immediate cause is the retardation of theprevailing, low-level offshore windby the SBC. The simulated boundary-layer development can be accounted for by using a simple 1D Lagrangian model of growth driven by the surface heatflux. This model is obtained as a straightforward modification ofan established MIBL analytic growth model.     

A moving grid finite-element model of the bulk properties of the atmospheric boundary layer

A moving-grid finite-element model has been developed to model numerically the vertically integrated properties of the atmospheric boundary layer (ABL) in one dimension. The model equations for mean wind velocity and potential temperature are combined with a surface energy budget and predictive equations for boundary-layer height to simulate both stable and unstable ABLs. The nodal position defining the top of the boundary layer is one of the model unknowns and is determined by boundary-layer dynamics. The finite-element method, being an integral method, has advantages of accurate representation of both bulk values and their vertical derivatives, the latter being essential properties of the nocturnal boundary layer. Compared with observations and results of other models, the present model predicts bulk properties very well while retaining a simple and economical form.Journal Paper No. J-12996 of the Iowa Agriculture and Home Economics Experiment Station, Ames, Iowa, Project No. 2779.     

A model of atmospheric boundary-layer flow above an isolated two-dimensional ‘hill’; an example of flow above ‘gentle topography’

A numerical model is developed for two-dimensional turbulent boundary-layer flow above gentle topography — defined as not giving rise to mean flow separation. Although the model is formulated in a framework of mixing length and turbulent energy equation models for the surface layer of the atmospheric boundary layer, it could be modified to include higher-order closure hypotheses and/or extended to model gentle topography for the planetary boundary layer or on the sea bed. Results are presented for flow above a specific shape of hill and the effects of surface roughness and hill height are investigated.     

土壤-大气系统中蒸发过程的数值模拟

本文利用一维土壤-大气耦合数值模式,对土面蒸发过程进行了模拟研究。模式较好反映了蒸发过程三阶段的趋势特征,并试验了不同土壤性质、太阳辐射强度和地转风速对蒸发过程影响。模式还揭示了蒸发过程中,土面热量平衡各分量相互转换过程。