Estimation of the Lagrangian Velocity Structure Function Constant C0 by Large-Eddy Simulation

The inertial subrange Kolmogorov constant C ₀, which determines the effective turbulent diffusion in velocity space, plays an important role in the Lagrangian modelling of pollutants. A wide range of values of the constant are found in the literature, most of them determined at low Reynolds number and/or under different assumptions. Here we estimate the constant C ₀ by tracking an ensemble of Lagrangian particles in a planetary boundary layer simulated with a large-eddy simulation model and analysing the Lagrangian velocity structure function in the inertial subrange. The advantage of this technique is that it easily allows Reynolds numbers to be achieved typical of convective turbulent flows. Our estimates of C ₀ is C ₀=4.3±0.3 consistent with values found in the literature     

From Planetary Boundaries to national fair shares of the global safe operating space — How can the scales be bridged?

The planetary boundaries framework proposes quantitative global limits to the anthropogenic perturbation of crucial Earth system processes, and thus marks out a planetary safe operating space for human activities. Yet, decisions regarding resource use and emissions are mostly made at less aggregated scales, by national and sub-national governments, businesses, and other local actors. To operationalize the planetary boundaries concept, the boundaries need to be translated into and aligned with targets that are relevant at these decision-making scales. In this paper, we develop a framework that addresses the biophysical, socio-economic, and ethical dimensions of bridging across scales, to provide a consistently applicable approach for translating the planetary boundaries into national-level fair shares of Earth’s safe operating space. We discuss our findings in the context of previous studies and their implications for future analyses and policymaking. In this way, we link the planetary boundaries framework to widely-applied operational and policy concepts for more robust strong sustainability decision-making.     

A Parameterisation of the Effective Roughness Length over Inhomogeneous, Flat Terrain

A parameterisation ofthe effective roughness length is presented for an arbitrary givenroughness distribution z₀(x,y) over flat terrainat neutral stratification.Beyond pure averaging, it takes into account the spatial structure of the distribution, especially the influence of length scales, and inflow direction.To allow for interactions between different rough areas, Boussinesq-approximated equations with a turbulence closure of first order are considered and solved using perturbation theory.As a result, the logarithm of the effective roughness length isrepresented as a sum over the product of the Fourier transformation of log z₀ and a so-called dynamic function, which describesthe response of the flow field to a single wavelength of z₀.Although the numerical expenditure is larger than for simple averaging formulae,this method could be used by large-scale models to calculateeffective roughness lengths in every close-to-surface grid cell.     

CFD simulation of airflow over a regular array of cubes. Part II: analysis of spatial average properties

In the first part of this study, results of a computational fluid dynamics simulation over an array of cubes have been validated against a set of wind-tunnel measurements. In Part II, such numerical results are used to investigate spatially-averaged properties of the flow and passive tracer dispersion that are of interest for high resolution urban mesoscale modelling (e.g. non resolved obstacle approaches). The results show that vertical profiles of mean horizontal wind are linear within the canopy and logarithmic above. The drag coefficient, derived from the numerical results using the classical formula for the drag force, is height dependent (it decreases with height). However, a modification of the formula is proposed (accounting for subgrid velocity scales) that makes the drag coefficient constant with height. Results also show that the dispersive fluxes are similar in magnitude to the turbulent fluxes, and that they play a very important role within the canopy. Vertical profiles of turbulent length scales (to be used in k–l closure schemes, where k is the turbulent kinetic energy and l a turbulent length scale) are also derived. Finally the distribution of the values around the mean over the reference volumes are analysed for wind and tracer concentrations.     

Subfilter-Scale Modelling Using Transport Equations: Large-Eddy Simulation of the Moderately Convective Atmospheric Boundary Layer

We perform large-eddy simulation (LES) of a moderately convective atmospheric boundary layer (ABL) using a prognostic subfilter-scale (SFS) model obtained by truncating the full conservation equations for the SFS stresses and fluxes. The truncated conservation equations contain production mechanisms that are absent in eddy-diffusivity closures and, thus, have the potential to better parametrize the SFS stresses and fluxes. To study the performance of the conservation-equation-based SFS closure, we compare LES results from the surface layer with observations from the Horizontal Array Turbulence Study (HATS) experiment. For comparison, we also show LES results obtained using an eddy-diffusivity closure. Following past studies, we plot various statistics versus the non-dimensional parameter, Λ _w /Δ, where Λ _w is the wavelength corresponding to the peak in the vertical velocity spectrum and Δ is the filter width. The LES runs are designed using different domain sizes, filter widths and surface fluxes, in order to replicate partly the conditions in the HATS experiment. Our results show that statistics from the different LES runs collapse reasonably and exhibit clear trends when plotted against Λ _w /Δ. The trends exhibited by the production terms in the modelled SFS conservation equations are qualitatively similar to those seen in the HATS data with the exception of SFS buoyant production, which is underpredicted. The dominant production terms in the modelled SFS stress and flux budgets obtained from LES are found to approach asymptotically constant values at low Λ _w /Δ. For the SFS stress budgets, we show that several of these asymptotes are in good agreement with their corresponding theoretical values in the limit Λ _w /Δ → 0. The modelled SFS conservation equations yield trends in the mean values and fluctuations of the SFS stresses and fluxes that agree better with the HATS data than do those obtained using an eddy-diffusivity closure. They, however, underpredict considerably the level of SFS anisotropy near the wall when compared to observations, which could be a consequence of the shortcomings in the model used for the pressure destruction terms. Finally, we address the computational cost incurred due to the use of additional prognostic equations.     

Velocity profiles,the resistance law and the dissipation rate of mean flow kinetic energy in a neutrally and stably stratified planetary boundary layer

The logarithmic + polynomial approximation is suggested for vertical profiles of velocity components in a planetary boundary layer (PBL) at neutral and stable stratification. The resistance law functions A and B are determined on the basis of this approximation, using integral relations derived from the momentum equations, the Monin-Obukhov asymptotic formula for the wind profile in a stably stratified near-surface layer and the known expressions for the PBL depth. This result gives a realistic and convenient method for calculating the surface friction velocity and direction and the total dissipation rate of mean flow kinetic energy in terms of geostrophic velocity, buoyancy flux at the surface, the roughness parameter and the Coriolis parameter. In the course of these derivations a review is given of current views on the main problems of the neutral and stable PBL.     

Topographic generation of intermediate-scale motions by wind-driven currents — A model for the circulation in the vicinity of the Indian-Antarctic ridge

When a broad ocean current encounters a large-scale topographic feature, standing Rossby wave patterns can be generated. Short Rossby waves with a scale L_i = √ Q/β (Q is the speed of the approaching flow; β is the meridional gradient of f) are generated east of the topography. If the zonal scale of the topography, L, is planetary, long standing Rossby waves can be generated west of the topography, when the current has a meridional component. The long waves focus the disturbance zonally and produce alternating regions of intensified or reduced zonal flow. The meridional scale that characterizes these zonal bands is the intermediate scales, L = L_i^2/3L^1/3. When the meridional topographic scale is comparable to L, the amplitude of the long-wave disturbance is dominant. Using multiple-scale methods to exploit the scale gap between the planetary, intermediate and Rossby wave scales, the topographically induced pressure and velocity fields due to a zonal ridge are obtained. When the planetary-scale flow field is directed poleward, a westward counterflow can occur along the poleward flank of the ridge. The meridional scales of these topographically induced flows are comparable to those observed along the Indian-Antarctic Ridge by Callahan (1971).     

Chemical perturbations in the planetary boundary layer and their relevance for chemistry transport modelling

The role of perturbations of reactive trace gas concentration distributions in turbulent flows in the planetary boundary layer (PBL) is discussed. The paper focuses on disturbances with larger spatial scales. Sequential nesting of a chemical transport model is applied to assess the effect of neglecting subgrid chemical perturbations on the formation and loss of ozone, NO _x , peroxyacetyl nitrate (PAN) and HNO₃ calculated with a highly complex chemical mechanism. The results point to characteristic differences regarding the process of mixing of chemically reactive species in the PBL and lower troposphere.     

2007/2008和2008/2009冬季平流层强、弱极涡事件对应的行星波活动的对比分析

在对逐日气象资料进行纬向谐波分析的基础上, 对比和讨论了2007/2008年冬季强极涡期间和2008/2009冬季弱极涡期间平流层和对流层不同波数的行星波的变化特征, 特别关注强极涡或弱极涡发生之后, 500 hPa 沿60°N和30°N行星波1波和2波振幅和位相的差异, 以及相应的500 hPa位势场的差异, 进而讨论为什么不同的平流层极涡异常会对东亚有不同的影响, 特别讨论为什么同一种极涡异常, 对我国南北方近地面气温的影响会不同。结果表明:平流层极涡发生异常时, 平流层行星波活动有明显的异常。随着极涡异常的下传, 对流层行星波的振幅和位相也有明显的变化, 而且, 对于不同的纬度带, 其变化又有不同, 表现为:2008年1月强极涡发生之后, 500 hPa行星波1波和2波的扰动都向南伸, 而2009年1月的弱极涡(SSW)期间和之后, 1波和2波的扰动都偏北; 在对流层, 强极涡和弱极涡发生之后不但行星波1波和2波的振幅有所差异, 其位相也有明显的不同。特别是, 其位相的差异还随纬度而变化。就同一年(或者说对于同是强极涡或者同是弱极涡)而言, 无论是1波还是2波, 在60°N和30°N附近的扰动相比, 几乎反位相。这样就使得它们的500 hPa 位势场也有明显不同:在东半球, 主要表现为乌拉尔高压和东亚大槽的强度和位置不同。2008年1月强极涡发生之后, 乌拉尔高压和东亚大槽东移, 不利于冷空气向欧亚大陆北部(包括我国北方)的输送, 使这些地区的温度偏高;而2009年1月弱极涡之后, 东亚大槽西退, 利于冷空气向欧亚大陆北部输送, 导致这些地区较冷。对于同一种极涡异常(如2008强极涡或者2009弱极涡)由于南方和北方行星波扰动的位相不同, 对南方和北方冷暖空气的输送也就不一样。所以同一种极涡异常对(我国)南北地区的温度影响是不同的。     

Turbulence in the Stable Boundary Layer at Higher Richardson Numbers

We present some algebraic and numerical simulations of the stable boundary layer. We also discuss the problem of the existence of a critical Richardson number (Ri), beyond which the turbulence is suppressed. We compare the results of a second-order algebraic model with those of a third-order numerical model and, to this purpose, numerical simulations of a wind-tunnel flow, which is characterized by various Richardson numbers, were performed. As far as the second-order model is concerned, solutions, for the Richardson number greater than any critical value, can be obtained by modifying the time scales of the second-order equation pressure correlation terms in order to account for a buoyancy damping factor. We show that using a third-order model allows the same results (no critical Richardson number) to be obtained without modifications to the time scales. It is suggested that the non-locality, accounted for by the third-order moments, could allow the turbulence to persist also for Ri > 1.     

Direct numerical simulation of a vigorously heated low Reynolds number convective boundary layer

Computations of the buoyantly unstable Ekman layer are performed at low Reynolds number. The results are obtained by directly solving the three-dimensional time-dependen Navier-Stokes equations with the Boussinesq buoyancy approximation, resolving all relevant scales of motion (no turbulence closure is needed). The flow is capped by a stable temperature inversion and heated from below at a rate that produces an inversion-height to Obukhov-length ratio −z_i/L_* = 32. Temperature and velocity variance profiles are found to agree well with those from an earlier vigorously heated under-resolved computation at higher Reynolds number, and with experimental data of Deardorff and Willis (Boundary-Layer Meteorol., 32: 205–236, 1985). Significant helicity is found in the layer, and helical convection patterns of the scale of the inversion height are observed.     

On the influence of stratospheric conditions on forced tropospheric waves in a steady‐state primitive equation model

Abstract

The structure of a forced planetary wave is computed by means of a linearized steady‐state primitive equation model on a hemisphere. The vertical velocity in pressure coordinates is specified at the lower boundary to simulate orographie forcing. The vertical finite differences are on equally spaced pressure levels with a moderately high vertical resolution. The upper boundary condition dp/dt = 0 is applied at p =0 in the model. Numerical experiments show that the tropospheric structure of forced planetary waves is sensitive to the stratospheric background conditions in the model.     

Land-Surface Heterogeneity Effects in the Planetary Boundary Layer

We investigate the cumulative added value of assimilating temperature, moisture, and wind observations in the three-dimensional non-hydrostatic Fifth-Generation Pennsylvania State University/National Center for Atmospheric Research Mesoscale Model MM5 and use these forecasts to analyze the relationship between surface forcing and planetary boundary-layer (PBL) depth. A data assimilation methodology focused on the surface and the PBL, previously tested in a one-dimensional version of MM5, is applied to 29 May, 6 June, and 7 June 2002 during the International $\hbox {H}_{2}\hbox {O}$ Project over the Southern Great Plains. Model-predicted PBL depth is evaluated against PBL depth diagnosed from data across 4,800 km of airborne lidar data (flight tracks 100–300 km long). The forecast with data assimilation verifies better against observations and is thus used to investigate the environmental conditions that govern PBL depth. The spatial structure in PBL depth is found to be most affected by spatial variations in surface buoyancy flux and capping inversion strength. The spatial scales of surface flux forcing reflected in the PBL depth are found through Fourier analysis and multiresolution decomposition. Correlations are ${<}0.50$ at scales of 64 km or less and increase at larger scales for 29 May and 6 June, but on 7 June low correlations are found at all scales, possibly due to greater within-PBL wind speeds, a stronger capping inversion on this day, and clouds. The results suggest a minimum scale, a function of wind speed, below which heterogeneity in surface buoyancy fluxes is not reflected directly in PBL depth.     

New Approaches in Two-equation Turbulence Modelling for Atmospheric Applications

The turbulence closure in atmospheric boundary-layer modelling utilizing Reynolds Averaged Navier–Stokes (RANS) equations at mesoscale as well as at local scale is lacking today a common approach. The standard kɛ model, although it has been successful for local scale problems especially in neutral conditions, is deficient for mesoscale flows without modifications. The k–ɛ model is re-examined and a new general approach in developing two-equation turbulence models is proposed with the aim of improving their reliability and consequently their range of applicability. This exercise has led to the replacement of the ɛ-transport equation by the transport equation for the turbulence inverse length scale (wavenumber). The present version of the model is restricted to neutrally stratified flows but applicable to both local scale and mesoscale flows. The model capabilities are demonstrated by application to a series of one-dimensional planetary boundary-layer problems and a two-dimensional flow over a square obstacle. For those applications, the present model gave considerably better results than the standard k–ɛ model.     

Sensitivity of PBL model predictions to model design and uncertainties in environmental inputs

A simple time-dependent one-dimensional model of the planetary boundary layer (PBL) is described and used to examine the degree to which model design decisions affect model output variables. The model's sensitivity to changes in the environmental conditions is also explored. Averages of the surface fluxes, near-ground wind speeds and other PBL properties from 48 h simulations are compared to control runs. The model-calculated surface fluxes are most sensitive, in decreasing order of importance, to the vertical grid spacing, the form of closure between the surface temperature and the atmosphere, the use of vertical diffusivity smoothing, the choice of maximum time step and choice of turbulence closure scheme. These fluxes are relatively insensitive to mixing-length scaling or choice of implicit time step weighting factor. Sensitivity to changes in soil type exceeds any of the design criteria tested. The modeled fluxes are moderately sensitive to small variations in the horizontal pressure gradient, to unsteadiness in the geostrophic wind and to variations in surface roughness. They are relatively insensitive to uncertainties in local vertical velocities and small (25%) variations applied separately to soil thermal diffusivity or heat capacity. The sensitivity of the average PBL depth (Z _i ) to model and environmental changes are similar to those of surface fluxes except thatZ _i is more sensitive to changes in mixing length, albedo and imposed vertical velocity then are the surface fluxes.     

Airborne measurements of turbulent trace gas fluxes and analysis of eddy structure in the convective boundary layer over complex terrain

Using the new high-frequency measurement equipment of the research aircraft DO 128, which is described in detail, turbulent vertical fluxes of ozone and nitric oxide have been calculated from data sampled during the ESCOMPTE program in the south of France. Based on airborne turbulence measurements, radiosonde data and surface energy balance measurements, the convective boundary layer (CBL) is examined under two different aspects. The analysis covers boundary-layer convection with respect to (i) the control of CBL depth by surface heating and synoptic scale influences, and (ii) the structure of convective plumes and their vertical transport of ozone and nitric oxides. The orographic structure of the terrain causes significant differences between planetary boundary layer (PBL) heights, which are found to exceed those of terrain height variations on average. A comparison of boundary-layer flux profiles as well as mean quantities over flat and complex terrain and also under different pollution situations and weather conditions shows relationships between vertical gradients and corresponding turbulent fluxes. Generally, NO_x transports are directed upward independent of the terrain, since primary emission sources are located near the ground. For ozone, negative fluxes are common in the lower CBL in accordance with the deposition of O₃ at the surface.The detailed structure of thermals, which largely carry out vertical transports in the boundary layer, are examined with a conditional sampling technique. Updrafts mostly contain warm, moist and NO_x loaded air, while the ozone transport by thermals alternates with the background ozone gradient. Evidence for handover processes of trace gases to the free atmosphere can be found in the case of existing gradients across the boundary-layer top. An analysis of the size of eddies suggests the possibility of some influence of the heterogeneous terrain in mountainous area on the length scales of eddies.     

Momentum Transfer and Turbulent Kinetic Energy Budgets within a Dense Model Canopy

Second-order closure models for the canopy sublayer (CSL) employ aset of closure schemes developed for `free-air' flow equations andthen add extra terms to account for canopy related processes. Muchof the current research thrust in CSL closure has focused on thesecanopy modifications. Instead of offering new closure formulationshere, we propose a new mixing length model that accounts for basicenergetic modes within the CSL. Detailed flume experiments withcylindrical rods in dense arrays to represent a rigid canopy areconducted to test the closure model. We show that when this lengthscale model is combined with standard second-order closureschemes, first and second moments, triple velocity correlations,the mean turbulent kinetic energy dissipation rate, and the wakeproduction are all well reproduced within the CSL provided thedrag coefficient (C_D) is well parameterized. The maintheoretical novelty here is the analytical linkage betweengradient-diffusion closure schemes for the triple velocitycorrelation and non-local momentum transfer via cumulant expansionmethods. We showed that second-order closure models reproducereasonably well the relative importance of ejections and sweeps onmomentum transfer despite their local closure approximations.Hence, it is demonstrated that for simple canopy morphology (e.g.,cylindrical rods) with well-defined length scales, standard closureschemes can reproduce key flow statistics without much revision.When all these results are taken together, it appears that thepredictive skills of second-order closure models are not limitedby closure formulations; rather, they are limited by our abilityto independently connect the drag coefficient and the effectivemixing length to the canopy roughness density. With rapidadvancements in laser altimetry, the canopy roughness densitydistribution will become available for many terrestrialecosystems. Quantifying the sheltering effect, the homogeneity andisotropy of the drag coefficient, and more importantly, thecanonical mixing length, for such variable roughness density isstill lacking.     

Numerical Modelling of the Turbulent Flow Developing Within and Over a 3-D Building Array, Part I: A High-Resolution Reynolds-Averaged Navier—Stokes Approach

A study of the neutrally-stratified flow within and over an array of three-dimensional buildings (cubes) was undertaken using simple Reynolds-averaged Navier—Stokes (RANS) flow models. These models consist of a general solution of the ensemble-averaged, steady-state, three-dimensional Navier—Stokes equations, where the k-ε turbulence model (k is turbulence kinetic energy and ε is viscous dissipation rate) has been used to close the system of equations. Two turbulence closure models were tested, namely, the standard and Kato—Launder k-ε models. The latter model is a modified k-ε model designed specifically to overcome the stagnation point anomaly in flows past a bluff body where the standard k-ε model overpredicts the production of turbulence kinetic energy near the stagnation point. Results of a detailed comparison between a wind-tunnel experiment and the RANS flow model predictions are presented. More specifically, vertical profiles of the predicted mean streamwise velocity, mean vertical velocity, and turbulence kinetic energy at a number of streamwise locations that extend from the impingement zone upstream of the array, through the array interior, to the exit region downstream of the array are presented and compared to those measured in the wind-tunnel experiment. Generally, the numerical predictions show good agreement for the mean flow velocities. The turbulence kinetic energy was underestimated by the two different closure models. After validation, the results of the high-resolution RANS flow model predictions were used to diagnose the dispersive stress, within and above the building array. The importance of dispersive stresses, which arise from point-to-point variations in the mean flow field, relative to the spatially-averaged Reynolds stresses are assessed for the building array.     

The exchange of fluid mass between quasi-geostrophic and ageostrophic motions during the reflection of Rossby waves from a coast. I. The case of an infinite rectilinear coast

When the problem of the reflection of spatially localized Rossby waves from a coast is treated using the quasigeostrophic (QG) approximation, the total fluid mass and the along-shore circulation calculated from the geostrophic height field are not conserved. To understand the correct mass balance and the degree to which the QG equations and boundary conditions may be in error, we analyze an initial-value problem for the Laplace tidal equations on a β-plane in the asymptotic limit 1, where is the ratio of the spatial scale of the motion to the Earth's radius.It is shown that there is a coupling between QG and O() fields. Physically, the coupling occurs by a peculiar adjustment process in the O() approximation in which fast gravity waves are permanently generated to build up a quasi-stationary edge Kelvin wave. Different temporal scales (large for O(1) Rossby waves and small for the O() gravity waves make comparable the contributions of the waves to the mass and circulation balance equations. However, QG analysis itself describes the reflection of Rossby waves correctly, but is incomplete, and for satisfactory balances one has to take into account the fields of both orders of the approximation.Applications of the results to closed basins, baroclinicity, and variable bottom topography are discussed. It is conjectured that an interaction of strong oceanic eddies with a coast (continental slope) may give rise to noticeable along-shore jet currents.     

The structure and propagation of stationary planetary wave packet in the barotropic atmosphere

Monthly or seasonally mean anomalies of large-scale atmospheric circulation are better represented by wave packets or their combination. Both qualitative and quantitative analyses of equations of wave packet dynamics, which are obtained by the use of WKB approximation, are very helpful for the understanding of structure, formation and propagation of stationary and quasi-stationary planetary wave packet patterns in the atmosphere. Indeed, these equations of wave packet dynamics can be directly solved by the method of characteristic lines, and the results can be simply and clearly interpreted by physical laws. In this paper, a quasi-geostrophic barotropic model is taken for simplicity, and the wave packets superimposed on several ideal profiles of the basic current and excited by some ideal forcings are investigated in order to make comparison of the accuracy of calculation with the analytical solution. It is revealed that (a) the rays of stationary planetary wave packet do not coincide with but go away from the great circle with significant difference if the shear of the basic zonal flow is not too small; (b) being superimposed on a westerly jet flow with positive shear (Uλ/y>0), the stationary wave packets excited by low-latitudinal forcing are first intensified during their northeastward propagation in the Northern Hemisphere, then reach their maximum of amplitude at some critical latitude, and after that weaken again; (c) the connected line of extremes (the positive and negative centres) of wave packet does not coincide with but crosses the ray by an angle, the larger the scale of external forcing, the larger the angle; and (d) the whole pattern of a trapped stationary wave packet is complicated by the interference between the incident and reflected waves.