A fluid experiment of large-scale topography effect on baroclinic wave flows

The effects of topography on baroclinic wave flows are studied experimentally in a thermally driven rotating annulus of fluid.Fourier analysis and complex principal component (CPC) analysis of the experimental data show that, due to topographic forcing, the flow is bimodal rather than a single mode. Under suitable imposed experimental parameters, near thermal Rossby number ROT = 0.1 and Taylor number Ta = 2.2 × 107, the large-scale topography produces low-frequency oscillation in the flow and rather long-lived flow pattern resembling blocking in the atmospheric cir-culation. The ‘blocking’ phenomenon is caused by the resonance of travelling waves and the quasi-stationary waves forced by topography.The large-scale topography transforms wavenumber-homogeneous flows into wavenumber-dispersed flows, and the dispersed flows possess lower wavenumbers.     

Semi-Implicit Scheme for the DWD Lokal-Modell

Summary  The fourth generation of numerical weather prediction (NWP) models is currently under development at the Deutscher Wetterdienst (DWD) consisting of a global grid point model (GME) and limited-area Lokal-Modell (LM). The nonhydrostatic fully compressible LM has been designed to meet high-resol ution regional fore-cast requirements at meso-β and meso-γ scales. The initial LM implementation is based on the NCAR/Penn State MM5 with the addition of a novel generalized terrain-following coordinate and rotated lat-lon grid. A fully 3D semi-implicit time-stepping scheme has been implemented by retaining the full buoyancy term instead of the approximate form found in MM5. In contrast with earlier schemes, mass-lumping is not applied to simplify the elliptic operator on an Arakawa-C/Lorenz grid. The resulting variable-coefficient elliptic problem is solved using a minimal residual Krylov iterative method with line relaxation preconditione rs. The new semi-implicit scheme is compared with a variant of the Klemp–Wilhelmson split-explicit scheme (horizontal explicit, vertical implicit) on the basis of computational efficiency and accuracy at resolutions ranging from 7 km to 400 m. Both idealized 3D mountain wave flows and naturally occuring flows are analyzed. Below the tropopause, the 3D semi-implicit scheme can be more efficient for low Mach number M ≪ 1 flows when the number of small time steps Δt _s of the split-explicit approach increases with the sound-speed Courant number. Revised October 7, 1999     

The Theoretical Model of Atmospheric Turbulence Spectrum in Surface Layer

It is shown that the slope of energy spectrum obtained from the velocity solution of Kdv-Burgers equation lies between -5/3 and -2 in the dilogarithmic coordinates paper. The spectrum is very close to one of Kolmogorov’s isotropic turbulence and Frisch’s intermittent turbulence in inertial region. In this paper, the Kdv-Burgers equation to describe atmospheric boundary layer turbulence is obtained. In the equation, the 1 / Re, corresponds to dissipative coefficient v, to dispersive coefficient β, then (v/ 2β)2 corresponds to .We prove that the wave number corresponding to maximum energy spectrum decreases with the decrease of stability (i.e., the increase of in eddy-containing region. And the spectrim amplitude decreases with the increase of (i.e., the decrease of stability). These results are consistent with actual turbulence spectrum of atmospheric surface layer from turbulence data.     

Local Similarity in the Stable Boundary Layer and Mixing-Length Approaches: Consistency of Concepts

In stably stratified flows vertical movement of eddies is limited by the fact that kinetic energy is converted into potential energy, leading to a buoyancy displacement scale z _B . Our new mixing-length concept for turbulent transport in the stable boundary layer follows a rigid-wall analogy, in the sense that we assume that the buoyancy length scale is similar to neutral length scaling. This implies that the buoyancy length scale is: ℓ _B  = κ _B z _B , with κ _B ≈ κ, the von Karman constant. With this concept it is shown that the physical relevance of the local scaling parameter z/Λ naturally appears, and that the α coefficient of the log-linear similarity functions is equal to c/κ ², where c is a constant close to unity. The predicted value α ≈ 1/κ ² = 6.25 lies within the range found in observational studies. Finally, it is shown that the traditionally used inverse linear interpolation between the mixing length in the neutral and buoyancy limits is inconsistent with the classical log-linear stability functions. As an alternative, a log-linear consistent interpolation method is proposed.     

Surface and atmospheric controls on the heating coefficient

Measurements of the surface radiation budgets for three surfaces—grass, soil and a cornfield—are used to evaluate the ‘heating coefficient’β, and its componentsβ _↑(=dL _↑/dR _n ) andβ _↓(=dL _↓/dR _n ). This resolution permits an analysis of the sensitivity of β to surface and atmospheric influences.β _↑ is shown, both theoretically and empirically, to be determined by surface properties. For grass and soil, the parameter functions as an index of surface desiccation.β _↓ values are large (even under clear conditions) and variable, accounting for part of the variance in β and the anomalously small and negative values reported in the literature.β _↓ values for cloudy conditions may be larger or smaller than those for clear skies. It is concluded that, unless a predictive procedure can be developed forβ _↓, the Monteith and Szeicz model is of limited use for the routine estimation of net radiation.     

Turbulent Flow at 190 m Height Above London During 2006–2008: A Climatology and the Applicability of Similarity Theory

Flow and turbulence above urban terrain is more complex than above rural terrain, due to the different momentum and heat transfer characteristics that are affected by the presence of buildings (e.g. pressure variations around buildings). The applicability of similarity theory (as developed over rural terrain) is tested using observations of flow from a sonic anemometer located at 190.3 m height in London, U.K. using about 6500 h of data. Turbulence statistics—dimensionless wind speed and temperature, standard deviations and correlation coefficients for momentum and heat transfer—were analysed in three ways. First, turbulence statistics were plotted as a function only of a local stability parameter z/Λ (where Λ is the local Obukhov length and z is the height above ground); the σ _i /u _* values (i = u, v, w) for neutral conditions are 2.3, 1.85 and 1.35 respectively, similar to canonical values. Second, analysis of urban mixed-layer formulations during daytime convective conditions over London was undertaken, showing that atmospheric turbulence at high altitude over large cities might not behave dissimilarly from that over rural terrain. Third, correlation coefficients for heat and momentum were analyzed with respect to local stability. The results give confidence in using the framework of local similarity for turbulence measured over London, and perhaps other cities. However, the following caveats for our data are worth noting: (i) the terrain is reasonably flat, (ii) building heights vary little over a large area, and (iii) the sensor height is above the mean roughness sublayer depth.     

Atmospheric Stability Effects on Wind Fields and Scalar Mixing Within and Just Above a Subalpine Forest in Sloping Terrain

Air temperature T _a , specific humidity q, CO₂ mole fraction χ _c , and three-dimensional winds were measured in mountainous terrain from five tall towers within a 1 km region encompassing a wide range of canopy densities. The measurements were sorted by a bulk Richardson number Ri _b . For stable conditions, we found vertical scalar differences developed over a “transition” region between 0.05 < Ri _b < 0.5. For strongly stable conditions (Ri _b > 1), the vertical scalar differences reached a maximum and remained fairly constant with increasing stability. The relationships q and χ _c have with Ri _b are explained by considering their sources and sinks. For winds, the strong momentum absorption in the upper canopy allows the canopy sublayer to be influenced by pressure gradient forces and terrain effects that lead to complex subcanopy flow patterns. At the dense-canopy sites, soil respiration coupled with wind-sheltering resulted in CO₂ near the ground being 5–7 μmol mol⁻¹ larger than aloft, even with strong above-canopy winds (near-neutral conditions). We found Ri _b -binning to be a useful tool for evaluating vertical scalar mixing; however, additional information (e.g., pressure gradients, detailed vegetation/topography, etc.) is needed to fully explain the subcanopy wind patterns. Implications of our results for CO₂ advection over heterogenous, complex terrain are discussed.     

Consistent Two-Equation Closure Modelling for Atmospheric Research: Buoyancy and Vegetation Implementations

A self-consistent two-equation closure treating buoyancy and plant drag effects has been developed, through consideration of the behaviour of the supplementary equation for the length-scale-determining variable in homogeneous turbulent flow. Being consistent with the canonical flow regimes of grid turbulence and wall-bounded flow, the closure is also valid for homogeneous shear flows commonly observed inside tall vegetative canopies and in non-neutral atmospheric conditions. Here we examine the most often used two-equation models, namely and E − ω (where is the dissipation rate of turbulent kinetic energy, E, and is the specific dissipation), comparing the suggested buoyancy-modified closure against Monin–Obukhov similarity theory. Assessment of the closure implementing both buoyancy and plant drag together has been done, comparing the results of the two models against each other. It has been found that the E − ω model gives a better reproduction of complex atmospheric boundary-layer flows, including less sensitivity to numerical artefacts, than does the model. Re-derivation of the equation from the ω equation, however, leads to the model implementation that produces results identical to the E − ω model. Overall, numerical results show that the closure performs well, opening new possibilities for application of such models to tasks related to the atmospheric boundary layer—where it is important to adequately account for the influences of both vegetation and atmospheric stability.     

A Study in Search of Interconnection between Surface Parameters and Surrounding Synoptic And Subsynoptic Features

The paper reveals that the variations in parameters like u*, the scaling velocity and θ*. The scaling tempera-ture during the various phases of monsoon might be linked with subsynoptic features. The rise in u* is mainly connected with the presence of lower tropospheric cyclonic vorticity over a subsynoptic scale of the site. However the variations in θ* is mainly linked with the various phases of monsoon and θ* shows a sharp rise in presence of low level convective cloud.Besides the correlation studies of u and u*, θv and θv* , θv-θv0 and θv* are undertaken. The correlation be?tween θv and θv* is poor. In other two cases correlations are good. Besides u/u* , has shown good coefficient of variation values within the ζ range.     

A K-Theory of Dispersion,Settling and Deposition in the Atmospheric Surface Layer

Dispersion estimates with a Gaussian plume model are often incorrect because of particle settling (β), deposition (γ) or the vertical gradient in diffusivity (K _v (z) = K ₀ + μz). These “non-Gaussian” effects, and the interaction between them, can be evaluated with a new Hankel/Fourier method. Due to the deepening of the plume downwind and reduced vertical concentration gradients, these effects become more important at greater distance from the source. They dominate when distance from the source exceeds L _β = K ₀ U/β ², L _γ  = K ₀ U/γ ² and L _μ = K ₀ U/μ ² respectively. In this case, the ratio β/μ plays a central role and when β/μ = 1/2 the effects of settling and K gradient exactly cancel. A general computational method and several specific closed form solutions are given, including a new dispersion relation for the case when all three non-Gaussian effects are strong. A more general result is that surface concentration scales as C(x) ~ γ ⁻² whenever deposition is strong. Categorization of dispersion problems using β/μ, L _γ and L _μ is proposed.     

On the Time Evolution of the Turbulent Kinetic Energy Spectrum for Decaying Turbulence in the Convective Boundary Layer

Our focus is the time evolution of the turbulent kinetic energy for decaying turbulence in the convective boundary layer. The theoretical model with buoyancy and inertial transfer terms has been extended by a source term due to mechanical energy and validated against large-eddy simulation data. The mechanical effects in a boundary layer of height z _i at a convective surface-layer height z = 0.05z _i are significant in the time evolution of the vertical component of the spectrum, i.e. they enhance the decay time scale by more than an order of magnitude. Our findings suggest that shear effects seem to feedback to eddies with smaller wavenumbers, preserving the original shape of the spectrum, and preventing the spectrum from shifting towards shorter wavelengths. This occurs in the case where thermal effects only are considered.     

Atmospheric turbidity of urban and desert areas of the Nile Basin in the aftermath of Mt. Pinatubo’s eruption

Summary  The Linke T _L , ?ngstr?m β and Unsworth-Monteith δ _a turbidity parameters are investigated for two sites in Egypt: Cairo, a densely populated urban area, and Aswan, an arid unpolluted area. These three turbidity parameters are calculated from broadband pyrheliometric measurements recorded hourly over the period 1992–96. Monthly averages of T _L , β and δ _a show relatively flat and identical seasonal variations with a marked main maxima during spring at both sites, due to Khamsin depressions coming from the Great Sahara. A secondary maximum is observed at Aswan in summer, due to dust haze which prevails during that season, and at Cairo in autumn, due to the northern extension of the Sudan monsoon trough, which is accompanied by small scale depressions with dust particles. Annual mean values of T _L , β and δ _a (5.59, 0.250 and 0.372, respectively) at Cairo are larger than at Aswan (3.89, 0.139 and 0.213, respectively). In the same way, the seasonal mean values of T _L , β and δ _a at Cairo are larger than at Aswan. More generally, the monthly and yearly average turbidity values are significantly larger in Cairo than in Aswan for the whole period 1992–96, which is attributable in part to the urbanization/industrialization effect of Cairo. An estimate of the corresponding overburden is obtained by comparison between the present data and older T _L data from 1922–27. It is also shown that turbidity over both sites is largest during 1992, just after the eruption of Mount Pinatubo in 1991. The dependence of β on some meteorological parameters such as wind speed and direction, precipitable water, relative humidity, temperature and visibility, is also analyzed. This reveals in particular that visibility is not a good predictor of turbidity at either site. Conversely, the wind direction and speed have a definite effect on turbidity, and consequently, largest turbidities occur when the wind carries aerosols from the main industrial particle source areas around Cairo. For any season of the year, the average turbidity at the latter site is larger than that at other big cities such as Athens, Rome, and Toronto, but is lower than at Dhahran, Saudi Arabia. Received February 3, 2000 Revised August 25, 2000     

Daily reference evapotranspiration estimates by the Penman-Monteith equation in Southern Italy. Constant vs. variable canopy resistance

Summary ?The performance of the Penman-Monteith (PM) equation to estimate daily reference evapotranspiration (ET_O) was investigated by attributing three distinct features to the canopy resistance (r _c): (i) r _c constant at 70 s m⁻¹ (Allen et al., 1998; FAO Irrigation and Drainage Paper n. 56), (ii) r _c variable as linear function of a critical resistance r _c, depending on weather variables and empirical parameters relating r _c to r ^* (Katerji and Perrier, 1983; Agronomie, 3[6]: 513–521) and (iii) r _c variable as a mechanistic function of weather variables only (Todorovic, 1999; J. Irrig. Drainage Eng., ASCE, 125[5]: 235–245). Daily weather and grass lysimeter data, measured for a period of seven years at Policoro (Southern Italy), were used. The results confirmed the relative robustness of the PM method with constant r _c while better estimates were obtained only when variable r _c was used. The mechanistic approach of Todorovic (1999) provided the best estimates, while the approach of Katerji and Perrier (1983), with empirically derived parameters, has shown to be not conservative enough to be extended to different locations without calibration. Received January 2, 2002; revised October 31, 2002; accepted December 7, 2002     

The radiative influence of aerosol effects on liquid-phase cumulus and stratiform clouds based on sensitivity studies with two climate models

Aerosol effects on warm (liquid-phase) cumulus cloud systems may have a strong radiative influence via suppression of precipitation in convective systems. A consequence of this suppression of precipitation is increased liquid water available for large-scale stratiform clouds, through detrainment, that in turn affect their precipitation efficiency. The nature of this influence on radiation, however, is dependent on both the treatment of convective condensate and the aerosol distribution. Here, we examine these issues with two climate models—CSIRO and GISS, which treat detrained condensate differently. Aerosol–cloud interactions in warm stratiform and cumulus clouds (via cloud droplet formation and autoconversion) are treated similarly in both models. The influence of aerosol–cumulus cloud interactions on precipitation and radiation are examined via simulations with present-day and pre-industrial aerosol emissions. Sensitivity tests are also conducted to examine changes to climate due to changes in cumulus cloud droplet number (N _c); the main connection between aerosols and cumulus cloud microphysics. Results indicate that the CSIRO GCM is quite sensitive to changes in aerosol concentrations such that an increase in aerosols increases N _c, cloud cover, total liquid water path (LWP) and reduces total precipitation and net cloud radiative forcings. On the other hand, the radiative fluxes in the GISS GCM appear to have minimal changes despite an increase in aerosols and N _c. These differences between the two models—reduced total LWP in the GISS GCM for increased aerosols, opposite to that seen in CSIRO—appear to be more sensitive to the detrainment of convective condensate, rather than to changes in N _c. If aerosols suppress convective precipitation as noted in some observationally based studies (but not currently treated in most climate models), the consequence of this change in LWP suggests that: (1) the aerosol indirect effect (calculated as changes to net cloud radiative forcing from anthropogenic aerosols) may be higher than previously calculated or (2) lower than previously calculated. Observational constrains on these results are difficult to obtain and hence, until realistic cumulus-scale updrafts are implemented in models, the logic of detraining non-precipitating condensate at appropriate levels based on updrafts and its effects on radiation, will remain an uncertainty.     

Large-Eddy Simulation of Flow and Pollutant Transports in and Above Two-Dimensional Idealized Street Canyons

A large-eddy simulation (LES) model, using the one-equation subgrid-scale (SGS) parametrization, was developed to study the flow and pollutant transport in and above urban street canyons. Three identical two-dimensional (2D) street canyons of unity aspect ratio, each consisting of a ground-level area source of constant pollutant concentration, are evenly aligned in a cross-flow in the streamwise direction x. The flow falls into the skimming flow regime. A larger computational domain is adopted to accurately resolve the turbulence above roof level and its influence on the flow characteristics in the street canyons. The LES calculated statistics of wind and pollutant transports agree well with other field, laboratory and modelling results available in the literature. The maximum wind velocity standard deviations σ _i in the streamwise (σ _u ), spanwise (σ _v ) and vertical (σ _w ) directions are located near the roof-level windward corners. Moreover, a second σ _w peak is found at z ≈ 1.5h (h is the building height) over the street canyons. Normalizing σ _i by the local friction velocity u _*, it is found that σ _u /u _* ≈ 1.8, σ _v /u _* ≈ 1.3 and σ _w /u _* ≈ 1.25 exhibiting rather uniform values in the urban roughness sublayer. Quadrant analysis of the vertical momentum flux u′′w′′ shows that, while the inward and outward interactions are small, the sweeps and ejections dominate the momentum transport over the street canyons. In the x direction, the two-point correlations of velocity R _v,x and R _w,x drop to zero at a separation larger than h but R _u,x (= 0.2) persists even at a separation of half the domain size. Partitioning the convective transfer coefficient Ω _T of pollutant into its removal and re-entry components, an increasing pollutant re-entrainment from 26.3 to 43.3% in the x direction is revealed, suggesting the impact of background pollutant on the air quality in street canyons.     

Convection in PIDCAP Part I: Evaluating LAM convection

Summary One goal of BALTEX is to quantify the energy and water cycle over the Baltic region through LAMs. The convective components of these cycles cannot be explicitly modelled and thus are parameterized. In this two-part paper we consider the total convective heat flux c as key sub-gridscale quantity. We use it, in the present Part I, for a model intercomparison. We have four European LAMs (EM, BM, REMO and HIRLAM) at our disposal and run them, both in the prognostic and in the diagnostic mode, for a 9-day target period (Aug/Sep 1995) located within the convectively active part of PIDCAP. We find in the prognostic mode that the vertical profile of c is in most cases characterized by boundary layer convection. The flux is strongest at the surface ( W/m²) and becomes weaker at higher levels. In three deep convection experiments we select columns with a maximum number of lightning events, a maximum precipitation rate and a maximum c. These events are rare (order of 1 % of all columns) but show free atmosphere fluxes of up to  W/m². This result is supported by the net condensation rate CON which is complementary to c. The differences between the LAMs in terms of c and CON are small but significant. From the analysed gridscale fields of the LAMs we further derive, in the diagnostic mode, profiles of c, CON and other convective parameters and compare them with their forecast equivalents. The comparison is best for the HIRLAM model. This suggests that the physical package of HIRLAM, due to its hybrid convection scheme STRACO, may be a bit more realistic than those of EM, BM and REMO. Conclusion of Part I is that the total convective heat flux c is useful to make the model physics comparable across LAMs. In Part II, the problem how to diagnose c from observed gridscale data will be reviewed. Received January 8, 2001 Revised April 26, 2001     

SHEBA flux–profile relationships in the stable atmospheric boundary layer

Measurements of atmospheric turbulence made during the Surface Heat Budget of the Arctic Ocean Experiment (SHEBA) are used to examine the profile stability functions of momentum, φ _m , and sensible heat, φ _h , in the stably stratified boundary layer over the Arctic pack ice. Turbulent fluxes and mean meteorological data that cover different surface conditions and a wide range of stability conditions were continuously measured and reported hourly at five levels on a 20-m main tower for 11 months. The comprehensive dataset collected during SHEBA allows studying φ _m and φ _h in detail and includes ample data for the very stable case. New parameterizations for φ _m (ζ) and φ _h (ζ) in stable conditions are proposed to describe the SHEBA data; these cover the entire range of the stability parameter ζ = z/L from neutral to very stable conditions, where L is the Obukhov length and z is the measurement height. In the limit of very strong stability, φ _m follows a ζ ^1/3 dependence, whereas φ _h initially increases with increasing ζ, reaches a maximum at ζ ≈ 10, and then tends to level off with increasing ζ. The effects of self-correlation, which occur in plots of φ _m and φ _h versus ζ, are reduced by using an independent bin-averaging method instead of conventional averaging.     

The stability of large—scale horizontal air motion in the non—linear basic zephyr flow under the effect of rossby parameter

The stability of large-scale horizontal motion in the atmosphere is discussed in this paper by using qualitative analysis theory of non-linear ordinary differential equations. Both the non-linear distribution of basic Zephyr flow and the variation of geostrophic vorticity along the latitude (f=f0 + βδy ) are all included in this paper's mathematical model so as to analogue the background field of large-scale horizontal air motion more really in the rotating reference frame of the earth. Some significant results are drawn out from this paper and the conclusions of Li(1986)'s and Wan et al.(1990)'s are extended widely.     

A generalized scaling for convective shear flows

During the last two decades, different scalings for convective boundary layer (CBL) turbulence have been proposed. For the shear-free regime, Deardorff (1970) introduced convective velocity and temperature scales based on the surface potential temperature flux,Q _s , the buoyancy parameter, , and the time-dependent boundary-layer depth,h. Wyngaard (1983) has proposed decomposition of turbulence into two components, bottom-up (b) and top-down (t), the former characterized byQ _s , the latter, by the potential temperature flux due to entrainment,Q _h . Sorbjan (1988) has devised height-dependent velocity and temperature scales for both b- and t-components of turbulence.Incorporating velocity shear, the well known similarity theory of Monin and Obukhov (1954) has been developed for the atmospheric surface layer. Zilitinkevich (1971, 1973) and Betchov and Yaglom (1971) have elaborated this theory with the aid of directional dimensional analysis for a particular case when different statistical moments of turbulence can be alternatively attributed as being of either convective or mechanical origin.In the present paper, we attempt to create a bridge between the two approaches pointed out above. A new scaling is proposed on the basis of, first, decomposition of statistical moments of turbulence into convective (c), mechanical (m) and covariance (c&m) contributions using directional dimensional analysis and, second, decomposition of these contributions into bottom-up and top-down components using height-dependent velocity and temperature scales. In addition to the statistical problem, the scaling suggests a new approach of determination of mean temperature and velocity profiles with the aid of the budget equations for the mean square fluctuations.Notation ATL alternative turbulence layer - CBL convective boundary layer - CML convective and mechanical layer - FCL free convection layer - MTL mechanical turbulence layer