An observational study of the structure of the nocturnal urban boundary layer

The formation mechanism of the nocturnal urban boundary layer (UBL), especially in the winter nighttime, was investigated based on the extensive field observations conducted during November 1984 in Sapporo, Japan. A strong, elevated inversion formed over the Sapporo urban area and the inversion base height was approximately twice the average building height. Velocity fluctuations _u, _w and Reynolds stress % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaaca% WG1bWaaWbaaSqabeaacaaIXaaaaGGaaOGae8hiaaIaam4DamaaCaaa% leqabaGaaGymaaaaaaaaaa!3A9C!\[\overline {u^1 w^1 } \] had nearly uniform profiles within the nocturnal UBL and decreased with height above the UBL. On the other hand, temperature fluctuations _t , and heat fluxes % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaaca% WG1bWaaWbaaSqabeaacaaIXaaaaGGaaOGae8hiaaIaeqiUde3aaWba% aSqabeaacaaIXaaaaaaaaaa!3B56!\[\overline {u^1 \theta ^1 } \] and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaaca% WG3bWaaWbaaSqabeaacaaIXaaaaGGaaOGae8hiaaIaeqiUde3aaWba% aSqabeaacaaIXaaaaaaaaaa!3B58!\[\overline {w^1 \theta ^1 } \] had peaks at the inversion base and small values within the nocturnal UBL. The turbulent kinetic energy budget showed that the turbulent transport term and shear generation from urban canopy elements are important in the nocturnal UBL development; the role of the buoyancy term is small. The turbulence data analysis and application of a simple advective model showed that the mechanism of UBL formation may be controlled by the downward transport of sensible heat from the elevated inversion caused by mechanically-generated turbulence.Nomenclature g accelaration due to gravity, m s^-2 - k turbulent kinetic energy, m² s^-1 - K _m eddy viscosity, m² s^-1 - L Monin-Obukhov lenght, m - p pressure, Kg m^-2 - U, V, W mean wind speed in the downwind, crosswind, and vertical directions, respectively, m s^-1 - u ¹, w ¹ wind speed fluctuation in the downwind and vertical direction, respectively, m s^-1 - u ₁ friction velocity, m s^-1 - % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaaca% WG1bWaaWbaaSqabeaacaaIXaaaaGGaaOGae8hiaaIaam4DamaaCaaa% leqabaGaaGymaaaaaaaaaa!3A9C!\[\overline {u^1 w^1 } \] momentum flux, m²s^-2 - % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaaca% WG1bWaaWbaaSqabeaacaaIXaaaaGGaaOGae8hiaaIaam4DamaaCaaa% leqabaGaaGymaaaaaaaaaa!3A9C!\[\overline {u^1 \theta^1 } \] sensible heat flux, m²s^-1°C - WD wind direction, deg - WS wind speed, m s^-1 - z altitude, m - Z _h inversion base height, m - Z _j wind maximum height, m - Z _t inversion top height, m - T _u-r heat island intensity, °C - temperature lapse rate at rural site, °C m^-1 - energy dissipation rate, m²s^-3 - ¹ Potential temperature fluctuation, °C - _* scaling temperature, (=-% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaaca% WG1bWaaWbaaSqabeaacaaIXaaaaGGaaOGae8hiaaIaeqiUde3aaWba% aSqabeaacaaIXaaaaaaaaaa!3B56!\[\overline {u^1 \theta ^1 } \]/u_*) °C - mean potential temperature fluctuation, K - density of air, Kgm^-3 - K von Kármán constant (=0.4) - _u, _v, _w standard deviation of wind speed in the downwind, crosswind, and vertical directions, respectively, m s^-1 - standard diviation of temperature, °C     

Etude experimentale d'une relation entre le coefficient de frottement et le facteur de rafales en regime de vent faible et au-dessus d'une surface d'eau

In usual aerodynamic bulk formulas, the drag coefficient C _d has been best estimated in the 5 to 16 m s^–1 range of mean wind velocity; a value of 1.3 × 10^–3 is often considered for operational use. However, in the 0 to 5 m s^–1 range of mean wind velocity, corresponding to meteorological conditions of very light wind, experimental results have not resulted in any convincing agreement between various authors (Hicks et al., 1974; Wu, 1969; Kondo and Fujinawa, 1972; Mitsuta, 1973; Brocks and Krugermeyer, 1970).In the present paper, the drag coefficient is experimentally determined in conditions of very light wind and limited fetch (about 250 m). Due to this limited fetch, we have to be cautious in the extrapolation of our results to other sites. Nevertheless, some of experimental results are worth describing, considering the paucity of data in light wind conditions.Mean value and standard deviation (respectively 1.84 × 10^–3 and 1.24 × 10^–3) are obtained from 70 runs of 10-min duration. Mean wind velocities observed at 2 m above water surface are found to lie between 1.2 and 3.6 m s^–1. Whereas this mean value is in fair agreement with C _{d 10} = 1.3 × 10^–3, usually given for the 5 to 16 m s^–1 range (Kraus, 1972), the above value for the standard deviation seems too large to be left without further analysis.A more exhaustive analysis of the 70 values obtained for C _d shows that it depends on a parameter characteristic of longitudinal fluctuations of the wind velocity. A similar idea was put forward earlier by Kraus (1972). Relations between the drag coefficient and wind fluctuations may be tentatively given by: % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa% aaleaacaWGKbGaaGOmaaqabaGccqGH9aqpdaqadaqaaiabgkHiTiaa% igdacaGGUaGaaGimaiaaiEdacqGHRaWkcaaIXaGaaGinaiaac6caca% aIZaGaaGinamaalaaabaGaeq4Wdm3aaSbaaSqaaiaadwhacaGGNaaa% beaaaOqaaaaaaiaawIcacaGLPaaaruqqYLwySbacfaGaa8hEaiaa-b% cacaaIXaGaaGimamaaCaaaleqabaGaeyOeI0IaaG4maaaakiaabcca% caqGGaGaaeiiaiaabccacaqGXaGaaeOlaiaabAdacaqGGaGaaeyBai% aabccacaqGZbWaaWbaaSqabeaacaqGTaGaaeymaaaakiabgsMiJkqa% dwhagaqeamaaBaaaleaacaaIYaaabeaakiabgsMiJkaaiodacaGGUa% GaaGOnaiaab2gacaqGGaGaae4CamaaCaaaleqabaGaaeylaiaabgda% aaaaaa!634E!\[C_{d2} = \left( { - 1.07 + 14.34\frac{{\sigma _{u'} }}{{}}} \right)x 10^{ - 3} {\text{ 1}}{\text{.6 m s}}^{{\text{ - 1}}} \leqslant \bar u_2 \leqslant 3.6{\text{m s}}^{{\text{ - 1}}} \] and % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa% aaleaacaWGKbGaaGOmaaqabaGccqGH9aqpdaqadaqaaiabgkHiTiaa% iodacaGGUaGaaGioaiaaiAdacqGHRaWkcaaIZaGaaiOlaiaaiodaca% aI2aGaam4raaGaayjkaiaawMcaaerbbjxAHXgaiuaacaWF4bGaa8hi% aiaaigdacaaIWaWaaWbaaSqabeaacqGHsislcaaIZaaaaOGaaeilaa% aa!4B42!\[C_{d2} = \left( { - 3.86 + 3.36G} \right)x 10^{ - 3} {\text{,}}\] where _u/\-u ₂ and G, respectively, represent the standard deviation of u normalized with \-u ₂ and the longitudinal gust factor quoted in Smith (1974).We have established a relationship between these fluctuation parameters and the stability as given by a bulk layer Richardson number (between 0 and 2 m). These relations are given by: % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSaaaeaacq% aHdpWCdaWgaaWcbaGaamyDaiaacEcaaeqaaaGcbaGabmyDayaaraWa% aSbaaSqaaiaaikdaaeqaaaaakiabg2da9iaaicdacaGGUaGaaGymai% aaikdacqGHRaWkcaaIZaGaaiOlaiaaiIdacaaI1aGaaeiiaiaabkfa% caqGPbWaaSbaaSqaaiaabcdacaqGTaGaaeOmaaqabaaaaa!4802!\[\frac{{\sigma _{u'} }}{{\bar u_2 }} = 0.12 + 3.85{\text{ Ri}}_{{\text{0 - 2}}} \] and G=1.35+14.56 Ri_0–2. The increase in gustiness with stability is in qualitative agreement with Goptarev (1957)'s experimental results.In spite of the high-level correlation between C _d and _u/\-u ₂(G) on the one hand and between _u/\-u ₂(G) and Ri_0–2on the other hand, we found a poor relationship between C _d and Ri_0–2. It is worth noting too that the trend observed here for C _d to increase with stability is in complete disagreement with the usual theoretical expectation for C _d to decrease with increasing layer stability above water.

---

---

E.R.A. du C.N.R.S. n 259.     

Characteristics of turbulent fluxes of sensible heat and momentum in the surface boundary layer during the Indian summer monsoon

Turbulent fluctuations of wind and temperature were measured using a three-component sonic anemometer at 8 m on a 30 m micro-meteorological tower erected at the Indian Institute of Technology (IIT) Kharagpur (22.3° N, 87.2° E), India, as part of the Monsoon Trough Boundary Layer Experiment (MONTBLEX). Diurnal and nocturnal variations of fluxes of sensible heat and momentum were estimated by the eddy correlation technique from 42 observations, each of 10 min duration during 6–8 July in the monsoon season of 1989. The estimated heat flux shows a diurnal trend while the momentum flux shows variability but no particular trend. The nocturnal heat flux changes sign intermittently.Fluctuations of vertical wind velocity _wand temperature when normalised with the respective scaling parameters u _*and _* are found to scale with Z/L in accordance with the Monin-Obukhov similarity hypothesis: % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqOXdy2aaS% baaSqaaiaadEhaaeqaaOGaamiEaiaacIcacaWGAbGaai4laiaadYea% caGGPaWaaWbaaSqabeaacaaIXaGaai4laiaaiodaaaaaaa!3FE8!\[\phi _w x(Z/L)^{1/3} \], % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqOXdy2aaS% baaSqaaiabeI7aXbqabaGccaWG4bGaaiikaiaadQfacaGGVaGaamit% aiaacMcadaahaaWcbeqaaiaaigdacaGGVaGaaG4maaaaaaa!40A2!\[\phi _\theta x(Z/L)^{1/3} \] during unstable conditions and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqOXdy2aaS% baaSqaaiaadEhaaeqaaOGaamiEaiaacIcacaWGAbGaai4laiaadYea% caGGPaaaaa!3D90!\[\phi _w x(Z/L)\], % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqOXdy2aaS% baaSqaaiabeI7aXbqabaGccaWG4bGaaiikaiaadQfacaGGVaGaamit% aiaacMcadaahaaWcbeqaaiabgkHiTiaaigdaaaaaaa!401F!\[\phi _\theta x(Z/L)^{ - 1} \] during stable conditions. Correlation coefficients for heat and momentum flux % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4SdC2aaS% baaSqaaiaadEhacqaH4oqCaeqaaaaa!3A71!\[\gamma _{w\theta } \] and _uwshow stability dependence. For unstable conditions, % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4SdC2aaS% baaSqaaiaadEhacqaH4oqCaeqaaaaa!3A71!\[\gamma _{w\theta } \] increases with increasing ¦Z/L¦ whereas _uwdecreases. During stable conditions, % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4SdC2aaS% baaSqaaiaadEhacqaH4oqCaeqaaaaa!3A71!\[\gamma _{w\theta } \] decreases with increasing Z/L while _uwdecreases very slowly from a value -0.36 to -0.37.     

An analytical-empirical method for determining the roughness length and zero-plane displacement

A method based on the principle of the Method of Weighted Residuals for evaluating the roughness-length (z ₀) and zero-plane displacement (d) is presented. This method not only can minimize errors involved during the calculation process but can also smooth and re-distribute the already minimized error in a most favourable manner via using appropriate weighting functions. With the inclusion of d in addition to z ₀, formulae for wind and temperature profiles in the surface layer are presented by:% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGvbGaeyypa0% ZaaSaaaeaacaWG1bWaaSbaaSqaaiaacQcaaeqaaaGcbaGaam4Aaaaa% daWadaqaaiGacYgacaGGUbWaaeWaaeaadaWcaaqaaiaadQhacqGHsi% slcaWGKbaabaGaamOEamaaBaaaleaacaaIWaaabeaaaaaakiaawIca% caGLPaaacqGHRaWkcqaHipqEaiaawUfacaGLDbaaaaa!43FC!\[U = \frac{{u_* }}{k}\left[ {\ln \left( {\frac{{z - d}}{{z_0 }}} \right) + \psi } \right]\]and% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH4oqCcqGHsi% slcqaH4oqCdaWgaaWcbaGaaGimaaqabaGccqGH9aqpcqaH4oqCdaWa% daqaaiGacYgacaGGUbWaaeWaaeaadaWcaaqaaiaadQhacqGHsislca% WGKbaabaGaamOEamaaBaaaleaacaaIWaaabeaaaaaakiaawIcacaGL% PaaacqGHRaWkcqaHipqEdaWgaaWcbaacbmGaa8hvaaqabaaakiaawU% facaGLDbaaaaa!485A!\[\theta - \theta _0 = \theta \left[ {\ln \left( {\frac{{z - d}}{{z_0 }}} \right) + \psi _T } \right]\]where and _T are the integrated diabetic influence functions' for velocity and temperature profiles respectively.Analytical expressions for both and _T as functions of wind shear or, implicitly in terms of the Richardson number have been derived.Presented at the 10th Annual Congress of the Canadian Meteorological Society, Quebec City, Canada, May 26–28, 1976.     

Review of some basic characteristics of the atmospheric surface layer

Some of the fundamental issues of surface layer meteorology are critically reviewed. For the von Karman constant (k), values covering the range from 0.32 to 0.65 have been reported. Most of the data are, however, found in a rather narrow range between 0.39 and 0.41. Plotting all available atmospheric data against the so-called roughness Reynolds number, % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeOuaiaabw% gadaWgaaWcbaGaaeimaaqabaGccqGH9aqpcaWG1bWaaSbaaSqaaiaa% cQcaaeqaaOGaamOEamaaBaaaleaacaaIWaaabeaakiaac+cacqaH9o% GBaaa!3FD0!\[{\rm{Re}}_{\rm{0}} = u_* z_0 /\nu \] or against the surface Rossby number, % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeOuaiaab+% gadaWgaaWcbaGaaeimaaqabaGccqGH9aqpcaWGhbGaai4laiaadAga% caWG6bWaaSbaaSqaaiaaicdaaeqaaaaa!3DF1!\[{\rm{Ro}}_{\rm{0}} = G/fz_0 \] gives no clear indication of systematic trend. It is concluded that k is indeed constant in atmospheric surface-layer flow and that its value is the same as that found for laboratory flows, i.e. about 0.40.Various published formulae for non-dimensional wind and temperature profiles, _m and _h respectively, are compared after adjusting the fluxes so as to give % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4Aaiabg2% da9iaaicdacaGGUaGaaGinaiaaicdaaaa!3AC6!\[k = 0.40\] and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaeWabeaaii% GacqWFgpGzdaWgaaWcbaGaamiAaaqabaGccaGGVaGae8NXdy2aaSba% aSqaaiaad2gaaeqaaaGccaGLOaGaayzkaaWaaSbaaSqaaiaadQhaca% GGVaGaamitaiabg2da9iaaicdaaeqaaOGaeyypa0JaaGimaiaac6ca% caaI5aGaaGynaaaa!4655!\[\left( {\phi _h /\phi _m } \right)_{z/L = 0} = 0.95\]. It is found that for % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaqWabeaaca% WG6bGaai4laiaadYeaaiaawEa7caGLiWoacqGHKjYOcaaIWaGaaiOl% aiaaiwdaaaa!3F72!\[\left| {z/L} \right| \le 0.5\] the various formulae agree to within 10–20%. For unstable stratification the various formulations for _h continue to agree within this degree of accuracy up to at least % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOEaiaac+% cacaWGmbGaeyisISRaeyOeI0IaaGOmaaaa!3BC9!\[z/L \approx - 2\]. For _m in very unstable conditions results are still conflicting. Several recent data sets agree that for unstable stratification % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeOuaiaabM% gacqGHijYUcaaIXaGaaiOlaiaaiwdacaWG6bGaai4laiaadYeaaaa!3E0D!\[{\rm{Ri}} \approx 1.5z/L\] up to at least % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeyOeI0Iaam% OEaiaac+cacaWGmbGaeyypa0JaaGimaiaac6cacaaI1aaaaa!3C8D!\[ - z/L = 0.5\] and possibly well beyond.For the Kolmogorov streamwise inertial subrange constant, _u , it is concluded from an extensive data set that % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySde2aaS% baaSqaaiaadwhaaeqaaOGaeyypa0JaaGimaiaac6cacaaI1aGaaGOm% aiabgglaXkaaicdacaGGUaGaaGimaiaaikdaaaa!4178!\[\alpha _u = 0.52 \pm 0.02\]. The corresponding constant for temperature is much more uncertain, its most probable value being, however, about 0.80, which is also the most likely value for the corresponding constant for humidity.The turbulence kinetic energy budget is reviewed. It is concluded that different data sets give conflicting results in important respects, particularly so in neutral conditions.It is demonstrated that the inertial-subrange method can give quite accurate estimates of the fluxes of momentum, sensible heat and water vapour from high frequency measurements of wind, temperature and specific humidity alone, provided apparent values of the corresponding Kolmogorov constants are used. For temperature and humidity, the corresponding values turn out to be equal to the true constants, so % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqOSdi2aaS% baaSqaaiaadgeaaeqaaOGaeyisISRaeqOSdiMaeyisISRaaGimaiaa% c6cacaaI4aGaaGimaaaa!4074!\[\beta _A \approx \beta \approx 0.80\]. For momentum, however, the apparent constant % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySde2aaS% baaSqaaiaadwhacaWGbbaabeaakiabgIKi7kaaicdacaGGUaGaaGOn% aiaaicdaaaa!3E18!\[\alpha _{uA} \approx 0.60\].Based on an invited paper presented at the EGS Workshop Instrumental and Methodical Problems of Land Surface Flux Measurements, Grenoble 22–26 April, 1994.     

Measuring C T 2 , C Q 2 , and C TQ in the unstable surface layer,and relations to the vertical fluxes of heat and moisture

The structure parameters of temperature (C _T ² ), humidity (C _Q ² ) and temperature-humidity (C _TQ ) were observed at a height of 4 m in the unstable surface layer using thin platinum wires and two Ly- hygrometers. Two ways of measuring structure parameters were employed: one using spaced sensors, the other using time-delayed observations at one location. It is found that the three structure parameters follow free-convection scaling down to -z/L 0.02. The scaling functions % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGMbaaaa!3346!\[f\]₁ (of C _T ² ), % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGMbaaaa!3346!\[f\]₂ (of C _TO ) and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGMbaaaa!3346!\[f\]₃ (of C _Q ² ) are found to be related through % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGMbaaaa!3346!\[f\]₂/% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGMbaaaa!3346!\[f\]₁ 0.69 and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGMbaaaa!3346!\[f\]₃/% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGMbaaaa!3346!\[f\]₁ 0.84. The usefulness of the structure parameters for inferring the fluxes of heat and water vapor, as well as the Bowen ratio, is demonstrated. The scatter is about 30% on either side of the mean.This work was done while the author was a visiting scientist at the Wave Propagation Laboratory, NOAA, ERL, Boulder, U.S.A. He received support from the Netherlands Minister for Science Policy and the U.S. Army Research Office.     

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.     

A framework for evaluating air quality models

This paper describes a framework to evaluate air quality model predictions against observations. We propose the following relationship between observations and predictions from an adequate model% MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabm4qayaaja% WaaSbaaSqaaiaaicdaaeqaamXvP5wqonvsaeHbfv3ySLgzaGqbaOGa% e8hkaGIaamiEamaaBaaaleaacaaIXaaabeaakiaacYcacaWG4bWaaS% baaSqaaiaaikdaaeqaaOGae8xkaKIaeyypa0Jabm4qayaajaWaaSba% aSqaaiaadchaaeqaaOGae8hkaGIaamiEamaaBaaaleaacaaIXaaabe% aakiab-LcaPiab-TcaRiabew7aLjab-HcaOiaadIhadaWgaaWcbaGa% aGOmaaqabaGccqWFPaqkaaa!4F93!\[\hat C_0 (x_1 ,x_2 ) = \hat C_p (x_1 ) + \varepsilon (x_2 )\],where x ₁ refers to the inputs used in the model prediction C _p(x ₁), and x ₂denotes unknown variables which affect the observed concentration C ₀. The hats associated with C _pand C ₀denote transformations to convert the residual to a white noise sequence which is normally distributed. In this paper we assume % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabm4qayaaja% GaeyyyIORaciiBaiaac6gacaWGdbaaaa!3B39!\[\hat C \equiv \ln C\].The standard deviation of determines the expected deviation between model prediction and observation. The purpose of model improvement is to make this deviation as small as possible.The formalism we have proposed is applied to the evaluation of two models developed by this author. We show how careful analysis of residuals can lead to improvements in the model. We have also estimated for each of the models.In the last part of the part of the paper we show how the statistics of can be used to interpret model predictions.     

Coherent eddies and turbulence in vegetation canopies: The mixing-layer analogy

This paper argues that the active turbulence and coherent motions near the top of a vegetation canopy are patterned on a plane mixing layer, because of instabilities associated with the characteristic strong inflection in the mean velocity profile. Mixing-layer turbulence, formed around the inflectional mean velocity profile which develops between two coflowing streams of different velocities, differs in several ways from turbulence in a surface layer. Through these differences, the mixing-layer analogy provides an explanation for many of the observed distinctive features of canopy turbulence. These include: (a) ratios between components of the Reynolds stress tensor; (b) the ratio K _H/K _M of the eddy diffusivities for heat and momentum; (c) the relative roles of ejections and sweeps; (d) the behaviour of the turbulent energy balance, particularly the major role of turbulent transport; and (e) the behaviour of the turbulent length scales of the active coherent motions (the dominant eddies responsible for vertical transfer near the top of the canopy). It is predicted that these length scales are controlled by the shear length scale % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamitamaaBa% aaleaacaWGtbaabeaakiabg2da9iaadwfacaGGOaGaamiAaiaacMca% caGGVaGabmyvayaafaGaaiikaiaadIgacaGGPaaaaa!3FD0!\[L_S = U(h)/U'(h)\] (where h is canopy height, U(z) is mean velocity as a function of height z, and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabmyvayaafa% Gaeyypa0JaaeizaiaadwfacaGGVaGaaeizaiaadQhaaaa!3C32!\[U' = {\rm{d}}U/{\rm{d}}z\]). In particular, the streamwise spacing of the dominant canopy eddies is _x = mL _s, with m = 8.1. These predictions are tested against many sets of field and wind-tunnel data. We propose a picture of canopy turbulence in which eddies associated with inflectional instabilities are modulated by larger-scale, inactive turbulence, which is quasi-horizontal on the scale of the canopy.     

The influence of anisotropy on stress estimation by the indirect dissipation method

A comparison of observations by different authors reveals that systematic differences exist between momentum fluxes measured directly, and momentum fluxes determined indirectly by the dissipation method. This discrepancy is attributed to systematic errors due to the indirect determination of energy dissipation from the presumed inertial subrange spectrum of the horizontal wind component. The discrepancy increases with increasing degree of anisotropy, indicated by the ratio (vertical wind spectrum): (horizontal wind spectrum) deviating from % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWcaaqaaiaais% daaeaacaaIZaaaaaaa!33E6!\[\frac{4}{3}\]The results support a value of 0.48 for Kolmogoroff's constant.     

Kolmogoroff constants at the 1976 ITCE

The high-frequency data from 12 sensors at the ITCE 1976^* are analysed to determine the Kolmogoroff constants for velocity, temperature and humidity fluctuation, _u , _T , and _q . The occurrence of aliasing in the spectral analysis in some cases together with the limited response of some sensors at the higher frequencies introduces some uncertainties into the analysis. The Soviet sonic anemometer, fine-wire thermometer and infrared hygrometer and the Australian infrared hygrometer provide the best information, namely that% MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9qq-f0-yqaqVeLsFr0-vr% 0-vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaHXoqyda% WgaaWcbaGaamyDaaqabaGccqGH9aqpcaaIWaGaaiOlaiaaiwdacaaI% 5aGaeyySaeRaaGimaiaac6cacaaIWaGaaGymaiaacYcacaqGGaGaae% iiaiaabccacaqGGaGaeqySde2aaSbaaSqaaiaadsfaaeqaaOGaeyyp% a0JaaGimaiaac6cacaaI2aGaaGioaiabgglaXkaaicdacaGGUaGaaG% imaiaaikdacaGGSaGaaeiiaiaabccacaqGGaGaaeiiaiabeg7aHnaa% BaaaleaacaWGXbaabeaakiabg2da9iaaicdacaGGUaGaaG4naiaaiA% dacqGHXcqScaaIWaGaaiOlaiaaicdacaaIZaaaaa!6248!\[\alpha _u = 0.59 \pm 0.01,{\text{ }}\alpha _T = 0.68 \pm 0.02,{\text{ }}\alpha _q = 0.76 \pm 0.03\]where the errors quoted refer solely to statistical errors. The other instruments provide general support to these values.The technique of using spectral density measurements to determine eddy fluxes is illustrated.International Turbulence Comparison Experiment.     

Stability functions correct at the free convection limit and consistent for for both the surface and Ekman layers

For the thermal stability function _h used to calculate heat and moisture fluxes in the surface layer, we choose a formulation which has the theoretically correct free convection limit % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeikaiabgk% HiTGqaciaa-PhacaqGVaGaamitaiaabMcadaahaaWcbeqaaiabgkHi% TiaaigdacaGGVaGaaG4maaaaaaa!3DFE!\[{\rm{(}} - z{\rm{/}}L{\rm{)}}^{ - 1/3} \]. We then use the experimental result that z/L Ri to deduce a formulation with an exponent -1/6 for the momentum stability function _m. This formulation also resolves the matching problem at the interface between the surface and Ekman layers. The proposed functions are found to remain reasonably close to another formulation that is well supported by observations and has exponents -1/2 for _h and -1/4 for _m. The intent of the proposals is mainly to clarify and simplify the parameterization of the convective boundary layer in present day atmospheric models, without significantly altering the results.     

Mesoscale nocturnal jetlike winds within the planetary boundary layer over a flat,open coast

Mesoscale nocturnal jetlike winds have been observed over a flat, open coast. They occur within the planetary boundary layer between 100 and 600 m. At times the wind shear may reach 15 m s^-1 per 100 m. Unlike the common low-level jet that occurs most often at the top of the nocturnal inversion and only with a wind from the southerly quadrant, this second kind of jet exists between nocturnal ground-based inversion layers formed by the cool pool, or mesohigh, and the elevated mesoscale inversion layer over the coast. It occurs mostly when light % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaiikaiabgs% MiJkaaiwdacqGHsislcaaI2aGaaeyBaiaabccacaqGZbWaaWbaaSqa% beaacqGHsislcaaIXaaaaOGaaiykaaaa!3FCF!\[( \leqslant 5 - 6{\text{m s}}^{ - 1} )\] geostrophic winds blow from land to sea and when the air temperature over adjacent seas is more than 5 °C warmer than that over the coast. This phenomenon may be explained by combined Venturi and gravity-wind effects existing in a region from just above the area a few kilometres offshore to 100–600 m in height approximately 40–50 km inland because this region is sandwiched between the aforementioned two inversion layers.     

Observations in the nocturnal boundary layer

Low-latitude observations of the stably-stratified planetary boundary layer (SBL) above rough terrain are compared to observations of the mid-latitude SBL mainly through the depth h and its dependence upon surface fluxes. This involves the quantity h/L and the similarity prediction h = (u _* L/f)^1/2.Mid-latitude observations are consistent with model calculations for nighttime-averaged quantities and their deviations, as functions of latitude and surface roughness, from the equilibrium values found at large t. The above applies to horizontally-homogeneous terrain.Low-latitude observations of % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafq4SdCMbae% baaaa!37AB!\[\bar \gamma \] and h/L are significantly smaller than mid-latitude values, apparently the result of katabatic flows at the site and not the differences in latitude. This is consistent with model calculations for non-zero slope terrain.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Modelling velocity spectra in the lower part of the planetary boundary layer

Principles used when constructing models for velocity spectra are reviewed. Based upon data from the Kansas and Minnesota experiments, simple spectral models are set up for all velocity components in stable air at low heights, and for the vertical spectrum in unstable air through a larger part of the planetary boundary layer. Knowledge of the variation with stability of the (reduced) frequency % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaaccaGae8NKby% kaaa!37B5!\[f\] _m for the spectral maximum is utilized in this modelling. Stable spectra may be normalized so that they adhere to one curve only, irrespective of stability, and unstable w-spectra may also be normalized to fit one curve.The problem of using filtered velocity variances when modelling spectra is discussed. A simplified procedure to provide a first estimate of the filter effect is given.In stable, horizontal velocity spectra, there is often a gap at low frequencies. Using dimensional considerations and the spectral model previously derived, an expression for the gap frequency is found.     

The source area influencing a measurement in the Planetary Boundary Layer: The “footprint” and the “distribution of contact distance”

This paper considers the ground area which affects the properties of fluid parcels observed at a given spot in the Planetary Boundary Layer (PBL). We examine two source-area functions; the footprint, giving the source area for a measurement of vertical flux: and the distribution of contact distance, the distance since a particle observed aloft last made contact with the surface. We explain why the distribution of contact distance extends vastly farther upwind than the footprint, and suggest for the extent of the footprint the inequalities: % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqaqpepeea0xe9qqVa0l% b9peea0lb9Lq-JfrVkFHe9peea0dXdarVe0Fb9pgea0xa9W8qr0-vr% 0-viWZqaceaabiGaciaacaqabeaadaqaaqaaaOabaeqabaGaamyvam% aalaaabaGaamiAaaqaaiabeo8aZnaaBaaaleaacaWGxbaabeaakiaa% cIcacaWGObGaaiykaaaacqGH8aapcaWG4bGaeyipaWJaamyvaKazaa% iadaGabaqaamaaDaaajqwaacqaaiaadIgacaGGVaGabmOEayaacaGa% aiilaiaabccacaGGVbGaaiiDaiaacIgacaGGLbGaaiOCaiaacEhaca% GGPbGaai4CaiaacwgaaeaacaWGubWaaSbaaKazcaiabaGaamitaaqa% baqcKfaGaiaacIcacaWGObGaaiykaiaabYcacaqGGaGaaeiAaiaabc% cacaGGHbGaaiOyaiaac+gacaGG2bGaaiyzaiaabccacaGGZbGaaiyD% aiaackhacaGGMbGaaiyyaiaacogacaGGLbGaeyOeI0IaaiiBaiaacg% gacaGG5bGaaiyzaiaackhaaaaajqgaacGaay5EaaaakeaaaeaacaGG% 8bGaamyEaiaacYhacqGH8aapcqaHdpWCdaWgaaWcbaGaamODaaqaba% GccaGGOaGaamiAaiaacMcadaWcaaqaaiaadIhaaeaacaWGvbaaaaaa% aa!7877!\[\begin{array}{l} U\frac{h}{{\sigma _W (h)}} < x < U\left\{ {_{h/\dot z,{\rm{ }}otherwise}^{T_L (h){\rm{, h }}above{\rm{ }}surface - layer} } \right. \\ \\ |y| < \sigma _v (h)\frac{x}{U} \\ \end{array}\] where U is the mean streamwise (x) velocity, h is the observation height, _L is the Lagrangian timescale, _v and _w are the standard deviations of the cross-stream horizontal (y) and vertical (z) velocity fluctuations, and is the Lagrangian Similarity prediction for the rate of rise of the centre of gravity of a puff released at ground.Simple analytical solutions for the contact-time and the footprint are derived, by treating the PBL as consisting of two sub-layers. The contact-time solutions agree very well with the predictions of a Lagrangian stochastic model, which we adopt in the absence of measurements as our best estimate of reality, but the footprint solution offers no improvement over the above inequality.     

Sea surface winds derived from cloud motion vectors over the tropical Atlantic

The relationship between satellite-derived low-level cloud motion, surface wind and geostrophic wind vectors is examined using GATE data. In the trades, surface wind speeds can be derived from cloud motion vectors by the linear relation: V = 0.62 V _s + 1.9 m s^–1 with a mean scatter of ±1.3 m s^–1. The correlation coefficient between surface and satellite wind speed is 0.25. Considering baroclinicity, i.e., the influence of the thermal wind, the correlation coefficient does not increase, because of the uncertainty of the thermal wind vectors. The ratios of surface to geostrophic wind speed and surface to satellite wind speed are 0.7 and 0.8, respectively, with a statistical uncertainty of ±0.3. Calculations of the ratio of surface to geostrophic wind speed on the basis of the resistance law yield V/V _g = 0.8 ± 0.2, in agreement with experimental results. The mean angle difference between the surface and the satellite wind vectors amounts to - 18 °, taking into account baroclinicity. This value is in good agreement with the mean ageostrophic angle - 25 °.     

Kinetics of the reaction of nitrogen dioxide with ozone

The kinetics of the reaction of NO₂ with O₃ have been investigated at 296 K, using UV absorption spectroscopy to monitor decay of NO₂ or O₃ and infrared laser absorption spectroscopy to monitor formation of the reaction product N₂O₅. The results both for the rate coefficient at 296 K (k ₁=3.5×10^-17 cm³ molecule^-1 s^-1) and the reaction stoichiometry (NO₂/O₃=1.85±0.09) are in good agreement with previous studies, confirming that the two step mechanism involving formation of symmetrical NO₃ as an intermediate is predominant.% MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeOtaiaab+% eadaWgaaWcbaGaaeOmaaqabaGccqGHRaWkcaqGpbWaaSbaaSqaaiaa% bodaaeqaaOWaa4ajaSqaaaqabOGaayPKHaGaaeOtaiaab+eadaWgaa% WcbaGaae4maaqabaGccqGHRaWkcaqGpbWaaSbaaSqaaiaabkdaaeqa% aaaa!41D7!\[{\text{NO}}_{\text{2}} + {\text{O}}_{\text{3}} \xrightarrow{{}}{\text{NO}}_{\text{3}} + {\text{O}}_{\text{2}} \]% MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeOtaiaab+% eadaWgaaWcbaGaae4maaqabaGccqGHRaWkcaqGobGaae4tamaaBaaa% leaacaqGYaaabeaakiabgUcaRiaab2eadaGdKaWcbaaabeGccaGLsg% cacaqGobWaaSbaaSqaaiaabkdaaeqaaOGaae4tamaaBaaaleaacaqG% 1aaabeaakiabgUcaRiaab2eaaaa!4464!\[{\text{NO}}_{\text{3}} + {\text{NO}}_{\text{2}} + {\text{M}}\xrightarrow{{}}{\text{N}}_{\text{2}} {\text{O}}_{\text{5}} + {\text{M}}\]A possible minor role for the unsymmetrical ONOO species is suggested to account for the lower-than-expected stoichiometry factor. The importance of this reaction in the oxidation of atmospheric NO₂ is discussed.     

An international turbulence comparison experiment (ITCE 1976)

Turbulence data for the International Turbulence Comparison Experiment (ITCE) held at Conargo, N.S.W. (35° 18′ S., 145° 10′ E.) during October, 1976 are analysed. The standard deviation % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaiikaiqado% hagaqbamaaCaaaleqabaGaaGOmaaaakiaacMcadaahaaWcbeqaaiaa% igdacaGGVaGaaGOmaaaaaaa!3B93!\[(s'^2 )^{1/2} \] and covariance % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaa0aaaeaace% WG3bGbauaaceWGZbGbauaaaaaaaa!3809!\[\overline {w's'} \] measured by a number of instruments and instrument arrays have been compared to assess their field performance and calibration accuracy. Satisfactory agreement, i.e. typically 5% for % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaiikamaana% aabaGabm4CayaafaWaaWbaaSqabeaacaaIYaaaaaaakiaacMcadaah% aaWcbeqaaiaaigdacaGGVaGaaGOmaaaaaaa!3BA4!\[(\overline {s'^2 } )^{1/2} \] (except in humidity) and of the order of 20% for % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaiikaiqado% hagaqbamaaCaaaleqabaGaaGOmaaaakiaacMcadaahaaWcbeqaaiaa% igdacaGGVaGaaGOmaaaaaaa!3B93!\[(s'^2 )^{1/2} \], was achieved, but only after consideration of:

1. Instrumental response at high frequencies.
2. Flow distortion induced by instruments and supporting structures.
3. Spatial separation of instruments used for covariance measurements.
4. Statistical errors associated with single point measurements over a finite averaging time, and with lateral separation of two sensor arrays being compared.

     

Boundary-layer dynamics of a developing vortex

The dynamics of the boundary layer of a developing vortex is studied. To keep the mathematical treatment simple, axisymmetric vortices in both solid and non-solid rotation are considered. The governing equations are subjected to a similarity transformation and the resulting non-linear parabolic partial differential equations are solved by a Galerkin technique for simplified initial conditions which permit non-zero tangential and vertical velocities at the lower and upper boundaries. The application of a time-dependent Taylor boundary condition at the lower extremity of the vortex causes the fluid to spin-down gradually. By doing so, the no-slip steady-state solutions are progressively approached. On the other hand the employment of the customary Taylor boundary condition made the vortex to spin down rapidly causing a large-amplitude oscillatory vertical velocity.Analyses of the linearized equations and of the numerical solutions of the non-linear problem indicate the presence in the radial and tangential fields of inertial oscillations of frequency % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaGOmaKaaaj% aad6gakmaaleaaleaaruqqYLwySbacfaGaa8xmaaqaaiaa-jdaaaac% ciGccqGFPoWvaaa!3D19!\[2n\tfrac{1}{2}\Omega \] (: angular velocity of the flow and n: a parameter such that n=1 signifies solid rotation and % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9qq-f0-yqaqVeLsFr0-vr% 0-vr0db8meaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaaIWaWefv% 3ySLgznfgDOjdaryqr1ngBPrginfgDObcv39gaiuaacqWFKjcHcaWG% Ubaaaa!461C!\[0 \leqslant n\] <1 non-solid rotation). Thus the frequency of inertial oscillations is reduced for non-rigid rotation.