Estimation of the Monin-Obukhov similarity functions from a spectral model

From measured one-dimensional spectra of velocity and temperature variance, the universal functions of the Monin-Obukhov similarity theory are calculated for the range –2 z/L + 2. The calculations show good agreement with observations with the exception of a range –1 z/L 0 in which the function _m , i.e., the nondimensional mean shear, is overestimated. This overestimation is shown to be caused by neglecting the spectral divergence of a vertical transport of turbulent kinetic energy. The integral of the spectral divergence over the entire wave number space is suggested to be negligibly small in comparison with production and dissipation of turbulent kinetic energy.Notation a,b,c contants (see Equations (–4)) - C_i constants i=u, v, w, (see Equation (5) - k_me,k_mT peak wave numbers of 3-d moel spectra of turbulent kinetic energy and of temperature variance, respectively - k_mi peak wave numbers of 1-d spectra of velocity components i=u, v, w and of temperature fluctuations i= - k_sb, k_c characteristics wave numbers of energy-feeding by mechanical effects being modified by mean buoyancy, and of convective energy feeding, respectively - L Monin-Obukhov length - % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% Gabeivayaaraaaaa!3C5B!\[{\rm{\bar T}}\] difference of mean temperature and mean potential temperature - T* Monin-Obukhov temperature scale - velocity of mean flow in positive x-direction - u* friction velocity - u, v, w components of velocity fluctuations - z height above ground - von Kármanán constant - temperature fluctuation - _m nondimensional mean shear - _H nondimensional mean temperature gradient - nondimensional rate of lolecular dissipation of turbulent kinetic energy - _D nondimensional divergence of vertical transports of turbulent linetic energy     

Further Studies Of Atmospheric Turbulence In Layers Near The Surface: Scaling The Tke Budget Above The Roughness Sublayer

The second of two experimental studies of the TKE budget conducted on sites of different roughness is described, and results are compared. The first took place within a shallow layer above a small field of mostly bare, cultivated soil; the second was carried out above a roughness sublayer of significant depth on an extensive plain of tall dry grass. Budget terms observed in the second study were scaled with a modified u which compensated for effects of an unusually large stress gradient and ensured that the _m functions would be collinear. By showing that the modification becomes negligible in smaller gradients, it is demonstrated that in normal conditions, budgets observed above significant roughness sublayers should be normalized by scaling in terms of the unreduced Reynolds stress at the sublayer's upper surface. This procedure is shown to be consistent with the expectation that TKE budgets in layers near the surface all scale in fundamentally the same way.Other findings include: (1) the fact that most _m functions previously reported are not quite collinear is attributed to a type of overspeeding known to affect three-cup anemometers; (2) revised _m functions, collinear and largely free of the effects of overspeeding, are determined from a well-established characteristic of the linear _m relation for the stable case; (3) data that define collinear _m functions can also be represented with single hyperbolic curves; (4) dissipation is found to be 10 to 15% too small to balance total TKE production in unstable and neutral conditions and to decrease with increasing z/L in thestable regime; and (5) new relations for based on the observed behaviour of the dissipation deficit provide an improved closure for the set of equations that express the budget terms as functions of _m and z/L.     

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.     

An Alternative Function For The Wind And Temperature Gradients In Unstable Surface Layers

The function ()=(1+|z/L|^2/3)^1/2,where z is the height, L the Obukhov length, and a constant,is proposed for the nondimensional wind speed and temperaturegradients (flux-profile relationships) in anunstable surface layer. This function agrees quite well withboth wind speed and temperature data,has the theoretically correct behaviour in convective conditions,and leads to simple results when integrated to produce the mean profiles.     

Turbulence in a Katabatic Flow

Turbulence measurements performed in a stable boundary layer over the sloping ice surface of the Vatnajökull in Iceland are described. The boundary layer, in which katabatic forces are stronger than the large-scale forces, has a structure that closely resembles that of a stable boundary layer overlying a flat land surface, although there are some important differences. In order to compare the two situations the set-up of the instruments on an ice cap in Iceland was reproduced on a flat grass surface at Cabauw, the Netherlands. Wind speed and temperature gradients were calculated and combined with flux measurements made with a sonic anemometer in order to obtain the local stability functions _m and _h as a function of the local stability parameter z/L. Unlike the situation at Cabauw, where _m was linear as a function of z/L, in the katabatically forced boundary layer, the dependence of _m on stability was found to be non-linear and related to the height of the wind maximum. Thermal stratification and the depth of the stable boundary layer however seem to be rather similar under these two different forcing conditions.Furthermore, measurements on the ice were used to construct the energy balance. These showed good agreement between observed melt and components contributing to the energy balance: net radiation (supplying 55% of the energy), sensible heat flux (30%) and latent heat flux (15%).Local sources and sinks in the turbulent kinetic energy budget are summed and indicate a reasonable balance in near-neutral conditions but not in more stable situations. The standard deviation of the velocity fluctuations _u, _v, and _w, can be scaled satisfactorily with the local friction velocity u_* and the standard deviation of the temperature fluctuation with the local temperature scale _*.     

Flux-profile Relationships for Wind Speed and Temperature in the Stable Atmospheric Boundary Layer

Wind and temperature profiles in the stable boundary layer were analyzed in the context of MoninObukhov similarity. The measurements were made on a 60-m tower in Kansas during October 1999 (CASES-99). Fluxprofile relationships, obtained from these measurements in their integral forms, were established for wind speed and temperature. Use of the integral forms eliminates the uncertainty and accuracy issues resulting from gradient computations. The corresponding stability functions, which were nearly the same for momentum and virtual sensible heat, were found to exhibit different features under weakly stable conditions compared to those under strongly stable conditions. The gradient stability functions were found to be linear, namely _m = 1+ 5.8 and _h = 1 + 5.4 up to a limit of the MoninObukhov stability parameter = 0.8; this is consistent with earlier findings. However, for stronger stabilities beyond a transition range, both functions were observed gradually to approach a constant, with a value of approximately 7. To link these two distinct regimes, a general but pliable functional form with only two parameters is proposed for the stability functions, covering the entire stability range from neutral to very stable conditions.     

Airflow above changes in surface heat flux,temperature and roughness; an extension to include the stable case

A numerical model of airflow above changes in surface roughness and thermal conditions is extended to include cases with stable thermal stratification within the internal boundary-layer. The model uses a mixing-length approach with empirical forms for M and H.Results are presented for some basic cases and an attempt is then made to compare results given by the model with the experimental results of Rider, Philip and Bradley. Tolerable agreement is achieved. The importance of roughness change and thermal stability effects in the diffusion of heat and moisture near a leading edge is emphasised.Notation A Refers to Taylor (1970) - B Businger-Dyer constant (= 16.0) in forms for _M and _H - C Constant in form for in stable case - c _p Specific heat at constant pressure - E Scaled absolute humidity - g Acceleration due to gravity - H Upward vertical heat flux - H ₀, H ₁ Surface heat fluxes for x <0, x0 - H _E Upward latent heat flux - k Von Kármán's constant (= 0.4) - K _H K _W Eddy transfer coefficients for heat and water vapour - L Monin-Obukhov length - L _H Latent heat of evaporation for water - m Ratio of roughness lengths ( = z ₁/z ₀) - RPB Refers to Rider et al. (1964) - RL* Non-dimensional parameter (see Equations (9), (20a), (22a), (24a)) - R* Net radiation less ground heat flux (see Equations (15), (16)) - T Scaled temperature - T ₁ Downstream scaled surface temperature - u ₀ u ₁(x) Surface friction velocities for x <0, x0 - U, W Horizontal and vertical mean velocities - x, z Horizontal and vertical co-ordinates - Z _i Local roughness length - z ₀, z _i Roughness lengths for x < 0, x 0 - Temperature - ₀, ₁ Surface temperatures for x<0, x0 - _E Non-dimensional absolute humidity gradient - _H Non-dimensional temperature gradient of heat flux - _M Non-dimensional wind shear - = _M = _H = _E an assumption used in stable conditions - Air density - Absolute humidity - _w Density of water - Kinematic shear stress - Logarithmic height scale (= ln(z+z ₁)/z ₁)     

Effects of heat and water vapor transport on eddy covariance measurement of CO2 fluxes

Flux densities of carbon dioxide were measured over an arid, vegetation-free surface by eddy covariance techniques and by a heat budget-profile method, in which CO₂ concentration gradients were specified in terms of mixing ratios. This method showed negligible fluxes of CO₂, consistent with the bareness of the experimental site, whereas the eddy covariance measurements indicated large downward fluxes of CO₂. These apparently conflicting observations are in quantitative agreement with the results of a recent theory which predicts that whenever there are vertical fluxes of sensible or latent heat, a mean vertical velocity is developed. This velocity causes a mean vertical convective mass flux (= _cw for CO₂, in standard notation). The eddy covariance technique neglects this mean convective flux and measures only the turbulent flux _c w. Thus, when the net flux of CO₂ is zero, the eddy covariance method indicates an apparent flux which is equal and opposite to the mean convective flux, i.e., _c w = – _c w. Corrections for the mean convective flux are particularly significant for CO₂ because _cw and _c w are often of similar magnitude. The correct measurement of the net CO₂ flux by eddy covariance techniques requires that the fluxes of sensible and latent heat be measured as well.     

On Trajectory Curvature as a Selection Criterion for Valid Lagrangian Stochastic Dispersion Models

Recently Wilson and Flesch (Boundary-Layer Meteorology, 84, 411-426, 1997) suggested that the average increment d z to the orientation = arctan(w/u) of the Lagrangian velocity-fluctuation vector can be used to distinguish the better Lagrangian stochastic models within the well-mixed class. Here it is demonstrated that the specification of d z constitutes neither a sufficient or universally applicable criterion to distinguish the better Lagrangian stochastic models within the well-mixed class. The hypothesis made by Wilson and Flesch that Lagrangian stochastic models with /P_E irrotational are zero-spin models, having d z=0, is proven     

Experiments On Buoyant Plume Dispersion In A Laboratory Convection Tank

Plume dispersion in the convective boundary layer (CBL) is investigated experimentally in a laboratory convection tank. The focusis on highly-buoyant plumes that loft near or become trapped in the CBL capping inversion and resistdownward mixing. Such plumes are defined by dimensionless buoyancy fluxes F_* 0.1, where F_* = F_b/(U w_* ² z_i), F_b is the stack buoyancy flux,U is the mean wind speed, w_* is the convective velocity scale, and z_i is the CBL depth. The aim is to obtain statistically-reliable mean (C) and root-mean-square (rms, _c) concentration fields as a function of F_* and the dimensionless distance X = w_*x/(U z_i), where x is the distance downstream of the source.The experiments reveal the following mainresults: (1) For 3 X 4and F_* 0.1, the crosswind-integrated concentration (CWIC) fields exhibit distinctly uniform profiles below z_i with a CWIC maximum aloft, in contrast to the nonuniform profiles obtained earlier by Willis and Deardorff. (2) The lateral dispersion (_y) variation with X is consistent with Taylor's theory for _* 0.1 and a buoyancy-enhanced dispersion, _y/z_i F_* ^1/3X^2/3, forF_* = 0.2 and 0.4. (3) The entrapment, the plume fraction above z_i, has a mean (E) that follows a systematic variationwith X and F_*, and a variability (_e/E) that is broad ( 0.3 to 2) near the source but subsides to 0.25 far downstream. (4) Vertical profiles of the concentration fluctuation intensity (c/C) are uniform for z < z_i and X > 1.5, but exhibit significant increases: (a) at the surface and close to the source (X 1.5), and(b) in the entrainment zone. (5) The cumulative distribution functions (CDFs) of the scaled concentration fluctuations (c/_c) separate into mixed-layer and entrainment-layer CDFs for X 2, with the mixed-layer group collapsing to a single distribution independent of z.These are the first experiments to obtain all components of the lateral and vertical dispersion parameters (rms meander, relative dispersion, total dispersion) for continuous buoyant releases in a convection tank. They also are the first tank experiments to demonstrate agreement with field observations of: (1) the scaled ground-level concentration along the plume centreline, and (2) the dimensionless lateral dispersion _y/z_i of buoyant plumes.     

A study of the sensitivity of bare soil evaporation schemes to soil surface wetness,using the coupled soil moisture and surface temperature prediction model,BARESOIL

Summary The performance of evaporation schemes with and approach and their combination within resistance representation of evaporation from bare soil surface is discussed. For this purpose nine schemes, based on different functions of or , on the ratio of the volumetric soil moisture content and its saturated value are used.The quality of the chosen schemes has been evaluated using the results of time integration by the coupled soil moisture and surface temperature prediction model, BARESOIL, using in situ data. A sensitivity analysis was made using two sets of data derived from the volumetric soil moisture content of the top soil layer. One with values below the wilting point (0.17 m³m^–3) and the second with values above 0.20m³m^–3. Data sets were obtained at the experimental site Rimski anevi, Yugoslavia, from the bare surface of a chernozem soil.With 4 Figures     

On the turbulence structure in the stable boundary layer over the Greenland ice sheet

Turbulence measurements from a 30 m tower in the stably stratifiedboundary layer over the Greenland ice sheet are analyzed. The observationsinclude profile and eddy-correlation measurements at various levels. Atfirst, the analysis of the turbulence data from the lowest level (2 m aboveground) shows that the linear form of the non-dimensional wind profile(m) is in good agreement with the observations for z/L <0.4, whereL represents the Obukhov length. A linear regression yieldsm=1+5.8z/L. The non-dimensional temperature profile (h) at the2m level shows no tendency to increase with increasing stability. The datafrom the upper levels of the tower are analyzed in terms of both localscaling and surface-layer scaling. The m and the h values show atendency to level off at large stability (z/>0.4) where represents the local Obukhov length. Hence, the linear form of the functions is no longer appropriate under such conditions. The bestcorrespondence to the data can be achieved when using the expression ofBeljaars and Holtslag for m and h. The vertical profiles of theturbulent fluxes, the wind velocity variances and temperature variance arealso determined. The momentum flux profile and the profiles of wind speedvariances are in general agreement with other observations if a welldeveloped low-level wind maximum occurs, and the height of this maximum isused as a height scale.     

Influence of heterogeneous land surfaces on surface energy and mass fluxes

Summary Land-surface heterogeneity affects surface energy fluxes. The magnitudes of selected land-surface influences are quantified by comparing observations with model simulations of the FIFE (First ISLSCP Field Experiment) domain. Several plausible heterogeneous and homogeneous initial and boundary conditions are examined, although soilmoisture variability is emphasized. It turns out that simple spatial averages of surface variation produced biased flux values. Simulated maximum latent-heat fluxes were approximately 30 to 40 W m^–2 higher, and air temperatures 0.4 °C lower (at noon), when computations were initialized with spatially averaged soil-moisture and leaf-area-index fields. The planetary boundary layer (PBL) height and turbulent exchanges were lower as well. It additionally was observed that (largely due to the nonlinear relationship between initial soil-moisture availability and the evapotranspiration rate), real latent-heat flux can be substantially less than simulated latent-heat flux using models initialized with spatially averaged soil-moisture fields. Differences between real and simulated fluxes also vary with the resolution at which real soil-moisture heterogeneity is discretized.With 8 Figures     

Thermal comfort of man in different urban environments

Summary On July 29, 1985, a hot summer day, biometeorological measurements were performed simultaneously in three different urban structures within the city of Munich and in the trunk space of a nearby tall spruce forest. Based on the results of these experiments the following thermophysiologically relevant biometeorological indices were calculated: Predicted mean vote, skin wettedness and physiologically equivalent temperature. These three indices are derived from different models for the human energy balance. They allow the assessment of the thermal components of the microclimates at the selected sites with regard to application in urban planning. The results quantitatively show the great heat stress in the urban structure street canyon, exposed to south, whereas in the trunk space of the tall spruce forest there is nearly an optimal climate even on hot summer days. Between these extremes the results for street canyon, exposed to north show a little higher heat load than for backyard with trees.

---

Zusammenfassung An einem heißen Sommertag, dem 29. Juli 1985, wurden in drei Stadtstrukturen in München und im Stammraum eines nahegelegenen Fichtenhochwaldes zeitgleich biometeorologische Messungen durchgeführt. Mit den Meßergebnissen wurden folgende thermophysiologisch relevante biometeorologische Indizes berechnet: Predicted mean vote, Hautbenetzungsgrad und physiologisch äquivalente Temperatur. Diese drei Indizes beruhen auf verschiedenen Modellen zur menschlichen Energiebilanz. Mit den drei Indizes wurden die thermischen Komponenten der Mikroklimate an den ausgewählten Meßplätzen im Hinblick auf Stadtplanungsaufgaben bewertet. Die Ergebnisse zeigen quantitativ die relativ große Hitzebelastung bei der Stadtstruktur Straßenschlucht, nach Süd exponiert, während im Stammraum des Fichtenhochwaldes selbst an heißen Sommertagen nahezu optimale Bedingungen herrschen. Zwischen diesen Extremen liegen die Ergebnisse für die anderen Meßplätze, wobei für Straßenschlucht, nach Nord exponiert die Wärmebelastung etwas höher als für Innenhof mit Bäumen ist.

---

---

With 6 Figures     

Bioaccumulation of Cobalt in Parmelia sulcata

Six locations across mainland Portugal were selected for exposing Parmelia sulcata, for a one-year period (8 months for one site), with simultaneous measurement of total (dry + wet) deposition (one-month periods). The exposed lichens and the total (dry + wet) deposition were analysed for cobalt contents by INAA (instrumental neutron activation analysis) and ICP-MS (inductively coupled plasma mass spectroscopy), respectively. The designated wet deposition was evaluated through the collected water volume; the designated dry deposition was assessed after the (dried) residual mass of the wet deposition. An excellent agreement between Co contents in exposed lichens and the cumulative (1) Co contents in the dry deposition, (2) dry deposition, and (3) wet deposition has been found for the locations with alternate drought and precipitation months, high dry deposition, and high Co contents in the latter. Continuous rainfall was found to hinder the Co accumulation in the lichen due to its release from the lichen and/or lower Co contents in the dry deposition. At three locations, P. sulcata Co contents, after subtraction of the background (before exposure), equalled or exceeded the Co contents in the cumulative dry deposition at the end of the exposure time. The optimal exposure period for this species likely depends on the exposure conditions.     

Turbulence mechanisms at an agricultural site

An extensive set of turbulence data from the 3- and 12-m heights taken over an agricultural site (Marsta, Sweden) are analyzed and compared with data from ideal sites.In unstable air, Monin-Obukhov similarity is found to be valid for the non-dimensional gradients of wind, _m , temperature, _h , and humidity, _e , for (only a few data), for _T /|T _*|,/ _E /|E _*| and for the non-dimensionalized inertial subrange spectra of temperature and humidity. Where comparison is possible, the unstable data also agree with those found in the Kansas study, with one remarkable exception, the inertial subrange constant of the temperature spectrum, ₁, being only 0.39, compared to the value 0.80 found at the Kansas site.On the stable side, most similarity predictions break down, with most of the data differing systematically from the corresponding Kansas results, the only exception being . The inertial subrange constants for temperature, ₁, and for humidity, ₁ are found to have the same values, 0.39 and 0.30, respectively, as they do on the unstable side. Remarkable similarity is found for the shape of the stable u- and - and e-spectra. In addition, this shape is found to be identical with that found in Kansas. The peak wavelength of the stable u-, and -spectra is found to be about four times larger than it is for the corresponding Kansas spectra. This is interpreted to be a result of the increased macro-roughness at the Marsta site as compared with that at the Kansas site. A possible explanation for the low ₁-value is discussed, suggesting that ₁ is not a universal constant, but instead dependent on the turbulent structure.     

Trace gas exchange at the air/water interface: Measurements of mass accommodation coefficients

A liquid jet of 90 m diameter and variable length has been utilized to determine absorption rates and, hence, mass accommodation coefficients , of atmospheric trace gases. The compounds investigated are HCl (0.01), HNO₃ (0.01), N₂O₅ (0.005), peroxyacetyl nitrate (>0.001), and HONO (0.005). It is concluded that the absorption of these trace gases by liquid atmospheric water is not significantly retarded by interfacial mass transport. The strengths and limitations of the liquid jet technique for measuring mass accommodation coefficients are explored.     

Similarity theory of diffusion and the observed vertical spread in the diabatic surface layer

The extension of Lagrangian similarity theory of diffusion to stratified flow is examined, to improve its prediction of the vertical spread of a passive substance. In the basic equation, where is the average height of a cluster of particles,u _* is the friction velocity andL is Monin-Obukhov length. It is shown theoretically, under the assumption of an equivalence between the diffusivities of heat and matter, that the unspecified function is the reciprocal of a more familiar meteorological parameter _n , the dimensionless temperature gradient. The universal constantb is found to be approximately equal to von Karman's constant for various stability conditions. The predicted effect of stability on vertical spread shows excellent agreement with that of the published data from the O'Neill experiments.     

Northern hemisphere atmospheric response to changes of atlantic ocean SST on decadal time scales: a GCM experiment

Analyses indicate that the Atlantic Ocean seasurface temperature (SST) was considerably colder at the beginning than in the middle of the century. In parallel, a systematic change in the North Atlantic sea-level pressure (SLP) pattern was observed. To find out whether the SST and SLP changes analyzed are consistent, which would indicate that the SST change was real and not an instrumental artifact, a response experiment with a low-resolution (T21) atmospheric GCM was performed. Two perpetual January simulations were conducted, which differ solely in the Atlantic Ocean (40° S-60° N) SST: the cold simulation utilizes the SSTs for the period 1904–1913; the warm simulation uses the SSTs for the period 1951–1960. Also, a control run with the model's standard SST somewhat between the cold and warm SST was made. For the response analysis, a rigorous statistical approach was taken. First, the null hypothesis of identical horizontal distributions was subjected to a multivariate significance test. Second, the level of recurrence was estimated. The multivariate statistical approaches are based on hierarchies of test models. We examined three different hierarchies: a scale-dependent hierarchy based on spherical harmonics (S), and two physically motivated ones, one based on the barotropic normal modes of the mean 300 hPa flow (B) and one based on the eigenmodes of the advection diffusion operator at 1000 hPa (A). The intercomparison of the cold and warm experiments indicates a signal in the geostrophic stream function that in the S-hierarchy is significantly nonzero and highly recurrent. In the A-hierarchy, the low level temperature field is identified as being significantly and recurrently affected by the altered SST distribution. The SLP signal is reasonably similar to the SLP change observed. Unexpectedly, the upper level stream-function signal does not appear to be significantly nonzero in the B-hierarchy. If, however, the pairs of experiments warm versus control and cold versus control are examined in the B-hierarchy, a highly significant and recurrent signal emerges. We conclude that the cold versus warm response is not a small disturbance that would allow the signal to be described by eigenmodes of the linear system. An analysis of the three-dimensional structure of the signal leads to the hypothesis that two different mechanisms are acting to modify the model's mean state. At low levels, local heating and advection are dominant, but at upper levels the extratropical signal is a remote responce to modifications of the tropical convection.This paper was presented at the International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 11–15 September 1989 under the auspices of the Meteorological Institute of the University of Hamburg and the Max Planck Institute for Meteorology. Guest Editor for these papers is Dr. L. Dilmenil.AWI Publication no. 254