Stability classification by acoustic remote sensing

Two different Doppler acoustic sounders have been operated at the Kernforschungszentrum Karlsruhe (KfK) since 1982. It has been investigated whether meteorological data from these sounders can be used for dispersion modeling and monitoring in the environment of pollutant-emitting plants. Data from the sounders and from a 200 m high meteorological tower have been sampled continuously for intercomparison.Two schemes of stability classification are presented. They are based on 30-min mean values of the following meteorological data measured by the acoustic sounders: (a) standard deviation σ_w of the vertical wind speed and horizontal wind speed u, at a height of 100 m; and (b) standard deviation σ_φ of the vertical wind direction at a height of 100 m and vertical profile of the backscattered amplitude A_w.The class limits applied in these schemes are determined by “statistical equivalence” with a standard classification scheme. This standard scheme is based on σ_φ, measured by a vector vane at the 100 m level of the tower. Statistical equivalence in this context means that the frequency distributions of the classes are approximately equal at the same site and during the same period.The reliability of these schemes is investigated and compared to the standard scheme by correlation analysis. Finally, the schemes are compared with other commonly applied classification methods.     

Vertical divergence of fogwater fluxes above a spruce forest

Two almost identical eddy covariance measurement setups were used to measure the fogwater fluxes to a forest ecosystem in the “Fichtelgebirge” mountains (Waldstein research site, 786 m a.s.l.) in Germany. During the first experiment, an intercomparison was carried out with both setups running simultaneously at the same measuring height on a meteorological tower, 12.5 m above the forest canopy. The results confirmed a close agreement of the turbulent fluxes between the two setups, and allowed to intercalibrate liquid water content (LWC) and gravitational fluxes. During the second experiment, the setups were mounted at a height of 12.5 and 3 m above the canopy, respectively. For the 22 fog events, a persistent negative flux divergence was observed with a greater downward flux at the upper level. To extrapolate the turbulent liquid water fluxes measured at height z to the canopy of height h_c, a conversion factor 1/[1+0.116(z−h_c)] was determined. For the fluxes of nonvolatile ions, no such correction is necessary since the net evaporation of the fog droplets appears to be the primary cause of the vertical flux divergence. Although the net evaporation reduces the liquid water flux reaching the canopy, it is not expected to change the absolute amount of ions dissolved in fogwater.     

The reflection and diffraction of internal waves from the junction of a slit and a half-space, with application to submarine canyons

We consider the three-dimensional reflection and diffraction properties of internal waves in a continuously stratified rotating fluid which are incident on the junction of a vertical slit and a half-space. This geometry is a model for submarine canyons on continental slopes in the ocean, where various physical phenomena embodying reflection and diffraction effects have been observed. Three types of incident wave are considered: (1) Kelvin waves in the slit (canyon); (2) Kelvin waves on the slope; and (3) plane internal waves incident from the half-space (ocean). These are scattered into Kelvin and Poincaré waves in the slit, a Kelvin wave on the slope and Poincaré waves in the half-space. Most of the discussion is centered around case (1). Various properties of the wave field are calculated for ranges of the parameters c/cot θ, γα and ƒ/ω where cot θ is the topographic slope, c is the internal wave ray slope, α is the canyon half-width, γ is the down-slope wave-number, ƒ is the Coriolis parameter and ω is the wave frequency. Analytical results are obtained for small γα and some approximate results for larger values of γα. The results show that significant wave trapping may occur in oceanic situations, and that submarine canyons may act as source regions for internal Kelvin waves on the continental slope.     

Stratified flow over three-dimensional topography

In order to investigate flows over topography in an atmospheric context, we have studied experimentally the wake structure of axi-symmetric Gaussian obstacles towed through a linearly stratified fluid. Three dimensionless parameters govern the flow dynamics: F, the Froude number based on the topography height h; Re, the Reynolds number and the aspect ratio r = h/L, where L is the topography horizontal scale. Two-dimensional (2-D), saturated lee wave (SLW) and three-dimensional (3-D) regimes, as defined in Chomaz et al. (1993), are found to be functions of F and r only (Fig. 1) as soon as Re is larger than Re_c ≈ 2000. For F < 0.7 the flow goes around the obstacle and the motion in the wake is quasi-two-dimensional. This 2-D layer is topped by a region affected by lee wave motions with amplitude increasing with r and F. For 0.7 < F < 1/r, the flow is entirely dominated by a lee wave of saturated amplitude which suppresses the separation of the boundary layer from the obstacle. Above the critical value 1/r, the lee wave amplitude decreases with F and a recirculating zone appears behind the obstacle. Simultaneously, coherent large-scale vortices start to be shed periodically from the wake at a Strouhal number which decreases as 1/F until it reaches its neutral asymptotic value.     

Vortex–ridge interaction in a rotating fluid

The evolution of barotropic vortices interacting with a topographic ridge on a f-plane is studied by means of laboratory experiments in a rotating tank and numerical simulations. The initial condition in all experiments is a cyclonic vortex created at a certain distance from the ridge. The results are presented in two main scenarios: (a) weak interactions, which occur at early stages of the experiments, when the vortex is far from the ridge, and thus weakly experiences the influence of the topography. In these situations, the vortex slowly drifts towards the ridge with a leftward inclination due to the ascending slope of the topography. Such a behaviour is similar to the “northwestern” motion of cyclones over a weak sloping bottom. The circular shape of the monopolar vortex is preserved. (b) Strong interactions, in which the vortex core reaches the ridge and presents a more complicated evolution. The cyclone “climbs” to the top of the topography and crosses to the other side. Once the vortex experiences the opposite slope, it moves backwards trying to return to the original side of the ridge. For strong enough vortices, this process may be repeated a number of times until the vortex is dissipated by viscous effects. During these interactions the shape of the vortex is strongly deformed and several filaments are produced. In some cases the vortex is cleaved in two parts when crossing the ridge, one at each side of it and moving in opposite directions.Weak and strong interactions are numerically simulated by using a quasi-two-dimensional model. The results confirm that the vortex behaviour is governed by stretching and squeezing effects associated with changes in depth over the ridge and, at latter stages, by Ekman damping due to the solid bottom. The main results observed during strong interactions on a f-plane are also found on preliminar topographic β-plane experiments.     

Spectral scaling in a high reynolds number laboratory boundary layer

Spectra and co-spectra of the streamwise (u) and normal or vertical (w) velocity fluctuations have been measured in the inner region of a large Reynolds number laboratory boundary layer over a rough wall. There is reasonable evidence of ak ₁ ^–1 range in theu spectrum (wherek ₁ is the streamwise wavenumber). Such a range results from an overlap between a spectral region dominated by largescale, inactive motion, which scales on the boundary-layer thickness, and a region dominated by smaller-scale, active motion which scales on the distance from the wall. Spectra ofw, anduw cospectra, scale in a manner consistent with the dominance by active motion. The present spectral data do not support local isotropy over the inertial subrange. A comparison between measuredw spectra and calculations based on isotropy indicates that the inertial subrange anisotropy is only slightly affected by the magnitude of the non-dimensional mean shear.     

Observational evidence for the Monin-Obukhov similarity under all stability conditions

Data collected in the surface layer in a northern suburban area of Nanjing from 15 November to 29 December 2007 were analyzed to examine the Monin-Obukhov similarity for describing the turbulent fluctu- ations of 3D winds under all stability conditions and to obtain the turbulence characteristics under different weather conditions. The results show that the dimensionless standard deviations of turbulent velocity com- ponents (σ u /u* , σ v /u* , σ w /u * ) and dimensionless turbulent kinetic energy (TKE) can be well described by "1/3" power law relationships under stable, neutral, and unstable conditions, with σ u /u * > σ v /u * > σ w /u* . Land use and land cover changes mainly impact dimensionless standard deviations of horizontal component fluctuations, but they have very little on those of the vertical component. The dimensionless standard devi- ations of wind components and dimensionless TKE are remarkably affected by different weather conditions; the deviations of horizontal wind component and dimensionless TKE present fog day > clear sky > overcast > cloudy; the trend of the vertical wind component is the reverse. The surface drag coefficient at a Nan- jing suburban measurement site during the observation period was obviously higher than at other reported plains and plateau areas, and was approximately one order larger in magnitude than the reported plains areas. Dimensionless standard deviation of temperature declined with increasing |z /L| with an approximate "-1/3" slope in unstable stratification and "-2/3" slope in stable stratification.     

A laboratory study of the mechanism of the oxygen airglow

Further laboratory studies of emission by O(¹ S) and by O₂ A ³ _u ⁺ ,A³ _u andc ¹ _u ^– in the oxygen afterglow lead to the conclusion that Barth's mechanism for the excitation of the auroral green line O ₂ ^* +O(³P=O₂+O(¹S)–(1) is correct and that levelsv=6 and 7 of O₂ A ³ _u ⁺ are Barth precursors. The value ofk ₁=7×10^–11 cm³ s^–1 deduced for these levels is shown to be in fair agreement with atmospheric measurements.     

Measurements of spray at an air-water interface

The spray content in the surface boundary layer above an air—water interface was determined by a series of measurements at various feteches and wind speeds in a laboratory facility. The droplet flux density N(z) can be described in terms of the scaling flux density N_* and von Karman constant K throguh the equation, N(z)/N_* = −(1/K) ln(z/z_0d) where z is height above the mean water level and z_0d is the droplet boundary layer thickness. N_* is given by a unique relationship in terms of the roughness Reynolds number u_*σ/ν where σ is the root-mean-square surface displacement. Spray inception occurred for u_* 0.3. The dominant mode of spray generation in the present and most other laboratory tests, as well as in available field data, appears to be bubble bursting.     

Air-sea bulk transfer coefficients in diabatic conditions

On the basis of recent data for the roughness Reynolds number of the sea surface, and using the Owen-Thomson theory on the transfers of heat and mass between a rough surface and the flow above it, the bulk transfer coefficients of the sea surface have been estimated. For a reference height of 10 m, the neutral-lapse transfer coefficient for water vapor is larger by only a few percent than that for sensible heat. When the wind speed at the 10-m height is u ₁₀>3 m s^–1, the coefficient for sensible heat C _H is larger by about 10% than that for momentum C _D . For u ₁₀<5 m s^–1, however, the value of C _D exceeds the value of C _H , and for u ₁₀=15 m s^–1 it is shown that C _H 0.8C _D . It may be also proposed that 10³ C _D =1.11 to 1.70, 10³ C _E =1.18 to 1.30, and 10³ C _H =1.15 to 1.26 for a range of u ₁₀=4 to 20 m s^–1. A plot of diabatic transfer coefficients versus wind speed is obtained by using a parameter of the sea-air temperature difference. For practical purposes, the coefficients are approximated by empirical formulae.     

On Transpiration and Soil Moisture Content Sensitivity to Soil Hydrophysical Data

Sensitivity of evapotranspiration E and root zone soil moisture content θ to the parameterization of soil water retention Ψ(θ) and soil water conductivity K(Ψ), as well as to the definition of field capacity soil moisture content, is investigated by comparing Psi1-PMSURF and Theta-PMSURF models. The core of PMSURF (Penman–Monteith Surface Fluxes) consists of a 3-layer soil moisture prediction module based on Richard’s equation in combination with the Penman–Monteith concept for estimating turbulent heat fluxes. Psi1- PMSURF and Theta-PMSURF differ only in the parameterization of the moisture availability function F_ma. In Psi1,F_ma is parameterized by using Ψ(θ) and K(Ψ) hydrophysical functions; in Theta, F_ma is parameterized by using hydrophysical parameters: the field capacity θ_f and wilting point θ_w soil moisture contents. Both Psi1 and Theta are based on using soil hydrophysical data, that is, there is no conceptual difference between them in the parameterization of E even if in Psi1F_ma depends on 12 parameters, while in Theta only on two soil/vegetation parameters. Sensitivity tests are performed using the Cabauw dataset. Three soil datasets are used: the vG (van Genuchten), CH/vG (Clapp and Hornberger/van Genuchten) and CH/PILPS (Clapp and Hornberger/Project for Intercomparison of Land-surface Parameterization Schemes) datasets. The vG dataset is used in van Genuchten’s parameterization, while in Clapp and Hornberger’s the CH/vG and CH/PILPS datasets are used. It is found that the consistency of soil hydrophysical data in the simulation of transpiration is quite important. The annual sum of E obtained by Psi1E^Psi1, differs from the annual sum of E obtained by Theta, E^Theta, because of the inconsistency between the fitting parameters of Ψ(θ) and K(Ψ) and the θ_f, and not because of the differencies in the parameterization of F_ma. Further, θ_f can be estimated not only on the basis of using soil hydrophysical functions (the θ_f so obtained is θ^Soil_f) but also on the basis of analysing the transpiration process (the θ_fso obtained is θ^tr_f). θ^tr_f values estimated from the condition E^Theta ≈ E^Psi1 are in acceptable accordance with the θ^Soil_f values proposed by Wösten and co-workers. The results are useful in optimizing the parameterization of transpiration in land-surface schemes.     

Vertical Structure of Selected Turbulence Characteristics above an Urban Canopy

Summary In this paper the results of an urban measurement campaign are presented. The experiment took place from July 1995 to February 1996 in Basel, Switzerland. A total of more than 2000 undisturbed 30-minute runs of simultaneous measurements of the fluctuations of the wind vector u′, v′, w′ and the sonic temperature θ _s ′ at three different heights (z=36, 50 and 76 m a.g.l.) are analysed with respect to the integral statistics and their spectral behaviour. Estimates of the zero plane displacement height d calculated by the temperature variance method yield a value of 22 m for the two lower levels, which corresponds to 0.92 h (the mean height of the roughness elements). At all three measurement heights the dimensionless standard deviation σ _w /u _* is systematically smaller than the Monin-Obukhov similarity function for the inertial sublayer, however, deviations are smaller compared to other urban turbulence studies. The σ_θ/θ_* values follow the inertial sublayer prediction very close for the two lowest levels, while at the uppermost level significant deviations are observed. Profiles of normalized velocity and temperature variances show a clear dependence on stability. The profile of friction velocity u _* is similar to the profiles reported in other urban studies with a maximum around z/h=2.1. Spectral characteristics of the wind components in general show a clear dependence on stability and dimensionless measurement height z/h with a shift of the spectral peak to lower frequencies as thermal stability changes from stable to unstable conditions and as z/h decreases. Velocity spectra follow the −2/3 slope in the inertial subrange region and the ratios of spectral energy densities S _w (f)/S _u (f) approach the value of 4/3 required for local isotropy in the inertial subrange. Velocity spectra and spectral peaks fit best to the well established surface layer spectra from Kaimal et al. (1972) at the uppermost level at z/h=3.2. Received September 26, 1997 Revised February 15, 1998     

Structure of the atmospheric surface layer over an industrialized equatorial area

This paper is written to report observations of the structure of the atmospheric surface layer over a coastal industrialized equatorial area. The observations were recorded at Prai Industrial Park, Penang (5° 22′ N, 100° 23′ E) a relatively simple terrain area during the south-west monsoon season in the period of three months using slow response systems. The limitations of the instruments used and its effects on the results are discussed. Wind turbulence and temperature were measured on a 10 m tower and analyzed using eddy correlation method and Monin–Obukhov similarity relations to obtain the normalized standard deviation of longitudinal (σ_u/u), lateral (σ_v/u) and vertical wind velocity fluctuations (σ_w/u) with respect to stability parameter z/L. From the results of the analysis, we found that most of turbulence is generated by shear or mechanical force. It was found that the average neutral value of σ_u/u is 2.35, 1.98 for σ_v/u and 1.47 for σ_w/u with a significantly lower than the proportionality to the power of 1/3 during unstable atmospheric conditions, and thus do not obey Monin–Obukhov similarity theory. It was observed that σ_u/u and σ_v/u values increase linearly in the range of 0 < z/L < 2 and fairly well correlated while σ_w/u does not.     

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.     

Large-Eddy Simulation of Flow and Pollutant Transports in and Above Two-Dimensional Idealized Street Canyons

A large-eddy simulation (LES) model, using the one-equation subgrid-scale (SGS) parametrization, was developed to study the flow and pollutant transport in and above urban street canyons. Three identical two-dimensional (2D) street canyons of unity aspect ratio, each consisting of a ground-level area source of constant pollutant concentration, are evenly aligned in a cross-flow in the streamwise direction x. The flow falls into the skimming flow regime. A larger computational domain is adopted to accurately resolve the turbulence above roof level and its influence on the flow characteristics in the street canyons. The LES calculated statistics of wind and pollutant transports agree well with other field, laboratory and modelling results available in the literature. The maximum wind velocity standard deviations σ _i in the streamwise (σ _u ), spanwise (σ _v ) and vertical (σ _w ) directions are located near the roof-level windward corners. Moreover, a second σ _w peak is found at z ≈ 1.5h (h is the building height) over the street canyons. Normalizing σ _i by the local friction velocity u _*, it is found that σ _u /u _* ≈ 1.8, σ _v /u _* ≈ 1.3 and σ _w /u _* ≈ 1.25 exhibiting rather uniform values in the urban roughness sublayer. Quadrant analysis of the vertical momentum flux u′′w′′ shows that, while the inward and outward interactions are small, the sweeps and ejections dominate the momentum transport over the street canyons. In the x direction, the two-point correlations of velocity R _v,x and R _w,x drop to zero at a separation larger than h but R _u,x (= 0.2) persists even at a separation of half the domain size. Partitioning the convective transfer coefficient Ω _T of pollutant into its removal and re-entry components, an increasing pollutant re-entrainment from 26.3 to 43.3% in the x direction is revealed, suggesting the impact of background pollutant on the air quality in street canyons.     

Bulk transfer coefficient over a snow surface

The drag coefficient C _D and the bulk transfer coefficient for sensible heat C _H over a flat snow surface were determined experimentally. Theoretical considerations reveal that C _D depends on the friction velocity u _* as well as on the geometrical roughness h of the snow surface. It is found that C _D increases with increasing u _* and/or h. The dependency of C _H on u _* and h is so small that it is possible to consider C _H as a constant for practical purposes: C _{H, 1} = 2.0 × 10^–3 for a reference height of 1 m. The bulk transfer coefficient for water vapor is estimated at C _{E, 1} = 2.1 × 10^–3 for a reference height of 1 m.     

A model of tidal flushing of an estuary by dipole formation

We discuss a laboratory model of the flushing and mixing between a bay and a sea connected by a narrow channel. Tidal changes in water elevation lead to periodic flow in and out of the bay through the narrow channel. Due to flow separation there is an asymmetry between the flow entering and leaving the channel on different phases of the tide. A consequence of this asymmetry is that every tidal cycle a fraction of fluid is not returned through the channel, so that a passive tracer in the bay is successively flushed away. We consider what happens when the creation of vorticity in the channel leads to the formation of a dipole that can propagate away from the returning flow of the later phase of the tide. When the dimensionless ratio of maximum channel velocity U, width of channel W and tidal period T is such that W/UT<0.13, we find that dipoles can propagate away from the source region without being drawn back into the sink. This process leads to a larger exchange of water between the estuary and sea, than the classical model of tidal flushing by the jet-sink asymmetry of [Stommel, H., Farmer, H. G., 1952. On the nature of estuarine circulation WHOI Tech. Rep. 52–88, 131].     

Laboratory modeling of topographic Rossby normal modes

The physical modeling of topographic Rossby normal modes carried out at the “Coriolis” Rotating Platform (Grenoble), is presented. The basic feature of the bottom topography is a linear slope of 4.3 m×2 m delimited by two lateral walls. Since the studied motions are essentially barotropic, homogeneous water was used. Unsheared currents were generated by a simple movement of a wavemaker located in front of the topographic barrier. The conservation of potential vorticity for the currents flowing onto the channel slope produced Rossby waves: reflections at the lateral boundaries then led to the formation of propagating barotropic Rossby normal modes, whose frequencies and spatial structures were selected by the physical system. The currents were measured through the correlation imaging velocimetry (CIV) method, which allowed an extremely detailed synoptic map of the horizontal velocities in an area (13 m²) including the slope to be obtained every 30 s.A variety of experiments were performed in order to provide a complete process study in which the effect of different channel lengths and rotation periods could be tested. Two different lengths of the linear slope, 4.3 and 3.3 m, and rotation periods ranging from 30 to 50 s were considered. The qualitative analysis of the 2D current patterns, and the good agreement found between the measured eigenperiods and the periods obtained by means of a simple analytical model, show that in all cases the first Rossby normal mode was generated. Moreover, numerical simulations based on the shallow-water equations, for a geometry and paddle movements that match closely the experimental setup, allow to calibrate the analytical model and provide useful information on a discrepancy found between experimental and analytical eigenperiods due to an oscillation of the normal mode trajectory.     

Small-particle concentration fluctuations at a coastal site

Aerosol size spectra (d=10 nm–10 μm) were measured with an electrical aerosol spectrometer (EAS) at Mace Head on the west coast of Ireland. Several small aerosol particle (diameter 10–32 nm) concentration bursts were observed during the measurement period. Relationships between the events, air mass trajectories, tide height, and meteorological parameters are examined. Series of bursts were observed when a spectral transformation due to subsequent particle growth from 10 to 56–100 nm can be identified in an Eulerian experiment. Particle growth rates of between 1 and 3 nm/h were determined. These bursts appear in cold and comparatively clean arctic or polar air masses with temperature and relative humidity fluctuations, and do not correlate with low tide in some cases. These episodes, similar to those frequently found in the continental boundary layer, are thought to occur over a wide area and, for clear detection, require stable airflow for a few days. Elevated small-particle concentration events are more common during low tide or shortly after, and are typically associated with low wind speeds. Here, the increased shore exposure during low tide is thought to influence the nucleation and the subsequent growth of these aerosol particles. The occurrences of the bursts are found to depend on local wind direction. The highest d=10–32 nm particle concentrations appeared for wind sectors furthest from the tidal regions when the wind direction was 150–160°(south-easterly). Most of the events occurred during daytime when solar irradiation is most intense.