Aerodynamic roughness over an inhomogeneous ground surface

The aerodynamic roughness parameter z ₀ over inhomogeneous ground surfaces, such as cities, rural towns and so on, is determined by analyzing the wind data at AMeDAS observatories in the Tohoku and Kanto districts of Japan, by making use of Rossby number similarity theory. It is found that the aerodynamic roughness parameter is proportional to the average size of the roughness elements.A practical method of estimating the aerodynamic roughness parameter over an extensive area with various inhomogeneities is developed. In this method, the Digital National Land Information data bank is employed. As an example, the roughness parameter distribution around Tsukuba Academic City is presented.     

Simultaneous measurements of turbulence over land and water

Turbulence and mean flow variables under unstable conditions are examined with special emphasis on the consequences of roughness and surface elevation change. An interpolation formula for _w ², between neutral and free convection, is shown to bring order to the data. The spectral distribution of vertical wind variance is found to be in good agreement with results over horizontal homogeneous terrain, both with respect to form and position. In particular, the length scale _m corresponding to the maximum of nS _w(n) is unchanged. Another turbulent length scale, (k/)z, is found to be substantially reduced in the disturbed zone of the internal boundary layer. To a first approximation, the flow-acceleration effect on the non-dimensional wind shear can be separated from the diabatic effect.     

Characteristics of secondary circulations over an inhomogeneous surface simulated with large-eddy simulation

Large-eddy simulation is used to study secondary circulations in the convective boundary layer modulated as a result of horizontally varying surface properties and surface heat fluxes over flat terrain. The presence of heat flux heterogeneity and its alignment with respect to geostrophic wind influences the formation, strength and orientation of organized thermals. Results show boundary-attached roll formation along heat flux maxima in the streamwise direction. The streamwise organization of the updrafts and downdrafts formed downwind of heterogeneities leads to counter-rotating secondary circulations in the crosswind plane. The distribution of resolved-scale pressure deviations shows large pressure gradients in the crosswind plane. Spanwise and vertical velocity variances and heat flux profiles depict considerable spatial variability compared to a homogeneous forest simulation. Secondary circulations are observed for various ambient wind scenarios parallel and perpendicular to heterogeneities. In the presence of increased wind speed, thermals emerging from the heat flux heterogeneity are elongated, and organize along and downwind of large-scale heterogeneity in the streamwise direction. Simulation with a reduced heat flux shows a shallower circulation with a lower aspect ratio. Point measurements of heat flux inside the roll circulation could be overestimated by up to 15–25% compared to a homogeneous case.     

Determination of evaporation from integrated profiles of humidity and temperature over an inhomogeneous surface

In this study a moderate-sized alfalfa field was downwind of a large dry region. Measured vertical profiles of temperature, humidity and wind at upwind and downwind sides of the field were used to calculate the short-term evaporation rate, as well as the contribution of horizontally transported or advected heat energy to the evaporation. The vertical profiles must be measured at least to the height at which air is modified by the new surface. In this case that height was as large as 18 m over a several hundred meter long traverse.Evaporation rates calculated by such an approach were in very close agreement to surface fluxes measured by an eddy correlation system near the surface. The difference between calculated and measured values averaged 9.5%. The reduction in sensible heat content of the air of the control volume was substantial during passage over the field. If all this energy was assumed to have been used in evaporation, then the advection of heat contributed from 35 to 86% of the total evaporation rate.It appears that for an inhomogeneous surface, knowledge of the distortion of some properties of local boundary layers can yield reliable estimates of local evaporation. The approach has little empiricism as it is based on simple conservation laws.     

The breakup of marine boundary-layer clouds over an inhomogeneous sea surface temperature field

This numerical study examines the breakup of marine atmospheric boundary-layer (MABL) clouds through various physical processes over an inhomogeneous sea surface temperature (SST) field. Three regimes are identified under which the cloud layer will break up. (A) advection of drier air into the MABL for the California case. (B) daytime absorption of solar radiation, occurring most easily over the cold water. (C) mesoscale fluctuations in the flow, producing holes in the cloud layer.The budget study of these three situations concludes that large-scale subsidence, solar radiation, local mesoscale advection, and inhomogeneous surface fluxes cannot be neglected in modeling cloud breakup. This study also confirms the belief that the mixing process alone induced by evaporative entrainment is generally insufficient to predict the breakup of the cloud layer.Sections of this paper are based on an extended abstract by the author and Dr. Steven Stage for the Ninth Symposium on Turbulence and Diffusion held at Riso, Denmark, 1990.     

Boundary layer structure over an inhomogeneous surface: Simulation with a non-hydrostatic mesoscale model

A three-dimensional, non-hydrostatic mesoscale model is used to study boundary-layer structure over an area characterized by the city of Copenhagen, the Øresund strait, and adjacent coastal farmland. Simulations are compared with data obtained on June 5, 1984 during the Øresund experiment.Under moderately strong wind conditions, a stable internal boundary layer (IBL) developed over the Øresund strait during the day. Near-surface winds decelerate over water due to diminished vertical momentum transfer.The turbulent kinetic energy field closely reflects the surface roughness distribution due to the imposed relatively strong wind forcing. TKE budgets over water, farmland and a city area are discussed by inspection of vertical profiles of the individual terms. The buoyancy term is used to indicate IBL heights because it changes sign at the boundary between different stability regimes. Measured and simulated dissipation rates show a decrease in the transition zone as the air travels over water and an abrupt increase when the IBL over a downwind city area is intersected. The top of the stable IBL is characterized by a minimum in the vertical TKE profile.     

Modelling of turbulence over low hills and varying surface roughness

Neutral surface layer flow over low hills and varying surface roughness is considered with emphasis on closure schemes in relation to the prediction of turbulence quantities. The equations are linearised, Fourier transformed in the two horizontal directions and solved by means of a finite difference method in the vertical. Three closure schemes are. employed, namely mixing length, E- and e-- closure where E, and indicate that differential equations are used for turbulent kinetic energy, dissipation rate and shear stress. Model calculations are compared with experimental data for the step in roughness problem and for the Askervein hill. The mean flow results turn out to be relatively insensitive to the closure scheme. The shear stress and the dimensionless shear, however, are much better predicted with the E- equations than with mixing length closure. In the outer layer of the hill problem, advection of shear stress becomes important. An equation for is needed here.     

Evaluation of an urban land surface scheme over a tropical suburban neighborhood

A two-step statistical downscaling method has been reviewed and adapted to simulate twenty-first-century climate projections for the Gulf of Fonseca (Central America, Pacific Coast) using Coupled Model Intercomparison Project (CMIP5) climate models. The downscaling methodology is adjusted after looking for good predictor fields for this area (where the geostrophic approximation fails and the real wind fields are the most applicable). The method’s performance for daily precipitation and maximum and minimum temperature is analysed and revealed suitable results for all variables. For instance, the method is able to simulate the characteristic cycle of the wet season for this area, which includes a mid-summer drought between two peaks. Future projections show a gradual temperature increase throughout the twenty-first century and a change in the features of the wet season (the first peak and mid-summer rainfall being reduced relative to the second peak, earlier onset of the wet season and a broader second peak).     

Surface fluxes estimation over agricultural areas. Comparisonof methods and the effects of land surface inhomogeneity

Summary  The Bowen Ratio-Energy Balance (BREB) and the aerodynamic method were used to estimate turbulent fluxes of sensible and latent heat flux over an irrigated agricultural area (IAA) and over two dry agricultural areas (DAA1 and DAA2). These turbulent fluxes were analysed and particular attention paid to two specific areas. First, a quantitative analysis of sensible and latent heat fluxes obtained by the BREB method was carried out, taking into account different soil type, vegetation and surface conditions. The results showed that in IAA latent heat flux was higher than sensible heat flux, except in summer months, while in DAA1 and DAA2, sensible heat flux was higher except in the months when the vegetation was at the stage of maximum development. Second, sensible and latent heat fluxes estimates from the BREB method were compared with those obtained from the aerodynamic method. In this comparison factors such as soil type, soil vegetation cover, homogeneity or inhomogeneity of terrain and mesoscale effects such as orography and wind patterns were taken into account. The results show that in conditions of light wind, the two methods only concur if the condition of horizontal homogeniety is fulfilled. The influence of inhomogeneity seems to decrease and agreement between methods improves, if the wind is stronger and the effects of meso and synoptic scales are predominant. Received May 18, 1999/Revised March 15, 2000     

Statistical description of anisotropic and inhomogeneous turbulence

Statistical functions such as correlation functions and spectra yield length scales, each of which describes a certain aspect of the large-scale structure of turbulence. In this study not only statistics of wind fluctuations measured at the same place have been analysed but also statistics of fluctuations at different heights. Thus a spatial picture of the energy-containing and the momentum-transporting eddies has been obtained and some features of the momentum transport are described.     

Characteristics of atmospheric turbulence in the surface layer over Antarctica

This paper presents meteorological measurements made during the antarctic summer period, on two 9 m and 3 m towers, on the rocky and ice shelf terrains of the Indian antarctic stations Maitri and Dakshin Gangotri, respectively. The measurements of fluctuations in temperature and wind speed made with relatively lesser precision instrumentation pertain to smaller wave numbers ~10^-2 m^-1 appropriate to outer scale L ₀ of the atmospheric turbulence spectrum. Autocorrelation analysis of the fluctuations in temperature and wind speed has been performed. A new autoregressive scheme has been developed to represent the computed autocorrelation functions by a Yule statistical model, and to estimate the correlation period T ₀ of the turbulent medium. Height profiles of outer scale L ₀ of turbulence may be given in terms of T ₀ and mean wind speed u. Further, the similarity theory of Monin-Obukhov has been used to compute height profiles of temperature structure parameter C _T ². At Maitri, values of L ₀ and C _T ² are higher between 03–22 h local time than between 22–03 h. Values of L ₀ and C _T ² are smaller over the ice shelf terrain of the Dakshin Gangotri station, compared to those over the rocky terrain of the Maitri station.     

Estimation of the average surface heat flux over an inhomogeneous terrain from the vertical velocity variance

An indirect method of estimating the surface heat flux from observations of vertical velocity variance at the lower mid-levels of the convective atmospheric boundary layer is described. Comparison of surface heat flux estimates with those from boundary-layer heating rates is good, and this method seems to be especially suitable for inhomogeneous terrain for which the surface-layer profile method cannot be used.     

Stability dependence of the drag and bulk transfer coefficients over a coastal sea surface

Bulk transfer coefficients were evaluated from eddy correlation flux measurements on a fixed pier during onshore winds. The mean values are C _D = 1.69 × 10^-3, C _H = 2.58 × 10^-3 and C _E = 1.51 × 10^-3. The drag coefficient, C _D, gradually increases with wind speed but C _H and C _E are independent of wind speed. According to theory and empirical formulas based on experimental results over flat grassland, the transfer coefficients should gradually increase with increasing instability. This is confirmed experimentally in the stable region in our case. However, the drag coefficient appears to decrease with increasing instability, which is against the theoretical result. A stability dependence is not clearly observed for C _H or C _E.     

Surface energy balance complexity in GCM land surface models

 A large number of land surface models (LSMs) have been designed for use in atmospheric general circulation models (GCMs) and GCM modellers therefore have a large number of options when selecting an LSM for their GCM. This study provides information to aid LSM design choices. A framework within which sensitivity to LSM design can be tested is presented and a series of experiments carried out to investigate how general aspects of surface energy balance parametrisation affect land-atmosphere evaporation. Firstly, it is shown that a combination of surface energy balance complexity and aerodynamic parametrisation can be used to explain the gross simulation differences obtained in the Project for Intercomparison of Land-surface Parameterization Schemes (PILPS). Secondly, a simple surface energy balance parametrisation with a constant surface resistance is found to be as appropriate as more complex method for simulating annual, monthly and seasonally averaged diurnal cycles of evaporation. However, complex aspects of surface energy balance parametrisation (canopy interception, bare ground evaporation and canopy resistance) are shown to contain substantial geographic and daily functionality that is not present in the simpler parametrisation. Received: 12 June 1998/Accepted: 24 December 1998     

Turbulence and diffusion over inhomogeneous terrain

This paper provides an overview of some aspects of atmospheric boundary-layer dispersion processes over homogeneous and complex terrain. Special emphasis is placed on a discussion of the boundarylayer scaling regimes over homogeneous terrain and the characteristics of the dispersion processes associated with each of these regimes. The paper points out that vertical concentration profiles usually deviate substantially from a Gaussian distribution. The mean flow and turbulence over a low hill is dealt with, and in the inner layer the turbulence levels are increased due to the mean flow speed-up. In the outer layer the turbulence is modified by the rapid distortion effect. In a middle layer the turbulence is reduced due to the effect of a hill-induced streamline curvature. The paper concludes that the flow perturbations introduced by large-scale hills and valleys invalidate the use of simple approximations for describing atmospheric dispersion processes, and that it is necessary to utilize the full set of equations of motion.     

A scheme for pixel-scale aerodynamic surface temperature over hilly land

Hilly-land satellite pixel-scale aerodynamic surface temperatures (AdST) are investigated using LAS (Large Aperture Scintillometer) and meteorological observations during 21-22 May 2001, indicating that the calculated temperatures are predominantly subject to estimated roughness lengths and, to a less extent, to estimated Bowen ratios, with errors to within 3.0 K between the AdST calculations and hilly radiometric surface temperatures retrieved from satellite data with the split window model. The errors depend heavily on the model used and the zenith angles and azimuth of the satellite and sun with respect to the observational site.     

A rapid-response 2-D drag anemometer for atmospheric turbulence measurements

A rapid response drag anemometer for measuring streamwise and lateral components of horizontal windspeed is described. Theory of operation, design and calibration are discussed with emphasis on the electronic preconditioning of signals and problems associated with using a mechanically resonant system as a sensor. Field comparisons showed half-hourly means and standard deviations of the streamwise component to be within 8% and 5% of respective values obtained from a 3-dimensional sonic anemometer. The lateral component from the drag anemometer was significantly more noisy than that from the 3-D sonic due to induced oscillations arising from vortex shedding. After mechanical and electronic filtering, half-hourly standard deviation comparisons agreed to within 6% for this component. Friction velocities obtained from the drag anemometer in combination with a 1-D sonic, agreed with measurements from the 3-D sonic anemometer to within 4% over a measured range of 0.05 to 1.2 m s^-1     

Turbulent variance characteristics of temperature and humidity over a non-uniform land surface for an agricultural ecosystem in China

This paper describes the application of the variance method for flux estimation over a mixed agricultural region in China. Eddy covariance and flux variance measurements were conducted in a near-surface layer over a non-uniform land surface in the central plain of China from 7 June to 20 July 2002. During this period, the mean canopy height was about 0.50 m. The study site consisted of grass (10% of area), beans (15%), corn (15%) and rice (60%). Under unstable conditions, the standard deviations of temperature and water vapor density (normalized by appropriate scaling parameters), observed by a single instrument, followed the Monin-Obukhov similarity theory. The similarity constants for heat (CT ) and water vapor (Cq) were 1.09 and 1.49, respectively. In comparison with direct measurements using eddy covariance techniques, the flux variance method, on average, underestimated sensible heat flux by 21% and latent heat flux by 24%, which may be attributed to the fact that the observed slight deviations (20% or 30% at most) of the similarity “constants” may be within the expected range of variation of a single instrument from the generally-valid relations.