The Entrainment Rate for Buoyant Plumes in a Crossflow

We consider large-eddy simulations (LES) of buoyant plumes from a circular source with initial buoyancy flux F ₀ released into a stratified environment with constant buoyancy frequency N and a uniform crossflow with velocity U. We make a systematic comparison of the LES results with the mathematical theory of plumes in a crossflow. We pay particular attention to the limits [(U)\tilde] << 1{\tilde{U}\ll1} and [(U)\tilde] >> 1{\tilde{U}\gg 1}, where [(U)\tilde]=U/(F₀ N)^1/4{\tilde{U}=U/(F_0 N)^{1/4}}, for which analytical results are possible. For [(U)\tilde] >> 1{\tilde{U}\gg 1}, the LES results show good agreement with the well-known two-thirds law for the rise in height of the plume. Sufficiently far above the source, the centreline vertical velocity of the LES plumes is consistent with the analytical z ^−1/3 and z ^−1/2 scalings for respectively [(U)\tilde] << 1{\tilde{U}\ll 1} and [(U)\tilde] >> 1{\tilde{U}\gg 1}. In the general case, where the entrainment is assumed to be the sum of the contributions from the horizontal and vertical velocity components, we find that the discrepancy between the LES data and numerical solutions of the plume equations is largest for [(U)\tilde]=O(1){\tilde{U}=O(1)}. We propose a modified additive entrainment assumption in which the contributions from the horizontal and vertical velocity components are not equally weighted. We test this against observations of the plume generated by the Buncefield fire in the U.K. in December 2005 and find that the results compare favourably. We also show that the oscillations of the plume as it settles down to its final rise height may be attenuated by the radiation of gravity waves. For [(U)\tilde] << 1{\tilde{U}\ll 1} the oscillations decay rapidly due to the transport of energy away from the plume by gravity waves. For ${\tilde{U}>rsim 1}${\tilde{U}>rsim 1} the gravity waves travel in the same direction and at the same speed as the flow. In this case, the oscillations of the plume do not decay greatly by radiation of gravity waves.     

A numerical diffusion model for continuous releases

The purpose of this study was to develop a diffusion model for a continuous point source which takes into account the increase of wind speed with height, and to compare this model with short-range diffusion experiments. The main problem was to find a good expression for the vertical diffusion coefficient. It turned out that good agreement between theory and experiment could only be obtained by introducing a settling speed W for the tracer combined with a conventional expression for the vertical diffusion (K(Z) = K ₀Z^1−p). An empirical relation was found between K ₀ and τ _vU and between W and bar σ _vU².     

Atmospheric inorganic nitrogen in marine aerosol and precipitation and its deposition to the North and South Pacific Oceans

Aerosol and rain samples were collected between 48°N and 55°S during the KH-08-2 and MR08-06 cruises conducted over the North and South Pacific Ocean in 2008 and 2009, to estimate dry and wet deposition fluxes of atmospheric inorganic nitrogen (N). Inorganic N in aerosols was composed of ~68% NH₄⁺ and ~32% NO₃^– (median values for all data), with ~81% and ~45% of each species being present on fine mode aerosol, respectively. Concentrations of NH₄⁺ and NO₃⁻ in rainwater ranged from 1.7–55 μmol L⁻¹ and 0.16–18 μmol L⁻¹, respectively, accounting for ~87% by NH₄⁺ and ~13% by NO₃⁻ of total inorganic N (median values for all data). A significant correlation (r = 0.74, p < 0.05, n = 10) between NH₄⁺ and methanesulfonic acid (MSA) was found in rainwater samples collected over the South Pacific, whereas no significant correlations were found between NH₄⁺ and MSA in rainwater collected over the subarctic (r = 0.42, p > 0.1, n = 6) and subtropical (r = 0.33, p > 0.5, n = 6) western North Pacific, suggesting that emissions of ammonia (NH₃) by marine biological activity from the ocean could become a significant source of NH₄⁺ over the South Pacific. While NO₃⁻ was the dominant inorganic N species in dry deposition, inorganic N supplied to surface waters by wet deposition was predominantly by NH₄⁺ (42–99% of the wet deposition fluxes for total inorganic N). We estimated mean total (dry + wet) deposition fluxes of atmospheric total inorganic N in the Pacific Ocean to be 32–64 μmol m⁻² d⁻¹, with 66–99% of this by wet deposition, indicating that wet deposition plays a more important role in the supply of atmospheric inorganic N than dry deposition.     

Temporal and Spatial Variations of the Aerodynamic Roughness Length in the Ablation Zone of the Greenland Ice Sheet

To understand the response of the Greenland ice sheet to climate change the so-called ablation zone is of particular importance, since it accommodates the yearly net surface ice loss. In numerical models and for data analysis, the bulk aerodynamic method is often used to calculate the turbulent surface fluxes, for which the aerodynamic roughness length (z ₀) is a key parameter. We present, for the first time, spatial and temporal variations of z ₀ in the ablation area of the Greenland ice sheet using year-round data from three automatic weather stations and one eddy-correlation mast. The temporal variation of z ₀ is found to be very high in the lower ablation area (factor 500) with, at the end of the summer melt, a maximum in spatial variation for the whole ablation area of a factor 1000. The variation in time matches the onset of the accumulation and ablation season as recovered by sonic height rangers. During winter, snow accumulation and redistribution by snow drift lead to a uniform value of z ₀≈ 10⁻⁴ m throughout the ablation area. At the beginning of summer, snow melt uncovers ice hummocks and z ₀ quickly increases well above 10⁻² m in the lower ablation area. At the end of summer melt, hummocky ice dominates the surface with z ₀ > 5  ×  10⁻³ m up to 60 km from the ice edge. At the same time, the area close to the equilibrium line (about 90 km from the ice edge) remains very smooth with z ₀ = 10⁻⁵ m. At the beginning of winter, we observed that single snow events have the potential to lower z ₀ for a very rough ice surface by a factor of 20 to 50. The total surface drag of the abundant small-scale ice hummocks apparently dominates over the less frequent large domes and deep gullies. The latter results are verified by studying the individual drag contributions of hummocks and domes with a drag partition model.     

Influence of the boundary layer height on the global air–sea surface fluxes

Results from large-eddy simulations and field measurements have previously shown that the velocity field is influenced by the boundary layer height, z _i , during close to neutral, slightly unstable, atmospheric stratification. During such conditions the non-dimensional wind profile, φ _m , has been found to be a function of both z/L and z _i /L. At constant z/L, φ _m decreases with decreasing boundary layer height. Since φ _m is directly related to the parameterizations of the air–sea surface fluxes, these results will have an influence when calculating the surface fluxes in weather and climate models. The global impact of this was estimated using re-analysis data from 1979 to 2001 and bulk parameterizations. The results show that the sum of the global latent and sensible mean heat fluxes increase by 0.77 W m⁻² or about 1% and the mean surface stress increase by 1.4 mN m⁻² or 1.8% when including the effects of the boundary layer height in the parameterizations. However, some regions show a larger response. The greatest impact is found over the tropical oceans between 30°S and 30°N. In this region the boundary layer height influences the non-dimensional wind profile during extended periods of time. In the mid Indian Ocean this results in an increase of the mean annual heat fluxes by 2.0 W m⁻² and an increase of the mean annual surface stress by 2.6 mN m⁻².     

An atmospheric dispersion model with continuous vertical variation of dispersion parameters

Summary A dispersion model is proposed to predict the continuous vertical variation of the dispersion parameters _y and _z in case of hot pollutant release to the atmosphere. In such a case, the plume rises far above the ground and is subject to varying levels of turbulence. The framework in this paper can be divided into three approaches: (1) determination of the eddy diffusivitiesK _y (z, _y ) andK _z (z, _z ) as functions of height above ground and plume dimensions, (2) determination of both the plume rise and its vertical velocity using a modified version of Brigg's formula, and (3) numerical solution of actual problems with buoyant plumes at each time step. The model results have been applied to a case of pollutant release from fire destruction of a chemical storehouse roof.With 15 Figures     

Dispersion of a Point-Source Release of a Passive Scalar Through an Urban-Like Array for Different Wind Directions

The dispersion of a point-source release of a passive scalar in a regular array of cubical, urban-like, obstacles is investigated by means of direct numerical simulations. The simulations are conducted under conditions of neutral stability and fully rough turbulent flow, at a roughness Reynolds number of Re _τ  = 500. The Navier–Stokes and scalar equations are integrated assuming a constant rate release from a point source close to the ground within the array. We focus on short-range dispersion, when most of the material is still within the building canopy. Mean and fluctuating concentrations are computed for three different pressure gradient directions (0°, 30°, 45°). The results agree well with available experimental data measured in a water channel for a flow angle of 0°. Profiles of mean concentration and the three-dimensional structure of the dispersion pattern are compared for the different forcing angles. A number of processes affecting the plume structure are identified and discussed, including: (i) advection or channelling of scalar down ‘streets’, (ii) lateral dispersion by turbulent fluctuations and topological dispersion induced by dividing streamlines around buildings, (iii) skewing of the plume due to flow turning with height, (iv) detrainment by turbulent dispersion or mean recirculation, (v) entrainment and release of scalar in building wakes, giving rise to ‘secondary sources’, (vi) plume meandering due to unsteady turbulent fluctuations. Finally, results on relative concentration fluctuations are presented and compared with the literature for point source dispersion over flat terrain and urban arrays.     

Water interaction with MgCl<Subscript>2</Subscript>×6H<Subscript>2</Subscript>O and NaCl surfaces: measurements of the uptake coefficient

The uptake of water vapor on MgCl₂×6H₂O and NaCl salt dry solid films was studied over the temperature range 240 to 340 K and at 1 Torr pressure of helium using a flow reactor coupled to a modulated molecular beam mass spectrometer. The H₂O to salt uptake data were obtained from the kinetics of H₂O loss on salt coated Pyrex rods. The following Arrhenius expression was obtained for the initial uptake coefficient of H₂O on MgCl₂×6H₂O films: γ ₀ (MgCl₂) = (6.5 ± 1.0) × 10⁻⁶ exp[(470 ± 40)/T] (calculated with specific BET surface area, quoted uncertainties are 1σ statistical). The rate of H₂O adsorption on NaCl was found to be much lower than on MgCl₂×6H₂O, and only an upper limit was determined for the corresponding uptake coefficient: γ (NaCl) ≤ 5.6 × 10⁻⁶ at T = 300 K. The results show that the rate of H₂O adsorption to salt surfaces is drastically dependent on the salt sample composition.     

The Fallacy of Drifting Snow

A common parametrization over snow-covered surfaces that are undergoing saltation is that the aerodynamic roughness length for wind speed (z ₀) scales as au_*²/g{\alpha u_\ast^2/g}, where u _* is the friction velocity, g is the acceleration of gravity, and α is an empirical constant. Data analyses seem to support this scaling: many published plots of z ₀ measured over snow demonstrate proportionality to u_*²{u_\ast^2 }. In fact, I show similar plots here that are based on two large eddy-covariance datasets: one collected over snow-covered Arctic sea ice; another collected over snow-covered Antarctic sea ice. But in these and in most such plots from the literature, the independent variable, u _*, was used to compute z ₀ in the first place; the plots thus suffer from fictitious correlation that causes z ₀ to unavoidably increase with u _* without any intervening physics. For these two datasets, when I plot z ₀ against u _* derived from a bulk flux algorithm—and thus minimize the fictitious correlation—z ₀ is independent of u _* in the drifting snow region, u _* ≥ 0.30 ms⁻¹. I conclude that the relation z₀ = au_*²/g{z_0 = \alpha u_\ast^2/g} when snow is drifting is a fallacy fostered by analyses that suffer from fictitious correlation.     

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.     

Analytical solutions for the Ekman layer

The PBL equation that governs the transition from the constant-stress surface layer to the geostrophic wind in a neutrally stratified atmosphere for which the eddy viscosityK(z) is assumed to vary smoothly from the surface-layer value U _*z (0.4,U _*=friction velocity,z=elevation) to the geostrophic asymptoteK _GU _*d forzd is solved through an expansion in fd/U _*1 (f=Coriolis parameter). The resulting solution is separated into Ekman's constant-K solution an inner component that reduces to the classical logarithmic form forzd and isO() relative to the Ekman component forzd. The approximationKU _*d is supported by the solution of Nee and Kovasznay's phenomenological transport equation forK(z), which yieldsK–U _*d exp(–z/d), where is an empirical constant for which observation implies, 1. The parametersA andB in Kazanskii and Monin's similarity relation forG/U _* (G=geostrophic velocity) are determined as functions of . The predicted values ofG/U _* and the turning angle are in agreement with the observed values for the Leipzig wind profile. The predicted value ofB based on the assumption of asymptotically constantK is 4.5, while that based on the Nee-Kovasznay model is 5.1; these compare with the observed value of 4.7 for the Leipzig profile. A thermal wind correction, an asymptotic solution for arbitraryK(z) and 1, and an exact (unrestricted ) solution forK(z)=U _*d[1–exp(–z/d)] are developed in appendices.     

Sulfur dioxide plume dispersion and subsequent absorption by coniferous forests: Field measurements and model calculations

A Forest SO₂ Absorption Model (ForSAM) was developed to simulate (1) SO₂ plume dispersion from an emission source, (2) subsequent SO₂ absorption by coniferous forests growing downwind from the source. There are three modules: (1) a buoyancy module, (2) a dispersion module, and (3) a foliar absorption module. These modules were used to calculate hourly abovecanopy SO₂ concentrations and in-canopy deposition velocities, as well as daily amounts of SO₂ absorbed by the forest canopy for downwind distances to 42 km. Model performance testing was done with meteorological data (including ambient SO₂ concentrations) collected at various locations downwind from a coal-burning power generator at Grand Lake in central New Brunswick, Canada. Annual SO₂ emissions from this facility amounted to about 30,000 tonnes. Calculated SO₂ concentrations were similar to those obtained in the field. Calculated SO₂ deposition velocities generally agreed with published values.Notation c air parcel cooling parameter (non-dimensional) - E foliar absorption quotient (non-dimensional) - f areal fraction of foliage free from water (non-dimensional) - f _w SO₂ content of air parcel - h height of the surface layer (m) - H height of the convective mixing layer (m) - H _stack stack height (m) - k time level - k _drag coefficient of drag on the air parcel (non-dimensional) - K _z eddy viscosity coefficient for SO₂ (m²·s^–1) - L Monin-Obukhov length scale (m) - L _A single-sided leaf area index (LAI) - n degree-of-sky cloudiness (non-dimensional) - N number of parcels released with every puff (non-dimensional) - PAR photosynthetically active radiation (W m^–2) - Q emission rate (kg s^–2) - r _b diffusive boundary-layer resistance (s m^–1) - r _c canopy resistance (s m^–1) - r _cuticle cuticular resistance (s m^–1) - r _m mesophyllic resistance (s m^–1) - r _s stomatal resistance (s m^–1) - r _exit smokestack exit radius (m) - R normally distributed random variable with mean of zero and variance of t (s) - u _* frictional velocity scale, (m s^–1) - v lateral wind vector (m s^–1) - v _d SO₂ dry deposition velocity (m s^–1) - VCD water vapour deficit (mb) - z _can mean tree height (m) - Z zenith position of the sun (deg) - environmental lapse rate (°C m^–1) - dry adiabatic lapse rate (0.00986°C m^–1) - von Kármán's constant (0.04) - _B vertical velocities initiated by buoyancy (m s^–1) - canopy extinction coefficient (non-dimensional) - ()a denotes ambient conditions - ()can denotes conditions at the top of the forest canopy - ()h denotes conditions at the top of the surface layer - ()H denotes conditions at the top of the mixed layer - ()s denotes conditions at the canopy surface - ()p denotes conditions of the air parcels     

Monin–Obukhov Similarity Functions of the Structure Parameter of Temperature and Turbulent Kinetic Energy Dissipation Rate in the Stable Boundary Layer

The Monin–Obukhov similarity theory (MOST) functions f_ε and f_T, of the dissipation rate of turbulent kinetic energy (TKE). ε, and the structure parameter of temperature, C_T², were determined for the stable atmospheric surface layer using data gathered in the context of CASES-99. These data cover a relatively wide stability range, i.e. ζ=z/L of up to 10, where z is the height and L the Obukhov length. The best fits were given by f_ε = 0.8 + 2.5ζ and f_T= 4.7[ 1+1.6(ζ)^2/3], which differ somewhat from previously published functions. ε was obtained from spectra of the longitudinal wind velocity using a time series model (ARMA) method instead of the traditional Fourier transform. The neutral limit f_ε =0.8 implies that there is an imbalance between TKE production and dissipation in the simplified TKE budget equation. Similarly, we found a production-dissipation imbalance for the temperature fluctuation budget equation. Correcting for the production-dissipation imbalance, the ‘standard’ MOST functions for dimensionless wind speed and temperature gradients (φ_m and φ_m) were determined from f_ε and f_T and compared with the φ_m and φ_h formulations of Businger and others. We found good agreement with the Beljaars and Holtslag [J. Appl. Meteorol. 30, 327–341 (1991)] relations. Lastly, the flux and gradient Richardson numbers are discussed also in terms of f_ε and f_T.     

An Evaluation of the Flux–Gradient Relationship in the Stable Boundary Layer

Data collected during the SHEBA and CASES-99 field programs are employed to examine the flux–gradient relationship for wind speed and temperature in the stably stratified boundary layer. The gradient-based and flux-based similarity functions are assessed in terms of the Richardson number Ri and the stability parameter z/Λ_*, z being height and Λ_* the local Obukhov length. The resulting functions are expressed in an analytical form, which is essentially unaffected by self-correlation, when thermal stratification is strong. Turbulence within the stably stratified boundary layer is classified into four regimes: “nearly-neutral” (0 < z/Λ_* < 0.02), “weakly-stable” (0.02 < z/Λ_* < 0.6), “very-stable” (0.6 < z/Λ_* < 50), and “extremely-stable” (z/Λ_* > 50). The flux-based similarity functions for gradients are constant in “nearly-neutral” conditions. In the “very-stable” regime, the dimensionless gradients are exponential, and proportional to (z/Λ_*)^3/5. The existence of scaling laws in “extremely-stable” conditions is doubtful. The Prandtl number Pr decreases from 0.9 in nearly-neutral conditions and to about 0.7 in the very-stable regime. The necessary condition for the presence of steady-state turbulence is Ri < 0.7.     

Scaling Properties of Temperature Spectra and Heat-Flux Cospectra in the Surface Friction Layer Beneath an Unstable Outer Layer

Temperature variance and temperature power spectra in the unstable surface layer have always presented a problem to the standard Monin-Obukhov similarity model. Recently that problem has intensified with the demonstration by Smedman et al. (2007, Q J Roy Meteorol Soc 133: 37–51) that temperature spectra and heat-flux cospectra can have two distinct peaks in slightly unstable conditions, and by McNaughton et al. (2007, Nonlinear Process Geophys 14: 257–271) who showed that the wavenumber of the peak of temperature spectra in a convective boundary layer (CBL), closely above the surface friction layer (SFL), can be sensitive to the CBL depth, z _i. Neither the two-peak form at slight instability nor the dependence of peak position on z _i at large instability is compatible with the Monin-Obukhov model. Here we examine the properties of temperature spectra and heat-flux cospectra from between these extremes, i.e. from within the unstable SFL, in two experiments. The analysis is based on McNaughton’s model of the turbulence structure in the SFL. According to this model, heat is transported through most of the SFL by sheet plumes, created by the action of impinging outer eddies. The smallest and most effective of these outer eddies have sizes that scale on SFL depth, z _s. The z _s-scale eddies and plumes are organised within the overall convection pattern in the CBL, and in turn they organise the motion of smaller eddies within the SFL, whose sizes scale on height, z. The main experimental results are: (1) the peak amplitudes of the temperature spectra in the SFL are collapsed with a scaling factor (z_sz)^1/3e_o^2/3{(z_{\rm s}z)^{1/3}\varepsilon_{\rm o}^{2/3}} divided by the square of the surface temperature flux, where e_o{\varepsilon_{\rm o}} is the dissipation rate of turbulent energy in the outer CBL (above the SFL); (2) the peak wavenumbers of the temperature spectra are collapsed with the mixed length scale (z _i z _s)^1/2; (3) the peak wavenumbers of the heat-flux cospectra are collapsed with the doubly-mixed length scale (z _i z _s)^1/4 z ^1/2; (4) for z/z _s < 0.03, the peak in the cospectrum is replaced by another peak at a wavenumber about a magnitude larger. This peak’s position scales on z; (5) all these findings are consistent with the observations of Smedman et al.     

A review of flux-profile relationships

Flux-profile relationships in the constant flux layer are reviewed. The preferred relationships are found to be those of Dyer and Hicks (1970), namely, _H = _W =(1–16(z/L))^–1/2, _M =(1–16(z/L))^–1/4 for the unstable region, and _H = _W = _M = 1+5(z/L) for the stable region.The carefully determined results of Businger et al. (1971) remain a difficulty which calls for considerable clarification.     

Turbulence decay in stratified and homogeneous marine layers

A spectral approach is applied to shear-induced turbulence in stratified layers. A system of spectral equations for stationary balance of turbulent energy and temperature variances was deduced in the vicinity of the local shear scale L_U = (ε/U_Z³)^1/2. At wavenumbers between the inertial-convective (k^−5/3) and wak turbulence (k⁻³) subranges, additional narrow spectral intervals—‘production’ subranges—may appear (E k⁻¹, E_T k⁻²). The upper boundary of these subranges is determined as L_U, and the lower boundaries as L_R (ε/U_ZN²)^1/2(χ/T_Z²). It is shown that the scale L_U is a unique spectral scale that is uniform up to a constant value for every hydrophysical field. It appears that the spectral scale L_U is equivalent to the Thorpe scale L_Th for the active turbulence model. Therefore, if turbulent patches are generated in a background of permanent mean shear, a linear relation between temperature and mass diffusivities exists. In spectral terms, the fossil turbulence model corresponds to the regime of the Boldgiano-Obukhov buoyancy subrange (E k^−11/5, E_T k^−7/5). During decay the buoyancy subrange is expanded to lower and higher wavenumbers. At lower wavenumbers the buoyancy subrange is bounded by L_** = 3(χ^1/2/N^1/2T_Z), which is equivalent to the Thorpe scale L_Th. In such a transition regime only, when the viscous dissipation rate is removed from the set of main turbulence parameters, the Thorpe scale does not correlate with the buoyancy scale L_N ε^1/2/N^3/2 and fossil turbulence is realized. Oceanic turbulence measurements in the equatorial Pacific near Baker Island confirm the main ideas of the active and fossil turbulence models.     

Organosulfates – A New Component of Humic-Like Substances in Atmospheric Aerosols?

Ion trap mass spectrometry (ITMS) was used to obtain further qualitative information about the chemical composition of humic-like substances (HULIS) in atmospheric particulate matter. Particles ≤10 μm (PM₁₀) were collected on quartz fiber filters for 24 h in the region of Basel (Switzerland) and extracted with water. HULIS were separated from inorganic salts by size exclusion chromatography (SEC) and detected by electrospray ionization in the negative ion mode (ESI(−)). Series of consecutive fragment ion spectra (MSⁿ) were recorded by ITMS. Full scan mass spectra of the extracts showed a mass distribution pattern characteristic for HULIS. Different molecular ions were selected from this pattern for further fragmentations. Among them the molecular ion m/z 299 was considered as representative and intensively studied. Many MS² and MS³ fragment spectra contained a fragment m/z 97 and a neutral loss of 80 u. Time-of-flight (TOF) MS and deuterium exchange experiments identified m/z 97 as hydrogen sulfate. MS² and MS³ fragment spectra supported the existence of sulfate covalently bound to HULIS. The fragmentation behavior of sulfated HULIS could be confirmed by model compounds.     

Critical Evaluation of Scalar Roughness Length Parametrizations Over a Melting Valley Glacier

We present a field investigation over a melting valley glacier on the Tibetan Plateau. In the ablation zone, aerodynamic roughness lengths (z _0M ) vary on the order of 10⁻⁴–10⁻² m, whose evolution corresponds to three melt phases with distinct surface cover and moisture exchange: snow (sublimation/evaporation), bare ice (deposition/condensation), and ice hummocks (sublimation/evaporation). Bowen-ratio similarity is validated in the stably stratified katabatic winds, which suggests a useful means for data quality check. A roughness sublayer is regarded as irrelevant to the present ablation season, because selected characteristics of scalar turbulence over smooth snow are quite similar to those over hummocky ice. We evaluate three parametrizations of the scalar roughness lengths (z _0T for temperature and z _0q for humidity), viz. key factors for the accurate estimation of sensible heat and latent heat fluxes using the bulk aerodynamic method. The first approach is based on surface-renewal models and has been widely applied in glaciated areas; the second has never received application over an ice/snow surface, despite its validity in (semi-)arid regions; the third, a derivative of the first, is proposed for use specifically over rough ice defined as z _0M > 10⁻³ m or so. This empirical z _0M threshold value is deemed of general relevance to glaciated areas (e.g. ice sheet/cap and valley/outlet glaciers), above which the first approach gives notably underestimated z _0T,q . The first and the third approaches tend to underestimate and overestimate turbulent heat/moisture exchange, respectively, frequently leading to relative errors higher than 30%. Comparatively, the second approach produces fairly low errors in energy flux estimates both in individual melt phases and over the whole ablation season; it thus emerges as a practically useful choice to parametrize z _0T,q in glaciated areas. Moreover, we find all three candidate parametrizations unable to predict diurnal variations in the excess resistances to humidity transfer, thus encouraging more efforts for improvement.     

Can Water Vapour Raman Lidar Resolve Profiles of Turbulent Variables in the Convective Boundary Layer?

High-resolution water vapour measurements made by the Atmospheric Radiation Measurement (ARM) Raman lidar operated at the Southern Great Plains Climate Research Facility site near Lamont, Oklahoma, U.S.A. are presented. Using a 2-h measurement period for the convective boundary layer (CBL) on 13 September 2005, with temporal and spatial resolutions of 10 s and 75 m, respectively, spectral and autocovariance analyses of water vapour mixing ratio time series are performed. It is demonstrated that the major part of the inertial subrange was detected and that the integral scale was significantly larger than the time resolution. Consequently, the major part of the turbulent fluctuations was resolved. Different methods to retrieve noise error profiles yield consistent results and compare well with noise profiles estimated using Poisson statistics of the Raman lidar signals. Integral scale, mixing-ratio variance, skewness, and kurtosis profiles were determined including error bars with respect to statistical and sampling errors. The integral scale ranges between 70 and 130 s at the top of the CBL. Within the CBL, up to the third order, noise errors are significantly smaller than sampling errors and the absolute values of turbulent variables, respectively. The mixing-ratio variance profile rises monotonically from ≈0.07 to ≈3.7 g² kg⁻² in the entrainment zone. The skewness is nearly zero up to 0.6 z/z _i , becomes −1 around 0.7–0.8 z/z _i , crosses zero at about 0.95 z/z _i , and reaches about 1.7 at 1.1 z/z _i (here, z is the height and z _i is the CBL depth). The noise errors are too large to derive fourth-order moments with sufficient accuracy. Consequently, to the best of our knowledge, the ARM Raman lidar is the first water vapour Raman lidar with demonstrated capability to retrieve profiles of turbulent variables up to the third order during daytime throughout the atmospheric CBL.