Turbulent kinetic energy budgets from a large-eddy simulation of airflow above and within a forest canopy

The output of a large-eddy simulation was used to study the terms ofthe turbulent kinetic energy (TKE) budget for the air layers above andwithin a forest. The computation created a three-dimensional,time-dependent simulation of the airflow, in which the lowest third ofthe domain was occupied by drag elements and heat sources to representthe forest. Shear production was a principal source of TKE in theupper canopy, diminishing gradually above tree-top height and moresharply with depth in the canopy. The transfer of energy to subgridscales (dissipation) was the main sink in the upper part of the domainbut diminished rapidly with depth in the canopy. Removal ofresolved-scale TKE due to canopy drag was extremely important,occurring primarily in the upper half of the forest where the foliagedensity was large. Turbulent transport showed a loss at the canopytop and a gain within the canopy. These general features have beenfound elsewhere but uncertainty remains concerning the effects ofpressure transport. In the present work, pressure was calculateddirectly, allowing us to compute the pressure diffusion term. Wellabove the canopy, pressure transport was smaller than, and opposite insign to, the turbulent transport term. Near the canopy top andbelow, pressure transport acted in concert with turbulent transport toexport TKE from the region immediately above and within the uppercrown, and to provide turbulent energy for the lower parts of theforest. In combination, the transport terms accounted for over half ofthe TKE loss near the canopy top, and in the lowest two-thirds of thecanopy the transport terms were the dominant source terms in thebudget. Moreover, the pressure transport was the largest source ofturbulent kinetic energy in the lowest levels of the canopy, beingparticularly strong under convective conditions. These resultsindicate that pressure transport is important in the plant canopyturbulent kinetic energy budget, especially in the lowest portion ofthe stand, where it acts as the major driving force for turbulentmotions.     

A subgrid-scale model for large-eddy simulation of planetary boundary-layer flows

A long-standing problem in large-eddy simulations (LES) of the planetary boundary layer (PBL) is that the mean wind and temperature profiles differ from the Monin-Obukhov similarity forms in the surface layer. This shortcoming of LES has been attributed to poor grid resolution and inadequate sub-grid-scale (SGS) modeling. We study this deficiency in PBL LES solutions calculated over a range of shear and buoyancy forcing conditions. The discrepancy from similarity forms becomes larger with increasing shear and smaller buoyancy forcing, and persists even with substantial horizontal grid refinement. With strong buoyancy forcing, however, the error is negligible.In order to achieve better agreement between LES and similarity forms in the surface layer, a two-part SGS eddy-viscosity model is proposed. The model preserves the usual SGS turbulent kinetic energy formulation for the SGS eddy viscosity, but it explicitly includes a contribution from the mean flow and a reduction of the contributions from the turbulent fluctuations near the surface. Solutions with the new model yield increased fluctuation amplitudes near the surface and better correspondence with similarity forms out to a distance of 0.1–0.2 times the PBL depth, i.e., a typical surface-layer depth. These results are also found to be independent of grid anisotropy. The new model is simple to implement and computationally inexpensive.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Characteristics of secondary circulations over an inhomogeneous surface simulated with large-eddy simulation

Large-eddy simulation is used to study secondary circulations in the convective boundary layer modulated as a result of horizontally varying surface properties and surface heat fluxes over flat terrain. The presence of heat flux heterogeneity and its alignment with respect to geostrophic wind influences the formation, strength and orientation of organized thermals. Results show boundary-attached roll formation along heat flux maxima in the streamwise direction. The streamwise organization of the updrafts and downdrafts formed downwind of heterogeneities leads to counter-rotating secondary circulations in the crosswind plane. The distribution of resolved-scale pressure deviations shows large pressure gradients in the crosswind plane. Spanwise and vertical velocity variances and heat flux profiles depict considerable spatial variability compared to a homogeneous forest simulation. Secondary circulations are observed for various ambient wind scenarios parallel and perpendicular to heterogeneities. In the presence of increased wind speed, thermals emerging from the heat flux heterogeneity are elongated, and organize along and downwind of large-scale heterogeneity in the streamwise direction. Simulation with a reduced heat flux shows a shallower circulation with a lower aspect ratio. Point measurements of heat flux inside the roll circulation could be overestimated by up to 15–25% compared to a homogeneous case.     

A grid nesting method for large-eddy simulation of planetary boundary-layer flows

A method for performing nested grid calculations with a large-eddy simulation code is described. A common numerical method is used for all meshes, and the grid architecture consists of a single outer or coarse grid, and nested or fine grids, which overlap in some common region. Inter-grid communication matches the velocity, pressure and potential temperature fields in the overlap region. Resolved and sub-grid scale (SGS) turbulent fluxes and kinetic energy on the fine grid are averaged to the coarse grid using a conservation rule equivalent to Germano's identity used to develop dynamic SGS models.Simulations of a slightly convective, strong shear planetary boundary layer were carried out with varying surface-layer resolutions. Grid refinements in the (x, y, z) directions of up to (5, 5, 2) times were employed. Two-way interaction solutions on the coarse and fine meshes are successfully matched in the overlap region on an instantaneous basis, and the turbulent motions on the fine grid blend smoothly into the coarse grid across the grid interface. With surface-layer grid nesting, significant increases in resolved eddy fluxes and variances are found. The energy-scale content of the vertical velocity, and hence vertical turbulent fluxes, appear to be most influenced by increased grid resolution. Vertical velocity spectra show that the dominant scale shifts towards higher wavenumbers (smaller scales) and the magnitude of the peak energy is increased by more than a factor of 3 with finer resolution. Outside of the nested region the average heat and momentum fluxes and spectra are slightly influenced by the fine resolution in the surface layer. From these results we conclude that fine resolution is required to resolve the details of the turbulent motions in the surface layer. At the same time, however, increased resolution in the surface layer does not appreciably alter the ensemble statistics of the resolved and SGS motions outside of the nested region.     

A new large-eddy simulation model for simulating air flow and warm clouds above highly complex terrain. Part I: The dry model

This paper presents the dry version of a new large-eddy simulation (LES) model, which is designed to simulate air flow and clouds above highly complex terrain. The model is three-dimensional and nonhydrostatic, and the governing equations are sound filtered by use of the anelastic approximation. A fractional step method is applied to solve the equations on a staggered Cartesian grid. Arbitrarily steep and complex orography can be accounted for through the method of viscous topography. The dynamical model core is validated by comparing the results for a spreading density current against a benchmark solution. The model accuracy is further assessed through the simulation of turbulent flow across a quasi two-dimensional ridge. The results are compared with wind-tunnel data. The method of viscous topography is not restricted to moderately sloped terrain. Compared to models using curvilinear grids, it allows this model to be applied to a much wider range of flows. This is illustrated through the simulation of an atmospheric boundary-layer flow over a surface mounted cube. The results show that the dry model version is able to accurately represent the complex flow in the vicinity of three-dimensional obstacles. It is concluded that the method of viscous topography was successfully implemented into a micrometeorological LES model. As will be shown in Part II, this allows the detailed study of clouds in highly complex terrain.     

Ecosystem stewardship: A resilience framework for arctic conservation

Ecosystem stewardship is a framework for actively shaping trajectories of ecological and social change to foster a more sustainable future for species, ecosystems, and society. We apply this framework to conservation challenges and opportunities in the Arctic, where the rapid pace of human-induced changes and their interactions force us now to consider a new relationship between people and the rest of nature. Biodiversity, which has traditionally been the target of conservation efforts, is increasingly affected by human impacts such as energy demand and industrial development that are motivated more by short-term profits than by concerns for societal consequences of long-term arctic biodiversity change. We posit that effective approaches to conservation must (a) foster both ecosystem resilience and human wellbeing, (b) integrate ecological and social processes across scales, and (c) take actions that shape the future rather than seeking only to restore the past. To this end, we identify progress through actions that have been or could be taken at local, national, and international scales to promote arctic resilience and conservation. A stewardship approach to conservation aims to prevent undesirable changes and prepares for adaptation to rapid and uncertain changes that cannot be avoided and for transformation to avoid or escape undesirable states. The greatest opportunity for arctic stewardship at the local scale may lie in building upon culturally engrained (often indigenous) respect for nature and reliance on local environment, empowering it through knowledge and power sharing with national regulatory frameworks. This, in turn, allows connection of drivers with impacts across scales and raises awareness of the value of human–environment relationships. At national and international scales stewardship provides rules for coordinated action to reconcile local and regional conservation actions with those that are motivated by constraints at the global level, to foster ecosystem integrity and human wellbeing in the face of transformative changes in environment, landscapes, species, and society.     

Negative water vapour skewness and dry tongues in the convective boundary layer: observations and large-eddy simulation budget analysis

This study focuses on the intrusion of dry air into the convective boundary layer (CBL) originating from the top of the CBL. Aircraft in-situ measurements from the IHOP_2002 field campaign indicate a prevalence of negative skewness of the water vapour distribution within the growing daytime CBL over land. This negative skewness is interpreted according to large-eddy simulations (LES) as the result of descending dry downdrafts originating from above the mixed layer. LES are used to determine the statistical properties of these intrusions: their size and thermodynamical characteristics. A conditional sampling analysis demonstrates their significance in the retrieval of moisture variances and fluxes. The rapid CBL growth explains why greater negative skewness is observed during the growing phase: the large amounts of dry air that are quickly incorporated into the CBL prevent a full homogenisation by turbulent mixing. The boundary-layer warming in this phase also plays a role in the acquisition of negative buoyancy for these dry tongues, and thus possibly explains their kinematics in the lower CBL. Budget analysis helps to identify the processes responsible for the negative skewness. This budget study underlines the main role of turbulent transport, which distributes the skewness produced at the top or the bottom of the CBL into the interior of the CBL. The dry tongues contribute significantly to this turbulent transport.     

Low-frequency variability of the arctic climate: the role of oceanic and atmospheric heat transport variations

Changes in meridional heat transports, carried either by the atmosphere (HTRA) or by the ocean (HTRO), have been proposed to explain the decadal to multidecadal climate variations in the Arctic. On the other hand, model simulations indicate that, at high northern latitudes, variations in HTRA and HTRO are strongly coupled and may even compensate each other. A multi-century control integration with the Max Planck Institute global atmosphere-ocean model is analyzed to investigate the relative role of the HTRO and HTRA variations in shaping the Arctic climate and the consequences of their possible compensation. In the simulation, ocean heat transport anomalies modulate sea ice cover and surface heat fluxes mainly in the Barents Sea/Kara Sea region and the atmosphere responds with a modified pressure field. In response to positive HTRO anomalies there are negative HTRA anomalies associated with an export of relatively warm air southward to Western Siberia and a reduced inflow of heat over Alaska and northern Canada. While the compensation mechanism is prominent in this model, its dominating role is not constant over long time scales. The presence or absence of the compensation is determined mainly by the atmospheric circulation in the Pacific sector of the Arctic where the two leading large-scale atmospheric circulation patterns determine the lateral fluxes with varying contributions. The degree of compensation also determines the heat available to modulate the large-scale Arctic climate. The combined effect of atmospheric and oceanic contributions has to be considered to explain decadal-scale warming or cooling trends.     

Spatial representativeness of single tower measurements and the imbalance problem with eddy-covariance fluxes: results of a large-eddy simulation study

A large-eddy simulation (LES) study is presented that investigates the spatial variability of temporal eddy covariance fluxes and the systematic underestimation of representative fluxes linked to them. It extends a prior numerical study by performing high resolution simulations that allow for virtual measurements down to 20 m in a convective boundary layer, so that conditions for small tower measurement sites can be analysed. It accounts for different convective regimes as the wind speed and the near-surface heat flux are varied. Moreover, it is the first LES imbalance study that extends to the stable boundary layer. It reveals shortcomings of single site measurements and the necessity of using horizontally-distributed observation networks. The imbalances in the convective case are attributed to a locally non-vanishing mean vertical advection due to turbulent organised structures (TOS). The strength of the TOS and thus the imbalance magnitude depends on height, the horizontal mean wind and the convection type. Contrary to the results of a prior study, TOS cannot generally be responsible for large energy imbalances: at low observation heights (corresponding to small towers and near-surface energy balance stations) the TOS related imbalances are generally about one order of magnitude smaller than those in field experiments. However, TOS may cause large imbalances at large towers not only in the case of cellular convection and low wind speeds, as found in the previous study, but also in the case of roll convection at large wind speeds. In the stably stratified boundary layer for all observation heights neither TOS nor significant imbalances are observed. Attempting to reduce imbalances in convective situations by applying the conventional linear detrending method increases the systematic flux underestimation. Thus, a new filter method is proposed.     

对流边界层发展受覆盖逆温影响的大涡模拟研究

本文运用大涡模拟方法研究了存在覆盖逆温时对流边界层顶部的夹卷过程特征,并通过敏感性试验着重分析了在此情况之下的夹卷速度参数化问题,以及风切变的作用。模拟结果表明,初始覆盖逆温改变了夹卷层的结构特征,使得夹卷层结构参数明显增大,风切变使这一效应略有增强,以往关于夹卷速度的参数化方案不再适用。分析研究表明,在有覆盖逆温的情况下,夹卷层结构参数是控制夹卷速度的关键因子,应该将夹卷层结构参数作为变量引入夹卷速度参数化方案,该方案能够很好地预报出受覆盖逆温影响的对流边界层的发展过程。     

Turbulent heat transfer from a sparsely vegetated surface: Two-component representation

The conventional calculation of heat fluxes from a vegetated surface, involving the coefficient of turbulent heat transfer, which increases logarithmically with surface roughness (commonly taken as about 0.12 of the plant height), appears inappropriate for highly structured surfaces such as desertscrub or open forest. An approach is developed here for computing sensible heat flux from sparsely vegetated surfaces, where the absorption of insolation and the transfer of absorbed heat to the atmosphere are calculated separately for the plants and for the soil. This approach is applied to a desert-scrub surface for which the turbulent transfer coefficient of sensible heat flux from the plants is much larger than that from the soil below, as shown by an analysis of plant, soil and air temperatures measured in an animal exclosure in the northern Sinai. The plant density is expressed as the sum of products (plant-height) x (plant-diameter) of plants per unit horizontal surface area (the dimensionless silhouette parameter of Lettau). The solar heat absorbed by the plants is assumed to be transferred immediately to the airflow. The effective turbulent transfer coefficientk _g-eff for sensible heat from the desert-scrub/soil surface computed under this assumption increases sharply with increasing solar zenith angle, as the plants absorb a greater fraction of the incoming irradiation. The surface absorptivity (the co-albedo) also increases sharply with increasing solar zenith angle, and thus the sensible heat flux from such complex surfaces (which include open forests) is a much broader function of time of day than when computed under constantk _g-eff and constant albedo assumptions. The major role that desert-fringe plants play in the genesis of convection and advection cannot be evaluated properly in the conventional calculations....you know not the way of the wind... Ecclesiastes, Ch. XI, verse 5     

不同下垫面湍流输送计算方法的研究

文章应用三种不同下垫面近地面层风速和温度梯度资料，代入风速和温度廓线公式，迭代求得感热通量及热量交换系数。确定Kh /Kw与稳定参数ξ的关系式。分别利用空气动力学方法和波文比能量平衡法计算潜热通量，并进行了比较。     

A climate change simulation starting from 1935

Due to restrictions in the available computing resources and a lack of suitable observational data, transient climate change experiments with global coupled ocean-atmosphere models have been started from an initial state at equilibrium with the present day forcing. The historical development of greenhouse gas forcing from the onset of industrialization until the present has therefore been neglected. Studies with simplified models have shown that this cold start error leads to a serious underestimation of the anthropogenic global warming. In the present study, a 150-year integration has been carried out with a global coupled ocean-atmosphere model starting from the greenhouse gas concentration observed in 1935, i.e., at an early time of industrialization. The model was forced with observed greenhouse gas concentrations up to 1985, and with the equivalent C0₂ concentrations stipulated in Scenario A (Business as Usual) of the Intergovernmental Panel on Climate Change from 1985 to 2085. The early starting date alleviates some of the cold start problems. The global mean near surface temperature change in 2085 is about 0.3 K (ca. 10%) higher in the early industrialization experiment than in an integration with the same model and identical Scenario A greenhouse gas forcing, but with a start date in 1985. Comparisons between the experiments with early and late start dates show considerable differences in the amplitude of the regional climate change patterns, particularly for sea level. The early industrialization experiment can be used to obtain a first estimate of the detection time for a greenhouse-gas-induced near-surface temperature signal. Detection time estimates are obtained using globally and zonally averaged data from the experiment and a long control run, as well as principal component time series describing the evolution of the dominant signal and noise modes. The latter approach yields the earliest detection time (in the decade 1990–2000) for the time-evolving near-surface temperature signal. For global-mean temperatures or for temperatures averaged between 45°N and 45°S, the signal detection times are in the decades 2015–2025 and 2005–2015, respectively. The reduction of the cold start error in the early industrialization experiment makes it possible to separate the near-surface temperature signal from the noise about one decade earlier than in the experiment starting in 1985. We stress that these detection times are only valid in the context of the coupled model's internally-generated natural variability, which possibly underestimates low frequency fluctuations and does not incorporate the variance associated with changes in external forcing factors, such as anthropogenic sulfate aerosols, solar variability or volcanic dust.     

Characteristics of Turbulent Airflow Deduced from Rapid Surface Thermal Fluctuations: An Infrared Surface Anemometer

The intermittent nature of turbulent airflow interacting with the surface is readily observable in fluctuations of the surface temperature resulting from the thermal imprints of eddies sweeping the surface. Rapid infrared thermography has recently been used to quantify characteristics of the near-surface turbulent airflow interacting with the evaporating surfaces. We aim to extend this technique by using single-point rapid infrared measurements to quantify properties of a turbulent flow, including surface exchange processes, with a view towards the development of an infrared surface anemometer. The parameters for the surface-eddy renewal (\(\alpha \) and \(\beta )\) are inferred from infrared measurements of a single-point on the surface of a heat plate placed in a wind tunnel with prescribed wind speeds and constant mean temperatures of the surface. Thermally-deduced parameters are in agreement with values obtained from standard three-dimensional ultrasonic anemometer measurements close to the plate surface (e.g., \(\alpha = 3\) and \(\beta = 1/26~\hbox {(ms)}^{-1}\) for the infrared, and \(\alpha = 3\) and \(\beta = 1/19~\hbox {(ms)}^{-1}\) for the sonic-anemometer measurements). The infrared-based turbulence parameters provide new insights into the role of surface temperature and buoyancy on the inherent characteristics of interacting eddies. The link between the eddy-spectrum shape parameter \(\alpha \) and the infrared window size representing the infrared field of view is investigated. The results resemble the effect of the sampling height above the ground in sonic anemometer measurements, which enables the detection of larger eddies with higher values of \(\alpha \). The physical basis and tests of the proposed method support the potential for remote quantification of the near-surface momentum field, as well as scalar-flux measurements in the immediate vicinity of the surface.     

Evaporation from rain drops on leaves in a cereal canopy: A simulation model

In crop canopies, the persistence of discrete water drops on leaves after rain is of particular importance to the epidemiology of certain foliar pathogens. A model is described which simulates the heat and water vapour fluxes in a plant canopy and includes evaporation from water drops on the leaves. Energy balance equations allow for heat conducted to drops from the adjacent leaf tissue. Preliminary field tests of model performance for winter wheat, which compare predicted and visually assessed leaf wetness, are encouraging.     

The effect of an arctic polynya on the Northern Hemisphere mean circulation and eddy regime: a numerical experiment

The feedback of an arctic polynya, which is a large ice-free zone within the sea ice, on the hemispheric climate is studied with the ECMWF T21 GCM. For this purpose a control and an anomaly integration, in which a polynya was introduced in the Kara Sea, are compared. As the GCM, like the real atmosphere, shows a high level of low frequency variability, the mean response to the changed boundary conditions is obscured by internal noise. The necessary significance analyses are thus performed to enhance the signal-tonoise ratio within the framework of an a priori chosen guess pattern and a multivariate test statistic. The sensible and latent heat fluxes increased above the polynya, which resulted in a warming of the lower troposphere above and near the polynya. No statistically significant local or global sea-level pressure changes are associated with this heating. However we find a significant change of hemispheric extent of the geopotential fields at 300 hPa, if we use as guess patterns the eigenmodes of the barotropic vorticity equation. The different mean flow field is accompanied by significant changes of the synoptic transient eddy field. We find a significant variation in the barotropic and baroclinic forcing of the mean flow by the eddies, a change in the location and intensity of the storm tracks and in the conversion between eddy available and eddy kinetic energy. The additional heat flux from the polynya results in a reduction of the meridional heat flux by the synoptic eddies on the western Atlantic.     

A multi-scale simulation of an extreme downslope windstorm over complex topography

Summary A severe localized windstorm, with near-surface winds > 60 ms⁻¹, occurred in an isolated valley within the Alpine mountains (> 1800 m) of central Norway on 31 January 1995. A multi-scale numerical simulation of the event was performed with the Naval Research Laboratory (NRL)’s Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS), configured with four nested grids telescoping down to 1-km horizontal resolution. The windstorm occurred in response to topographic blocking and deformation of a lower-tropospheric warm front and attendant jet (> 35 ms⁻¹ at 2 km). The key findings are: i) mountain wave resonance and amplification arising from the interaction of the surface-based front and jet with complex orography, ii) sensitivity of the wave response to differential diabatic heating (vertical) gradients above the front, and iii) trapped response within the layer of large frontal stratification in the lower troposphere and subsequent amplification consistent with the theoretically-established two-layer windstorm analogue of Durran (1986). Received September 29, 1999 Revised December 30, 1999