[1] In the present study, the Penn State/NCAR Mesoscale Model (MM5) was modified by considering recent advances in the urban boundary layer. In particular, the modifications were carried out in two directions: (1) With respect to the thermal properties of an urban surface the surface energy balance was modified by taking into account the anthropogenic heat released in urban areas and the urban heat storage term to account for urban/building mass effects, including hysteresis; and (2) the surface stress and fluxes of heat and momentum were modified following recent advances in the atmospheric boundary layer over rough surfaces under unstable conditions. The whole process was supplemented by detailed information on land use cover, derived from satellite image analysis. The modifications were applied to the high-resolution nonlocal medium-range forecast planetary boundary layer parameterization scheme, based on work by Troen and Mahrt (1986). The improvements seen with the modified model, after comparison with available measurements of temperature and fluxes, refer to (1) the strengthening of the nocturnal urban heat island; (2) the changes in the temperature, which proved to be favorable through the whole diurnal cycle, resulting in decreasing the temperature amplitude wave; (3) the decrease of turbulence and fluxes during the daytime; and (4) the diffusion coefficient and potential temperature profiles that are reduced during daytime and are increased at the lower levels during the night and thus affect accordingly the mixing height.

SUMMARYClose to the surface large coherent eddies consisting of plumes and downdraughts cause convergent winds blowing towards the plume axes, which in turn cause wind shears and generation of turbulence. This mechanism strongly enhances the convective heat/mass transfer at the surface and, in contrast to the classical formulation, implies an important role of the surface roughness. In this context we introduce the stability-dependence of the roughness length. The latter is important over very rough surfaces, when the height of the roughness elements becomes comparable with the large-eddy Monin-Obukhov length. A consistent theoretical model covering convective regimes over all types of natural surfaces, from the smooth still sea to the very rough city of Athens, is developed; it is also comprehensively validated against data from measurements at different sites and also through the convective boundary layer. Good correspondence between model results, field observations and large-eddy simulation is achieved over a wide range of surface roughness lengths and convective boundary-layer heights.

No abstract

SUMMARYThe concept of the minimum friction velocity is studied using three different methods of scalar averaging for the calculation of the stresses. Particular emphasis is given to the extraction of the in uence of a non-zero ambient wind shear observed in eld measurements. Data from three different experimental sites in Athens with high roughness values are analysed in order to provide information concerning the dependence of the dimensionless minimum friction velocity on the dimensionless roughness length. Data from the BOREX-95 experiment have also been re-analysed according to the methodologies presented in this study. The results are compared to the large-eddy simulations that are considered to be a reference study on shear-free convection.

A simple scheme of forcing turbulence away from decay was introduced by Lundgren some time ago, the "linear forcing," which amounts to a force term that is linear in the velocity field with a constant coefficient. The evolution of linearly forced turbulence toward a stationary final state, as indicated by direct numerical simulations (DNS), is examined from a theoretical point of view based on symmetry arguments. In order to follow closely the DNS, the flow is assumed to live in a cubic domain with periodic boundary conditions. The simplicity of the linear forcing scheme allows one to rewrite the problem as one of decaying turbulence with a decreasing viscosity. Its late-time behavior can then be studied by scaling symmetry considerations. The evolution of the system in the description of "decaying" turbulence can be understood as the gradual symmetry breaking of a larger approximate symmetry to a smaller symmetry that is exact at late times. The latter symmetry implies a stationary state: In the original description all correlators are constant in time, while, in the "decaying" turbulence description, that state possesses constant Reynolds number and integral length scale. The finiteness of the domain is intimately related to the evolution of the system to a stationary state at late times: In linear forcing there is no other large scale than the domain size, therefore, it is the only scale available to set the magnitude of the necessarily constant integral length scale in the stationary state. A high degree of local isotropy is implied by the late-time exact symmetry, the symmetries of the domain itself, and the solenoidal nature of the velocity field. The fluctuations observed in the DNS for all quantities in the stationary state can be associated with deviations from isotropy that is necessarily broken at the large scale by the finiteness of the domain. Indeed, to strengthen this conclusion somewhat, self-preserving isotropic turbulence models are used to study evolution from a direct dynamical point of view. Simultaneously, the naturalness of the Taylor microscale as a self-similarity scale in this system is emphasized. In this context the stationary state emerges as a stable fixed point. We also note that self-preservation seems to be the reason behind a noted similarity of the third-order structure function between the linearly forced and freely decaying turbulence, where, again, the finiteness of the domain plays a significant role.

In the present study, an attempt is made to assess the atmospheric boundary-layer (ABL) depth over an urban area, as derived from different ABL schemes employed by the mesoscale model MM5. Furthermore, the relationship of the mixing height, as depicted by the measurements, to the calculated ABL depth or other features of the ABL structure, is also examined. In particular, the diurnal evolution of ABL depth is examined over the greater Athens area, employing four different ABL schemes plus a modified version, whereby urban features are considered. Measurements for two selected days, when convective conditions prevailed and a strong sea-breeze cell developed, were used for comparison. It was found that the calculated eddy viscosity profile seems to better indicate the mixing height in both cases, where either a deep convective boundary layer develops, or a more confined internal boundary layer is formed. For the urban scheme, the incorporation of both anthropogenic and storage heat release provides promising results for urban applications.

Thermal stratification and its seasonal variations in Kouris Dam in Cyprus were simulated, and the impact of five different water withdrawal schemes was studied, using the 2-D, laterally averaged CE-QUAL-W2 reservoir model. Based on the model simulations, it was found that the thermal stratification of the reservoir is significant for most of the year. Most importantly, a complete mixing of the water column, triggered by seasonal variations in meteorological conditions, occurs in late-January. Predicted thermal stratification and water temperature profiles in the reservoir are noticeably affected by water withdrawal schemes. It was found that deep-water withdrawals tend to facilitate heat transfer in the water column and deepen the water mixing layer (epilimnion), especially from September to the following January. These study results suggest that it is prudent for Kouris Dam to integrate selective water withdrawal schemes into reservoir management by using the water withdrawal effects on thermal stratification for different water quality management strategies.

No abstract