The turbulent wind flow over an embankment

Under neutral conditions and with low winds, profiles of mean and turbulent wind components have been measured at various points across an embankment with aspect ratio 0.3. These measurements have been compared with and related to those of undisturbed flow in a horizontal homogeneous area on the windward side. The speed-up ratio, the turbulent and mean kinetic energy and the turbulent shear stress are examined. It is found that the flow stagnates on the windward side, accelerates above the crest, and separates behind the crest. The results show a remarkable dependence on the angle of attack. With an angle smaller than 90 °, the influence of the embankment on the mean wind field is reduced but is increased on the turbulent part, as lateral gustiness components are amplified. With the incoming flow normal to the embankment, maximum turbulence is found on the top of the ridge near the surface but at greater heights farther downwind. The same is true for the shear stress, but only for oblique flow, whereas for normal flow a minimum is found above the crest and a maximum on the windward side. Therefore, with varying angle of attack the embankment acts in different ways on mean wind, turbulent kinetic energy, and turbulent stress. Although the winds were low, all effects are clearly evident in the data.     

Turbulence in an almond orchard: Vertical variations in turbulent statistics

Three-dimensional wind velocity components were measured above and within a uniform almond orchard. Turbulent statistics associated with the turbulent flow inside the canopy are examined in detail. Turbulence in an almond orchard is characterized by relatively high turbulent intensities and large skewness and kurtosis values. These results indicate that the frequency distribution of wind velocity components is non-Gaussian. Conditional sampling of the turbulent measurements show that large, infrequent sweeps provide the predominant mechanism for tangential momentum stress in the canopy crown. Deep inside the canopy, a secondary wind maximum and small, but positive, tangential momentum stresses are observed.     

Observations of fluxes over heterogeneous surfaces

This study analyzes data collected from repeated aircraft runs 30 m over alternating regions of irrigated and dry nonirrigated surfaces, each region on the order of 10 km across, during the California Ozone Deposition Experiment (CODE). After studying the scale dependence of the flow, the variables and their fluxes are decomposed into means for sublegs defined in terms of irrigated and nonirrigated regions and deviations from such subleg means. Since the repeated runs were flown over the same track, compositing the eight flight legs for each of the two days allows partial isolation of the influences of surface heterogeneity and transient mesoscale motions.A variance analysis is carried out to quantify the relative importance of surface heterogeneity and transient mesoscale motions on the variability of the turbulence fluxes. The momentum and ozone fluxes are more influenced by transient mesoscale motions while fluxes of heat, moisture and carbon dioxide are more influenced by surface heterogeneity. The momentum field is also influenced by a quasi-stationary mesoscale front and larger scale velocity gradients.For the present case, the mesoscale modulation of the turbulent flux is numerically more important than the direct mesoscale flux. This spatial modulation of the turbulent fluxes leads to extra Reynolds terms which act to reduce the area-averaged turbulent momentum flux and enhance the area-averaged turbulent heat flux.     

Measurements of the turbulent heat flux over the sea

The turbulent heat flux was measured with two instruments simultaneously over the Baltic Sea by means of the eddy-correlation method. In one observational period, a small but noticeable divergence in heat flux was found, which may be explained by the advection of colder air. The parameterization of heat flux by the bulk method leads to a value for C _Hof 1 × 10^–3.     

Numerical simulation of turbulent convective flow over wavy terrain

By means of a large-eddy simulation, the convective boundary layer is investigated for flows over wavy terrain. The lower surface varies sinusoidally in the downstream direction while remaining constant in the other. Several cases are considered with amplitude up to 0.15H and wavelength ofH to 8H, whereH is the mean fluid-layer height. At the lower surface, the vertical heat flux is prescribed to be constant and the momentum flux is determined locally from the Monin-Obukhov relationship with a roughness lengthz _o=10^–4 H. The mean wind is varied between zero and 5w _*, wherew _* is the convective velocity scale. After rather long times, the flow structure shows horizontal scales up to 4H, with a pattern similar to that over flat surfaces at corresponding shear friction. Weak mean wind destroys regular spatial structures induced by the surface undulation at zero mean wind. The surface heating suppresses mean-flow recirculation-regions even for steep surface waves. Short surface waves cause strong drag due to hydrostatic and dynamic pressure forces in addition to frictional drag. The pressure drag increases slowly with the mean velocity, and strongly with /H. The turbulence variances increase mainly in the lower half of the mixed layer forU/w _*>2.     

Simplified expressions for adjusting higher-order turbulent statistics obtained from open path gas analyzers

When density fluctuations of scalars such as CO₂ are measured with open-path gas analyzers, the measured vertical turbulent flux must be adjusted to take into account fluctuations induced by ‘external effects’ such as temperature and water vapour. These adjustments are needed to separate the effects of surface fluxes responsible for ‘natural’ fluctuations in CO₂ concentration from these external effects. Analogous to vertical fluxes, simplified expressions for separating the ‘external effects’ from higher-order scalar density turbulence statistics are derived. The level of complexity in terms of input to these expressions are analogous to that of the Webb–Pearman–Leuning (WPL), and are shown to be consistent with the conservation of dry air. It is demonstrated that both higher-order turbulent moments such as the scalar variances, the mixed velocity-scalar covariances, and the two-scalar covariance require significant adjustments due to ‘external effects’. The impact of these adjustments on the turbulent CO₂ spectra, probability density function, and dimensionless similarity functions derived from flux-variance relationships are also discussed.     

The influence of buoyancy on third-order turbulent velocity statistics within a deciduous forest

The influence of atmospheric stability on the behaviour of the third moment of flow velocities observed inside a deciduous forest canopy is examined. Results suggest that buoyancy plays a dominant role in dictating the magnitude of gusts observed inside tall vegetation. Furthermore, an examination of the turbulence recorded throughout leaf fall inside the same forest indicates that larger velocity skewnesses are observed inside a canopy in full leaf than inside a sparse canopy. The behaviour of the measured terms in the non-dimensionalized rate equation of the third moment of canopy flow velocities is also examined. Turbulent diffusion and turbulence gradient interaction terms are largest in stable conditions in the upper canopy layer while these are most important in unstable conditions in the lower canopy layer. In all stability regimes, the turbulent diffusion term is the main source of skewness. The turbulence gradient interaction term, the residual and buoyant production terms all contribute to destroy skewness in stable conditions.     

A wind tunnel study of turbulent flow over model hills

Detailed wind tunnel measurements have been made of mean flow and turbulence over a two-dimensional ridge and a circular hill, both having cosine-squared cross-section and maximum slope about 15 °. The measurements were made in an artificially thickened neutrally stratified boundary layer, and have been compared with results from linear models and rapid distortion theory as appropriate.Our study shows that linear theory gives generally good predictions of the mean flow on the upwind side of the hills, and especially of the flow speedup at the hill top, but that the turbulence is less well predicted. In particular, the measurements show a major increase in the vertical component of turbulence and in the shear stress on the upwind slope of both the two- and three-dimensional hills which is not predicted by either equilibrium or isotropic rapid-distortion theories, although this may be partly due to the effect of streamline curvature. Rapid-distortion theory is successful only in describing the streamwise component of turbulence in the outer region of the flow, while in the upper part of the inner region of the flow, the turbulence measurements show disagreement with both the equilibrium and the rapid-distortion theories. Our experiments also confirm that the equilibrium region is a very thin layer close to the surface, while above this region and below the outer region, there is a transitional region where all terms in the stress equation are important.The measurements over the three-dimensional hill suggest that the mean flow and turbulence are broadly similar to those over the two-dimensional ridge, but with reduced perturbation amplitudes. The major differences between the two cases are found on the upwind slope and in the wake where, respectively, horizontal divergence and convergence of the three-dimensional flow are most pronounced.     

The onset of separation in neutral,turbulent flow over hills

The onset of separation in turbulent, neutrally stratified, boundary-layer flow over hills is considered. Since the flows are fully turbulent, the occurrence of intermittent separation, in the sense of any reversal of near surface flow, will depend strongly on the detailed structure and behaviour of the turbulent eddies. Very little is known about such intermittent separation and the phenomenon cannot be studied with numerical models employing standard turbulence closures; eddy-resolving models are required. Therefore, here, as elsewhere in the literature, the arguably less physically significant process of mean flow separation is studied. Numerical simulations of flow over idealised two- and three-dimensional hills are examined in detail to determine the lowest slope, _crit, for which the mean flow separates.Previous work has identified this critical slope as that required to produce a zero surface stress somewhere over the hill. This criterion, when a mixing-length turbulence closure is applied, reduces to requiring the near-surface vertical velocity shear to vanish at some point on the hill's surface. By applying results from a recent linear analysis for the flow perturbations to this condition, a new expression for _crit is obtained. The expression is approximate but its relative simplicity makes it practically applicable without the need for use of a computer or for detailed mapping of the hill. The approach suggested differs from previous ones in that it applies linear results to a non-linear expression for the surface stress. In the past, a linear expression for the surface stress has been used. The proposed expression for _crit leads to critical angles that are about twice previous predictions. It is shown that the present expression gives good agreement with the numerical results presented here, as well as with other numerical and experimental results. It is also consistent with atmospheric observations.     

A wind tunnel study of turbulent flow over a two-dimensional ridge

We present a wind-tunnel simulation of adiabatic atmospheric flow normal to a rough, two-dimensional ridge. The data are analyzed in physical streamline coordinates, which are described in some detail. The mean velocity speed-up on the hill top is adequately predicted by existing formulae while the behaviour of the wake flow fits into a pattern that emerges from other wind-tunnel experiments. The turbulent stresses evolve in response to the extra strain rates induced by the hill, streamline curvature and acceleration: % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWG1bWaaWbaaS% qabeaaceaIYaGbaebaaaaaaa!3456!\[u^{\bar 2}\]is coupled strongly to acceleration while % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaadaqdaaqaaiaadw% hacaWG3baaaaaa!3462!\[\overline {uw}\]and % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWG3bWaaWbaaS% qabeaaceaIYaGbaebaaaaaaa!3458!\[w^{\bar 2}\]follow curvature. These differing responses lead to significant phase differences between the changes in the component stresses as the hill is traversed. An analogous response is seen in the components of turbulent stress divergence, which are computed as part of streamwise momentum budgets. Only very close to the surface is turbulent stress divergence comparable to the inertial and pressure terms in the momentum budget; over most of the flow regime, the mean flow response is approximately inviscid. Finally, we compare our results with data from other wind tunnel models and from real hills.     

Acoustic tomography and conventional meteorological measurements over heterogeneous surfaces

Summary Conventional micro-meteorological measuring methods are not particularly suitable for the investigation of the energy exchange under heterogeneous surface conditions. To consider the influence of the different surface properties, area-covered and spatially averaged meteorological measurements in combination with highly resolved simulations are necessary. In this context, the method of acoustic travel time tomography is introduced to provide information about the horizontal temperature and wind field.To check the applicability of the tomographic method for this problem and to provide data for the model initialisation a field experiment was carried out.The tomographic system and conventional meteorological equipment were utilised inside the investigation area with an extension of several hundred meters (300×700m²), which was arranged over areas with different surface properties: grassland and bare soil.The results of the field experiment show that differences between the measuring systems (in-situ and remote sensing) exist and the heterogeneity of the underlying surface is visible in the near surface temperature and wind field at a scale, which can be resolved with highly resolved numerical models. Depending on the incoming solar radiation and the local advection regional distinctions in the air temperature and wind field as well as in the vertical sensible heat fluxes were observed.The investigations demonstrate that the sensitivity of the acoustic tomography is sufficient to verify gradients in the meteorological fields even when the horizontal differences in sensible heat flux are small.     

Parameterizing turbulent exchange over sea ice: the ice station weddell results

A 4-month deployment on Ice Station Weddell (ISW) in the western Weddell Sea yielded over 2000 h of nearly continuous surface-level meteorological data, including eddy-covariance measurements of the turbulent surface fluxes of momentum, and sensible and latent heat. Those data lead to a new parameterization for the roughness length for wind speed, z₀, for snow-covered sea ice that combines three regimes: an aerodynamically smooth regime, a high-wind saltation regime, and an intermediate regime between these two extremes where the macroscale or `permanent' roughness of the snow and ice determines z₀. Roughness lengths for temperature, z_T, computed from this data set corroborate the theoretical model that Andreas published in 1987. Roughness lengths for humidity,z_Q, do not support this model as conclusively but are all, on average, within an order of magnitude of its predictions. Only rarely arez_Tand z_Q equal to z₀. These parameterizations have implications for models that treat the atmosphere-ice-ocean system.     

Bulk transfer velocity to and from natural and artificial surfaces

A wide set of published mass and heat transfer data is reviewed in terms of the stochastic renewal theory. A simple model is then proposed describing the transfer in terms of the fluid bulk properties. Two expressions are presented and they enable us to evaluate the bulk transfer velocity as a function of the turbulence of the flow for smooth and fully-rough surfaces in large ranges of Schmidt (or Prandtl) and Reynolds numbers. The transfer over transitional surfaces can also be evaluated by using a simple criterion for the choice of the right expression. The proposed relationships agree satisfactorily with data from laboratory and atmospheric measurements, for both solid and liquid surfaces.     

An estimate of day-time turbulent diffusivity over complex terrain from standard weather data

Summary A method for the calculation of atmospheric turbulent diffusivity over complex terrain during day-time is presented, which may improve the predictions based on diagnostic meteorological models. The proposed procedure takes into account the geographic location of the area (latitude and longitude), the time of the day, the inclination and exposition of the surface, the soil type and the cloud cover. These data are used to compute the amount of solar heat flux contributing to the heating of the air mass above the ground level, and, consequently, the atmospheric turbulence. The model accounts for the effect of shadows generated by mountain profiles, which determine a differential heating at the valley floor and induce spatial and temporal variations of turbulent diffusivity. Model calibration has been performed through ground data collected during a field campaign in the Adige valley in the surrounding of the town of Trento.     

Regional climate simulations over Vietnam using the WRF model

Time-irreversible symmetry is a fundamental property of nonlinear time series. Time-irreversible behaviors of mean temperature measured on 182 stations over China from 1960 to 2012 are analyzed by directed horizontal visibility graph (DHVG for short), and significance of results has been estimated by Monte Carlo simulations. It is found that dominated time irreversibility emerges in nearly all daily temperature anomaly variations over China. Further studies indicate that these time-irreversible behaviors result from asymmetric distributions of persistent daily temperature increments and decrements, and this kind of symmetry can be quantified by distributions of consecutive daily mean temperature increasing or decreasing steps. At the same time, the findings above have been confirmed by artificially generated time series with given value of multiscale asymmetry.     

Regional modeling of decadal rainfall variability over the Sahel

A regional climate model is used to investigate the mechanism of interdecadal rainfall variability, specifically the drought of the 1970s and 1980s, in the Sahel region of Africa. The model is the National Center for Environmental Prediction’s (NCEPs) Regional Spectral Model (RSM97), with a horizontal resolution of approximately equivalent to a grid spacing of 50 km, nested within the ECHAM4.5 atmospheric general circulation model (AGCM), which in turn was forced by observed sea surface temperature (SST). Simulations for the July–September season of the individual years 1955 and 1986 produced wet conditions in 1955 and dry conditions in 1986 in the Sahel, as observed. Additional July–September simulations were run forced by SSTs averaged for each month over the periods 1950–1959 and the 1978–1987. These simulations yielded wet conditions in the 1950–1959 case and dry conditions in the 1978–1987 case, confirming the role of SST forcing in decadal variability in particular. To test the hypothesis that the SST influences Sahel rainfall via stabilization of the tropospheric sounding, simulations were performed in which the temperature field from the AGCM was artificially modified before it was used to force the regional model. We modified the original 1955 ECHAM4.5 temperature profiles by adding a horizontally uniform, vertically varying temperature increase, taken from the 1986–1955 tropical mean warming in either the AGCM or the NCEP/National Center for Atmospheric Research Reanalysis. When compared to the 1955 simulations without the added tropospheric warming, these simulations show a drying in the Sahel similar to that in the 1986–1955 difference and to the decadal difference between the 1980s and 1950s. This suggests that the tropospheric warming may have been, at least in part, the agent by which the SST increases led to the Sahel drought of the 1970s and 1980s.     

Regional variations in fluctuations of seasonal rainfall over Nigeria

Summary Previous work on rainfall variations over Nigeria has concentrated on country-wide averages, which approach has tended to mask the regional contrasts in the country's rainfall. In this paper, variations of rainy season rainfall over the Southern, Middle Belt, and Northern regions of Nigeria as well as the country as a whole are examined over a 72-year secular period (1916–1987). The extent and nature of nonrandom changes, such as fluctuations, trend and persistence, are investigated.The trend analysis showed a tendency towards decreasing seasonal rainfall totals in all the regions, though only those for the Northern region and the country as a whole were significant. No significant persistence was however evident in any of the regions.Power spectrum analysis revealed the occurrence of significant oscillations with time periods of 2.53 to 2.67 and 3.69 to 4.36 years only in the Middle Belt, and 8.00 to 9.60 years in all the regions. The climatological identity of the Middle Belt and the implications of its rainfall oscillations for the large scale agricultural projects planned for the region are pointed out.With 3 Figures     

Regional Climate Change and Uncertainty Analysis based on Four Regional Climate Model Simulations over China

Four sets of climate change simulations at grid spacing of 50 km were conducted over East Asia with two regional climate models driven at the lateral boundaries by two global models for the period 1981–2050. The focus of the study was on the ensemble projection of climate change in the mid-21 st century(2031–50) over China. Validation of each simulation and the ensemble average showed good performances of the models overall, as well as advantages of the ensemble in reproducing present day(1981–2000) December–February(DJF), June–August(JJA), and annual(ANN) mean temperature and precipitation. Significant warming was projected for the mid-21 st century, with larger values of temperature increase found in the northern part of China and in the cold seasons. The ensemble average changes of precipitation in DJF, JJA, and ANN were determined, and the uncertainties of the projected changes analyzed based on the consistencies of the simulations. It was concluded that the largest uncertainties in precipitation projection are in eastern China during the summer season(monsoon precipitation).     

Comparison of areally averaged turbulent fluxes over non-homogeneous terrain: Results from the EFEDA-field experiment

In June 1991 the EFEDA-field experiment (ECHIVAL Field Experiment in a Desertification-Threatened Area) was carried out in the Spanish province Castilla-La Mancha, to improve the understanding of the interactions between the soil, the vegetation and the atmosphere.Here results of energy balance studies at the Barrax site are given, one of the three intensively studied experimental sites within Castilla-La Mancha. This area is characterized by a large fraction of irrigated fields (40%) while the remaining 60% was fallow land at the end of June 1991. The energy balances over these two characteristic land-use classes totally differ. While for the irrigated fields the evapotranspiration is dominant, for the non-irrigated fields the sensible and the soil heat fluxes dominate and the latent heat flux is nearly negligible.In order to achieve areally averaged turbulent fluxes, surface, SODAR and aircraft data have been used. Comparing the surface fluxes from all three facilities, it can be found that:The sensible heat flux estimation from SODAR (_w-method) gives reliable results when a calibration of _w is done with another independent system (e.g. aircraft).Aircraft measurements in conjunction with energy budget methods yield surface fluxes of sensible heat, which are about 20% lower than the areally averaged values calculated by the surface measurements. The areally averaged latent heat fluxes from aircraft and surface measurements agree better than the sensible heat fluxes.