Etude de la couche limite convective sahelienne en presence de brumes seches (Expérience ECLATS)

The ECLATS experiment was conducted in order to investigate the influence of radiative processes on the dynamics of the atmospheric boundary layer during its diurnal evolution. This experiment was carried out over Niger, near Niamey, by measuring continuously the energy balance at ground level and by using an instrumented aircraft for turbulence, radiative fluxes and aerosol measurements in the boundary layer during dusty conditions (brumes sèches). This paper is restricted to an analysis of the turbulent structure in the homogeneous and stationary convective boundary layer. The turbulence moments for kinetic energy and the spectral characteristics of the vertical velocity are discussed. These results are compared with a set of data obtained for clear convective boundary layers. The differences observed are quite important and seem, at least in part, due to radiative processes (infrared radiative divergence in the surface layer and absorption of solar radiation in the boundary layer).     

Instabilities of the tidally induced bottom boundary layer in the rotating frame and their mixing effect

To investigate the stability of the bottom boundary layer induced by tidal flow (oscillating flow) in a rotating frame, numerical experiments have been carried out with a two-dimensional non-hydrostatic model. Under homogeneous conditions three types of instability are found depending on the temporal Rossby number Ro_t, the ratio of the inertial and tidal periods. When Ro_t < 0.9 (subinertial range), the Ekman type I instability occurs because the effect of rotation is dominant though the flow becomes more stable than the steady Ekman flow with increasing Ro_t. When Ro_t > 1.1 (superinertial range), the Stokes layer instability is excited as in the absence of rotation. When 0.9 < Ro_t < 1.1 (near-inertial range), the Ekman type I or type II instability appears as in the steady Ekman layer. Being much thickened (100 m), the boundary layer becomes unstable even if tidal flow is weak (5 cm/s). The large vertical scale enhances the contribution of the Coriolis effect to destabilization, so that the type II instability tends to appear when Ro_t > 1.0. However, when Ro_t < 1.0, the type I instability rather than the type II instability appears because the downward phase change of tidal flow acts to suppress the latter. To evaluate the mixing effect of these instabilities, some experiments have been executed under a weak stratification peculiar to polar oceans (the buoyancy frequency N²  10⁻⁶ s⁻²). Strong mixing occurs in the subinertial and near-inertial ranges such that tracer is well mixed in the boundary layer and an apparent diffusivity there is evaluated at 150–300 cm²/s. This suggests that effective mixing due to these instabilities may play an important role in determining the properties of dense shelf water in the polar regions.     

Impact of Secular Climate Change on the Thermal Structure of a Large Temperate Central European Lake

Strong climate-related secular trends are apparent in a 52-yr long (1947–1998) uninterrupted series of monthly temperature profiles fromLake Zurich, a large, deep (136 m), temperate lake on the Swiss Plateau. Decadal mean water temperatures have undergone a secular increase at all depths, reflecting the high degree of regional warming that occurred in the European Alpine area during the 20th century. From the 1950s to the 1990s, high warming rates ( 0.24 K per decade) in the uppermost 20 m of the lake (i.e., the epi/metalimnion) combined with lower warming rates ( 0.13 K per decade) below 20 m (i.e., in the hypolimnion), have resulted in a20% increase in thermal stability and a consequent extension of 2–3 weeksin the stratification period. In common with many other parts of the world, 20th-century climate change on the Swiss Plateau has involved a steep secular increase in daily minimum (nighttime) air temperatures, but not in daily maximum (daytime) air temperatures. With respect to both secular change and decadal-scale variability, the temporal structure of the temperature of the surface mixed layer of Lake Zurich faithfully reflects that of the regional daily minimum air temperature, but not that of the daily maximum. The processes responsible for longer-term changes in the temperature structure of the lake therefore act during the night, presumably by suppressing nighttime convective cooling of the surface mixed layer. Application of a one-box heat exchange model suggests that the observed secular changes in thermal structure are due to shifts in the nighttime rate of emission of infrared radiation from the atmosphere and in the nighttime rates of latent and sensible heat exchange at the air-water interface. The increase in hypolimnetic temperatures is mainly a result of the increased prevalence of warm winters in Europe.     

Heat fluxes and roll circulations over the western Gulf Stream during an intense cold-air outbreak

Turbulence and heat fluxes in the marine atmospheric boundary layer (MABL) for the roll vortex regime, observed during the Genesis of Atlantic Lows Experiment (GALE) over the western Gulf Stream, have been studied. The spectral analysis suggests that cloud streets (roll vortices) are vertically organized convection in the MABL having the same roll scale for both the cloud layer and subcloud layer, and that the roll spacing is about three times the MABL depth. The roll circulations contribute significantly to the sensible (temperature) and latent heat (moisture) fluxes with importance increasing upward. Near the MABL top, these fluxes are primarily due to roll vortices which transfer both sensible heat and moisture upward in the lower half of the convective MABL. Near the MABL top, the roll circulations transfer sensible heat downward and moisture upward in the clear thermal-street region, but roll vortices influenced by evaporative cooling can transfer sensible heat upward and moisture downward in the cloud-street region. Near the cloud-top, the upward buoyancy flux due to evaporative cooling is highly related to the roll circulations near the inversion.For the lower half of the MABL, the normalized temperature flux decreases upward more rapidly than the humidity flux, which is mainly because potential temperature () increases slightly upward while humidity (q) decreases slightly upward above the unstable surface layer. The gradient production (associated with the gradient) is a source for the temperature flux in the unstable surface layer but changes to a sink in the mixed layer, while the gradient production (associated with the q gradient) acts as a source for the humidity flux in both the unstable surface and mixed layers. The results suggest that the entrainment at the MABL top might affect the budgets of temperature and humidity fluxes in the lower MABL, but not in the unstable surface layer.Caelum Research Corporation, Silver Spring, MD, 20901, U.S.A.     

A numerical study of the influence of an air temperature-inversion layer and a seawater density-jump layer on the structure of interacting boundary layers

The effects of an air-temperature inversion in the atmosphere and a seawater density jump in the ocean on the structure of the atmospheric and oceanic boundary layers are studied by use of a coupled model. The numerical model consists of a closed system of equations for velocities, turbulent kinetic energy, turbulent exchange coefficient, local turbulent length scale, and stratification expressions for both air and sea boundary layers. The effects of the temperature inversion and the density jump are incorporated into the equations of turbulent kinetic energy of the atmosphere and ocean by a parameterization. A series of numerical experiments was conducted to determine the effects of various strengths of the inversion layer and surface heat fluxes in the atmosphere and of the density-jump layer in the ocean on the structure of the interacting boundary layers.The numerical results show that the temperature inversion in the atmosphere and density jump in the ocean have strong influences on turbulent structure [especially on the turbulent exchange coefficient (TEC) and turbulent kinetic energy (TKE)] and on air-sea interaction characteristics. Maxima of TKE and TEC strongly decrease with increasing strength of the inversion layer, and they disappear for strong inversions in the atmosphere. Certain strengths (density differences between the upper and the lower layers) of the density-jump layer in the ocean (₂ 0.1 g/cm³) produce double maxima in TEC-profiles and TKE-profiles in the ocean. The magnitudes of air-sea interaction characteristics such as geostrophic drag coefficient, and surface drift current increase with increasing strength of the density-jump layer in the ocean. The density-jump layer plays the role of a barrier that limits vertical mixing in the ocean. The numerical results agree well with available observed data and accepted quantitive understanding of the influences of a temperature inversion layer and a density-jump layer on the interacting atmospheric and oceanic boundary layers.     

A model of airflow above changes in surface heat flux,temperature and roughness for neutral and unstable conditions

A numerical model of airflow in the lowest 50–100 m of the atmosphere above changes in surface roughness and temperature or heat flux has been developed based on boundary layer approximations, the Businger-Dyer hypotheses for the non-dimensional wind shear and heat flux and a mixing length hypothesis.Results have been obtained for several situations, in particular, airflow with neutral upstream conditions encountering a step change in surface temperature or heat flux with no roughness change. In these cases large increases in shear stress at the outer edge of the internal boundary layer are predicted. The case of unstable upstream flow encountering a step change to zero heat flux is also considered.Two situations that may be encountered near the shores of the Great Lakes are considered.Notation B Businger-Dyer constant (= 16.0) in form for _M, _H - c _p Specific heat at constant pressure - g Acceleration due to gravity - H Upward vertical heat flux - H ₀ , H ₁ Surface heat fluxes for x < 0, x 0 - k von Kármán's constant ( = 0.4) - l Mixing length - L Monin-Obukhov length - L ₀ Upstream value of L - m Ratio of roughness lengths (= z ₁/z ₀) - RL ^* Non-dimensional parameter, see Equations (20, 22 and 24) - RL ₁ ^* Same as RL ^* but with z ₁ scaling (= mRL ^*) - T Scaled temperature - T ₀ (z) Upstream temperature profile - u ₀, u ₁(x) Surface friction velocities for x < 0, x 0 - U, W Horizontal and vertical mean velocities - U ₀ (z) Upstream velocity profile - x, z Horizontal and vertical coordinates - z _i Local roughness length     

The turbulent heat flux from arctic leads

The turbulent transfer of heat from Arctic leads in winter is one of the largest terms in the Arctic heat budget. Results from the AIDJEX Lead Experiment (ALEX) suggest that the sensible component of this turbulent heat flux can be predicted from bulk quantities. Both the exponential relation N = 0.14R _x ^0.72 and the linear relation N = 1.6 × 10^–3 R _x+ 1400 fit our data well. In these, N is the Nusselt number formed with the integrated surface heat flux, and R _x is the Reynolds number based on fetch across the lead. Because of the similarity between heat and moisture transfer, these equations also predict the latent heat flux. Over leads in winter, the sensible heat flux is two to four times larger than the latent heat flux.The internal boundary layer (IBL) that develops when cold air encounters the relatively warm lead is most evident in the modified downwind temperature profiles. The height of this boundary layer, , depends on the fetch, x, on the surface roughness of the lead, z ₀ and on both downwind and upwind stability. A tentative, empirical model for boundary layer growth is % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWcaaqaaiabes% 7aKbqaaiaadQhadaWgaaWcbaGaaGimaaqabaaaaOGaeyypa0JaeqOS% di2aaeWaaeaacqGHsisldaWcaaqaaiaadQhadaWgaaWcbaGaaGimaa% qabaaakeaacaWGmbaaaaGaayjkaiaawMcaamaaCaaaleqabaGaaGim% aiaac6cacaaI4aaaaOWaaeWaaeaadaWcaaqaaiaadIhaaeaacaWG6b% WaaSbaaSqaaiaaicdaaeqaaaaaaOGaayjkaiaawMcaamaaCaaaleqa% baGaaGimaiaac6cacaaI0aaaaaaa!472D!\[\frac{\delta }{{z_0 }} = \beta \left( { - \frac{{z_0 }}{L}} \right)^{0.8} \left( {\frac{x}{{z_0 }}} \right)^{0.4} \] where L is the Obukhov length based on the values of the momentum and sensible heat fluxes at the surface of the lead, and is a constant reflecting upwind stability.Velocity profiles over leads are also affected by the surface nonhomogeneity. Besides being warmer than the upwind ice, the surface of the lead is usually somewhat rougher. The velocity profiles therefore tend to decelerate near the surface, accelerate in the mid-region of the IBL because of the intense mixing driven by the upward heat flux, and rejoin the upwind profiles above the boundary layer. The profiles thus have distinctly different shapes for stable and unstable upwind conditions.     

Sodar Detected Top-Down Convection in a Nocturnal Cloud-Topped Boundary Layer: A Case Study

Nocturnal convection, originating in a well-mixed marine cloud-topped boundary layer, advected onshore, was observed using a Doppler sodar on the Tyrrhenian coast in Italy. The horizontal and vertical dimensions of the downdrafts were evaluated. The oscillation frequency triggered by the downdrafts at the inversion layer, derived from the harmonic analysis of the sodar measured vertical velocity (w), is compared with the Brunt-Vaisala frequency, obtained from the rawinsonde temperature profile. A similarity function for the ²_w vertical profile was used to fit the sodar experimental data and to retrieve the depth of the mixing layer and the sensible heat flux at the top of the cloud layer. The results are in agreement with the convection layer depth observed in the sodar echoes facsimile record, and with the energy budget evaluated at the top of the cloud layer using the rawinsonde profiles.     

Mixed-layer heat advection and entrainment during the sea breeze

A model is developed to simulate the potential temperature and the height of the mixed layer under advection conditions. It includes analytic expressions for the effects of mixed-layer conditions upwind of the interface between two different surfaces on the development of the mixed layer downwind from the interface. Model performance is evaluated against tethersonde data obtained on two summer days during sea breeze flow in Vancouver, Canada. It is found that the mixed-layer height and temperature over the ocean has a small but noticeable effect on the development of the mixed layer observed 10 km inland from the coast. For these two clear days, the subsidence velocity at the inversion base capping the mixed layer is estimated to be about 30 mm s^–1 from late morning to late afternoon. When the effects of subsidence are included in the model, the mixed-layer height is considerably underpredicted, while the prediction for the mean potential temperature in the mixed layer is considerably improved. Good predictions for both height and temperature can be obtained when values for the heat entrainment ratio,c, 0.44 and 0.68 for these two days respectively for the period from 1000 to 1300 LAT, were used. These values are estimated using an equation including the additional effects on heat entrainment due to the mechanical mixing caused by wind shear at the top of the mixed layer and surface friction. The contribution of wind shear to entrainment was equal to, or greater than, that from buoyant convection resulting from the surface heat flux. Strong wind shear occurred near the top of the mixed layer between the lower level inland flow and the return flow aloft in the sea breeze circulation.Symbols c entrainment parameter for sensible heat - c _p specific heat of air at constant pressure, 1010 J kg^–1 K^–1 - d ₁ the thickness of velocity shear at the mixed-layer top, m - Q _H surface sensible heat flux, W m^–2 - u _m mean mixed-layer wind speed, m s^–1 - u _* friction velocity at the surface, m s^–1 - w subsidence velocity, m s^–1 - W subsidence warming,^oC s^–1 - w _e entrainment velocity, m s^–1 - w _* convection velocity in the mixed layer, m s^–1 - x downwind horizontal distance from the water-land interface, m - y dummy variable forx, m - Z height above the surface, m - Z _i height of capping inversion, m - Z _m mixed-layer depth, i.e.,Z _i–Z_s, m - Z _s height of the surface layer, m - lapse rate of potential temperature aboveZ _i, K m^–1 - potential temperature step atZ _i, K - u _h velocity step change at the mixed-layer top - _m mean mixed-layer potential temperature, K     

Surface renewal analysis for sensible and latent heat flux density

High frequency temperature measurements were recorded at five heights and surface renewal (SR) analysis was used to estimate sensible heat flux density (H) over 0.1 m tall grass. Traces of the temperature data showed ramp-like structures, and the mean amplitude and duration of these ramps were used to calculate H using structure functions. Data were compared with H values measured with a sonic anemometer. Latent heat flux density (E) was calculated using an energy balance and the results were compared with E computed from the sonic anemometer data. SR analysis provided good estimates of H for data recorded at all heights but the canopy top and at the highest measurement level, which was above the fully adjusted boundary layer.     

The structure of the stably stratified internal boundary layer in offshore flow over the sea

Observations obtained mainly from a research aircraft are presented of the mean and turbulent structure of the stably stratified internal boundary layer (IBL) over the sea formed by warm air advection from land to sea. The potential temperature and humidity fields reveal the vertical extent of the IBL, for fetches out to several hundred of kilometres, geostrophic winds of 20–25 m s^–1, and potential temperature differences between undisturbed continental air and the sea surface of 7 to 17 K. The dependence of IBL depth on these external parameters is discussed in the context of the numerical results of Garratt (1987), and some discrepancies are noted.Wind observations show the development of a low-level wind maximum (wind component normal to the coast) and rotation of the wind to smaller cross-isobar flow angles. Potential temperature () profiles within the IBL reveal quite a different structure to that found in the nocturnal boundary layer (NBL) over land. Over the sea, profiles have large positive curvature with vertical gradients increasing monotonically with height; this reflects the dominance of turbulent cooling within the layer. The behaviour is consistent with known behaviour in the NBL over land where curvature becomes negative (vertical gradients of decreasing with height) as radiative cooling becomes dominant.Turbulent properties are discussed in terms of non-dimensional quantities, normalised by the surface friction velocity, as functions of normalised height using the IBL depth. Vertical profiles of these and the normalised wavelength of the spectral maximum agree well with known results for the stable boundary layer over land (Caughey et al., 1979).     

Growth Of The Summer Daytime Convective Boundary Layer At Anand

The heights of the daytime convective boundary layer (CBL), computed by a one-dimensional model for a bare soil surface at a semi-arid station,Anand, during the dry and hot summer month of May 1997, are presented. As input, the model requires surface heat flux, friction velocity and air temperature as functions of time. Temperature data at the one-metre level from a tower and sonic anemometer data at 9.5 m collected during the period 13–17 May 1997 in the Land Surface Processes Experiment (LASPEX-97) are used to compute hourly values of surface heat flux, friction velocity and Obukhov length following the operational method suggested by Holtslag and Van Ulden [J. Climate Appl. Meteorol. 22,517–529 (1983)]. The model has been tested with different values for the potential temperature gradient ( ) above the inversion. The model-estimated CBL heights comparefavourably with observed heights obtained from radiosonde ascents.     

The dependence of the measured cool skin of the ocean on wind stress and total heat flux

The temperature drop T between the ocean surface and the 5-cm depth was recorded during GATE, Phase III. With measured values of the total heat flux Q and an assumption about the thickness of the viscous boundary layer of the ocean, the wind-speed dependence of the factor of proportionality between T and Q is determined. This factor depends on the deviations of the thickness of the conductive layer from the thickness of the viscous layer and possibly partially on the wind stress. A further assumption about the thickness of the conductive layer leads to a wind-speed dependence of the ratio between total wind stress and its wave supporting part of it. This ratio increases from a value 1.5 at 1 m s^–1 to 9 at 10 m s^–1, which is in agreement with existing estimates.     

On Monin–Obukhov Similarity In The Stable Atmospheric Boundary Layer

Atmospheric measurements from several field experiments have been combined to develop a better understanding of the turbulence structure of the stable atmospheric boundary layer. Fast response wind velocity and temperature data have been recorded using 3-dimensional sonic anemometers, placed at severalheights (1 m to 4.3 m) above the ground. The measurements wereused to calculate the standard deviations of the three components of the windvelocity, temperature, turbulent kinetic energy (TKE) dissipation andtemperature variance dissipation. These data were normalized and plottedaccording to Monin–Obukhov similarity theory. The non-dimensional turbulencestatistics have been computed, in part, to investigate the generalapplicability of the concept of z-less stratification for stable conditions. From the analysis of a data set covering almost five orders ofmagnitude in the stability parameter = z/L (from near-neutral tovery stable atmospheric stability), it was found that this concept does nothold in general. It was only for the non-dimensional standard deviation oftemperature and the average dissipation rate of turbulent kinetic energythat z-less behaviour has been found. The other variables studied here(non-dimensional standard deviations of u, v, and w velocity components and dissipation of temperature variance) did not follow the concept of z-less stratification for the very stable atmospheric boundary layer. An imbalance between production and dissipation of TKE was found for the near-neutral limit approached from the stable regime, which matches with previous results for near-neutral stability approached from the unstable regime.     

Flux–Profile Relationships Over a Fetch Limited Beech Forest

The influence of an internal boundary layer and a roughness sublayer on flux–profile relationships for momentum and sensible heat have been investigated for a closed beech forest canopy with limited fetch conditions. The influence was quantified by derivation of local scaling functions for sensible heat flux and momentum (_h and _m) and analysed as a function of atmospheric stability and fetch. For heat, the influences of the roughness sublayer and the internal boundary layer were in agreement with previous studies. For momentum, the strong vertical gradient of the flow just above the canopy top for some wind sectors led to an increase in _m, a feature that has not previously been observed. For a fetch of 500 m over the beech forest during neutral atmospheric conditions, there is no height range at the site where profiles can be expected to be logarithmic with respect to the local surface. The different influence of the roughness sublayer on _h and _m is reflected in the aerodynamic resistance for the site. The aerodynamic resistance for sensible heat is considerably smaller than the corresponding value for momentum.     

Similarity scaling applied to sodar observations of the convective boundary layer above an irregular hill

Nine profiles of the temperature structure parameter C _T ² and the standard deviation of vertical velocity fluctuations ( _w) in the convective boundary layer (CBL) were obtained with a monostatic Doppler sodar during the second intensive field campaign of the First ISLSCP Field Experiment in 1987. The results were analyzed by using local similarity theory. Local similarity curves depend on four parameters: the height of the mixed layer (z _i ), the depth of the interfacial layer (), and the temperature fluxes at the top of the mixed layer (Q _i ) and the surface (Q _o). Values of these parameters were inferred from sodar data by using the similarity curve for C _T ² and observations at three points in its profile. The effects of entrainment processes on the profiles of C _T ² and _wnear the top of the CBL appeared to be described well by local similarity theory. Inferred estimates of surface temperature flux, however, were underestimated in comparison to fluxes measured by eddy correlation. The measured values of _wappeared to be slightly smaller than estimates based on available parmeterizations. These discrepancies might have been caused by experimental error or, more likely, by the distortion of turbulence structure above the site by flow over the nonuniform terrain at the observation site.     

Model of convective heat exchange due to isolated thermals in the atmospheric boundary layer

A numerical model of convective heat transfer due to isolated thermals in the atmospheric boundary layer is used to describe the temperature profile transformation in undisturbed conditions as a result of intensive dry free convection. Based on some assumptions, the heat transfer Equation (2) is transformed to the form (14) in which the coefficients and the function F are expressed by (d/dz)(ln ) and by parameters of thermals. Equation (14) has been solved numerically with the help of Equation (15) obtained from the statics equation because of Equation (8). The size distribution function f(z, r, t) of the thermals is discrete (Table I), according to Vulf'son (1961). On Figures 1 and 2 are plotted successive temperature profiles for a ground inversion, transformed due to free convection and turbulence (Figures 1a and 2a), and due to turbulence only (Figures 1b and 2b). The profiles are computed from Equation 14 (Figures 1a and 2a) and Equation 16 (Figures 1b and 2b) for k _z= 1 m² s^–1 (Figure 1) and k _z= 10 m² s^–1 (Figure 2). On Figure 3 the real temperature profiles in Sofia for June 22nd 1976 are compared with the profiles computed using the real initial profile for 4.30 h local time. Good qualitative agreement can be seen.     

A bulk similarity approach in the atmospheric boundary layer using radiometric skin temperature to determine regional surface fluxes

Profiles of wind velocity and temperature in the outer region of the atmospheric boundary layer (ABL) were used together with surface temperature measurements, to determine regional shear stress and sensible heat flux by means of transfer parameterizations on the basis of bulk similarity. The profiles were measured by means of radiosondes and the surface temperatures by infrared radiation thermometry over hilly prairie terrain in northeastern Kansas during the First ISLSCP Field Experiment (FIFE). In the analysis, the needed similarity functions were determined and tested; the main scaling variables used for the ABL were h _i , the height of the convectively mixed layer, and V _a and _a, the wind speed and potential temperature averaged over the mixed layer. Good agreement (r = 0.80) was obtained between values of friction velocity u _* determined by this ABL bulk similarity approach and those obtained by Monin-Obukhov similarity in the surface sublayer. Similarly, values of surface flux of sensible heat H determined by this method compared well (r = 0.90) with the regional means measured at six ground stations. The corresponding regional evaporation values, determined with the energy budget equation, also compared favorably (r = 0.94).     

Numerical studies of heat island circulations

A two dimensional model has been set up to investigate the circulation induced by an urban heat island in the absence of synoptic winds. The boundary conditions need to be formulated carefully and due to difficulties arising here, we restrict our attention to cases of initially stable thermal stratification. Heat island circulations are allowed to develop from rest and prior to the appearance of the final symmetric double cell pattern, a transitional multi-cell pattern is observed in some cases. The influence on the steady state circulation of various parameters is studied, among which are eddy transfer coefficients, the heat island intensity, the initial temperature stratification and the heat island size. Some results are presented for a case in which differential surface cooling beneath an initially stable atmosphere produces a circulation and an unstable layer capped by an elevated inversion over the city. It is hoped that this case is vaguely representative of the night-time heat island with no geostrophic wind.Notation c_p Specific heat at constant pressure - g Acceleration due to gravity - H Top of integration region - K_z Vertical eddy transfer coefficient - K_x, K_xH, K_xm Horizontal eddy transfer coefficients for heat and momentum - l ixing length - p Pressure - p₀ Reference surface pressure (1000 mb) - P_H (x, t) Pressure at z = H - R Specific gas constant for dry air - t Time - u, w Horizontal and vertical velocities - x, z Horizontal and vertical coordinates - x₁, x₂ Positions of discontinuities in surface temperature field (see Figure 2) - x_a Heat island half-width - x_b Boundary of integration region - Parameter in formula for eddy coefficients (variable-K case) = 18.0 - _s Intensity of heat island - Potential temperature field - Reference absolute temperature (variable-K case) - _r Reference temperature (° C) - _s Surface temperature - Q Air density     

A generalized scaling for convective shear flows

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